Tese de doutoramento. Engenharia Mecânica. Faculdade de Engenharia. Universidade do Porto. 2007

O presente trabalho tem como escopo validar qualitativa e quantitativamente o ensaio de swirl em cabeçotes de motores Diesel realizado de maneira contínua, comparando-o com o padrão estático amplamente aplicado na indústria no controle de qualidade deste componente. A modificação metodológica do ensaio permite ensaios muito mais velozes do que os atuais, podendo, com todas as verificações necessárias, ser extensível a uma linha de montagem, permitindo controle de qualidade diretamente após o processo de usinagem do fundido ou antes de montar o cabeçote no bloco do motor. O presente trabalho apresenta a metodologia aplicada para o tratamento do problema bem como os resultados obtidos até o momento, indicando também os que ainda serão obtidos. No escopo do projeto, apresenta-se de maneira sucinta todo o procedimento utilizado para as simulações estáticas e transiente dos métodos de medição a partir do método do volume finitos de um código comercial, destacando a metodologia para determinação das malhas, seleção dos modelos e esquemas de discretização utilizados, destacando por fim uma comparação entre diferentes aberturas de válvula do ensaio estático e os primeiros resultados da simulação com abertura de válvula ao longo do tempo de simulação. Do ponto de vista experimental, descreveu-se o equipamento utilizado bem como as adaptações realizadas, principalmente no que tangem a automação do controle de válvula, no escopo do presente projeto, apresentando por fim alguns resultados do ensaio contínuo na bancada em comparação com o estático. Os ensaios e simulações numéricas estática e transientes apresentaram coerência em seus resultados. As simulações, quando comparadas entre si e com os resultados do PIV, mostram que os perfis de escoamento são bastante parecidos em todas as situações, o que permite concluir que existe a possibilidade de se realizar o ensaio de swirl de forma dinâmica para esta configuração específica. / The present work aims to validate qualitatively and quantitatively the swirl test in diesel engine heads performed continuously, comparing it with the static standard widely applied in industry in the quality control of this component. The methodological modification of the test allows for tests much faster than the current ones and can be extensible to an assembly line, allowing quality control directly after the cast machining process or before mounting the engine block. The present work presents the applied methodology for the treatment of the problem as well as the results obtained so far, also indicating the ones that will be obtained. In the scope of the project, the whole procedure used for the static and transient simulations of the measurement methods based on the finite volume method of a commercial code is presented, highlighting the methodology for determining the meshes, selection of models and schemes of discretization. At last, a comparison between different valve openings of the static test and the first results of the simulation with valve opening during the simulation time will be presented. From the experimental point of view, the equipment used was described, as well as the adaptations made, mainly in what concerns the automation of valve control, within the scope of the present project. Some results of the continuous test in the bench are presented and compared with the static ones. The static and dynamic tests and simulations presented good coherence among their results. The simulations, when compared one with each other and to the PIV test results, showed great similarities among the in-cylinder flow velocity profiles, what might lead to the conclusion that it is possible to test cylinder heads dynamically, at least with the presented configuration.

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.

Neste trabalho é proposta uma metodologia numérica como ferramenta para avaliação de ambientes térmicos com manequins. A simulação de CFD do ambiente térmico em simulador comercial é integrada à simulação do sistema térmico do corpo humano realizada em código acadêmico. As soluções dos fluxos de calor e temperaturas nas peles são retroalimentadas e a transferência de informações é realizada via arquivo. A geometria do ambiente térmico é simplificada para minimizar os efeitos de problemas com a malha computacional na simulação de CFD, permitindo uma melhor análise do método interativo proposto. O manequim é separado em 15 segmentos cilíndricos representando cabeça, pescoço, tronco, braços, antebraços, mãos, coxas, pernas e pés. Cada segmento é subdividido em quatro quadrantes, totalizando 60 zonas, para capturar assimetrias térmicas e aerodinâmicas. Foi conduzido estudo de validação geométrica de manequim virtual pela comparação de resultados dos coeficientes térmicos obtidos com CFD e resultados obtidos de ensaios experimentais da literatura. A qualidade da malha e o tratamento de parede são discutidos. Os resultados tornam evidente que uma geometria simplificada do manequim é suficiente para estudos e avaliações de ambiente térmico e de conforto térmico quando se utiliza técnicas numéricas de CFD. Os resultados a partir da integração dos simuladores mostram que o método numérico pode ser instável nos segmentos com baixo metabolismo e baixa vazão de sangue, como nos pés e mãos. Pretende-se introduzir na metodologia proposta algum mecanismo que identifique automaticamente este fenômeno, para evitar a divergência do método e tornar a ferramenta mais robusta. / In this work is proposed a new numerical methodology as a tool for thermal comfort evaluation. This method promotes the interaction of the thermal environment simulation and the thermal system of the human body simulation. The commercial CFD simulator FLUENT R and an academician code for human body simulation are used. The solutions are fed back and the transfer is made by file. The geometry of the room is simplified to minimize the effects of problems with the computational mesh in the CFD simulation, allowing a better analysis of the proposed interactive method. The dummy is separated into 15 cylindrical segments representing head, neck, trunk, arms, forearms, hands, thighs, legs and feet. Each segment is subdivided into four quadrants, totaling 60 zones, to capture asymmetries in the heat flux field and temperature field. Was conducted a geometric validation of virtual dummy by comparing results of heat transfer coefficients from literature and CFD simulation. The mesh quality and near wall treatment are discussed. The results show that a simplified geometry of the dummy is sufficient for thermal environment studies and evaluations in CFD simulations. The results from coupled simulations show that the numerical method can be unstable in the segments with low metabolism and low blood flow, as the feet and hands. So, its intended to introduce mechanisms in the methodology to automatically identify this phenomenon and to avoid the divergence of the method to make more robust this methodology.

