Uma Simulação computacional do passeio aleatório simples

Em 1828 foi observado um fenômeno no microscópio em que se visualizava minúsculos grãos de pólen mergulhados em um líquido em repouso que mexiam-se de forma aleatória, desenhando um movimento desordenado. A questão era compreender este movimento. Após cerca de 80 anos, Einstein (1905) desenvolveu uma formulação matemática para explicar este fenômeno, tratado por movimento Browniano, teoria cada vez mais desenvolvida em muitas das áreas do conhecimento, inclusive recentemente em modelagem computacional. Objetiva-se pontuar os pressupostos básicos inerentes ao passeio aleatório simples considerando experimentos com e sem problema de valor de contorno para melhor compreensão ao no uso de algoritmos aplicados a problemas computacionais. Foram explicitadas as ferramentas necessárias para aplicação de modelos de simulação do passeio aleatório simples nas três primeiras dimensões do espaço. O interesse foi direcionado tanto para o passeio aleatório simples como para possíveis aplicações para o problema da ruína do jogador e a disseminação de vírus em rede de computadores. Foram desenvolvidos algoritmos do passeio aleatório simples unidimensional sem e com o problema do valor de contorno na plataforma R. Similarmente, implementados para os espaços bidimensionais e tridimensionais,possibilitando futuras aplicações para o problema da disseminação de vírus em rede de computadores e como motivação ao estudo da Equação do Calor, embora necessita um maior embasamento em conceitos da Física e Probabilidade para dar continuidade a tal aplicação.

Análise das deformações dos aterros que compõem a obra do Arco Metropolitano do Rio de Janeiro.

O presente trabalho tem como objetivo a análise numérica do comportamento dos aterros instrumentados que compõem a obra do Arco Metropolitano do Rio de Janeiro. Os resultados da instrumentação de campo serão retroanalisados, juntamente com dados disponíveis na literatura, buscando-se a definição de parâmetros confiáveis, representativos do comportamento da argila compressível da região em estudo. O Arco Metropolitano do Rio de Janeiro é uma rodovia projetada que servirá como ligação entre as principais rodovias que cortam o município do Rio de Janeiro. Devido a presença de grandes espessuras de solo compressível em alguns trechos da região, cortados pelo traçado da rodovia, instrumentos de monitoramento, como placas de recalque e inclinômetros, foram utilizados para avaliar os deslocamentos verticais e horizontais dos aterros durante o processo construtivo. Para este trabalho foram selecionados trechos de aterros representativos, devido à magnitude dos recalques, qualidade do resultado da instrumentação e diferentes métodos construtivos. A partir da análise dos parâmetros de projeto e dos parâmetros encontrados na literatura, procede-se à simulação numérica do processo construtivo dos aterros selecionados com o programa PLAXIS, de elementos finitos, através de modelagem bidimensional. Os resultados numéricos são confrontados com a instrumentação de campo (placas de recalque e inclinômetros) e com os resultados de previsões teóricas (teoria de adensamento unidimensional). Os resultados comprovaram que a modelagem numérica mostrou-se uma ferramenta adequada para a previsão dos recalques totais, tempos de adensamento e ganho de resistência ao longo do tempo. A retroanálise do comportamento de aterros sobre solos moles permite a reavaliação das premissas de projeto, uma vez que as limitações das teorias de análise e a dificuldade na seleção de parâmetros, muitas vezes acarretam em estimativas de recalque incoerentes com as observações de campo.

Impacto de melhorias em infraestrutura de transportes sobre variáveis econômicas no Brasil: uma abordagem com o modelo GTAP

