Custo-efetividade da Terapia Tripla para pacientes adultos portadores do genótipo 1 da hepatite C crônica.

As ações de prevenção, diagnóstico e tratamento da hepatite C crônica integram as agendas das políticas de saúde do Brasil e do mundo, pois se trata de uma doença com grande número de acometidos, com alto custo tratamento e que ocasiona graves desfechos e incapacidade, o que acaba por onerar seu custo social. Os protocolos clínicos e diretrizes terapêuticas demonstram os esforços de inúmeras entidades no combate da hepatite C, pois informam aos profissionais de saúde, pacientes e familiares e cidadãos em geral, qual seria a melhor forma, comprovada cientificamente, de se proceder frente a uma infecção desta natureza. Realizouse uma análise de custoefetividade, sob a perspectiva do SUS, das estratégias: tratamento e retratamento com a terapia dupla, tratamento com a terapia dupla e retratamento com a terapia tripla e tratamento com a terapia tripla. Através de modelo de simulação baseado em cadeias Markov foi criada uma coorte hipotética de 1000 indivíduos adultos, acima de 40 anos, de ambos os sexos, sem distinção declasse socioeconômica, com diagnóstico confirmado para hepatite C crônica, monoinfectados pelo genótipo 1 do VHC e com ausência de comorbidades. A simulação foi iniciada com todos os indivíduos portando a forma mais branda da doença, tida como a classificação histológica F0 ou F1 segundo a escala Metavir. Os resultados demonstram que as duas opções, ou seja, a terapia dupla/tripla e a terapia tripla estão abaixo do limiar de aceitabilidade para incorporação de tecnologia proposto pela OMS (2012) que é de 72.195 (R$/QALY) (IBGE, 2013; WHO, 2012). Ambas são custoefetivas, visto que o ICER da terapia dupla/tripla em relação alinha de base foi de 7.186,3 (R$/QALY) e o da terapia tripla foi de 59.053,8 (R$/QALY). Entretanto o custo incremental de terapia tripla em relação à dupla/tripla foi de 31.029 e a efetividade incremental foi de 0,52. Em geral, quando as intervenções analisadas encontramse abaixo do limiar, sugerese a adoção do esquema de maior efetividade. A terapia tripla, apesar de ter apresentado uma efetividade um pouco acima da terapia dupla/tripla, apresentou custo muito superior. Assim, como seria coerente a adoção de uma ou da outra para utilização no SUS, visto que este sistema apresenta recursos limitados, indicase a realização de um estudo de impacto orçamentário para obterse mais um dado de embasamento da decisão e assim poder apoiar o protocolo brasileiro existente ou sugerir a confecção de novo documento.

Avaliação de modelos de permeabilidade em meios porosos não consolidados

As simulações computacionais tem sido amplamente empregadas no estudo do escoamento darciano e não-darciano em meios porosos consolidados e não-consolidados. Neste trabalho, através de uma nova formulação para a equação de Forchheimer, foram identificadas duas novas propriedades denominados fator de comportamento do fluido, que atua paralelamente a permeabilidade, e permeabilidade equivalente global, resultado da relação anterior. Este comportamento foi estudado e validado através da implementação de um aparato experimental e um código computacional baseado no modelo de regressão-linear que, além disso, demonstrou que o escoamento, ainda que em regime não darciano, comporta-se linearmente como a equação de Darcy, ainda que o coeficiente angular desta diminuiu de acordo com a faixa do número de Reynolds atingida, sendo esta dependente do tipo de leito empregado. Ainda neste trabalho, foi implementado o método de otimização R2W para estimar os parâmetros da equação de Kozeny-Carman a partir de dados experimentais obtidos por Dias et al, a fim de simular o escoamento darciano em meios porosos. Por fim, foi alcançada excelente concordância entre os dados simulados pelo método R2W / equação de Kozeny-Carman e os dados reais.

Experimentos numéricos para deposição de parafinas com modelo Multisólido

Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo.

