Método multigrid algébrico: reutilização das estruturas multigrid no transporte de contaminantes

A necessidade de obter solução de grandes sistemas lineares resultantes de processos de discretização de equações diferenciais parciais provenientes da modelagem de diferentes fenômenos físicos conduz à busca de técnicas numéricas escaláveis. Métodos multigrid são classificados como algoritmos escaláveis.Um estimador de erros deve estar associado à solução numérica do problema discreto de modo a propiciar a adequada avaliação da solução obtida pelo processo de aproximação. Nesse contexto, a presente tese caracteriza-se pela proposta de reutilização das estruturas matriciais hierárquicas de operadores de transferência e restrição dos métodos multigrid algébricos para acelerar o tempo de solução dos sistemas lineares associados à equação do transporte de contaminantes em meio poroso saturado. Adicionalmente, caracteriza-se pela implementação das estimativas residuais para os problemas que envolvem dados constantes ou não constantes, os regimes de pequena ou grande advecção e pela proposta de utilização das estimativas residuais associadas ao termo de fonte e à condição inicial para construir procedimentos adaptativos para os dados do problema. O desenvolvimento dos códigos do método de elementos finitos, do estimador residual e dos procedimentos adaptativos foram baseados no projeto FEniCS, utilizando a linguagem de programação PYTHONR e desenvolvidos na plataforma Eclipse. A implementação dos métodos multigrid algébricos com reutilização considera a biblioteca PyAMG. Baseado na reutilização das estruturas hierárquicas, os métodos multigrid com reutilização com parâmetro fixo e automática são propostos, e esses conceitos são estendidos para os métodos iterativos não-estacionários tais como GMRES e BICGSTAB. Os resultados numéricos mostraram que o estimador residual captura o comportamento do erro real da solução numérica, e fornece algoritmos adaptativos para os dados cuja malha retornada produz uma solução numérica similar à uma malha uniforme com mais elementos. Adicionalmente, os métodos com reutilização são mais rápidos que os métodos que não empregam o processo de reutilização de estruturas. Além disso, a eficiência dos métodos com reutilização também pode ser observada na solução do problema auxiliar, o qual é necessário para obtenção das estimativas residuais para o regime de grande advecção. Esses resultados englobam tanto os métodos multigrid algébricos do tipo SA quanto os métodos pré-condicionados por métodos multigrid algébrico SA, e envolvem o transporte de contaminantes em regime de pequena e grande advecção, malhas estruturadas e não estruturadas, problemas bidimensionais, problemas tridimensionais e domínios com diferentes escalas.

A nova família de distribuições odd log-logística: teoria e aplicações

Neste trabalho, foi proposta uma nova família de distribuições, a qual permite modelar dados de sobrevivência quando a função de risco tem formas unimodal e U (banheira). Ainda, foram consideradas as modificações das distribuições Weibull, Fréchet, half-normal generalizada, log-logística e lognormal. Tomando dados não-censurados e censurados, considerou-se os estimadores de máxima verossimilhança para o modelo proposto, a fim de verificar a flexibilidade da nova família. Além disso, um modelo de regressão locação-escala foi utilizado para verificar a influência de covariáveis nos tempos de sobrevida. Adicionalmente, conduziu-se uma análise de resíduos baseada nos resíduos deviance modificada. Estudos de simulação, utilizando-se de diferentes atribuições dos parâmetros, porcentagens de censura e tamanhos amostrais, foram conduzidos com o objetivo de verificar a distribuição empírica dos resíduos tipo martingale e deviance modificada. Para detectar observações influentes, foram utilizadas medidas de influência local, que são medidas de diagnóstico baseadas em pequenas perturbações nos dados ou no modelo proposto. Podem ocorrer situações em que a suposição de independência entre os tempos de falha e censura não seja válida. Assim, outro objetivo desse trabalho é considerar o mecanismo de censura informativa, baseado na verossimilhança marginal, considerando a distribuição log-odd log-logística Weibull na modelagem. Por fim, as metodologias descritas são aplicadas a conjuntos de dados reais.

Propriedades de continuação única para soluções de equações de Schrödinger com ponto de interação

Neste trabalho, estudamos propriedades de continuação única para as soluções da equação tipo Schrödinger com um ponto interação centrado em x=0, \\partial_tu=i(\\Delta_Z+V)u, onde V=V(x,t) é uma função de valor real e -\\Delta_Z é o operador escrito formalmente como \\[-\\Delta_Z=-\\frac\\frac{d^2}{dx^2}+Z\\delta_0,\\] sendo \\delta_0 a delta de Dirac centrada em zero e Z qualquer número real. Logo, usamos estes resultados para ver o possível fenômeno de concentração das soluções, que explodem, da equação de tipo Schrödinger não linear com um ponto de interação em x=0, \\[\\partial_tu=i(\\Delta_Zu+|u|^u),\\] com ho>5. Também, mostramos que para certas condições sobre o potencial dependente do tempo V, a equação linear em cima tem soluções não triviais.

Avaliação quantitativa dos efeitos da levodopa e da estimulação do núcleo subtalâmico sobre o equilíbrio em pacientes com doença de Parkinson

INTRODUÇÃO: Os efeitos da levodopa (LD) e da estimulação cerebral profunda (ECP) de núcleo subtalâmico (STN) sobre o equilíbrio e sintomas axiais são até o momento controversos. OBJETIVOS: Avaliar quantitativamente os efeitos da ECP de STN e da LD sobre o equilíbrio estático em pacientes com DP operados, em comparação com a LD em pacientes não operados. MÉTODOS: Trinta e um pacientes submetidos a ECP de STN entre 3 meses e 1 ano e meio antes da avaliação e 26 controles portadores de DP não operados, estágios Hoehn e Yahr 2 a 4 foram avaliados usando UPDRS para avaliação clínica e plataforma de força para avaliar oscilações posturais. O primeiro grupo foi avaliado com ECP e sem medicação, com ECP e com medicação e sem ECP e sem medicação. O segundo grupo foi avaliado com e sem medicação. Cada paciente foi avaliado com os olhos abertos e fechados. O deslocamento do centro de pressão anteroposterior, laterolateral, a área, velocidade e deslocamento total linear foram medidos pela plataforma de força. Os dados paramétricos foram comparados usando o teste t de Student e os dados não-paramétricos foram comparados pelo teste de Kruskal-Wallis. A avaliação clínica consistiu na parte 3 da escala UPDRS e na escala Hoehn e Yahr. Nível de significância estatística considerada foi p=0,05. RESULTADOS: Os pacientes não operados oscilaram mais quando sob efeito da levodopa do que sem medicação. No grupo operado, a maior oscilação é no grupo com ECP desligada e sem medicação. Tende a reduzir sob efeito da ECP apresenta redução significativa sob efeito simultâneo de ECP e levodopa. CONCLUSÃO: A associação da ECP de NST com medicação tem impacto positivo sobre o controle postural. O efeito da ECP de NST reverte o efeito negativo da levodopa sobre as oscilações observadas em pacientes não operados