Neste Trabalho, modelamos a redução de custos de transportes decorrentes de hipotéticas melhorias na qualidade da infraestrutura de transportes terrestres da economia brasileira. A liberação de recursos gerada por essas melhorias foram traduzidas, em nosso modelo, como ganhos de produtividade para dez setores selecionados. A escolha dessa modelagem possui a vantagem de deixar que os agentes econômicos decidam o que fazer após o choque na produtividade total dos fatores, sem que estes recursos estejam comprometidos com quaisquer incentivos de natureza específica, como seria o caso de um subsídio à compra de determinada matéria-prima, por exemplo. O comportamento do setor, em consequência da realização do experimento, fica condicionado às hipóteses iniciais do modelo referentes aos parâmetros de decisão da firma, das famílias e do governo. Para os cálculos dos impactos sobre a produtividade, utilizamos as matrizes de insumoproduto estimadas por Martinez(2013). Os resultados mostraram ganhos expressivos para o produto interno bruto, a Balança Comercial e o volume de produção setorial, dentre outras variáveis analisadas.

Impacto de melhorias em infraestrutura de transportes sobre variáveis econômicas no Brasil: uma abordagem com o modelo GTAP

Neste Trabalho, modelamos a redução de custos de transportes decorrentes de hipotéticas melhorias na qualidade da infraestrutura de transportes terrestres da economia brasileira. A liberação de recursos gerada por essas melhorias foram traduzidas, em nosso modelo, como ganhos de produtividade para dez setores selecionados. A escolha dessa modelagem possui a vantagem de deixar que os agentes econômicos decidam o que fazer após o choque na produtividade total dos fatores, sem que estes recursos estejam comprometidos com quaisquer incentivos de natureza específica, como seria o caso de um subsídio à compra de determinada matéria-prima, por exemplo. O comportamento do setor, em consequência da realização do experimento, fica condicionado às hipóteses iniciais do modelo referentes aos parâmetros de decisão da firma, das famílias e do governo. Para os cálculos dos impactos sobre a produtividade, utilizamos as matrizes de insumoproduto estimadas por Martinez(2013). Os resultados mostraram ganhos expressivos para o produto interno bruto, a Balança Comercial e o volume de produção setorial, dentre outras variáveis analisadas.

Estudo do equilíbrio líquido-líquido e da estabilidade de sistemas envolvendo biodiesel, álcoois e água

O estudo da estabilidade de emulsões de biocombustíveis destaca-se, dentre outros aspectos, pela necessidade de identificação das características destas. Esta necessidade reforça-se pela dificuldade de detecção de biodiesel e suas misturas quando estão fora de especificação, além de se aplicar aos processos de separação e purificação do biodiesel após a transesterificação. O entendimento deste problema poderá ser alcançado a partir de estudos de estabilidade associados aos modelos termodinâmicos de equilíbrio líquido-líquido (modelo NRTL) em sistemas ternários, que serão a base da modelagem proposta neste trabalho. O levantamento de dados envolve técnicas de observação direta em titulação, e ensaios de estabilidade e tensiometria, com algumas alterações das técnicas clássicas. Desta forma, o objetivo principal do trabalho é propor uma maneira de identificar e caracterizar as interfaces das emulsões formadas e estudar a estabilidade das fases envolvidas. Além disso, estuda-se a introdução de alcoóis como agentes tensoativos para viabilizar a utilização de misturas de combustíveis, mesmo na região de mais de uma fase. Para o desenvolvimento do trabalho, emulsões envolvendo os constituintes básicos serão combinadas entre si em diferentes proporções. A etapa inicial constitui-se de um mapeamento das misturas e proporções que formam as emulsões, baseadas nas curvas de equilíbrio. Em seguida, faz-se o estudo de caracterização e estabilidade destes sistemas. Com estas informações é possível mapear o sentido de incremento de tensão na interface e instabilidade da separação de fases das emulsões em relação aos diagramas ternários correspondentes. Por fim, propõe-se uma estratégia para estimar os parâmetros termodinâmicos do NRTL para as espécies presentes nos sistemas estudados e prever o comportamento dos mesmos e de sistemas similares. Este modelo matemático proposto tem como objetivo complementar a técnica experimental adotada para construção das curvas de equilíbrio, sendo bem sucedido

