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. / The element helium, found mainly in natural gas reserves, condenses at temperature of 4.2K, and is the unique known substance that remains in liquid to absolute zero. In the liquid phase, the helium presents still another phase change in 2.19K, where passes of common liquid to superfluous liquid, with almost zero viscosity. These properties give the helium important applications. One of the major applications is as a coolant in superconductors, such as in the particle accelerator LHC, which is being built in the French border with Switzerland, in magnetic resonance devices, artificial satellites, etc.. In this paper, we present two mathematical models for heat transfer in liquid helium. The first model, considering only macroscopic movements, is derived based on constitutive laws of Fourier and Gorter-Mellink. The second model, based on techniques of Fremond, includes microscopic movements and can be seen as a regularization of the first model. Both models are governed by highly nonlinear differential equations resulting from the nonlinearity of the law of Gorter-Mellink and change of phase. Both models can be considered special cases of the Stefan problem in two phases, with phase one of the heat flux is governed by non-linear equation of the problem known as p-Laplacian, with p = 4/3. We also presented techniques to efficiently solve the problem of p-Laplacian, both for large values of p, p>> 2, and for values of p close to 1, which are major numerical challenges. Are proposed two simple iterative methods, one based on the method of quasi-Newton, with the relaxation term and the other by the Helmholtz decomposition, creating a system of equations whose matrices are constant, which reduces significantly the computational cost. Numerical experiments are conducted to test the efficiency of numerical models proposed and the algorithms developed for solving systems of nonlinear algebraic equations arising from approximations by finite elements. Are also presented results of studies of convergence, showing rates of optimal or near optimal convergence, comparable to that of interpolates. For the problem with phase change, due to the discontinuity of the gradient of temperature on the interface separating the two phases of liquid helium, the rate of convergence is not optimal. Using adaptive mesh, it is also great rates to the problem with change of phase. Using experimental data found in literature, for the parameters of thermal conductivity, density and specific heat, temperature dependent, are also presented for validation testing of the model and examples of possible applications. In tests for validating the model, compared to the numerical solution of the mathematical model with experimental results for the temperature found in literature.

Submitted by Erika Demachki (erikademachki@gmail.com) on 2015-01-16T17:13:32Z No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Tese - Marcos Roberto Marcial - 2014.pdf: 1577179 bytes, checksum: ac1649c996b2193bad6b704f05eca30c (MD5) / Approved for entry into archive by Erika Demachki (erikademachki@gmail.com) on 2015-01-16T17:40:15Z (GMT) No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Tese - Marcos Roberto Marcial - 2014.pdf: 1577179 bytes, checksum: ac1649c996b2193bad6b704f05eca30c (MD5) / Made available in DSpace on 2015-01-16T17:40:15Z (GMT). No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Tese - Marcos Roberto Marcial - 2014.pdf: 1577179 bytes, checksum: ac1649c996b2193bad6b704f05eca30c (MD5) Previous issue date: 2014-06-23 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / In this work, we study the topological structure of the solution set for a class of problems −Δpu = λ f (u)+μg(u)|∇u|p+Ψ(x) in Ω, u > 0 in Ω, u = 0 on ∂Ω, where Ω ⊂ IRN is a bounded domain with ∂Ω smooth, p, λ, μ are constants with p > 1, λ ≥ 0, μ ∈ IR and f ,g : (0,∞)→IR Ψ : Ω→IR are continuous functions. We will use Variational and Topological Methods, which includes minimization of energy functional and building connected components of solutions in a sense that we will define. Also we will employ arguments about the theory of regularity for p-Laplacian operator, approach arguments , maximum principles, results about sub and supersolutions and also arguments including monotonic type operators. / Neste trabalho estudamos a estrutura topológica do conjunto de soluções da classe de problemas −Δpu = λ f (u)+μg(u)|∇u|p+Ψ(x) em Ω, u > 0 em Ω, u = 0 sobre ∂Ω, onde Ω⊂IRN é um domínio limitado com fronteira ∂Ω regular, p, λ, μ são constantes com p > 1, λ ≥ 0, μ ∈ IR e f ,g : (0,∞)→IR, Ψ : Ω→IR são funções contínuas. Utilizamos Métodos Variacionais e Topológicos, que incluem minimização de funcionais energia e construção de componentes conexas de soluções em um sentido que definiremos. Empregamos também argumentos sobre a teoria da regularidade para o operador p- Laplaciano, argumentos de aproximação, bem como princípios de máximo, resultados sobre sub e supersoluções e também argumentos com operadores tipo monotônico.

[pt] Esta tese tem o intuito do estudo da existência, não existência e estimativas a priori de soluções não negativas de alguns tipos de problemas elípticos degenerados coercivos e não coercivos com um termo adicional dependendo do gradiente. Dentre outras coisas, obtemos condições integrais generalizadas tipo Keller-Osserman para a existência e não existência de soluções. Também mostramos que condições adicionais e diferentes são necessárias quando p é maior ou igual à 2 ou p é menor ou igual à 2, devido ao caráter degenerado do operador. As estimativas a priori são obtidas para super-soluções e soluções de EDPs elípticas superlineares o sistemas de tais tipos de equações em forma divergente com diferentes operadores e não linearidades. Além do mais, obtemos extensões até a fronteira de algumas desigualdades de Harnack fracas e lemas quantitativos de Hopf para operadores elípticos como o p-Laplaciano. / [en] This thesis concerns the study of existence, nonexistence and a priori estimates of nonnegative solutions of some types of degenerate coercive and non coercive elliptic problems involving an additional term which depends on the gradient. Among other things, we obtain generalized integral conditions of Keller-Osserman type for the existence and nonexistence of solutions. Also, we show that different conditions are needed when p is higher or equal to 2 or p is less than or equal to 2, due to the degeneracy of the operator. The uniform a priori estimates are obtained for supersolutions and solutions of superlinear elliptic PDE or systems of such PDE in divergence form that can contain different operators and nonlinearities. We also give full boundary extensions to some half Harnack inequalities and quantitative Hopf lemmas, for degenerate elliptic operators like the p-Laplacian.