Desenvolvimento de ferramenta computacional para metrologia com microtomografia computadorizada

Esta tese apresentada uma proposta de desenvolvimento de uma ferramenta computacional para metrologia com microtomografia computadorizada que possa ser implantada em sistemas de microtomógrafos convencionais. O estudo concentra-se nas diferentes técnicas de detecção de borda utilizadas em processamento de imagens digitais.Para compreender a viabilidade do desenvolvimento da ferramenta optou-se por utilizar o Matlab 2010a. A ferramenta computacional proposta é capaz de medir objetos circulares e retangulares. As medidas podem ser horizontais ou circulares, podendo ser realizada várias medidas de uma mesma imagem, uma medida de várias imagens ou várias medidas de várias imagens. As técnicas processamento de imagens digitais implementadas são a limiarização global com escolha do threshold manualmente baseado no histograma da imagem ou automaticamente pelo método de Otsu, os filtros de passa-alta no domínio do espaço Sobel, Prewitt, Roberts, LoG e Canny e medida entre os picos mais externos da 1 e 2 derivada da imagem. Os resultados foram validados através de comparação com os resultados de teste realizados pelo Laboratório de Ensaios Mecânicos e Metrologia (LEMec) do Intstituto Politécnico do Rio de Janeiro (IPRJ), Universidade do Estado do Rio de Janeiro (UERJ), Nova Friburdo- RJ e pelo Serviço Nacional da Indústria Nova Friburgo (SENAI/NF). Os resultados obtidos pela ferramenta computacional foram equivalentes aos obtidos com os instrumentos de medição utilizados, demonstrando à viabilidade de utilização da ferramenta computacional a metrologia.

Proposta de um otimizador para rastreamento do ponto de máxima potência de um conjunto de painéis solares conectados em um sistema híbrido de energia.

Este trabalho apresenta o projeto e os algoritmos de controle, de um sistema de geração de energia híbrido. Este sistema é formado por conversores de potência conectados em Back-to-Back associados a um arranjo solar fotovoltaico, que por sua vez é conectado no lado CC dos conversores. Em relação ao sistema de geração fotovoltaico, a contribuição consiste no desenvolvimento de cinco algoritmos para determinar o ponto de máxima potência (MPP) do arranjo fotovoltaico. O primeiro algoritmo consiste em uma versão modificada do algoritmo de Perturbar e Observar (PeO); o segundo algoritmo proposto é baseado no método do gradiente (MG); e o terceiro é baseado na otimização do MG (MGO). Porém, são desenvolvidos algoritmos híbridos que combinam rede neural com o método PeO, e rede neural com o algoritmo MGO. O sistema foi desenvolvido e simulado utilizando o Matlab/Simulink, e os resultados de simulação são apresentados com objetivo da avaliar o comportamento do sistema e a resposta dos diferentes algoritmos. Esta resposta foi avaliada para condições transitórias e de regime permanente, considerando diferentes requisitos de consumo na carga, irradiância e temperatura.

Metodologia de análise do abastecimento de veículos elétricos no contexto das redes inteligentes de energia.

Este trabalho apresenta um modelo de simulação para avaliação da recarga de veículos elétricos em redes de distribuição de energia, com foco na adequação da capacidade dos transformadores. Com base em medições efetuadas, foi possível o tratamento dos dados para a modelagem de curvas de recarga de veículos elétricos e a análise do acréscimo da demanda da recarga em transformadores de distribuição. O modelo de simulação permite avaliar os efeitos das variáveis relacionadas à distância percorrida e ao horário da conexão dos veículos elétricos à rede, tanto de forma determinística quanto de forma aleatória. Estes efeitos foram comparados e permitiram concluir a importância da utilização deste método mais realista para calcular a capacidade dos transformadores, visando suprir esta nova demanda. Além disso, o trabalho apresenta indicadores dos benefícios ambientais e energéticos da penetração dos veículos elétricos no que diz respeito ao gasto com reabastecimento e às emissões evitadas, respectivamente, ao se utilizar veículos elétricos à bateria em substituição aos tradicionais veículos à combustão interna.