Estudo do comportamento estrutural e análise de conforto humano de edifícios de concreto armado.

Atualmente, os projetos de edifícios altos necessitam cada vez mais de sistemas estruturais simples, que agilizem sua montagem, reduzindo os custos e promovendo maior flexibilidade de utilização para os espaços construídos. Com essa finalidade, estruturas com poucas vigas vêm sendo muito utilizadas. Entretanto, o sistema estrutural com poucas vigas pode ocasionar dois tipos de problemas, relacionados entre si, a saber: diminuição do sistema de contraventamento da edificação e vibrações excessivas. Portanto, é fundamental, nesses casos, a verificação da estabilidade global da estrutura, utilizando índices de sensibilidade além de outros parâmetros de projeto, como também, o desenvolvimento de um estudo minucioso acerca do conforto humano da edificação. Assim sendo, neste trabalho de pesquisa foram investigados quatro modelos estruturais de edifícios altos de concreto armado, com base no estudo da variação entre o número de pavimentos e a quantidade de vigas existentes em cada modelo, objetivando-se verificar quais os efeitos que tais variações podem vir a gerar sobre a estabilidade global e, bem como, sobre o conforto humano dos sistemas estruturais investigados. A modelagem numérica dos edifícios em estudo foi realizada através do emprego do programa ANSYS e, para tal, foram utilizadas técnicas básicas de discretização, por meio do método dos elementos finitos. As conclusões alcançadas ao longo da investigação versam acerca do estudo da resposta estrutural estática e dinâmica dos edifícios, no que diz respeito as variações dos valores dos parâmetros de instabilidade, dos valores dos deslocamentos e esforços, e, bem como, dos níveis de conforto humano de cada modelo estrutural analisado.

Arcabouço crustal profundo da parte Centro-Norte da margem de Angola: modelo de afinamento e contato de crostas

Este trabalho tem como objetivo a identificação de feições que permitam (1) a construção do arcabouço crustal profundo e da porção superior do manto em parte da Costa de Angola, (2) a comparação deste arcabouço com o embasamento aflorante e (3) a tentativa de adequar estes resultados aos modelos de ruptura continental já publicados. Para alcançar estes objetivos foi feita a interpretação de cinco linhas sísmicas de reflexão profundas (25 Km de profundidade) na costa de Angola, nas Bacias de Kwanza e Baixo Congo, adquiridas pela ION-GXT. As feições identificadas na sísmica de reflexão auxiliaram na determinação dos limites da crosta continental superior e inferior, no reconhecimento das camadas que compõem a crosta oceânica e na identificação da Descontinuidade de Mohorovicic (que marca o limite entre crosta e manto). A interpretação sísmica associada a dados da literatura (que proporcionaram valores de densidade para os pacotes identificados na interpretação sísmica) permitiram a realização de uma modelagem gravimétrica que foi comparada ao dado gravimétrico adquirido durante a aquisição sísmica. A modelagem gravimétrica serviu para validar a interpretação sísmica, atuando como um controle de qualidade para a interpretação. Caso a anomalia gravimétrica gerada pela modelagem não estivesse de acordo com a anomalia medida, a interpretação sísmica era revista na tentativa de um melhor ajuste entre o resultado modelado e o medido. Este ajuste, no entanto, sempre foi feito honrando os refletores que estavam bem marcados na sísmica. Somado a isto, ainda foi utilizado o dado magnético adquirido no campo, no auxilio da interpretação. O arcabouço crustal obtido com a utilização deste método permitiu a comparação dos resultados da interpretação com os modelos de evolução de margens passivas existentes na literatura, mostrando muitos pontos em comum aos modelos que defendem a possibilidade de ocorrência de manto exumado em margens passivas pobres em magmatismo. A interpretação final destes dados mostrou um domínio proximal marcado por uma crosta continental espessa porém pouco afinada em contato com um domínio distal marcado por uma crosta continental muito afinada (crosta hiper-estirada) e, em direção ao centro do oceano, uma região em que ocorre a exumação do manto. A passagem do domínio proximal para o distal ocorre de forma rápida em uma região denominada Zona de Estrangulamento. À oeste do manto exumado é possível identificar a crosta oceânica. O cruzamento dos resultados obtidos neste trabalho com dados do embasamento aflorante no continente africano sugerem um controle do deste nos valores finais de afinamento da crosta continental sob a bacia e nas regiões de manto exumado. Trabalhos recentes realizados na costa de Angola e do Brasil mostram feições semelhantes às identificadas nesta dissertação.

