Este trabalho lida com o escoamento incompressível ao redor de pares de cilindros rígidos e imóveis, posicionados de forma alinhada em relação a uma corrente uniforme. São dois os objetivos desta pesquisa: o primeiro é estabelecer relações de causalidade entre características físicas do escoamento e as mudanças observadas nas forças e no campo fluido com as variações do número de Reynolds (Re) e do espaçamento entre os corpos; e o segundo consiste em compreender a influência mútua entre tridimensionalidades e interferência no escoamento. Utilizou-se o método de elementos espectrais para realizar simulações bi e tridimensionais do escoamento. Os espaçamentos entre centros (lcc) analisados vão de 1,5 a 8 diâmetros, e eles são comparados com o caso de um cilindro isolado. A faixa de Re vai de 160 a 320, compreendendo a transição na esteira. O foco foi dado nas instabilidades de pequena escala (modos A e B). Dados referentes ao número de Strouhal, coeficiente de arrasto médio, flutuação do coeficiente de sustentação e correlação axial são apresentados. Com auxílio de visualizações do escoamento, são propostos mecanismos que explicam o fenômeno de interferência, refletido no comportamento das curvas. Os resultados mostraram que simulações bidimensionais são insuficientes para se prever a combinação (Re, lcc) de inversão do arrasto. Verificou-se também que, quando o espaçamento é menor do que o crítico, o processo de transição na esteira se dá de forma diferente do observado para um cilindro isolado. / This work deals with the incompressible flow around pairs of rigid and immovable circular cylinders in tandem arrangements. There are two goals in this research: the first one is to find causality relationships between physical characteristics of the flow and the changes that are observed in the forces and in the flow field with the variation of the Reynolds number (Re) and the distance between the bodies; and the second one is to comprehend the mutual influence between three-dimensional structures and interference. The spectral element method was employed to carry out two- and three-dimensional simulations of the flow. The centre-to-centre distance (lcc) of the investigated configurations varies between 1.5 and 8 diameters, and they are compared to the isolated cylinder case. The Re range goes from 160 to 320, covering the transition in the wake. We focused in the small scale instabilities (modes A and B). Data of Strouhal number, mean drag coefficient, RMS of the lift coefficient and axial correlation are presented. With aid of flow visualizations, we propose mechanisms to explain the interference phenomenon, which is reflected in the behaviour of the graphics. The results show that two-dimensional simulations are not sufficient to predict the (Re, lcc) pair correspondent to the drag inversion point. We also verified that, in the cases where lcc is lower than the critical spacing, the transition in the wake happens in a way different from the one observed in the flow around a single cylinder.