Este trabalho tem como objetivo desenvolver ferramenta computacional para simulação e análise do fenômeno de penetração e cravação de estacas torpedos em solo marítimo. A abordagem será baseada no método Moving Particle Semi-Implicit (MPS). Por se tratar de um método de partícula, sem malha, o mesmo apresenta grande flexibilidade na modelagem de problemas de interação fluido-sólido com fragmentação ou junção de superfície livre e grandes deslocamentos ou deformações dos sólidos, fenômenos esses presentes no impacto e cravação da estaca no solo marítimo. Para isso, dois desafios foram elencados: a modelagem dos solos como fluidos não-newtonianos e a determinação da força de arrasto viscosa na superfície de sólidos. A modelagem do fluido não-newtoniano foi feita considerando os modelos de Power Law, Bingham e Herschel-Bulkley. O cálculo da força de arrasto viscosa foi avaliado determinando-se o gradiente da velocidade do fluido na direção normal à parede com base na regressão polinomial. Por simplicidade, foi considerada a hipótese de que a variação da velocidade na direção tangencial da parede é muito menor se comparada a variação da mesma na direção do vetor normal. O método implementado, assim como o escoamento de fluidos não-newtonianos, foi validado por meio de comparação entre o resultado obtido de simulações com geometrias pré-definidas e as respostas analíticas para tais casos. Como exemplo de aplicação da ferramenta computacional desenvolvida, um caso simplificado de cravação das estacas torpedos foi simulado avaliando-se o seu deslocamento dentro do solo e os esforços cisalhantes a ela submetidas. / This work aims to develop computational tools to simulate and analysis the torpedo anchor penetration in marine soil. The approach will be based on the Moving Particle Semi-Implicit (MPS) method. Because it is a meshless method, it is extremely flexible to model fluid-solid interaction with fragmentation or junction of free surface and large displacements or deformations of solids, phenomena presented at the torpedo anchor impact. Two challenges were listed: the modeling of soils as non-Newtonian fluids and the determination of the viscous drag on the solids surface. The modeling of non-Newtonian fluid was done based on the Power Law, Bingham and Herschel-Bulkley models. The calculation of the viscous drag was evaluated by determining the velocity gradient in the normal direction of the wall based on polynomial regression considering the fluid particles near the solid wall. In this work, for sake of simplicity, the hypothesis that the velocity variation in the tangential direction of the wall is much smaller compared to its variation in the normal direction is adopted. The proposed technique, as well as the flow of non-Newtonian fluids, were validated comparing the results obtained in flow simulations with predefined geometries with the expected analytical responses. As an example of the application of the computational tools developed, a simplified case of torpedo penetration was simulated by evaluating its displacement and the shear stresses submitted to it.

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.

O transiente hidráulico na adução de Pequenas Centrais Hidrelétricas (PCH) é um importante tema de estudo para a correta escolha e dimensionamento das estruturas de amortecimento, que protegem o circuito adutor da variação abrupta da pressão interna causada pela variação na velocidade média do escoamento. Este trabalho propõe a utilização da hidrodinâmica computacional (CFD) para avaliar a oscilação de massa ou nível d\'água no interior da estrutura de proteção denominada Chaminé de Equilíbrio e comparar esses resultados com a solução teórica analítica, com os critérios de dimensionamento definidos em norma técnica e com os dados de medições experimentais. O problema foi modelado no software Ansys CFX que utiliza o método dos volumes finitos (MVF) para resolver os campos de velocidade e pressão do escoamento, tanto em regime permanente, quanto em regime transiente. Aplicando-se o modelo CFD proposto ao circuito hidráulico didático, os resultados foram validados com a formulação teórica e com as medições experimentais do nível d\'água na chaminé. O modelo foi aplicado num estudo de caso que consiste num circuito hidráulico de PCH. Os resultados da oscilação de massa no interior da chaminé para este estudo de caso confirmam o emprego da fluidodinâmica computacional para avaliar tais escoamentos em regime transitório. / The hydraulic transient in the circuit of Small Hydroelectric Power Plants (SHP) is an important study subject for the correct choice and design of damping structure which protects the conduit from an abrupt variation of the internal pressure caused by the changes in the average flow velocity. This work proposes the use of computational fluid dynamics (CFD) to evaluate the mass oscillation or water level inside the protection structure called Surge Tank and compares these results with theory, design criteria defined by the technical standard, and experimental data. The problem was modeling in Ansys CFX software that uses the finite volume method (FVM) to solve velocity and pressure flow fields, in both steady state and transient flow. The proposed model used in a hydraulic circuit led to results validation with experimental measurements of water level surge tank. For the proposed case study, the response of the numerical CFD simulation confirms the use of computational fluid dynamics to evaluate such transient flows.