Avaliação vestibular em indivíduos com cinetose

A cinetose tem elevada prevalência mundial, sendo mais frequente na infância e no sexo feminino, mas também pode acometer adultos. Resulta de um conflito vestíbulo-visual que incide durante a locomoção em diversos meios de transporte, como carro, ônibus, avião e barco. A forma mais dramática ocorre no transporte marítimo. Ocasiona náusea, vômito, sudorese, aumento da salivação, redução do apetite, hipotensão e mal-estar. Geralmente é produzida por estímulo vestibular, mas também pode ser induzida por estímulo visual. Tanto as acelerações lineares quanto angulares geram cinetose se persistirem por longo período em indivíduos susceptíveis. Recentemente ocorre maior interesse nessa afecção devido ao crescente uso de tecnologia de simulação de vôo e direção automobilística. O objetivo do estudo foi analisar as alterações vestibulares nos indivíduos adultos com cinetose. Trata-se de um estudo prospectivo, tipo série de casos. Os pacientes do ambulatório do Hospital Universitário Pedro Ernesto-UERJ foram avaliados através de anamnese geral e dirigida e exame fisico geral e otorrinolaringológico. Aqueles com histórico de doença otológica foram excluídos. Posteriormente realizaram audiometria e testes vestibulares com o registro gráfico dos nistagmos através da vectoeletronistagmografia. Nos resultados das provas calóricas encontramos 3,33% de pacientes com alteração de predomínio direcional, 6,67% com alteração de predomínio labiríntico, 3,33% com hiperreflexia esquerda, 3,33% com hiporreflexia direita e 3,33% com hiporreflexia esquerda. Encontrados algumas variações na análise dos movimentos sacádicos, nistagmo optocinético e rastreio pendular. Os resultados foram de encontro aos achados da literatura. Diante dos achados, foi observado que o exame otoneurológico com registro gráfico da cinetose mostrou-se muito importante para a avaliação dos pacientes com essa afecção. O estudo traz benefícios por contribuir para o melhor entendimento da cinetose e estimular novas pesquisas na área.

Mineração de dados para previsão de eventos extremos de precipitação

No presente trabalho foram desenvolvidos modelos de classificação aplicados à mineração de dados climáticos para a previsão de eventos extremos de precipitação com uma hora de antecedência. Mais especificamente, foram utilizados dados observacionais registrados pela estação meteorológica de superfície localizada no Instituto Politécnico da Universidade do Estado do Rio de Janeiro em Nova Friburgo RJ, durante o período de 2008 a 2012. A partir desses dados foi aplicado o processo de Descoberta de Conhecimento em Banco de Dados (KDD Knowledge Discovery in Databases), composto das etapas de preparação, mineração e pós processamento dos dados. Com base no uso de algoritmos de Redes Neurais Artificiais e Árvores de Decisão para a extração de padrões que indicassem um acúmulo de precipitação maior que 10 mm na hora posterior à medição das variáveis climáticas, pôde-se notar que a utilização da observação meteorológica de micro escala para previsões de curto prazo é suscetível a altas taxas de alarmes falsos (falsos positivos). Para contornar este problema, foram utilizados dados históricos de previsões realizadas pelo Modelo Eta com resolução de 15 km, disponibilizados pelo Centro de Previsão de Tempo e Estudos Climáticos do Instituto Nacional de Pesquisas Espaciais CPTEC/INPE. De posse desses dados, foi possível calcular os índices de instabilidade relacionados à formação de situação convectiva severa na região de Nova Friburgo e então armazená-los de maneira estruturada em um banco de dados, realizando a união entre os registros de micro e meso escala. Os resultados demonstraram que a união entre as bases de dados foi de extrema importância para a redução dos índices de falsos positivos, sendo essa uma importante contribuição aos estudos meteorológicos realizados em estações meteorológicas de superfície. Por fim, o modelo com maior precisão foi utilizado para o desenvolvimento de um sistema de alertas em tempo real, que verifica, para a região estudada, a possibilidade de chuva maior que 10 mm na próxima hora.