Modelos de agrupamento e classificação para os bairros da cidade do Rio de Janeiro sob a ótica da Inteligência Computacional: Lógica Fuzzy, Máquinas de Vetores Suporte e Algoritmos Genéticos

A partir de 2011, ocorreram e ainda ocorrerão eventos de grande repercussão para a cidade do Rio de Janeiro, como a conferência Rio+20 das Nações Unidas e eventos esportivos de grande importância mundial (Copa do Mundo de Futebol, Olimpíadas e Paraolimpíadas). Estes acontecimentos possibilitam a atração de recursos financeiros para a cidade, assim como a geração de empregos, melhorias de infraestrutura e valorização imobiliária, tanto territorial quanto predial. Ao optar por um imóvel residencial em determinado bairro, não se avalia apenas o imóvel, mas também as facilidades urbanas disponíveis na localidade. Neste contexto, foi possível definir uma interpretação qualitativa linguística inerente aos bairros da cidade do Rio de Janeiro, integrando-se três técnicas de Inteligência Computacional para a avaliação de benefícios: Lógica Fuzzy, Máquina de Vetores Suporte e Algoritmos Genéticos. A base de dados foi construída com informações da web e institutos governamentais, evidenciando o custo de imóveis residenciais, benefícios e fragilidades dos bairros da cidade. Implementou-se inicialmente a Lógica Fuzzy como um modelo não supervisionado de agrupamento através das Regras Elipsoidais pelo Princípio de Extensão com o uso da Distância de Mahalanobis, configurando-se de forma inferencial os grupos de designação linguística (Bom, Regular e Ruim) de acordo com doze características urbanas. A partir desta discriminação, foi tangível o uso da Máquina de Vetores Suporte integrado aos Algoritmos Genéticos como um método supervisionado, com o fim de buscar/selecionar o menor subconjunto das variáveis presentes no agrupamento que melhor classifique os bairros (Princípio da Parcimônia). A análise das taxas de erro possibilitou a escolha do melhor modelo de classificação com redução do espaço de variáveis, resultando em um subconjunto que contém informações sobre: IDH, quantidade de linhas de ônibus, instituições de ensino, valor m médio, espaços ao ar livre, locais de entretenimento e crimes. A modelagem que combinou as três técnicas de Inteligência Computacional hierarquizou os bairros do Rio de Janeiro com taxas de erros aceitáveis, colaborando na tomada de decisão para a compra e venda de imóveis residenciais. Quando se trata de transporte público na cidade em questão, foi possível perceber que a malha rodoviária ainda é a prioritária

Aspectos de segurança, persistência e nomeação em agentes de recursos para Internet das coisas

