Modelagem Automática de Escoamentos em Meios Porosos via Método dos Elementos finitos

• Autor(a) principal: Luna, Bruno Gustavo Borges
• Data de Publicação: 2012
• Tipo de documento: Dissertação
• Idioma: por
• Título da fonte: Repositório Institucional da UFPE
• Texto Completo: https://repositorio.ufpe.br/handle/123456789/11459
• topic: Escoamentos em meios porosos; Métodos dos elementos finitos; Modelagem automática

A simulação de escoamentos multifásicos em meios porosos impõe vários desafios de ordem numérica devido a uma série de fatores, como os meios altamente anisotrópicos e heterogêneos tratados neste tipo de análise, as Equações Diferencias Parciais (EDP) de natureza acoplada elíptica-hiperbólica que descrevem o fenômeno, entre outros. Uma vez definidas as formulações matemáticas e numéricas a serem utilizadas para modelar adequadamente o escoamento, encontra-se outra dificuldade na codificação destes métodos, já que, usualmente, despende-se um tempo considerável para o desenvolvimento de programas de computador que implementem formulações para casos gerais ou complexos. Este trabalho apresenta a implementação de um software criado utilizando a linguagem Python e a ferramenta computacional FEniCS para a geração automática de código de baixo-nível em C++ aplicado na solução numérica de escoamentos mono- e bifásicos em meios porosos usando o Método dos Elementos Finitos (MEF). Foram testados o MEF de Galerkin e o Método dos Elementos Finitos Mistos (MEFM) para a solução da equação da pressão (pressão e velocidade no caso do MEFM) e o MEF com estabilização via Streamline Upwind Petrov Galerkin (SUPG) e operador de captura de choque para a equação da saturação. Para a solução do sistema de equações lineares provenientes do MEF de Galerkin foram utilizadas técnicas de aceleração de convergência via Método Multigrid Algébrico (AMG). Os métodos descritos neste trabalho são gerais o suficiente para lidar com problemas tridimensionais, heterogêneos e anisotrópicos. Exemplos são apresentados e resultados discutidos para problemas uni- e bidimensionais com domínios homogêneos e heterogêneos com tensores de permeabilidade iso- e anisotrópicos. A comparação favorável dos resultados obtidos com soluções analíticas e referências da literatura demonstra o potencial da ferramenta desenvolvida para a simulação de escoamentos em meios porosos.