Este trabalho tem como objetivo desenvolver ferramenta computacional para simulação e análise do fenômeno de penetração e cravação de estacas torpedos em solo marítimo. A abordagem será baseada no método Moving Particle Semi-Implicit (MPS). Por se tratar de um método de partícula, sem malha, o mesmo apresenta grande flexibilidade na modelagem de problemas de interação fluido-sólido com fragmentação ou junção de superfície livre e grandes deslocamentos ou deformações dos sólidos, fenômenos esses presentes no impacto e cravação da estaca no solo marítimo. Para isso, dois desafios foram elencados: a modelagem dos solos como fluidos não-newtonianos e a determinação da força de arrasto viscosa na superfície de sólidos. A modelagem do fluido não-newtoniano foi feita considerando os modelos de Power Law, Bingham e Herschel-Bulkley. O cálculo da força de arrasto viscosa foi avaliado determinando-se o gradiente da velocidade do fluido na direção normal à parede com base na regressão polinomial. Por simplicidade, foi considerada a hipótese de que a variação da velocidade na direção tangencial da parede é muito menor se comparada a variação da mesma na direção do vetor normal. O método implementado, assim como o escoamento de fluidos não-newtonianos, foi validado por meio de comparação entre o resultado obtido de simulações com geometrias pré-definidas e as respostas analíticas para tais casos. Como exemplo de aplicação da ferramenta computacional desenvolvida, um caso simplificado de cravação das estacas torpedos foi simulado avaliando-se o seu deslocamento dentro do solo e os esforços cisalhantes a ela submetidas. / This work aims to develop computational tools to simulate and analysis the torpedo anchor penetration in marine soil. The approach will be based on the Moving Particle Semi-Implicit (MPS) method. Because it is a meshless method, it is extremely flexible to model fluid-solid interaction with fragmentation or junction of free surface and large displacements or deformations of solids, phenomena presented at the torpedo anchor impact. Two challenges were listed: the modeling of soils as non-Newtonian fluids and the determination of the viscous drag on the solids surface. The modeling of non-Newtonian fluid was done based on the Power Law, Bingham and Herschel-Bulkley models. The calculation of the viscous drag was evaluated by determining the velocity gradient in the normal direction of the wall based on polynomial regression considering the fluid particles near the solid wall. In this work, for sake of simplicity, the hypothesis that the velocity variation in the tangential direction of the wall is much smaller compared to its variation in the normal direction is adopted. The proposed technique, as well as the flow of non-Newtonian fluids, were validated comparing the results obtained in flow simulations with predefined geometries with the expected analytical responses. As an example of the application of the computational tools developed, a simplified case of torpedo penetration was simulated by evaluating its displacement and the shear stresses submitted to it.

O transiente hidráulico na adução de Pequenas Centrais Hidrelétricas (PCH) é um importante tema de estudo para a correta escolha e dimensionamento das estruturas de amortecimento, que protegem o circuito adutor da variação abrupta da pressão interna causada pela variação na velocidade média do escoamento. Este trabalho propõe a utilização da hidrodinâmica computacional (CFD) para avaliar a oscilação de massa ou nível d\'água no interior da estrutura de proteção denominada Chaminé de Equilíbrio e comparar esses resultados com a solução teórica analítica, com os critérios de dimensionamento definidos em norma técnica e com os dados de medições experimentais. O problema foi modelado no software Ansys CFX que utiliza o método dos volumes finitos (MVF) para resolver os campos de velocidade e pressão do escoamento, tanto em regime permanente, quanto em regime transiente. Aplicando-se o modelo CFD proposto ao circuito hidráulico didático, os resultados foram validados com a formulação teórica e com as medições experimentais do nível d\'água na chaminé. O modelo foi aplicado num estudo de caso que consiste num circuito hidráulico de PCH. Os resultados da oscilação de massa no interior da chaminé para este estudo de caso confirmam o emprego da fluidodinâmica computacional para avaliar tais escoamentos em regime transitório. / The hydraulic transient in the circuit of Small Hydroelectric Power Plants (SHP) is an important study subject for the correct choice and design of damping structure which protects the conduit from an abrupt variation of the internal pressure caused by the changes in the average flow velocity. This work proposes the use of computational fluid dynamics (CFD) to evaluate the mass oscillation or water level inside the protection structure called Surge Tank and compares these results with theory, design criteria defined by the technical standard, and experimental data. The problem was modeling in Ansys CFX software that uses the finite volume method (FVM) to solve velocity and pressure flow fields, in both steady state and transient flow. The proposed model used in a hydraulic circuit led to results validation with experimental measurements of water level surge tank. For the proposed case study, the response of the numerical CFD simulation confirms the use of computational fluid dynamics to evaluate such transient flows.

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.