Modelos para acontecimentos múltiplos

Ao longo dos últimos anos tem existido um interesse crescente em estudar o tempo até à ocorrência de vários acontecimentos, os quais podem ser observados mais do que uma vez para um mesmo indivíduo. Os dados respeitantes a acontecimentos múltiplos têm como principal característica o facto de se registar mais do que um tempo de vida para cada indivíduo, o que inviabiliza a aplicação direta do modelo de regressão de Cox. Assim, surgiu a necessidade de desenvolver novas extensões deste modelo, sendo que as mais utilizadas na prática foram sugeridas por: Prentice, Williams e Peterson (PWP); Andersen e Gill (AG); Wei, Lin e Weissfeld (WLW); e Lee, Wei e Amato (LWA). Um dos maiores obstáculos na aplicação destes modelos é a forte possibilidade de existir correlação intraindivíduos. Neste ponto, os quatro modelos referidos anteriormente são classificados como modelos marginais, uma vez que o vetor de parâmetros de regressão é estimado com base no ajustamento de um modelo que ignora a correlação entre acontecimentos. Para compensar esse facto, nestes modelos é usado um estimador robusto da matriz de covariância, o qual permite efetuar a correção necessária na estimativa da variância usual. Após realizar uma descrição detalhada de cada modelo marginal procedeu-se à respetiva implementação através do software estatístico R. Para o efeito, recorreu-se à simulação de dados relativos a acontecimentos múltiplos do mesma natureza, isto é, a acontecimentos recorrentes. Os resultados obtidos permitiram realçar e confirmar as características dos modelos estudados.

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Projeto Expande com balanço de uma experiência de programação a partir das necessidades.

O estudo que se segue é uma proposta de reflexão sobre a programação de saúde para a população usuária do SUS, em especial a assistência oncológica, as sucessivas estratégias político-administrativas e medidas de soluções a serem ofertadas à população portadora de câncer. Percebe-se o desequilíbrio entre a operacionalidade real e a operacionalidade suposta ideal, para dar conta dos casos novos de câncer (CNC)/ano, a atenção à saúde concentrada nos grandes centros urbanos, reflexo de um estímulo a investimentos em saúde direcionados aos centros urbanos, localidades de maior crescimento econômico e social. Objetivando o balanço de uma experiência de programação a partir das necessidades, tomamos como norte o Projeto Expande, para a análise da real assistência em oncologia no país, que se traduz no conjunto: oferta de serviço necessidade de tratamento específico para a população demanda, identificando déficit e necessidade. Nesse sentido, com base na programação e nas informações de produtividade SUS/Brasil, disponíveis no INCA , optamos por acompanhar os dados referentes ao ano 2008, quantificados e analisados, tendo como ponto de partida a estimativa de CNC/ano 2008, parâmetros assistenciais estimativa de necessidade máxima). Este estudo, com esse modelo, objetivou levantar subsídios que supostamente contribuam para o aprimoramento da Política de Atenção ao Câncer, em especial uma Política de Expansão Oncológica para o país.

Datu-trafikoaren monitorizazio-sistemetarako aplikatutako ereduen ebaluazio-erreminta baten diseinua eta garapena

[EU]Lan honetan, software diseinu bat sortu nahi da, zeinaren bidez datu-trafikoen monitorizazio sistemak aztertzeko garatu den eredu matematiko bat ebaluatuko den. Eredu horrentzat interfaze bat egin beharko da, eta interfaze horrek esandako eredua ebaluatzeko softwareaz gain, software gehiago biltzeko ahalmena eduki beharko du. Horrela, ikertzaileek Trafikoa Monitorizatzeko Sistemak aztertzeko sortzen diren eredu matematikoak sistema bakarra erabiliz ebaluatu ahalko dute ahalik eta modu errazenean.

Um método de matriz resposta com esquema iterativo de inversão parcial por região para problemas unidimensionais de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas

Um método de matriz resposta (RM) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas de transporte de nêutrons monoenergéticos e com fonte fixa, em geometria unidimensional na formulação de ordenadas discretas (SN). O método RM com esquema iterativo de inversão parcial por região (RBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações SN, afora os erros de arredondamento da aritmética finita computacional. Desenvolvemos um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em qualquer ponto do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Estudo do efeito das cargas de nutrientes no reservatório APM-Manso: uma abordagem através da modelagem numérica tridimensional.

Os gestores de recursos hídricos estão encarregados da gestão de longo prazo, a regulação e a proteção dos recursos hídricos. No entanto, reconhece-se que a estes gestores devem levar em conta a multiplicidade de usos dos recursos hídricos que são apresentadas pelas partes interessadas, tais como agricultores, fornecedores de água e grupos de ambientalistas. Dada a complexidade do sistema hidrológico, o desenvolvimento e a utilização de modelos matemáticos são muitas das vezes necessários. Neste contexto a modelagem ambiental é frequentemente realizada para avaliar os impactos da degradação do ecossistema devido à ação humana. Essa aplicação orientada a investigações proporciona um importante meio pelo qual os cientistas podem interagir e influenciar nas políticas a nível local, regional, nacional e internacional. No Mato Grosso, durante a implantação da hidroelétrica de Aproveitamento Múltiplo de Manso foram adotadas medidas de mitigação dos impactos socioeconômicos causados. Essas medidas geram uma tendência de aumento populacional associado a uma mudança das características sócio-econômicas, para toda a região do entorno do Reservatório, o que agrava o problema de poluição por nutrientes e denota que existe uma necessidade proeminente de estudos do impacto que estas cargas causariam no ecossistema do reservatório. Utilizando o modelo hidrodinâmico e termodinâmico tridimensional ELCOM, associado ao modelo biogeoquímico Caedym, este trabalho tem a finalidade de criar uma modelagem representativa das cargas dos principais nutrientes causadores da eutrofização cultural, sendo eles: a amônia (NH4), o nitrato (NO3) e o Ortofosfato (PO4), com a finalidade de estudar os efeitos das dinâmicas espaciais e temporais destas cargas no estado trófico deste reservatório no em torno dos pontos de lançamento de esgoto e na sua totalidade.

Simulações de problemas inversos com aplicações em engenharia nuclear usando técnicas de transporte de partículas neutras monoenergéticas na formulação unidimensional de ordenadas discretas

Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese.

Estudio quimiométrico de aceites lubricantes marinos

El proyecto realizado se basa la puesta en marcha de distintos procedimientos quimiométricos aplicados al control de los lubricantes presentes en sistemas de propulsión, para dos ensayos muy demandados por la importancia de las decisiones relacionadas con su análisis, el contenido de insolubles y el BN. Los sistemas de mantenimiento industrial se pueden dividir en varias etapas bajo la denominación de primera, segunda y tercera generación. La primera generación respondía a “reparar” en caso de avería. Una segunda generación se enfocaba a las “revisiones cíclicas” y una tercera generación iba enfocada al “análisis de las causas” y los efectos de los fallos. Es en esta tercera generación en donde tuvo el punto de arranque mi proyecto, habiendo enfocado mi atención en “optimizar el análisis de las causas” en función de la adaptación de las técnicas analíticas al siglo XXI. La herramienta usada para la puesta en marcha y la optimización de métodos alternativos a las normas usadas en el control y seguimiento de los lubricantes marinos en motores diesel para el análisis del contenido de insolubles y para la determinación del BN es la quimiometría, es decir, modelos matemáticos desarrollados a través de programas informáticos. Para el análisis de las muestras y la obtención de los modelos se ha utilizado la espectroscopía infrarroja mediante FTIR.

Redução de dimensionalidade em equilíbrio de fases fluidas

Atualmente, existem modelos matemáticos capazes de preverem acuradamente as relações entre propriedades de estado; e esta tarefa é extremamente importante no contexto da Engenharia Química, uma vez que estes modelos podem ser empregados para avaliar a performance de processos químicos. Ademais, eles são de fundamental importância para a simulação de reservatórios de petróleo e processos de separação. Estes modelos são conhecidos como equações de estado, e podem ser usados em problemas de equilíbrios de fases, principalmente em equilíbrios líquido-vapor. Recentemente, um teorema matemático foi formulado (Teorema de Redução), fornecendo as condições para a redução de dimensionalidade de problemas de equilíbrios de fases para misturas multicomponentes descritas por equações de estado cúbicas e regras de mistura e combinação clássicas. Este teorema mostra como para uma classe bem definidade de modelos termodinâmicos (equações de estado cúbicas e regras de mistura clássicas), pode-se reduzir a dimensão de vários problemas de equilíbrios de fases. Este método é muito vantajoso para misturas com muitos componentes, promovendo uma redução significativa no tempo de computação e produzindo resultados acurados. Neste trabalho, apresentamos alguns experimentos numéricos com misturas-testes usando a técnica de redução para obter pressões de ponto de orvalho sob especificação de temperaturas.

Modelagem computacional aplicada na estimação de parâmetros cinéticos em processos químicos e cromatográficos

A estimação de parâmetros cinéticos em processos químicos e cromatográficos utilizando técnicas numéricas assistidas por computadores tem conduzido para melhoria da eficiência e o favorecimento da compreensão das fenomenologias envolvidas nos mesmos. Na primeira parte deste trabalho será realizada a modelagem computacional do processo de produção de biodiesel via esterificação, sendo que, o método de otimização estocástica Random Restricted Window (R2W) será correlacionado com os dados experimentais da produção de biodiesel a partir da esterificação do ácido láurico com etanol anidro na presença do catalisador ácido nióbico (Nb2O5). Na segunda parte do mesmo será realizada a modelagem computacional do processo de cromatografia de adsorção (batch process) onde serão correlacionados os dados provenientes dos modelos cinéticos de HASHIM, CHASE e IKM2 com os dados experimentais da adsorção de amoxicilina com quitosana, e também serão correlacionados os dados experimentais da adsorção de Bovine Serum Albumin (BSA) com Streamline DEAE com os dados provenientes de uma nova aplicação do método R2W mediante a implementação de um modelo cinético reversível. Ademais, as constantes cinéticas para cada processo supracitado serão estimadas levando em consideração o valor mínimo da função resíduos quadrados.

Coste Óptimo para el tratamiento de limpieza de aguas residuales

[ES]El presente trabajo es un estudio sobre un problema de optimización matemático. En él, recabamos información de la Optimización Matemática como tal, y se analizan los modelos de optimización más relevantes, además de plantear y resolver un caso práctico. Primero, distinguiremos dos ramas principales: los métodos clásicos, en donde destacamos las denominadas Optimización Lineal y No Lineal; y los métodos metaheurísticos. Además, clasificaremos y plantearemos los modelos matemáticos en función de la presencia o no de restricciones. En la segunda parte del trabajo, con la finalidad de encuadrar el caso práctico, hablaremos sobre el problema, común a muchas zonas fuertemente industrializadas, de eliminación de las aguas residuales: sus antecedentes, situación actual y nuevas tendencias. Por último, plantearemos un problema práctico en este campo, lo modelizaremos matemáticamente y lo resolveremos con la ayuda de la herramienta Excel Solver. Al final del trabajo, expondremos las conclusiones derivadas del análisis del problema.

On some new mathematical models for infective diseases: analysis, equilibrium, positivity and vaccination controls

196 p.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Determinação de uma fonte de contaminantes no estuário do Rio Macaé através da técnica de problemas inversos

Um grande desafio da atualidade é a preservação dos recursos hídricos, bem como o correto manejo dos mesmos, frente à expansão das cidades e às atividades humanas. A qualidade de um corpo hídrico é usualmente avaliada através da análise de parâmetros biológicos, físicos e químicos. O comportamento de tais parâmetros pode convenientemente ser simulado através de modelos matemáticos e computacionais, que surgem assim como uma ferramenta bastante útil, por sua capacidade de geração de cenários que possam embasar, por exemplo, tomadas de decisão. Nesta tese são discutidas técnicas de estimação da localização e intensidade de uma fonte de contaminante conservativo, hipoteticamente lançado na região predominantemente fluvial de um estuário. O lançamento aqui considerado se dá de forma pontual e contínua e a região enfocada compreendeu o estuário do Rio Macaé, localizado na costa norte do Rio de Janeiro. O trabalho compreende a solução de um problema direto, que consiste no transporte bidimensional (integrado na vertical) desse contaminante hipotético, bem como a aplicação de técnicas de problemas inversos. Para a solução do transporte do contaminante, aqui modelada pela versão 2D horizontal da equação de advecção-difusão, foram utilizados como métodos de discretização o Método de Elementos Finitos e o Método de Diferenças Finitas. Para o problema hidrodinâmico foram utilizados dados de uma solução já desenvolvida para estuário do Rio Macaé. Analisada a malha de acordo com o método de discretização, foram definidas a geometria do estuário e os parâmetros hidrodinâmicos e de transporte. Para a estimação dos parâmetros propostos foi utilizada a técnica de problemas inversos, com o uso dos métodos Luus-Jaakola, Algoritmo de Colisão de Partículas e Otimização por Colônia de Formigas para a estimação da localização e do método Seção Áurea para a estimação do parâmetro de intensidade da fonte. Para a definição de uma fonte, com o objetivo de propor um cenário experimental idealizado e de coleta de dados de amostragem, foi realizada a análise de sensibilidade quanto aos parâmetros a serem estimados. Como os dados de amostragem de concentração foram sintéticos, o problema inverso foi resolvido utilizando-os com e sem ruído, esse introduzido de forma artificial e aleatória. Sem o uso de ruído, os três métodos mostraram-se igualmente eficientes, com uma estimação precisa em 95% das execuções. Já com o uso de dados de amostragem com ruídos de 5%, o método Luus-Jaakola mostrou-se mais eficiente em esforço e custo computacional, embora todos tenham estimado precisamente a fonte em 80% das execuções. Considerando os resultados alcançados neste trabalho tem-se que é possível estimar uma fonte de constituintes, quanto à sua localização e intensidade, através da técnica de problemas inversos. Além disso, os métodos aplicados mostraram-se eficientes na estimação de tais parâmetros, com estimações precisas para a maioria de suas execuções. Sendo assim, o estudo do comportamento de contaminantes, e principalmente da identificação de fontes externas, torna-se uma importante ferramenta para a gestão dos recursos hídricos, possibilitando, inclusive, a identificação de possíveis responsáveis por passivos ambientais.

Identificação de danos em estruturas bi-dimensionais via matriz de flexibilidade baseado em um modelo de dano contínuo

No presente trabalho, o modelo de identificação de danos apresentado por Stutz et al. (2005) é utilizado. A contribuição do presente trabalho consiste em avaliar alguns pontos da identificação de danos em vigas e, em seguida, expandir o modelo para identificar danos estruturais em placas. Uma avaliação do comportamento das frequências naturais e da matriz de flexibilidade para uma viga de Euler-Bernoulli simplesmente apoiada, na presença de danos simulados pelo parâmetro de coesão é realizada. Essa análise, permite também o conhecimento das regiões onde há maior sensibilidade ao dano, ajudando a traçar estratégias para melhorar a identificação de danos em regiões que sofrem poucas alterações na presença de falhas estruturais. Comparou-se o comportamento dos dois primeiros modos de vibração da viga simplesmente apoiada na presença de um dano estrutural, com os dois primeiros modos de vibração da estrutura intacta e corrompidos por ruído. Diversos métodos de localização de danos e de otimização são avaliados na tentativa de identificar os danos simulados através do campo de danos proposto por Stutz et al. (2005) na presença de dados ruidosos. Após a apresentação de resultados da identificação de danos obtidos para uma viga de Euler-Bernoulli simplesmente apoiada, uma análise do comportamento das frequências naturais e da matriz de flexibilidade de uma viga de Euler- Bernoulli engastada na presença de danos é apresentada, assim como os resultados de identificação de danos considerando-se diversos cenários e níveis de ruído. Uma importante contribuição do presente trabalho consiste em propor um método de identificação de danos via matriz de flexibilidade onde o campo de defeitos para a placa de Kirchoff é modelado via MEF. Uma análise do comportamento da matriz de flexibilidade devido à presença de danos na placa é apresentada, assim como os resultados numéricos da identificação de danos estruturais com e sem a presença de dados ruidosos. Com a finalidade de reduzir o custo computacional na identificação de danos em estruturas complexas, uma hibridização entre o método de otimização por enxame de particulas (PSO, do inglês, Particle Swarm Optimization) e o método de otimização Levenberg-Marquardt é proposta. Resultados numéricos da hibridização para uma estrutura do tipo placa são apresentados.