Este trabalho aborda o escoamento ao redor de cilindros flexíveis com grande razão de aspecto, representando os riser utilizados na indústria petrolífera. O escoamento foi dividido em seções bidimensionais ao longo do comprimento do riser. O método dos vórtices discretos foi utilizado para estimativa dos coeficientes hidrodinâmicos nas seções correspondentes e o acoplamento se dá através da dinâmica estrutural. Esta dinâmica é resolvida pelo método dos elementos finitos implementado no código Anflex (Mourelle et al. (2001)). Processamento paralelo é utilizado para acelerar o desempenho do método numérico. Um esquema master-slave (\"mestre-escravo\") utilizando MPI (Message Passing Interface1) é utilizado para que seja explorado o paralelismo da modelagem. As seções hidrodinâmicas são igualmente divididas ao longo dos nós de um cluster de computadores utilizado nos cálculos. Cada nó do cluster (núcleo de processamento) resolve o escoamento nas seções hidrodinâmicas necessárias. As forças hidrodinâmicas atuando nas seções correspondentes são enviadas ao processo mestre que também é responsável pela resolução da dinâmica estrutural. Uma das contribuições do presente trabalho é a possibilidade de análise de risers em regime pós-crítico, através da imposição do ponto de separação da camada-limite baseado em resultados experimentais, e da modelagem de casos com risers em tandem. O Método dos Vórtices Discretos (MVD) está detalhado de forma a facilitar a implementação computacional, levando em consideração a possibilidade de estudo de interferência entre esteiras de múltiplos corpos. Como primeira aproximação, foram simulados, casos de cilindros rígidos fixos e em base elástica para os quais foram examinadas as amplitudes máximas de oscilação e os coeficientes de arrasto e sustentação. Posteriormente, o MVD foi integrado ao Anflex dando origem ao AnflexCFD, com capacidade de processamento paralelo da parte correspondente ao CFD. Os resultados obtidos numericamente foram comparados a resultados experimentais realizados em laboratório. A metodologia de utilização do MVD para cálculo dos coeficientes hidrodinâmicos é baseada na comparação com resultados experimentais obtidos em laboratório para cilindros rígidos fixos e em base elástica para os quais se observa uma aderência satisfatória com relação aos resultados obtidos numericamente. O resultado do presente trabalho foi a criação de uma ferramenta alternativa no projeto de risers de grande razão de aspecto, como uma forma de avaliação do comportamento dinâmico da estrutura submetida a diferentes perfis de corrente e, consequentemente, a estimativa do tempo de sua vida útil. A característica puramente Lagrangiana do MVD possibilita a simulação de centenas de ciclos de oscilação da estrutura sendo que para os métodos Eulerianos tradicionais de CFD, como o método dos volumes finitos e o método dos elementos finitos, a simulação necessita de um tempo computacional inviável para projeto. A aderência do modelo adotado em relação aos resultados experimentais, aliada a eficiência do MVD, é um indicativo da vantagem da utilização de métodos Lagrangianos na modelagem do escoamento ao redor de cilindros flexíveis com grande razão de aspecto, em relação aos métodos Eulerianos tradicionais. A metodologia apresentada neste trabalho pode ser aplicada utilizando outros métodos de CFD, lembrando que isto pode implicar num aumento significativo do tempo computacional necessário para uma análise. / In this work the dynamic response of a high aspect ratio flexible cylinder due to vortex shedding is numerically investigated. The model is divided in two-dimensional sections along the riser length. The discrete vortex method (DVM) is employed for the assessment of the hydrodynamic forces acting on these two-dimensional sections. The hydrodynamic sections are solved independently, and the coupling among the sections is taken into account by the solution of the structure in the time domain by the finite element method implemented in Anflex code (Mourelle et al. (2001)). Parallel processing is employed to improve the performance of the method. A master-slave approach via MPI (Message Passing Interface) is used to exploit the parallelism of the present code. The riser sections are equally divided among the nodes of the cluster. Each node solves the hydrodynamic sections assigned to it. The forces acting on the sections are then passed to the master processor, which is responsible for the calculation of the displacement of the whole structure. One of the main contributions of the present work is the possibility of simulating the flow around flexible cylinders in the pos-critical regime and around bundle of risers.

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

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

Este trabalho aborda o escoamento ao redor de cilindros flexíveis com grande razão de aspecto, representando os riser utilizados na indústria petrolífera. O escoamento foi dividido em seções bidimensionais ao longo do comprimento do riser. O método dos vórtices discretos foi utilizado para estimativa dos coeficientes hidrodinâmicos nas seções correspondentes e o acoplamento se dá através da dinâmica estrutural. Esta dinâmica é resolvida pelo método dos elementos finitos implementado no código Anflex (Mourelle et al. (2001)). Processamento paralelo é utilizado para acelerar o desempenho do método numérico. Um esquema master-slave (\"mestre-escravo\") utilizando MPI (Message Passing Interface1) é utilizado para que seja explorado o paralelismo da modelagem. As seções hidrodinâmicas são igualmente divididas ao longo dos nós de um cluster de computadores utilizado nos cálculos. Cada nó do cluster (núcleo de processamento) resolve o escoamento nas seções hidrodinâmicas necessárias. As forças hidrodinâmicas atuando nas seções correspondentes são enviadas ao processo mestre que também é responsável pela resolução da dinâmica estrutural. Uma das contribuições do presente trabalho é a possibilidade de análise de risers em regime pós-crítico, através da imposição do ponto de separação da camada-limite baseado em resultados experimentais, e da modelagem de casos com risers em tandem. O Método dos Vórtices Discretos (MVD) está detalhado de forma a facilitar a implementação computacional, levando em consideração a possibilidade de estudo de interferência entre esteiras de múltiplos corpos. Como primeira aproximação, foram simulados, casos de cilindros rígidos fixos e em base elástica para os quais foram examinadas as amplitudes máximas de oscilação e os coeficientes de arrasto e sustentação. Posteriormente, o MVD foi integrado ao Anflex dando origem ao AnflexCFD, com capacidade de processamento paralelo da parte correspondente ao CFD. Os resultados obtidos numericamente foram comparados a resultados experimentais realizados em laboratório. A metodologia de utilização do MVD para cálculo dos coeficientes hidrodinâmicos é baseada na comparação com resultados experimentais obtidos em laboratório para cilindros rígidos fixos e em base elástica para os quais se observa uma aderência satisfatória com relação aos resultados obtidos numericamente. O resultado do presente trabalho foi a criação de uma ferramenta alternativa no projeto de risers de grande razão de aspecto, como uma forma de avaliação do comportamento dinâmico da estrutura submetida a diferentes perfis de corrente e, consequentemente, a estimativa do tempo de sua vida útil. A característica puramente Lagrangiana do MVD possibilita a simulação de centenas de ciclos de oscilação da estrutura sendo que para os métodos Eulerianos tradicionais de CFD, como o método dos volumes finitos e o método dos elementos finitos, a simulação necessita de um tempo computacional inviável para projeto. A aderência do modelo adotado em relação aos resultados experimentais, aliada a eficiência do MVD, é um indicativo da vantagem da utilização de métodos Lagrangianos na modelagem do escoamento ao redor de cilindros flexíveis com grande razão de aspecto, em relação aos métodos Eulerianos tradicionais. A metodologia apresentada neste trabalho pode ser aplicada utilizando outros métodos de CFD, lembrando que isto pode implicar num aumento significativo do tempo computacional necessário para uma análise. / In this work the dynamic response of a high aspect ratio flexible cylinder due to vortex shedding is numerically investigated. The model is divided in two-dimensional sections along the riser length. The discrete vortex method (DVM) is employed for the assessment of the hydrodynamic forces acting on these two-dimensional sections. The hydrodynamic sections are solved independently, and the coupling among the sections is taken into account by the solution of the structure in the time domain by the finite element method implemented in Anflex code (Mourelle et al. (2001)). Parallel processing is employed to improve the performance of the method. A master-slave approach via MPI (Message Passing Interface) is used to exploit the parallelism of the present code. The riser sections are equally divided among the nodes of the cluster. Each node solves the hydrodynamic sections assigned to it. The forces acting on the sections are then passed to the master processor, which is responsible for the calculation of the displacement of the whole structure. One of the main contributions of the present work is the possibility of simulating the flow around flexible cylinders in the pos-critical regime and around bundle of risers.

Este trabalho lida com o escoamento incompressível ao redor de pares de cilindros rígidos e imóveis, posicionados de forma alinhada em relação a uma corrente uniforme. São dois os objetivos desta pesquisa: o primeiro é estabelecer relações de causalidade entre características físicas do escoamento e as mudanças observadas nas forças e no campo fluido com as variações do número de Reynolds (Re) e do espaçamento entre os corpos; e o segundo consiste em compreender a influência mútua entre tridimensionalidades e interferência no escoamento. Utilizou-se o método de elementos espectrais para realizar simulações bi e tridimensionais do escoamento. Os espaçamentos entre centros (lcc) analisados vão de 1,5 a 8 diâmetros, e eles são comparados com o caso de um cilindro isolado. A faixa de Re vai de 160 a 320, compreendendo a transição na esteira. O foco foi dado nas instabilidades de pequena escala (modos A e B). Dados referentes ao número de Strouhal, coeficiente de arrasto médio, flutuação do coeficiente de sustentação e correlação axial são apresentados. Com auxílio de visualizações do escoamento, são propostos mecanismos que explicam o fenômeno de interferência, refletido no comportamento das curvas. Os resultados mostraram que simulações bidimensionais são insuficientes para se prever a combinação (Re, lcc) de inversão do arrasto. Verificou-se também que, quando o espaçamento é menor do que o crítico, o processo de transição na esteira se dá de forma diferente do observado para um cilindro isolado. / This work deals with the incompressible flow around pairs of rigid and immovable circular cylinders in tandem arrangements. There are two goals in this research: the first one is to find causality relationships between physical characteristics of the flow and the changes that are observed in the forces and in the flow field with the variation of the Reynolds number (Re) and the distance between the bodies; and the second one is to comprehend the mutual influence between three-dimensional structures and interference. The spectral element method was employed to carry out two- and three-dimensional simulations of the flow. The centre-to-centre distance (lcc) of the investigated configurations varies between 1.5 and 8 diameters, and they are compared to the isolated cylinder case. The Re range goes from 160 to 320, covering the transition in the wake. We focused in the small scale instabilities (modes A and B). Data of Strouhal number, mean drag coefficient, RMS of the lift coefficient and axial correlation are presented. With aid of flow visualizations, we propose mechanisms to explain the interference phenomenon, which is reflected in the behaviour of the graphics. The results show that two-dimensional simulations are not sufficient to predict the (Re, lcc) pair correspondent to the drag inversion point. We also verified that, in the cases where lcc is lower than the critical spacing, the transition in the wake happens in a way different from the one observed in the flow around a single cylinder.

Este trabalho tem por objetivo estudar através de simulações numéricas a interação fluido-estrutura de um escoamento ao redor de um cilindro livre para vibrar na direção transversal do escoamento, com o objetivo de compreender melhor o fenômeno muito recorrente em estruturas oceânicas, principalmente os dutos de extração de petróleo, chamados risers. A condição a ser estudada é a de escoamento com número de Reynolds igual a 10000, em um cilindro com parâmetro de massa igual a 3,3 e de amortecimento igual a 0,0026. Nestas condições, diversos modelos serão criados para representar toda a faixa de velocidades reduzidas em que ocorre o fenômeno de sincronização. Estes parâmetros foram escolhidos devido à presença de dados experimentais para o comparativo entre as simulações e os dados experimentais, que será apresentado ao final do trabalho. Resultados como a média do coeficiente de arrasto e o RMS do coeficiente de sustentação também serão comparados com dados experimentais disponíveis na literatura e uma análise dos perfis de escoamento na região da esteira é feita para a identificação de padrões recorrentes nestes tipos de caso. / The objective of this work is to study through numerical simulations the fluid-structure interaction in a cylinder, which is free to oscillate in the cross flow direction, aiming to better comprehend the phenomena that frequently appears in offshore structures, mainly in pipelines for oil extraction, called risers, The studied condition is a flow with Reynolds 10000 around a cylinder with mass ratio of 3.3 and damping factor of 0.0026. Under these conditions, multiples models are created to represent the reduced velocities range where the synchronization phenomena happen. Those parameters were chosen due to the availability of experimental data in the literature, which allows the comparison between the simulations and experiments that is presented in the end of the work. Results such as the drag coefficient mean and the lift coefficient RMS are also compared with the experimental data available in the literature and an analysis of the flow in the wake region is done with the intention of identifying the patterns found in this type of case.

Este trabalho tem por objetivo estudar através de simulações numéricas a interação fluido-estrutura de um escoamento ao redor de um cilindro livre para vibrar na direção transversal do escoamento, com o objetivo de compreender melhor o fenômeno muito recorrente em estruturas oceânicas, principalmente os dutos de extração de petróleo, chamados risers. A condição a ser estudada é a de escoamento com número de Reynolds igual a 10000, em um cilindro com parâmetro de massa igual a 3,3 e de amortecimento igual a 0,0026. Nestas condições, diversos modelos serão criados para representar toda a faixa de velocidades reduzidas em que ocorre o fenômeno de sincronização. Estes parâmetros foram escolhidos devido à presença de dados experimentais para o comparativo entre as simulações e os dados experimentais, que será apresentado ao final do trabalho. Resultados como a média do coeficiente de arrasto e o RMS do coeficiente de sustentação também serão comparados com dados experimentais disponíveis na literatura e uma análise dos perfis de escoamento na região da esteira é feita para a identificação de padrões recorrentes nestes tipos de caso. / The objective of this work is to study through numerical simulations the fluid-structure interaction in a cylinder, which is free to oscillate in the cross flow direction, aiming to better comprehend the phenomena that frequently appears in offshore structures, mainly in pipelines for oil extraction, called risers, The studied condition is a flow with Reynolds 10000 around a cylinder with mass ratio of 3.3 and damping factor of 0.0026. Under these conditions, multiples models are created to represent the reduced velocities range where the synchronization phenomena happen. Those parameters were chosen due to the availability of experimental data in the literature, which allows the comparison between the simulations and experiments that is presented in the end of the work. Results such as the drag coefficient mean and the lift coefficient RMS are also compared with the experimental data available in the literature and an analysis of the flow in the wake region is done with the intention of identifying the patterns found in this type of case.

O fenômeno de vibração induzida por vórtices em travessas de pré-distribuidores de turbinas hidráulicas tem sido estudado nos últimos anos e várias soluções têm sido adotadas para minimizar interferências na estrutura que podem causar fratura por fadiga. O princípio básico das modificações é alterar o perfil da travessa de modo que as frequências de emissão de vórtices não coincidam com as frequências naturais da estrutura. Este trabalho tem como objetivo avaliar através de uma série de simulações computacionais um perfil mais adequado para um pré-distribuidor de turbina Francis. Essas simulações envolvem o cálculo do escoamento ao redor da travessa e da vibração induzida por vórtices nele presentes, bem como uma técnica que combina as análises dinâmicas com uma otimização paramétrica. Para isso, foi utilizado um código comercial de CFD, ANSYS Fluent e o cálculo da resposta estrutural e seu acoplamento com as equações do escoamento foi feito através de uma UDF (User Defined Function). Para validar a metodologia, a resposta estrutural de um corpo prismático sobre base elástica foi calculada e comparada a dados previamente publicados na literatura. Por fim, um código desenvolvido controla a análise fluido-estrutural e passa as variáveis para o otimizador Mode Frontier, que trabalha para encontrar a estrutura mais eficiente variando-se os parâmetros pré-determinados da geometria da peça. A metodologia desenvolvida tem a vantagem de ajudar no projeto de tais componentes sem depender excessivamente de métodos experimentais ou regras empíricas. Dessa forma, torna possível modificar perfis existentes ou desenvolver perfis novos baseado nos melhores critérios de manufatura. / Vortex induced vibration phenomena in hydraulic turbines stay vanes have been studied in the last years and several solutions have been adopted in order to minimize interferences that can cause fatigue in the structure. The basic principle of all modifications is to change the stay vane profile so the natural vortex shedding frequency is different from the natural frequencies of the structure. This work presents a detailed computational simulation of a Francis turbine stay vane whose main objective is to find out a more suitable profile these components should assume. These simulations involve the calculation of the flow around the vanes and the associated vortex induced vibration in the structure in addition to a technique that combines the dynamic analysis with a parametric optimization In order to do that, a commercial CFD code, ANSYS Fluent, was adopted and the calculation of the structural response and its coupling with the flow equations was done with User Defined Functions. Validation of the methodology was made by comparing the structural response of an elastically-mounted prismatic body immersed in uniform flow with previously published data. Finally, a developed code controls the FSI analysis and provides information about the vibrations to the Mode Frontier optimizer, responsible to address the problem and determine the set of parameters that lead to the most efficient structure. The methodology developed has the advantage of helping the design of such components without depending excessively on experimental methods or empirical rules. Also, it allows either modifying existing profiles or choosing the best shape for new ones based on the best manufacturing criteria.

O fenômeno de vibração induzida por vórtices em travessas de pré-distribuidores de turbinas hidráulicas tem sido estudado nos últimos anos e várias soluções têm sido adotadas para minimizar interferências na estrutura que podem causar fratura por fadiga. O princípio básico das modificações é alterar o perfil da travessa de modo que as frequências de emissão de vórtices não coincidam com as frequências naturais da estrutura. Este trabalho tem como objetivo avaliar através de uma série de simulações computacionais um perfil mais adequado para um pré-distribuidor de turbina Francis. Essas simulações envolvem o cálculo do escoamento ao redor da travessa e da vibração induzida por vórtices nele presentes, bem como uma técnica que combina as análises dinâmicas com uma otimização paramétrica. Para isso, foi utilizado um código comercial de CFD, ANSYS Fluent e o cálculo da resposta estrutural e seu acoplamento com as equações do escoamento foi feito através de uma UDF (User Defined Function). Para validar a metodologia, a resposta estrutural de um corpo prismático sobre base elástica foi calculada e comparada a dados previamente publicados na literatura. Por fim, um código desenvolvido controla a análise fluido-estrutural e passa as variáveis para o otimizador Mode Frontier, que trabalha para encontrar a estrutura mais eficiente variando-se os parâmetros pré-determinados da geometria da peça. A metodologia desenvolvida tem a vantagem de ajudar no projeto de tais componentes sem depender excessivamente de métodos experimentais ou regras empíricas. Dessa forma, torna possível modificar perfis existentes ou desenvolver perfis novos baseado nos melhores critérios de manufatura. / Vortex induced vibration phenomena in hydraulic turbines stay vanes have been studied in the last years and several solutions have been adopted in order to minimize interferences that can cause fatigue in the structure. The basic principle of all modifications is to change the stay vane profile so the natural vortex shedding frequency is different from the natural frequencies of the structure. This work presents a detailed computational simulation of a Francis turbine stay vane whose main objective is to find out a more suitable profile these components should assume. These simulations involve the calculation of the flow around the vanes and the associated vortex induced vibration in the structure in addition to a technique that combines the dynamic analysis with a parametric optimization In order to do that, a commercial CFD code, ANSYS Fluent, was adopted and the calculation of the structural response and its coupling with the flow equations was done with User Defined Functions. Validation of the methodology was made by comparing the structural response of an elastically-mounted prismatic body immersed in uniform flow with previously published data. Finally, a developed code controls the FSI analysis and provides information about the vibrations to the Mode Frontier optimizer, responsible to address the problem and determine the set of parameters that lead to the most efficient structure. The methodology developed has the advantage of helping the design of such components without depending excessively on experimental methods or empirical rules. Also, it allows either modifying existing profiles or choosing the best shape for new ones based on the best manufacturing criteria.