Estudo experimental do fenômeno de vibração induzida por vórtices em cilindro rígido livre para oscilar com dois graus de liberdade.

O fenômeno de vibração induzida por vórtices (VIV) é um problema fundamental dentro da Mecânica dos Fluidos e um exemplo importante de interação fluido-estrutura. Esta tese investiga fenômeno de VIV quando um cilindro rígido, submetido a escoamento uniforme, está livre para oscilar na direção transversal e alinhada com a corrente incidente. A tese foi estruturada ao redor de sete perguntas relacionadas ao fenômeno de VIV: 1) O fenômeno e os resultados experimentais são repetitivos? 2) Como ocorre a transição entre ramos de resposta? 3) Qual é o papel da inércia da estrutura oscilante? 4) Qual é o papel de sua rigidez? 5) Quais são as frequências naturais mais importantes da estrutura? 6) Quais padrões de esteira se desenvolvem para VIV com dois graus de liberdade? 7) Quais são os efeitos do movimento na direção alinhada com a corrente no processo de formação e desprendimento de vórtices? O fenômeno de VIV é estudado de maneira experimental em uma base elástica pendular capaz de oscilar com o mesmo momento de inércia e frequência natural nas duas direções. Os experimentos de VIV foram realizados em canal de água recirculante e com diferentes condições de inércia e rigidez. A técnica de velocimetria por imagem de partículas foi usada e permitiu identificar diferentes padrões de esteira de vórtices. Verificou-se que o VIV é repetitivo a nível de amplitudes médias e frequências dominantes. A transição dos ramos pode ocorrer de maneira intermitente ou com histerese. Os parâmetros de inércia e rigidez da estrutura são capazes de mudar o regime de oscilação e, para algumas condições, suprimir as vibrações alinhadas com a corrente. Dentre os padrões de esteira observados, um deles não havia sido relatado na literatura e é definido nesta tese. O novo modo de emissão apresenta dois vórtices com circulação oposta e elevada intensidade emitidos por ciclo. A influência da direção alinhada com o escoamento está relacionada a dois efeitos: a velocidade relativa entre o cilindro e o fluido, responsável pelo aumento da circulação dos vórtices na esteira, e o ângulo de fase do movimento nas direções alinhada e transversal, capaz de mudar o processo de formação dos vórtices.

Desenvolvimento da modelagem de turbulência e interação fluido-estrutura para as vibrações induzidas por vórtices de cilindro rígido.

Esta tese apresenta o desenvolvimento e aplicação de modelos de turbulência, transição laminar-turbulenta e de interações fluido-estrutura ao escoamento externo em cilindro rígido estacionário e em vibrações induzidas por vórtices. Tais desenvolvimentos foram realizados no código ReFRESCO, baseado em técnicas de dinâmica de fluidos computacional (CFD). Realizou-se um estudo quanto ao desempenho do modelo k- SST em extensa faixa de números de Reynolds, segundo o qual se identificaram as deficiências de modelagem para este escoamento. A modelagem adaptativa das escalas (SAS) e o modelo de transição por correlações locais (LCTM), ambos combinados ao SST, melhoraram a aderência aos resultados experimentais para este escoamento, em uma contribuição original deste trabalho. A aplicação de técnicas de verificação e validação possibilitou a estimação de incertezas e erros para os modelos e números de Reynolds e também de identificada como outra contribuição deste trabalho. A combinação da modelagem em SST, SAS e LCTM com movimentos impostos de realizada para números de Reynolds moderados, diferentes frequências e amplitudes de vibração, algo que poucas publicações abordam em detalhes. Com relação aos movimentos livres, este trabalho traz contribuições com a aplicação dos modelos SST e SAS ao estudo de vibrações induzidas por vórtices em dois graus de liberdade, baixa razão de massa e números de Reynolds moderados, mais altos do que normalmente observados na literatura. Por fim, a investigação da importância relativa de efeitos da turbulência aos casos de movimentos livres e impostos, com relação ao caso de cilindro estacionário, comprovou a conjetura formulada na parte inicial deste trabalho, no que tange à escolha do modelo de turbulência em determinadas aplicações. Tal escolha mostrou-se menos decisiva no caso do cilindro em movimento imposto e ainda menos nos movimentos livres, em comparação ao caso estacionário, uma vez que a resposta em movimentos do corpo filtra grande parte dos efeitos turbulentos de ordem superior. Esta observação mostra-se relevante, uma vez que pode permitir simplificações na modelagem e aplicação de ferramentas de CFD em uma classe importante de projetos de engenharia.

