IMPLEMENTAÇÃO DO MÉTODO DAS CARACTERÍSTICAS NA MODELAGEM DE PROBLEMAS DE CONVECÇÃO NATURAL EM CAVIDADES CILÍNDRICAS - PDF
Download "IMPLEMENTAÇÃO DO MÉTODO DAS CARACTERÍSTICAS NA MODELAGEM DE PROBLEMAS DE CONVECÇÃO NATURAL EM CAVIDADES CILÍNDRICAS"
Madalena Quintão Marinho
1 Campus de Bauru IMPLEMENTAÇÃO DO MÉTODO DAS CARACTERÍSTICAS NA MODELAGEM DE PROBLEMAS DE CONVECÇÃO NATURAL EM CAVIDADES CILÍNDRICAS Alice Jordam Bauru 2010
2 Campus de Bauru IMPLEMENTAÇÃO DO MÉTODO DAS CARACTERÍSTICAS NA MODELAGEM DE PROBLEMAS DE CONVECÇÃO NATURAL EM CAVIDADES CILÍNDRICAS Alice Jordam Orientador: Prof. Dr. Vicente Luiz Scalon Dissertação apresentada à Faculdade de Engenharia da UNESP Campus de Bauru, como parte dos requisitos para obtenção do título de Mestre em Engenharia Mecânica. Bauru 2010
3 Jordam, Alice. Implementação do método das características na modelagem de problemas de convecção natural em cavidades cilíndricas / Alice Jordam, f. il. Orientador: Vicente Luiz Scalon Dissertação (Mestrado) Universidade Estadual Paulista. Faculdade de Engenharia, Bauru, CBS. 2. Convecção natural. 3. Método dos elementos finitos. 4. Cavidade cilíndrica. I. Universidade Estadual Paulista. Faculdade de Engenharia. II. Título.
5 i AGRADECIMENTOS Agradeço primeiramente à Deus, pela vida, por me propiciar viver tamanha experiência e realizar um sonho. Agradeço imensamente ao meu orientador, Prof. Dr. Vicente Luiz Scalon, que de forma crucial me auxiliou em todos os aspectos deste trabalho. Agradeço aos meus familiares e amigos que me apoiaram e me incentivaram. Agradeço também à CAPES pelo apoio financeiro concedido em forma de bolsa de estudo. Agradeço a todos, professores e funcionários da Unesp, que contribuíram, de certa forma, para o desenvolvimento deste trabalho.
6 ii A gravidade explica os movimentos dos planetas, mas não pode explicar quem colocou os planetas em movimento. Deus governa todas as coisas e sabe tudo que é ou que pode ser feito. (Isaac Newton)
7 iii Jordam, A. Implementação do método das características na modelagem de problemas de convecção natural em cavidades cilíndricas. 85f. Dissertação (Mestrado em Engenharia Mecânica) Faculdade de Engenharia, Universidade Estadual Paulista, RESUMO A fluidodinâmica computacional (CFD) tem sido utilizada, estudada e implementada ao longo das duas últimas décadas na solução dos mais diversos problemas de engenharia. O princípio básico desta ciência é a aplicação de métodos numéricos em problemas que envolvam mecânica dos fluidos. Nesse contexto, este trabalho utiliza essa técnica para analisar o comportamento de um fluido incompressível, que se encontra numa cavidade cilíndrica fechada onde as faces inferior e superior são adiabáticas e as superfícies laterais se encontram em diferentes temperaturas. Os perfis de velocidade e temperatura resultantes - ocasionados pela convecção natural- serão avaliados em todo o domínio do problema. Existe uma série de técnicas para a solução de problemas envolvendo escoamentos, sendo as mais comuns as que se utilizam do Esquema de Passo Fracionado proposto por Chorin no final da década de 60. Dentre as diversas soluções que se utilizam desta técnica, este trabalho optou pelo uso do método das características e do algoritmo CBS de solução proposto por Zienkiewicz e Codina (1995). Para a implementação do algoritmo de solução do problema proposto foi realizada uma discretização geral através do método dos elementos finitos usando-se de uma malha formada por elementos bilineares. A solução foi obtida a partir de um ambiente matemático adequado, o GNU-Octave (2008). Os resultados foram analisados para diferentes razões de curvatura, números de Rayleigh e métodos de solução, sendo plotados para as suas diversas variáveis buscando descrever o comportamento do fenômeno. PALAVRAS-CHAVE: CBS, Convecção Natural, Método dos Elementos Finitos, Cavidade Cilíndrica.
8 iv Jordam, A. Implementation of the characteristics method in the modeling of free convection problems in cylindrical cavity. 85f. Dissertação (Mestrado em Engenharia Mecânica) Faculdade de Engenharia, Universidade Estadual Paulista, ABSTRACT The Computational Fluid Dynamics (CFD) has been used, studied and performed through the last two decades to solve the series of problems in Engineering. The most basic aim of this science is the appliance of numerical methods in cases that envolve fluid mechanics. In this context, this work uses this technic to analyse the behaviour of an incompressible fluid, which is found in a closed cylindrical cavity, a place where the inferior and superior surfaces are adiabatic and the lateral surfaces are shown in different temperatures. The resultant profiles of speed and temperature induced by the free convection are going to be appraised in all the dominion of the problem. There is a set of technics to solve the problems which involve the drainage, but the most usual are those which use the technic Fractional Step Method offered by Chorin in the final of 60s. Among the several solutions that are solved through this technic, this research used the characteristics method and of the CBS algorithm, offered by Zienkiewicz e Codina (1995). For the implementation of the algorithm, it was realized a general discretization through the finite elements method, making use of a loop formed by bilinear elements. The resolution was obtained from an adequated mathematical ambient, the GNU-Octave (2008). The results were analysed for different curvature ratios, Rayleigh numbers and methods of solution, being plotted for its different variables searching to describe the behaviour of the phenomenon. KEYWORDS: CBS, Free Convection, Finite Elements Method, Cylindrical Cavity.
9 v SUMÁRIO 1 INTRODUÇÃO Objetivos REVISÃO BIBLIOGRÁFICA METODOLOGIA O método das características aplicado às equações de transporte Estimativa da velocidade ao longo da linha de característica Aproximações das diversas grandezas no ponto característico Aplicação do método das características para sistemas radiais Determinação da forma geral para a equação de transporte Aplicação às equações de balanço de quantidade de movimento Cálculo das matrizes elementares Tratamento do termo fonte Q A equação da pressão A equação da correção das velocidades Compressibilidade artificial Incremento de tempo Método semi-implícito Ausência do termo característico A equação de conservação de energia Procedimento geral de solução RESULTADOS E DISCUSSÕES Validação do modelo Analisando um problema no plano cartesiano Análise dos perfis de velocidades Análise do problema para os esquemas explícito, semi-implícito, com e sem o uso dos termos característicos Resultados para cavidades cilíndricas anulares com diferentes razões de curvatura Resultados para cavidades cilíndricas com diferentes números de Rayleigh CONCLUSÕES E PROPOSTAS DE TRABALHOS FUTUROS REFERÊNCIAS BIBLIOGRÁFICAS...67
10 vi LISTA DE FIGURAS Figura 1: Problema proposto...4 Figura 2: Diagrama do algoritmo CBS de solução...9 Figura 3: Condições de contorno para o problema proposto na Figura Figura 4: Procedimento de determinação da linha pseudo-característica de Galerkin...14 Figura 5: Fluxograma geral do processo de solução...34 Figura 6: Malha utilizada na solução do problema...37 Figura 7: Geometria do problema proposto por Schneider e Straub (1992)...39 Figura 8: Resultado para o perfil de velocidades obtido por Schneider e Straub (1992)...39 Figura 9: Resultado para o perfil de velocidades obtido neste trabalho...40 Figura 10: Esquema mostrando as linhas onde estão localizados os resultados para a montagem dos gráficos comparativos...41 Figura 11: Comportamento da velocidade v em relação ao seu eixo central para Ra=10³...42 Figura 12: Comportamento da velocidade u em relação ao seu eixo central para Ra=10³...43 Figura 13: Comportamento da velocidade u em relação ao seu eixo central para Ra= Figura 14: Comportamento da velocidade v em relação ao seu eixo central para Ra= Figura 15: Perfil de temperatura para Ra=10³...46 Figura 16: Perfil de temperatura para Ra= Figura 17: Perfil de velocidade e função corrente para Ra=10³...48 Figura 18: Perfil de velocidade e função corrente para Ra= Figura 19: Desvio da pressão em função do número de iterações para Ra=10³...50 Figura 20: Desvio da pressão em função do número de iterações para Ra= Figura 21: Tempo computacional gasto em função do número de iterações para Ra= Figura 22: Tempo computacional gasto em função do número de iterações para Ra= Figura 23: Comportamento do número de Nusselt em função da parede aquecida para os diversos esquemas de solução para Ra= Figura 24: Comportamento do número de Nusselt em função da parede aquecida para os diversos esquemas de solução para Ra= Figura 25: Comportamento do Nu na superfície quente em função da altura para as diferentes razões de curvatura...55 Figura 26: Comportamento do perfil de temperaturas em função da razão de curvatura...56
11 vii Figura 27: Comportamento do perfil de velocidades e função corrente em função da razão de curvatura...57 Figura 28: Comportamento da velocidade u em relação ao seu eixo central para os números de Rayleigh considerados...59 Figura 29: Comportamento da velocidade v em relação ao seu eixo central para os números de Rayleigh considerados...60 Figura 30: Perfil de temperatura em função do número de Rayleigh...61 Figura 31: Perfil de velocidade e função corrente em função do número de Rayleigh...63 Figura 32: Comportamento do Nu na superfície quente em função da altura para os números de Rayleigh considerados...64 Figura 33: Desvio da pressão em função do número de iterações para os números de Rayleigh considerados...64
12 viii LISTA DE TABELAS Tabela 1: Comparação de alguns resultados, utilizando malha 21x21, com outros disponíveis na literatura...38
13 ix LISTA DE ABREVIATURAS E SIGLAS AC Artificial Compressibility Compressibilidade Artificial CBS Characteristic Based Split Algoritmo utilizando o método das características
14 x SIMBOLOGIA Latinos B matriz da derivada da função de aproximação elementar c velocidade do som adimensional C matriz auxiliar Fo número de Fourier [Fo= t / H 2 ] F s matriz relacionada às condições de contorno g aceleração da gravidade [m/s 2 ] G matriz auxiliar h menor dimensão do elemento H comprimento característico H matriz auxiliar H * altura da cavidade cilíndrica k constante difusiva genérica K u K k M n N matriz auxiliar matriz auxiliar matriz auxiliar vetor unitário na direção normal à superfície matriz da função de aproximação elementar Nu número de Nusselt [Nu=h H / k ] p pressão adimensional P pressão [ Pa ] P matriz auxiliar Pr número de Prandtl [Pr= / ] Q termo fonte genérico Ra número de Rayleigh [Ra=Pr g T q T f H 3 / 2 ] Re * R e * R i t número de Reynolds [Re= u H / ] raio externo raio interno tempo adimensional
15 xi T temperatura [K ] T f temperatura mínima de referência (fria) [K ] T q temperatura máxima de referência (quente) [K ] u componente adimensional da velocidade na direção radial U velocidade adimensional na direção radial v componente adimensional da velocidade na direção vertical V velocidade adimensional na direção vertical r coordenada adimensional na direção radial z coordenada adimensional na direção vertical Gregos difusividade térmica fator de compressibilidade artificial coeficiente de expansão volumétrica [1/ K ] comprimento da linha de corrente durante o instante de tempo T diferença de temperatura de referência variação da propriedade que o acompanha propriedade genérica representação da superfície de controle razão de curvatura viscosidade cinemática [m 2 /s ], 1, 2 temperatura adimensional fatores, entre 0 e 1, que determinam o esquema de discretização no tempo ρ massa específica [kg/ m 3 ] representação do volume de controle (domínio) Sobrescritos n instante atual n 1 instante posterior n instante intermediário - valor médio + valor dimensional * variável auxiliar ' mudança de uma variável independente para trabalhar coincidente com a linha
16 xii corrente Subscritos r, z i, j, k ref u p direção de atuação notação tensorial valor referencial associado ao campo de velocidades associado ao campo de pressões
17 1 1 INTRODUÇÃO As simulações computacionais são alvos de estudos há décadas e cada vez mais vêm chamando a atenção dos pesquisadores. Este grande interesse se deve ao fato de que, por ser um método de relativa simplicidade e ter grande versatilidade, este diminui significativamente os trabalhos de laboratórios. A evolução da velocidade e capacidade de armazenamento dos computadores nas últimas três décadas representou um grande estímulo para a aplicação de novas técnicas numéricas de solução. Além disto, contribui muito também o fato de que a diminuição dos ensaios experimentais reduz os custos dos testes de maneira substancial. Pode-se dizer, de uma maneira geral, que as simulações de soluções para problemas científicos, deve ser capaz, através de técnicas matemáticas e computacionais, imitar um processo ou operação do mundo real, em outras palavras, deve abranger desde a análise dos fenômenos, o desenvolvimento de modelos matemáticos para sua descrição até a elaboração de códigos computacionais para obtenção das soluções desejadas. É uma área em expansão, que tem uma vasta aplicação, como por exemplo, em desenvolvimento de produtos industriais, pesquisas científicas básicas e aplicadas, simulações e previsões temporais e espaciais de fenômenos, dentre outras. Existem vários métodos estudados na modelagem computacional, podendose destacar os que têm aplicações diretas na solução de problemas típicos de engenharias, tais como: métodos dos elementos finitos, métodos dos elementos de contorno, métodos dos volumes finitos, métodos das diferenças finitas, métodos integrais e variacionais e métodos auto adaptativos. Além destes, há ainda técnicas computacionais, com o propósito de aprimorar e, principalmente, aumentar a velocidade de solução dos modelos anteriormente descritos, como por exemplo, redes e malhas computacionais, computação vetorial e paralela
18 2 aplicada, entre outras. Por fim tem-se ainda, para auxiliar nas etapas de estabelecimento do domínio computacional e visualização das soluções, princípios e, muitas vezes, programas que permitem etapas de pré e pós-processamento gráfico e otimização, sistemas de orientação espacial e modelagem. A solução por métodos numéricos de um problema é uma aproximação da solução de equações diferenciais, tanto ordinárias como parciais. Pode-se descrever simplificadamente o método da seguinte forma, as equações diferenciais parciais são resolvidas primeiramente através da discretização da equação, trazendo para dentro de um subespaço dimensional finito (isto pode ser feito através do método de elementos finitos, método das diferenças finitas ou método dos volumes finitos). Ao fazer isto, muitas vezes é necessário a solução de sistemas lineares ou não lineares de equações, que serão descritos no parágrafo à seguir. A justificativa teórica desses métodos envolvem teoremas de análise funcional que reduzem o problema para solução de equações algébricas. No presente estudo esta análise foi realizada utilizando-se o método dos elementos finitos. Os sistemas lineares ou não lineares de equações, representam outra importante interface entre os métodos numéricos e a análise matemática. Existem duas alternativas de soluções: o método direto e o iterativo. No método direto, a solução ocorre em um número finito de passos e fornece uma resposta precisa se utilizados com precisão aritmética infinita. No método iterativo não existe um número especificado de passos a serem calculados, começa com uma estimativa inicial que posteriormente com sucessivas aproximações converge para a resposta do problema. Um critério de convergência é especificado para decidir o quanto a solução está suficientemente precisa. Em geral, nos procedimentos de solução numérica os métodos iterativos são mais comumente utilizados que os métodos diretos, mesmo porque, são os indicados em sistemas de equações não-lineares. Cabe ressaltar, entretanto, que a análise numérica é anterior em muitos séculos ao uso dos computadores. Há anos atrás, por exemplo, já era utilizada a interpolação linear. Vários grandes matemáticos do passado já estavam preocupados com a análise numérica que pode-se refletir hoje, como resultado, nomes de importantes algoritmos como o método de Newton, interpolação polinomial de Lagrange, eliminação Gaussiana, o método de Euler entre outros. Foram produzidos muitos livros com fórmulas e tabelas contendo dados, tais como, pontos de interpolação e coeficientes de funções que eram e ainda são utilizados atualmente.
19 3 1.1 Objetivos Conhecidas as principais ferramentas disponíveis e seguindo a tendência de opção pelos modelos numéricos, este trabalho objetiva implementar e aplicar um algoritmo numérico para a simulação dos efeitos da convecção natural num fluido incompressível, situado no interior de uma cavidade cilíndrica, onde duas de suas faces encontram-se a diferentes temperaturas e as demais são adiabáticas. A convecção natural consiste na transferência de calor dentro de um fluído através de movimentos do próprio fluido, ocorrendo como consequência de diferenças na densidade ocasionada pela distribuição de temperaturas associadas à força gravitacional local e forças centrífugas. Seu estudo é de considerável importância em problemas de engenharia, uma vez que esta configuração pode ser tomada como uma representação simplificada de muitas situações práticas em diversas áreas da engenharia (nuclear, solar, mecânica). Pode-se citar como exemplo, a geometria típica de sistemas coletor-armazenador e o processo de convecção natural de água no interior de tanques. Muitos algoritmos e programas podem ser utilizados para implementar e executar a análise numérica nos computadores. Neste trabalho optou-se por implementar as rotinas de teste num ambiente matemático onde as ferramentas de análise e solução se encontram pré-implementadas. Para tanto, será utilizado o programa GNU-Octave (2008), que é um programa livre, aberto e com distribuição gratuita e pode ser executado nos sistemas operacionais Linux e Windows. Para o pós processamento e elaboração dos gráficos, embora pudesse ser feito diretamente pelo programa GNU-Octave, será utilizada uma outra ferramenta que melhor se adequa às necessidades de geração de isocurvas. Esta ferramenta pode ser encontrada em Kelley e Galbraith, GRI (2008) e também é livre, aberto, com distribuição gratuita e pode ser executado nos sistema operacionais Linux e Windows. O problema a ser solucionado consiste numa cavidade cilíndrica fechada na qual há diferenças de temperaturas em suas paredes opostas, a superfície interna encontra-se quente e a externa fria, as paredes superior e inferior são adiabáticas, como pode ser observado na Figura 1. Pretende-se analisar o comportamento da convecção natural para diferentes razões de curvatura. Define-se a razão de curvatura como a diferença entre raio externo R e * e raio interno R i * sobre o raio interno = r * / R i *. A razão geométrica
20 4 r * / H *, representada pela diferença de raios sobre a altura da cavidade permanecerá constante e igual a 1. Para todos os casos será analisado também o comportamento do número de Nusselt (Nu), na superfície radial interna (superfície aquecida). Os perfis de velocidades e temperaturas resultantes dos diferentes métodos de solução também serão alvos de análise, bem como o comportamento do problema para diferentes números de Rayleigh. Parede Adiabática u=0 e v=0 n =0 z Parede Fria u=0 v=0 =0 R i * R e * Parede Quente u=0 v=0 =1 Figura 1: Problema proposto.
21 5 2 REVISÃO BIBLIOGRÁFICA Há uma vasta diversidade de trabalhos que abordam os métodos numéricos. Por serem versáteis, flexíveis e de uma significativa simplicidade vários autores utilizam-se de problemas de engenharias para testá-los. O problema clássico de cavidade recirculante para o teste do funcionamento de modelos numéricos foi solucionado por Ghia, Ghia e Shin (1982), utilizando uma formulação baseada nas equações para escoamentos bidimensionais e incompressíveis de Navier-Stokes. O trabalho apresentou, ainda, uma grande variedade de informações dos testes realizados, com os resultados avaliados para diversos números de Reynolds. Uma série de outros resultados para cavidade recirculante e com diferentes geometrias também podem ser vistos no trabalho de Erturk (2008). Um outro problema, que recentemente foi sugerido como um referencial para fluxos incompressíveis, foi solucionado por Nithiarasu e Liu (2005), onde utiliza-se o algoritmo CBS para resolver tanto fluxos estáveis quanto instáveis dentro de uma dupla cavidade não-retangular para diferentes números de Reynolds. Ambas as malhas, estruturadas e não estruturadas foram empregadas no presente estudo para certificar-se de que os resultados obtidos estariam próximos da realidade. A conclusão sobre os resultados observados é que, quando utiliza-se Reynolds variando de 50 a 1.000, o fluxo encontra-se num regime estacionário, para Reynolds entre e 6.000, o fluxo adquiri um estado de transitoriedade e afasta-se do estado citado acima. E para Reynolds igual ou maior que o fluxo se aproxima de um regime turbulento. Em relação ao método utilizado, Codina et al (2006) compararam alguns métodos numéricos para a solução das equações de Navier-Stokes. O primeiro é baseado no algoritmo CBS-compressibilidade artificial e os outros dois abordam particularmente o
22 6 método SGS (subgrid scale). O estudo pôs em destaque que estes mecanismos de estabilização são na verdade muito semelhantes. O CBS é melhor em alguns casos, enquanto o regime SGS tem desempenho melhor em outros. Em relação a utilização do computador, o tempo é sem dúvida menor ao utilizar-se o CBS. Nesta mesma linha, Massarotti et al. (2006) apresentaram em seu trabalho uma comparação dos métodos explícitos e semi-implícitos do algoritmo CBScompressibilidade artificial. O estudo mostrou que os dois métodos apresentaram bons resultados. Todas as soluções para o estado estacionário se mostraram praticamente idênticas, especialmente para problemas de convecção natural. Entretanto, a aplicação do método explícito em problemas desta natureza apresenta alguns problemas de convergência e forças de flutuações extremamente altas. Na solução de problemas transientes nota-se algumas diferenças com relação ao incremento de tempo utilizado em cada situação. Embora o algoritmo CBS já tenha sido utilizado com sucesso em problemas de escoamento turbulento e com elevados valores do número de Reynolds, modelos de turbulência ainda não haviam sido testados. Percebendo-se esta lacuna, Nithiarasu e Liu (2006) apresentaram um algoritmo utilizando compressibilidade artificial e o algoritmo CBS para fluxos contínuos e com instabilidades de turbulência que pode ser aplicado a fluidos incompressíveis. Testes realizados pelos autores demonstraram o sucesso da aplicação do modelo k, aumentando assim, a gama de uso de variações do algoritmo CBS. A questão conservativa também foi abordada no trabalho de Nithiarasu e Zienkiewicz (2006). Neste estudo, foi realizado uma análise do método explícito para fluidos incompressíveis utilizando o algoritmo CBS, empregando tanto a forma conservativa como a não-conservativa. A estabilização, convergência e aspectos de conservação do presente método foram discutidas. Um procedimento para eliminação dos erros de primeira ordem no tempo foi proposto. Os resultados foram comparados com outros métodos tendo apresentado uma boa concordância. Escoamentos sem viscosidade e compressíveis utilizando como base o algoritmo CBS também foram analisados por Thomas e Nithiarasu (2005). Os resultados obtidos com este trabalho demonstraram que o tamanho do elemento e a sua uniformização são aspectos importantes na avaliação de sua precisão. Através deste estudo, os autores demonstraram que algoritmos que produzem maior refinamento, na direção da linha aerodinâmica, melhoram a acuracidade da solução na região supersônica do escoamento.
23 7 Ampliando a aplicação do algoritmo CBS, Rojek, Oñate e Taylor (2006) fizeram uma adaptação do esquema, antes aplicado somente na mecânica dos fluidos, para o uso em problemas que envolvam mecânica dos sólidos. Vários exemplos foram resolvidos, obtendo resultados bastante semelhantes aos utilizados com outros métodos de estabilização. Confirmando a vasta aplicação do algoritmo CBS, Kulkarni et al. (2006), utilizaram um modelo numérico bidimensional para simular o comportamento do fluido nos moldes de enchimento no processo de encapsulamento de um circuito integrado. As equações que regem o problema foram resolvidas usando o método das características baseado no algoritmo CBS em regime semi-implícito, juntamente com o método de elementos finitos. A simulação foi realizada para uma determinada geometria de um flip-chip. Os resultados obtidos apresentaram boa concordância com os numéricos e experimentais disponíveis, tal situação se enquadra para baixos números de Reynolds, que é o caso do processo de encapsulamento, devido ao alto índice de viscosidade de resina. Assim, foi possível demonstrar que a aplicação do presente modelo numérico mostrou-se uma alternativa eficaz. Fazendo uma compilação dos trabalhos até então publicados, Zienkiewicz e Taylor (2000) apresentaram um livro com a descrição completa da metodologia CBS. Utilizando-se da discretização de Galerkin e do método de elementos finitos, o livro mostra os passos da resolução, as formulações utilizadas e as avaliações do método, entre outras informações. Essa metodologia, a partir deste seu maior detalhamento, serviu como base de desenvolvimento para diversos autores no mundo inteiro, que passaram a trabalhar e divulgar este novo modelo. Nesse mesmo período, destaca-se ainda, outro trabalho de Nithiarasu e Zienkiewicz (2000) onde são apresentadas as técnicas de estabilização com os passos de tempo. Este trabalho tinha o seu principal foco nos escoamentos incompressíveis e na aceleração convectiva, utilizando o algoritmo CBS. Os modelos foram novamente testados para o caso da cavidade recirculante e, assim, foi possível uma avaliação dos resultados obtidos com a nova técnica. Mudando alguns aspectos da discretização usual apresentada por Zienkiewicz e Taylor (2000) com a utilização da forma fraca da equação em alguns outros pontos, Nithiarasu (2003) apresenta uma proposta variante do algoritmo CBS tradicional. As alterações propostas, de acordo com os resultados apresentados, mostraram que a qualidade do resultados não foi afetada e que a facilidade de aplicação teve um ganho. A eliminação da
24 8 necessidade de determinadas matrizes e/ou procedimentos trouxe um ganho para a sua aplicação, principalmente em problemas dinâmicos mais complexos onde o tempo de processamento torna-se mais significativo. Em relação a outros métodos de solução, Thomas e Nithiarasu (2007) comparam dois métodos: o método LCG (locally conservative Galerkin) e o método GG (Global Galerkin). Ao final do trabalho, pode-se concluir que ambos os métodos são praticamente idênticos, exceto pelo fato do LCG ter taxas de convergência maiores. Ampliando o método LCG, Thomas, Nithiarasu e Bevan (2007) resolveram as equações para fluido incompressível de Navier-Stokes. Além da implementação do método, dois exemplos de referência foram apresentados para demonstrar a veracidade do método. Voltando ao algoritmo CBS, as suas diversas aplicações e pequenas atualizações levaram à publicação de um novo trabalho escrito por Nithiarasu, Codina e Zienkiewicz (2006) onde eles apresentam uma visão abrangente de várias técnicas da dinâmica dos fluidos aplicadas no algoritmo, ou seja, neste trabalho encontram-se suas principais dificuldades de implementação, seus aspectos mais importantes, enfim todas as considerações necessárias para a aplicação do sistema CBS em resolução de problemas são exaustivamente discutidas. Uma grande variedade de trabalhos poderia ainda ser citada, abrangendo a utilização do algoritmo CBS e AC (Artifical Compressibility) em diversas situações. Entretanto, foi dada uma maior ênfase aos trabalhos relacionados ao desenvolvimento da sua metodologia e à soluções de escoamento de fluido incompressível utilizando o algoritmo CBS associado à compressibilidade artificial. Cabe destacar entre as referências apresentadas, os trabalhos de Zienkiewicz e Taylor (2000) e Nithiarasu (2003) que propõem as bases dos modelos utilizados e serão citados de forma recorrente durante o texto. Para melhor entender o algoritmo CBS, a Figura 2 ilustra um resumo do método, baseado no trabalho de Zienkiewicz e Taylor (2000). Pode-se observar que a solução por tal algoritmo consiste em uma solução fracionada, sendo que: primeiramente calcula-se as velocidades sem inserir a pressão, no segundo passo resolve-se a equação da pressão com os valores das velocidades calculados anteriormente, após corrige-se as velocidades. Com estas novas velocidades repete-se o processo até que o mesmo seja considerado convergente para a solução do problema.
25 9 INÍCIO Rotinas Preliminares Verificação dos dados Passos Passo 1: Velocidades Erro Fazer mudanças Passo 2: Pressão Passo 3: Correção das velocidades Não Energia Condições de Contorno Sim Calcula a temperatura Regime estacionário Não Dados de saída FIM Figura 2: Diagrama do algoritmo CBS de solução.
26 10 Cabe destacar ainda, que a grande maioria dos trabalhos citados utiliza-se da solução através de elementos finitos por elementos bidimensionais triangulares. Para o trabalho proposto a solução será avaliada com outro tipo de elemento: os elementos bidimensionais quadrangulares ou bilineares. Já em relação à convecção natural em cavidades cilíndricas, salienta-se que poucos artigos foram encontrados, e dentre esses, nenhum conincide com o problema proposto neste trabalho, no entanto pode-se destacar os trabalhos de Kurian, Varma e Kannan (2009), onde estuda-se os efeitos da convecção natural em um cilindro vertical com diversas inclinações. As superfícies inferior e superior estão a diferentes temperaturas e as demais são adiabáticas. Os resultados foram analisados para Prandtl fixado em 0,71, Rayleigh variando de 10 3 a 3,1 x10 4 e ângulos de inclinação variando de 0 a 180. Relações empíricas entre a velocidade máxima e número de Nusselt em função do ângulo de inclinação e número de Rayleigh foram analisadas neste trabalho. Preocupando-se com os efeitos da inclinação do cilindro sobre o comportamento do fluido, Soong et al. (1996) relataram um estudo numérico de convecção natural numa cavidade cilíndrica aquecida e inclinada. Os efeitos foram estudados para diversos números de Rayleigh e razões de curvatura. Seguindo esta mesma linha, mas não somente com a intenção de examinar as razões de curvatura, Ivancic et al. (1999) analisaram também as condições de contorno num cilindro vertical aquecido por baixo, com parede lateral sob duas condições: adiabática e condutora. Os resultados foram analisados para diferentes números de Rayleigh e Prandtl. Um outro trabalho que se destaca é o de Schneider e Straub (1992), onde estuda-se o comportamento de fluidos sob os aspectos da convecção natural em um cilindro com diversas inclinações. Foram realizadas diversas análises sobre o ângulo de inclinação e velocidade máxima resultante em função do número de Nusselt para diferentes números de Rayleigh e Prandtl. Conforme mencionado na introdução, o estudo realizado no presente trabalho, se propõe a avaliar, baseado no algoritmo CBS e na metodologia de elementos finitos, o comportamento de um fluido incompressível que se encontra numa cavidade cilíndrica, onde duas de sua faces estão a diferentes temperaturas e as demais são adiabáticas.