A realização da Internet das Coisas (Internet of Things, IoT) requer a integração e interação de dispositivos e serviços com protocolos de comunicação heterogêneos. Os dados gerados pelos dispositivos precisam ser analisados e interpretados em concordância com um modelo de dados em comum, o que pode ser solucionado com o uso de tecnologias de modelagem semântica, processamento, raciocínio e persistência de dados. A computação ciente de contexto possui soluções para estes desafios com mecanismos que associam os dados de contexto com dados coletados pelos dispositivos. Entretanto, a IoT precisa ir além da computação ciente de contexto, sendo simultaneamente necessário soluções para aspectos de segurança, privacidade e escalabilidade. Para integração destas tecnologias é necessário o suporte de uma infraestrutura, que pode ser implementada como um middleware. No entanto, uma solução centralizada de integração de dispositivos heterogêneos pode afetar escalabilidade. Assim esta integração é delegada para agentes de software, que são responsáveis por integrar os dispositivos e serviços, encapsulando as especificidades das suas interfaces e protocolos de comunicação. Neste trabalho são explorados os aspectos de segurança, persistência e nomeação para agentes de recursos. Para este fim foi desenvolvido o ContQuest, um framework, que facilita a integração de novos recursos e o desenvolvimento de aplicações cientes de contexto para a IoT, através de uma arquitetura de serviços e um modelo de dados. O ContQuest inclui soluções consistentes para os aspectos de persistência, segurança e controle de acesso tanto para os serviços de middleware, como para os Agentes de Recursos, que encapsulam dispositivos e serviços, e aplicações-clientes. O ContQuest utiliza OWL para a modelagem dos recursos e inclui um mecanismo de geração de identificadores únicos universais nas ontologias. Um protótipo do ContQuest foi desenvolvido e validado com a integração de três Agentes de Recurso para dispositivos reais: um dispositivo Arduino, um leitor de RFID e uma rede de sensores. Foi também realizado um experimento para avaliação de desempenho dos componentes do sistema, em que se observou o impacto do mecanismo de segurança proposto no desempenho do protótipo. Os resultados da validação e do desempenho são satisfatórios

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.

Análise da evolução termomecânica da Bacia do Amazonas

Os dados geológicos e geofísicos escolhidos para o tema de estudo pertencem a Bacia do Amazonas, na região centro-norte do Brasil. A Bacia do Amazonas é uma bacia intracratônica com cerca de 500.000 km. A mesma está limitada ao norte pelo Escudo das Guianas e ao sul pelo Escudo Brasileiro. O limite oeste com a Bacia do Solimões é marcado pelo Arco de Purus, ao passo que o Arco de Gurupá constitui seu limite leste. Possui características inerentes a uma bacia intracratônica paleozóica, com uma longa história evolutiva, marcada por discordâncias expressivas e com uma cunha sedimentar relativamente rasa se comparada às bacias cretáceas brasileiras, o que levanta controvérsia a respeito da suficiência do soterramento para a eficiência de geração de hidrocarboneto. Podem ser reconhecidas nos 5000 m do preenchimento sedimentar da Bacia do Amazonas, duas seqüências de primeira ordem: uma paleozóica, intrudida por diques e soleiras de diabásio, na passagem do Triássico para o Jurássico, e uma mesozóica-cenozóica que representam um aspecto importante na evolução térmica da matéria orgânica que ocorre na primeira seqüência. Com relação à exploração de petróleo, apesar do fomento exploratório ocorrido nos últimos anos, a bacia ainda é considerada pouco explorada sendo sua maior reserva a da província de Urucu. Um dos fatores que dificultam bastante a exploração desta bacia assim como a bacia do Solimões a oeste é o acesso restrito, pois estão situadas em áreas remotas e florestadas, de difícil acesso, com muitas reservas indígenas e florestais, o que causa restrições logísticas, operacionais e legais. O efeito térmico das intrusões ígneas é considerado como o responsável pelo acréscimo de calor necessário à maturação da matéria orgânica e conseqüente geração de hidrocarbonetos. Este trabalho contribui com a reconstrução da história térmica desta bacia a partir da modelagem das variáveis termais e da história de soterramento. Para isso, foram utilizados modelos consagrados na literatura, que permitem, de forma simples, a estimativa do fluxo térmico através do embasamento e da seqüência sedimentar. Na análise da influência de intrusões ígneas na estrutura térmica da bacia, o modelo bidimensional desenvolvido pelo método de diferenças finitas se mostrou apropriado. Utilizou-se o fluxo térmico basal calculado nas condições de contorno da modelagem da influência térmica das ígneas. Como resultado obteve-se a estruturação térmica da bacia e a historia maturação de suas rochas geradoras

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.

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.