Estudo de fluxos numéricos para modelos magneto-hidrodinâmicos utilizando o código Flash

A teoria magneto-hidrodinâmicos permite a estruturação de modelos computacionais, designados modelos MHDs, que são uma extensão da dinâmica dos fluidos para lidar com fluidos eletricamente carregados, tais como os plasmas, em que se precisa considerar os efeitos de forças eletromagnéticas. Tais modelos são especialmente úteis quando o movimento exato de uma partícula não é de interesse, sendo que as equações descrevem as evoluções de quantidades macroscópicas. Várias formas de modelos MHD têm sido amplamente utilizadas na Física Espacial para descrever muitos tipos diferentes de fenômenos de plasma, tais como reconexão magnética e interações de ventos estelares com diferentes objetos celestiais. Neste trabalho, o objetivo é analisar o comportamento de diversos fluxos numéricos em uma discretização de volumes finitos de um modelo numérico de MHD usando um esquema de malha entrelaçada sem separação direcional considerando alguns casos testes. Para as simulações, utiliza-se o código Flash, desenvolvido pela Universidade de Chicago, por ser um código de amplo interesse nas simulações astrofísicas e de fenômenos no espaço próximo à Terra. A metodologia consiste na inclusão de um fluxo numérico, permitindo melhoria com respeito ao esquema HLL.

Um estudo de caso envolvendo a retroanálise de tensões residuais de cravação.

As tensões residuais de cravação ocorrem quando a estaca e o solo não são totalmente descomprimidos após dissipação completa da energia transferida à estaca durante sua instalação. Este assunto tem merecido a atenção de vários pesquisadores, uma vez que a presença das tensões residuais na estaca causa uma alteração no seu comportamento quando carregada. A presente dissertação reuniu estudos anteriores que estabeleceram os principais fatores que influenciam as tensões residuais de cravação. Procurou-se executar uma série de simulações objetivando melhor entender o desenvolvimento das tensões residuais com a porcentagem de resistência de ponta, em continuidade aos estudos de Costa (1994). Costa (1994) observou que a porcentagem de carga na ponta é um fator relevante na avaliação das tensões residuais de cravação. Uma análise paramétrica efetuada confirmou estudos anteriores, verificando que a razão entre a carga residual na ponta em relação à capacidade de carga global cresce, a medida que a porcentagem de carga na ponta aumenta, chegando a um valor máximo para, em seguida, diminuir. Este comportamento é similar ao da curva de compactação do solo, quando, ao aumentar a umidade, o peso específico seco aumenta, até um valor máximo, correspondente à umidade ótima, para em seguida reduzir. Verificouse, ainda, semelhantemente ao aumento da energia de compactação, que vai transladando a curva para cima e para à esquerda do gráfico umidade x peso específico seco, que o aumento do comprimento da estaca apresenta um comportamento similar. O aumento do comprimento leva a curva para cima e para a esquerda. Finalmente, selecionou-se um caso de obra com condições propícias para o desenvolvimento de altas tensões residuais de cravação. Trata-se de um caso de estacas metálicas longas, embutidas em solo residual jovem, de elevada resistência, bem documentado com extensa instrumentação, por ocasião da instalação. A retroanálise de cinco estacas do banco de dados mostrou que as cargas residuais previstas, num programa de simulação de cravação, se aproximaram muito dos valores experimentais. Foi observado também que a profundidade do ponto neutro previsto e medido apresentou uma excelente concordância. O resultado mais relevante desta pesquisa foi quando os valores de porcentagem de ponta foram introduzidos no eixo das abscissas e os valores da razão entre a carga residual na ponta e a capacidade de carga global, no eixo das ordenadas, e se observou o aspecto da curva semelhante à de compactação. A dissertação ilustra ser possível e simples a previsão das tensões residuais de cravação através de uma análise pela equação da onda.