Método para análise da interação fluido-estrutura em travessas do pré-distribuidor de turbinas hidráulicas.

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.

On the development of advanced techniques for mixed-elastohydrodynamic lubrification modelling of journal and sliding bearing systems.

The present thesis is focused on the development of a thorough mathematical modelling and computational solution framework aimed at the numerical simulation of journal and sliding bearing systems operating under a wide range of lubrication regimes (mixed, elastohydrodynamic and full film lubrication regimes) and working conditions (static, quasi-static and transient conditions). The fluid flow effects have been considered in terms of the Isothermal Generalized Equation of the Mechanics of the Viscous Thin Films (Reynolds equation), along with the massconserving p-Ø Elrod-Adams cavitation model that accordingly ensures the so-called JFO complementary boundary conditions for fluid film rupture. The variation of the lubricant rheological properties due to the viscous-pressure (Barus and Roelands equations), viscous-shear-thinning (Eyring and Carreau-Yasuda equations) and density-pressure (Dowson-Higginson equation) relationships have also been taken into account in the overall modelling. Generic models have been derived for the aforementioned bearing components in order to enable their applications in general multibody dynamic systems (MDS), and by including the effects of angular misalignments, superficial geometric defects (form/waviness deviations, EHL deformations, etc.) and axial motion. The bearing exibility (conformal EHL) has been incorporated by means of FEM model reduction (or condensation) techniques. The macroscopic in fluence of the mixedlubrication phenomena have been included into the modelling by the stochastic Patir and Cheng average ow model and the Greenwood-Williamson/Greenwood-Tripp formulations for rough contacts. Furthermore, a deterministic mixed-lubrication model with inter-asperity cavitation has also been proposed for full-scale simulations in the microscopic (roughness) level. According to the extensive mathematical modelling background established, three significant contributions have been accomplished. Firstly, a general numerical solution for the Reynolds lubrication equation with the mass-conserving p - Ø cavitation model has been developed based on the hybridtype Element-Based Finite Volume Method (EbFVM). This new solution scheme allows solving lubrication problems with complex geometries to be discretized by unstructured grids. The numerical method was validated in agreement with several example cases from the literature, and further used in numerical experiments to explore its exibility in coping with irregular meshes for reducing the number of nodes required in the solution of textured sliding bearings. Secondly, novel robust partitioned techniques, namely: Fixed Point Gauss-Seidel Method (PGMF), Point Gauss-Seidel Method with Aitken Acceleration (PGMA) and Interface Quasi-Newton Method with Inverse Jacobian from Least-Squares approximation (IQN-ILS), commonly adopted for solving uid-structure interaction problems have been introduced in the context of tribological simulations, particularly for the coupled calculation of dynamic conformal EHL contacts. The performance of such partitioned methods was evaluated according to simulations of dynamically loaded connecting-rod big-end bearings of both heavy-duty and high-speed engines. Finally, the proposed deterministic mixed-lubrication modelling was applied to investigate the in fluence of the cylinder liner wear after a 100h dynamometer engine test on the hydrodynamic pressure generation and friction of Twin-Land Oil Control Rings.

Confiabilidade aplicada ao problema de interação estaca-solo.

Este trabalho busca aplicar técnicas de confiabilidade ao problema de grupo de estacas utilizadas como fundação de estruturas correntes. Para isso, lança-se mão de um modelo tridimensional de interação estaca-solo onde estão presentes o Método dos Elementos de Contorno (MEC) e o método dos Elementos Finitos (MEF) que atuam de forma acoplada. O MEC, com as soluções fundamentais de Mindlin (meio semi-infinito, homogêneo, isotrópico e elástico-linear é utiliza), é utilizado para modelar o solo. Já o MEF é utilizado para modelar as estacas. Definido o modelo de funcionamento estrutural das estacas, parte-se para a aplicação de métodos trazidos da confiabilidade estrutural para avaliação da adequabilidade em relação aos estados limite de serviço e estados limites últimos. Os métodos de confiabilidade utilizados foram o Método de Monte Carlo, o método FOSM (First-Order Second-Moment) e o método FORM (First-Order Reliability Method).