Modelagem matemática e métodos numéricos para simulação da condução do calor no hélio líquido

• Ano de defesa: 2008
• Autor(a) principal: Senger, Erasmo
• Tipo de documento: Tese
• Tipo de acesso: Acesso aberto
• Idioma: por
• Instituição de defesa: Laboratório Nacional de Computação Científica; Serviço de Análise e Apoio a Formação de Recursos Humanos; BR; LNCC; Programa de Pós-Graduação em Modelagem Computacional
• Palavras-chave em Português: Métodos dos elementos finitos; Hélio líquido; Teorias não-lineares 4. p-laplaciano; Simulação numérica; Finite element method; Liquid helium; Phase change Nonlinear theories; p-laplacian; Numerical simulation; CNPQ::CIENCIAS EXATAS E DA TERRA::CIENCIA DA COMPUTACAO::TEORIA DA COMPUTACAO::COMPUTABILIDADE E MODELOS DE COMPUTACAO
• Link de acesso: https://tede.lncc.br/handle/tede/108

O elemento hélio, encontrado principalmente em reservas de gás natural, entra em condensação à temperatura de 4,2K, e é a única substância conhecida que permanece no estado líquido até o zero absoluto. Na fase liquida, o hélio apresenta ainda, em K, outra mudança de fase, onde passa de líquido comum à superfluido, com viscosidade praticamente nula. Estas propriedades conferem ao hélio importantes aplicações. Uma hdas principais aplicações é como agente refrigerante em supercondutores, como por exemplo, no acelerador de partículas LHC, que está sendo construído na fronteira da França com a Suíça, em aparelhos de ressonância magnética, satélites artificiais, etc. Neste trabalho, são apresentados dois modelos matemáticos para a transferência de calor no hélio líquido. O primeiro modelo, considerando apenas movimentos macroscópicos, é derivado com base nas leis constitutivas de Fourier e de Gorter-Mellink. O segundo modelo, baseado nas técnicas de Fremond, inclui movimentos microscópicos e pode ser visto como uma regularização do primeiro modelo. Os dois modelos são governados por equações diferenciais fortemente não lineares resultantes da não linearidade da lei de Gorter-Mellink e da mudança de fase. Ambos os modelos podem ser considerados casos particulares do problema de Stefan de duas fases, sendo que em uma das fases o fluxo de calor é governado pela equação não-linear do problema conhecido como p-laplaciano, com p=4/3. São também apresentadas técnicas para resolver de forma eficiente o problema do p-laplaciano, tanto para valores grandes de p, p>>2, quanto para valores de p próximos à 1, que constituem importantes desafios numéricos. Para tanto são propostos dois métodos iterativos simples, um baseado no método de quase-Newton, com termo de relaxação e, outro através da decomposição de Helmholtz, gerando um sistema de equações cujas matrizes são constantes, o que diminui significativamente o custo computacional. Experimentos numéricos são realizados para testar a eficiência dos modelos numéricos propostos bem como dos algoritmos desenvolvidos para resolver os sistemas de equações algébricas não lineares resultantes das aproximações por elementos finitos. São apresentados resultados de estudos de convergência, mostrando taxas de convergência ótimas ou quase ótimas, comparáveis às das interpolantes. Para o problema com mudança de fase, devido à descontinuidade do gradiente da temperatura sobre a interface que separa as duas fases do hélio líquido, as taxas de convergência não são ótimas. Usando malhas adaptativas, consegue-se taxas ótimas também para o problema com mudança de fase. Usando dados experimentais, encontrados na literatura, para os parâmetros de condutividade térmica, densidade e calor específico, dependentes da temperatura, são também apresentados testes de validação do modelo e exemplos de possíveis aplicações. Nos testes de validação do modelo, compara-se a solução numérica do modelo matemático com resultados experimentais para a temperatura, encontrados na literatura.