Análise dos problemas direto e inverso em um modelo de evolução populacional através de transformações integrais e inferência bayesiana

Modelos de evolução populacional são há muito tempo assunto de grande relevância, principalmente quando a população de estudo é composta por vetores de doenças. Tal importância se deve ao fato de existirem milhares de doenças que são propagadas por espécies específicas e conhecer como tais populações se comportam é vital quando pretende-se criar políticas públicas para controlar a sua proliferação. Este trabalho descreve um problema de evolução populacional difusivo com armadilhas locais e tempo de reprodução atrasado, o problema direto descreve a densidade de uma população uma vez conhecidos os parâmetros do modelo onde sua solução é obtida por meio da técnica de transformada integral generalizada, uma técnica numérico-analítica. Porém a solução do problema direto, por si só, não permite a simulação computacional de uma população em uma aplicação prática, uma vez que os parâmetros do modelo variam de população para população e precisam, portanto, ter seus valores conhecidos. Com o objetivo de possibilitar esta caracterização, o presente trabalho propõe a formulação e solução do problema inverso, estimando os parâmetros do modelo a partir de dados da população utilizando para tal tarefa dois métodos Bayesianos.

Análise numérica do efeito de sobrecargas assimétricas em estacas: estudo de pontes nas obras do Arco Metropolitano.

O presente trabalho aborda o efeito de sobrecargas assimétricas em estacas, através do estudo de caso de encontros da Ponte do Rio Capivari nas obras do Arco Metropolitano. Neste caso específico as estacas foram construídas previamente a execução de aterros em terra armada e para minimizar o efeito de Tschebotarioff o solo foi reforçado com colunas de brita. Com embasamento na revisão bibliográfica apresentada foi realizada a análise pelo Método dos Elementos Finitos. Esta análise, realizada pelo programa computacional Plaxis, teve como foco principal a obtenção dos deslocamentos e momentos fletores nas estacas para os dois casos analisados: sem colunas de brita e com solo homogêneo equivalente (com colunas de brita). Foi possível verificar a influência da introdução das colunas de brita na redução dos momentos fletores e deslocamentos horizontais nas estacas dos encontros. Realizou-se ainda uma simulação em que as estacas seriam construídas após a realização dos aterros em terra armada onde pode-se constatar que o efeito de sobrecargas assimétricas seria mitigado. Foi efetuada também a comparação entre os deslocamentos provenientes de dados coletados da instrumentação de campo (inclinômetros) e os obtidos pelas análises numéricas, estando os mesmos compatíveis entre si, demonstrando que a metodologia adotada para simulação das colunas de brita no Método dos Elementos finitos foi adequada.

Condições de contorno tipo albedo para cálculos globais de reatores nucleares com o modelo de dois grupos de energia na formulação de transporte SN em geometria bidimensional cartesiana

Os eventos de fissão nuclear, resultados da interação dos nêutrons com os núcleos dos átomos do meio hospedeiro multiplicativo, não estão presentes em algumas regiões dos reatores nucleares, e.g., moderador, refletor, e meios estruturais. Nesses domínios espaciais não há geração de potência nuclear térmica e, além disso, comprometem a eficiência computacional dos cálculos globais de reatores nucleares. Propomos nesta tese uma estratégia visando a aumentar a eficiência computacional dessas simulações eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas (baffle e refletor) em torno do núcleo ativo. Apresentamos algumas modelagens e discutimos a eficiência da aplicação dessas condições de contorno aproximadas tipo albedo para uma e duas regiões nãomultiplicativas, na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. A denominação Albedo, palavra de origem latina para alvura, foi originalmente definida como a fração da luz incidente que é refletida difusamente por uma superfície. Esta denominação latina permaneceu como o termo científico usual em astronomia e, nesta tese, este conceito é estendido para reflexão de nêutrons. Estas condições de contorno tipo albedo SN não-convencional substituem aproximadamente as regiões de baffle e refletor no em torno do núcleo ativo do reator, desprezando os termos de fuga transversal no interior dessas regiões. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta tese, são exatas. Por eficiência computacional entende-se a análise da precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos considerando dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência.

