Análise de discretizações e interpolações em malhas icosaédricas e aplicações em modelos de transporte semi-lagrangianos

• Autor(a) principal: Pedro da Silva Peixoto
• Data de Publicação: 2013
• Tipo de documento: Tese
• Idioma: por
• Título da fonte: Biblioteca Digital de Teses e Dissertações da USP
• Texto Completo: https://doi.org/10.11606/T.45.2013.tde-26062013-174032

A esfera é utilizada como domínio computacional na modelagem de diversos fenômenos físicos, como em previsão numérica do tempo. Sua discretização pode ser feita de diversas formas, sendo comum o uso de malha regulares em latitude/longitude. Recentemente, também para melhor uso de computação paralela, há uma tendência ao uso de malhas mais isotrópicas, dentre as quais a icosaédrica. Apesar de já existirem modelos atmosféricos que usam malhas icosaédricas, não há consenso sobre as metodologias mais adequadas a esse tipo de malha. Nos propusemos, portanto, a estudar em detalhe diversos fatores envolvidos no desenvolvimento de modelos atmosféricos globais usando malhas geodésicas icosaédricas. A discretização usual por volumes finitos para divergente de um campo vetorial utiliza como base o Teorema da Divergência e a regra do ponto médio nas arestas das células computacionais. A distribuição do erro obtida com esse método apresenta uma forte relação com características geométricas da malha. Definimos o conceito geométrico de alinhamento de células computacionais e desenvolvemos uma teoria que serve de base para explicar interferências de malha na discretização usual do divergente. Destacamos os impactos de certas relações de alinhamento das células na ordem da discretização do método. A teoria desenvolvida se aplica a qualquer malha geodésica e também pode ser usada para os operadores rotacional e laplaciano. Investigamos diversos métodos de interpolação na esfera adequados a malhas icosaédricas, e abordamos o problema de interpolação e reconstrução vetorial na esfera em malhas deslocadas. Usamos métodos alternativos de reconstrução vetorial aos usados na literatura, em particular, desenvolvemos um método híbrido de baixo custo e boa precisão. Por fim, utilizamos as técnicas de discretização, interpolação e reconstrução vetorial analisadas em um método semi-lagrangiano para o transporte na esfera em malhas geodésicas icosaédricas. Realizamos experimentos computacionais de transporte, incluindo testes de deformações na distribuição do campo transportado, que mostraram a adequação da metodologia para uso em modelos atmosféricos globais. A plataforma computacional desenvolvida nesta tese, incluindo geração de malhas, interpolações, reconstruções vetoriais e um modelo de transporte, fornece uma base para o futuro desenvolvimento de um modelo atmosférico global em malhas icosaédricas.