Otimização global topográfica aplicada a problemas de otimização com restrições de igualdade e desigualdade

Os métodos de otimização que adotam condições de otimalidade de primeira e/ou segunda ordem são eficientes e normalmente esses métodos iterativos são desenvolvidos e analisados através da análise matemática do espaço euclidiano n-dimensional, o qual tem caráter local. Esses métodos levam a algoritmos iterativos que são usados para o cálculo de minimizadores globais de uma função não linear, principalmente não-convexas e multimodais, dependendo da posição dos pontos de partida. Método de Otimização Global Topográfico é um algoritmo de agrupamento, o qual é fundamentado nos conceitos elementares da teoria dos grafos, com a finalidade de gerar bons pontos de partida para os métodos de busca local, com base nos pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem como objetivo a aplicação do método de Otimização Global Topográfica junto com um método robusto e eficaz de direções viáveis por pontos-interiores a problemas de otimização que tem restrições de igualdade e/ou desigualdade lineares e/ou não lineares, que constituem conjuntos viáveis com interiores não vazios. Para cada um destes problemas, é representado também um hiper-retângulo compreendendo cada conjunto viável, onde os pontos amostrais são gerados.

Estimativa de fontes externas uniformes de nêutrons em regiões combustíveis para uma distribuição prescrita de potência

Projetos de reatores nucleares foram classificados em quatro gerações (Gen) pelo Departamento de Energia dos Estados Unidos da América (DOE), quando o DOE introduziu o conceito de reatores de geração IV (Gen IV). Reatores Gen IV são um conjunto de projetos de reator nuclear, em sua maioria teóricos, atualmente sendo pesquisados. Entre os projetos Gen IV, incluem-se os projetos dos ADS (Accelerator Driven Systems), que são sistemas subcríticos estabilizados por fontes externas estacionárias de nêutrons. Estas fontes externas de nêutrons são normalmente geradas a partir da colisão de prótons com alta energia contra os núcleos de metais pesados presentes no núcleo do reator, fenômeno que é conhecido na literatura como spallation, e os prótons são acelerados num acelerador de partículas que é alimentado com parte da energia gerada pelo reator. A criticalidade de um sistema mantido por reações de fissão em cadeia depende do balanço entre a produção de nêutrons por fissão e a remoção por fuga pelos contornos e absorção de nêutrons. Um sistema está subcrítico quando a remoção por fuga e absorção ultrapassa a produção por fissão e, portanto, tende ao desligamento. Entretanto, qualquer sistema subcrítico pode ser estabilizado pela inclusão de fontes estacionárias de nêutrons em seu interior. O objetivo central deste trabalho é determinar as intensidades dessas fontes uniformes e isotrópicas de nêutrons, que se deve inserir em todas as regiões combustíveis do sistema, para que o mesmo estabilize-se gerando uma distribuição prescrita de potência elétrica. Diante do exposto, foi desenvolvido neste trabalho um aplicativo computacional em linguagem Java que estima as intensidades dessas fontes estacionárias de nêutrons, que devem ser inseridas em cada região combustível para que estabilizem o sistema subcrítico com uma dada distribuição de potência definida pelo usuário. Para atingir este objetivo, o modelo matemático adotado foi a equação unidimensional de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas (SN) e o convencional método de malha fina diamond difference (DD) foi utilizado para resolver numericamente os problemas SN físicos e adjuntos. Resultados numéricos para dois problemas-modelos típicos são apresentados para ilustrar a acurácia e eficiência da metodologia proposta.