Problemas inversos em processos difusivos com retenção

Um Estudo para a solução numérica do modelo de difusão com retenção, proposta por Bevilacqua et al. (2011), é apresentado, bem como uma formulação implícita para o problema inverso para a estimativa dos parâmetros envolvidos na formulação matemática do modelo. Através de um estudo minucioso da análise de sensibilidade e do cálculo do coeficiente de correlação de Pearson, são identificadas as chances de se obter sucesso na solução do problema inverso através do método determinístico de Levenberg-Marquardt e dos métodos estocásticos Algoritmo de Colisão de Partículas (Particle Collision Algorithm - PCA) e Evolução Diferencial (Differential Evolution - DE). São apresentados os resultados obtidos através destes três métodos de otimização para três casos de conjunto de parâmetros. Foi observada uma forte correlação entre dois destes três parâmetros, o que dificultou a estimativa simultânea dos mesmos. Porém, foi obtido sucesso nas estimativas individuais de cada parâmetro. Foram obtidos bons resultados para os fatores que multiplicam os termos diferenciais da equação que modela o fenômeno de difusão com retenção.

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

Teoria de controle ótimo com aplicações a sistemas biológicos

Neste trabalho apresentamos as etapas para a utilização do método da Programação Dinâmica, ou Princípio de Otimização de Bellman, para aplicações de controle ótimo. Investigamos a noção de funções de controle de Lyapunov (FCL) e sua relação com a estabilidade de sistemas autônomos com controle. Uma função de controle de Lyapunov deverá satisfazer a equação de Hamilton-Jacobi-Bellman (H-J-B). Usando esse fato, se uma função de controle de Lyapunov é conhecida, será então possível determinar a lei de realimentação ótima; isto é, a lei de controle que torna o sistema globalmente assintóticamente controlável a um estado de equilíbrio. Como aplicação, apresentamos uma modelagem matemática adequada a um problema de controle ótimo de certos sistemas biológicos. Este trabalho conta também com um breve histórico sobre o desenvolvimento da Teoria de Controle de forma a ilustrar a importância, o progresso e a aplicação das técnicas de controle em diferentes áreas ao longo do tempo.

Al Berto e a poética do fogo

O fogo, do latim focus, conota uma multiplicidade de imagens levadas sempre para a intensidade, para o extremo, porque por si só, o fogo é exuberante, fascinante, encantador, deriva-se em chamas, calor, brasa, paixão, erotismo, lume, incêndio e, na arte, uma imaginação ardente. A partir dessas considerações, pretende-se mostrar que a obra do poeta português Al Berto, Horto de Incêndio, relaciona-se com a imagem poética do fogo, desdobrando-se ora em sexualidade, por meio da erotização da palavra, ora em morte, por intermédio do uso dessa imagem como profecia do seu próprio fim, antecipado nas páginas que se queimam no jardim de incêndio. Para tanto, apóia-se a presente pesquisa nos livros do filósofo francês Gaston Bachelard, acerca do fogo, suas simbologias e significações na arte poética

Métodos espectrais aplicados à relatividade numérica: determinação dos dados iniciais

Neste trabalho aplicamos métodos espectrais para a determinação da configuração inicial de três espaços-tempos contendo buracos negros. Para isto apresentamos primeiro a foliação do espaço-tempo em hipersuperfícies tridimensionais espaciais parametrizadas pela função temporal t. Este processo é chamado de decomposição 3+1 [2] [5]. O resultado deste processo são dois conjuntos de equações classificadas em equações de vínculo e evolução [4]. As equações de vínculo podem ser divididas em vínculos Hamiltoniano e dos momentos. Para a obtenção dos dados iniciais dos problemas estudados aqui, apenas a equação de vínculo Hamiltoniano será resolvida numericamente, pois as equações de vínculo dos momentos possuem solução analítica nestes casos. Uma pequena descrição dos métodos espectrais é apresentada, destacando-se os método de Galerkin, método pseudoespectral ou de colocação e método de Tau, que são empregados na resolução das equações de vínculo Hamiltoniano dos problemas estudados. Verificamos que os resultados obtidos neste trabalho superam aqueles produzidos por Kidder e Finn [15], devido a uma escolha diferente das funções de base, que aqui satisfazem uma das condições de contorno.

Comparação do Gasto energético de repouso obtido pela calorimetria indireta com as equações de predição em pacientes idosos em tratamento crônico de hemodiálise

Assim como na população geral, as necessidades energéticas diárias dos pacientes em tratamento crônico de hemodiálise (HD) podem ser calculadas multiplicando-se o gasto energético de repouso (GER) pelo nível de atividade física. Até o momento, não há estudos que avaliaram se as equações de predição são precisas para se estimar o GER de idosos em HD. O objetivo do presente estudo foi avaliar a concordância entre o GER obtido pela calorimetria indireta e as equações de predição de Harris&Benedict, Schofield e a proposta pelo documento da Organização Mundial de Saúde de 1985 (FAO 1985) nos pacientes idosos em HD. Tratou-se de um estudo transversal, onde foi avaliado o GER de 57 pacientes idosos não institucionalizados (> 60anos) em tratamento crônico de HD mensurado pela calorimetria indireta e comparado com as equações de predição de Harris&Benedict, Schofield e FAO 1985.A concordância entre o GER medido e as equações foi realizada pelo coeficiente de correlação intraclasse e pela análise de Bland-Altman. Neste estudo pode-se observar que o GER estimado pelas 3 equações foi significantemente maior do que o obtido pela calorimetria indireta. Um grau de reprodutibilidade moderado foi observado entre a calorimetria indireta e as equações. A superestimação foi o principal erro observado, sendo presente na metade dos pacientes. A subestimação foi vista em aproximadamente em 10 % dos pacientes. Com base nesses achados podemos concluir que as 3 equações tiveram uma performance similar ao estimar o GER. E estas podem ser utilizadas para calcular o GER de idosos em HD, na medida em que os nutricionistas reconheçam seus possíveis erros, principalmente quando as equações de predição subestimam o GER medido.

Representação social da confiança

Esta tese focaliza como tema principal a representação social da confiança. Qual a importância da confiança para a sociedade contemporânea, para as organizações de produção e para as relações interpessoais. A fundamentação teórica é a Teoria das Representações Sociais tomada em sua equação: representação = figura/significado. A metodologia da pesquisa consistiu na realização de grupos focais e na expressão imagética e verbal da confiança. A análise dos textos visuais conjuntamente com os verbais possibilitou uma percepção bastante acurada desta representação social.

Estudo do comportamento em fluência de juntas metálicas coladas

Juntas coladas têm se mostrado como a principal alternativa de união entre componentes fabricados em materiais distintos ou não, tendo como vantagem principal a ausência de aporte de calor e melhor distribuição de tensões. A literatura mostra que estudos anteriores foram realizados com o objetivo de caracterizar diversas propriedades associados às juntas, porém, pouco ainda se conhece sobre o comportamento dessas uniões relacionado à fluência. Os ensaios de fluência têm como principal desvantagem a baixa produtividade de resultados visto que, na maioria dos casos, os equipamentos aptos a realizarem este tipo de teste ensaiam apenas um corpo de prova por vez. Neste caso, o tratamento estatístico é baseado em uma pequena quantidade de resultados, reduzindo a confiabilidade na predição do comportamento em fluência. O Laboratório de Adesão e Aderência (LAA/UERJ/IPRJ) desenvolveu o equipamento pneumático de fluência (EPF) capaz de realizar dez ensaios simultâneos com parâmetros distintos. Este trabalho investiga o comportamento de juntas de cisalhamento simples (SLJ) compostas por substratos metálicos colados com adesivo epóxi e poliuretano. As juntas foram fabricadas em conformidade a norma ASTM D 1002, seguindo procedimento de preparo RT-003/08 TMC/CENPES/PETROBRAS. Os resultados mostram que, para projetos onde esteja prevista utilização das juntas durante longos períodos de tempo sob carregamento, conhecer apenas suas tensões médias de ruptura não é o suficiente para garantir a segurança do empreendimento. É proposto um modelo inicial do comportamento sob fluência de juntas coladas.

Estímulos fiscais e a interação entre as políticas monetária e fiscal no Brasil

Este trabalho estima, utilizando dados trimestrais de 1999 a 2011, o impacto dinâmico de um estímulo fiscal no Brasil sobre as principais variáveis macroeconômicas Brasileiras. Na estimativa dos impactos permitiu-se que as expectativas dos agentes econômicas fossem afetadas pela existência e probabilidade de alternância de regimes (foram detectados dois regimes) na política monetária do país. Os parâmetros da regra da política monetária, nos dois regimes detectados, foram estimados através de um modelo - composto apenas pela equação da regra da política monetária - que permite uma mudança de regime Markoviana. Os parâmetros do único regime encontrado para a política fiscal foram estimados por um modelo Vetorial de Correção de Erros (Vector Error Correction Model - VEC), composto apenas pelas variáveis pertencentes à regra da política fiscal. Os parâmetros estimados, para os diversos regimes das políticas monetária e fiscal, foram utilizados como auxiliares na calibragem de um modelo de equilíbrio geral estocástico dinâmico (MEGED), com mudanças de regime, com rigidez nominal de preços e concorrência monopolística (como em Davig e Leeper (2011)). Após a calibragem do MEGED os impactos dinâmicos de um estímulo fiscal foram obtidos através de uma rotina numérica (desenvolvida por Davig e Leeper (2006)) que permite obter o equilíbrio dinâmico do modelo resolvendo um sistema de equações de diferenças de primeira ordem expectacionais dinâmicas não lineares. Obtivemos que a política fiscal foi passiva durante todo o período analisado e que a política monetária foi sempre ativa, porém sendo em determinados momentos menos ativa. Em geral, em ambas as combinações de regimes, um choque não antecipado dos gastos do governo leva ao aumento do hiato do produto, aumento dos juros reais, redução do consumo privado e (em contradição com o resultado convencional) redução da taxa de inflação.

Caracterização geoquímica e efeito térmico de intrusão de diabásio na matéria orgânica da Formação Irati, poço FP-12-SP, bacia do Paraná

O estudo visou à caracterização geoquímica orgânica detalhada (carbono orgânico total, resíduo insolúvel, enxofre, pirólise Rock-Eval e biomarcadores) da Formação Irati através do furo pioneiro FP-12-SP localizado próximo da região de Anhembi (SP), e à análise das alterações destes parâmetros devido ao efeito térmico de intrusões de diabásio. A partir destes dados definiu-se 8 (oito) unidades quimioestratigráficas, determinou se o potencial gerador de cada unidade, o tipo de querogênio e a assinatura geoquímica através de parâmetros dos biomarcadores saturados e, por último, indicou-se as alterações causadas pelo efeito térmico de intrusivas nos sedimentos. O Membro Taquaral possui baixo teor de carbono orgânico total, querogênio tipo III/IV, com alta proporção de alcanos lineares e cicloalcanos, predominância dos esteranos C27 em relação ao C29, com indicações de ambiente marinho com salinidade normal. O Membro Assistência é caracterizado por apresentar valores de carbono orgânico total relativamente alto com aumento para o topo, intercalação de níveis siliciclásticos com carbonáticos, potencial gerador que chega a excelente e índice de hidrogênio predominante para óleo e condensado, baixa proporção de n-alcanos em relação aos cicloalcanos e alcanos ramificados, aparecimento do gamacerano, indicando ambiente marinho hipersalino com afogamento para o topo. O efeito térmico é visto a partir de 3 metros abaixo e 10 metros acima da soleira com maturação da matéria orgânica nos 9,5m mais próximos da soleira. A propagação de calor acima da soleira foi maior do que abaixo devido à diferença de condutividade térmica entre o carbonato e o folhelho, onde a máxima restrição com o aumento da salinidade na Unidade Quimioestratigráfica E, indicada pela presença de gamacerano, pode ter amplificado a passagem de calor para o topo.

Uma dedução dos critérios de multicriticalidade e o cálculo de pontos tricríticos via otimização global

Uma dedução dos critérios de multicriticalidade para o cálculo de pontos críticos de qualquer ordem representa a formalização de ideias utilizadas para calcular pontos críticos e tricríticos e ainda amplia tais ideias. De posse desta dedução pode-se compreender os critérios de tricriticalidade e, com isso, através de uma abordagem via problema de otimização global pode-se fazer o cálculo de pontos tricríticos utilizando um método numérico adequado de otimização global. Para evitar um excesso de custo computacional com rotinas numéricas utilizou-se aproximações na forma de diferenças finitas dos termos que compõem a função objetivo. Para simular a relação P v - T optou-se pela equação de estado cúbica de Peng-Robinson e pela regra clássica de fluidos de van der Vaals, para modelagem do problema também se calculou os tensores de ordem 2, 3, 4 e 5 da função do teste de estabilidade. Os resultados obtidos foram comparados com dados experimentais e por resultados obtidos com outros autores que utilizaram métodos numéricos, equação de estado ou abordagem diferente das utilizadas neste trabalho.

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.

Copolímeros à base de estireno e anidrido maleico: caracterização e estudo viscosimétrico

Copolímeros à base de estireno e anidrido maleico (SMA) são materiais sintéticos comercialmente disponíveis, obtidos pela reação dos dois monômeros citados em diferentes proporções, resultando em materiais versáteis, e disponíveis em diferentes graus de massas e porcentagens molares de anidrido maleico. São considerados polímeros funcionais devido à reatividade do grupamento anidrido maleico presente na cadeia polimérica. Por este motivo, esses materiais possuem vasta gama de aplicações, e elevada importância em diversas áreas, principalmente por terem baixa toxicidade, boa resistência térmica e boa estabilidade dimensional. Dessa forma, para melhor aplicação desses copolímeros, é muito importante o conhecimento dos parâmetros relativos ao seu comportamento em solução. A viscosimetria, em especial, é um método simples, útil e apropriado para fornecer essas informações. Os parâmetros viscosimétricos podem ser matematicamente calculados por extrapolação gráfica, entretanto a geração dos dados experimentais é mais demorada. Em contrapartida, é possível que a determinação experimental seja feita de forma mais rápida, por um único ponto, procedimento esse que desperta tanto o interesse acadêmico quanto o industrial. Neste trabalho, foram empregados os dois métodos de cálculo, utilizando solventes puros, misturas de solventes e três amostras de copolímeros à base de SMA. As determinações foram conduzidas a 40C. Os copolímeros utilizados possuiam teores de anidrido maleico de 50%, 45% e 40%, sendo os dois últimos esterificados com butil-metil-éster e sec-butil-metil-éster, respectivamente. Os solventes utilizados foram: N-metil-pirrolidona (NMP), tetrahidrofurano (THF) e suas respectivas misturas 1:1 com metil-etil-cetona (MEK), ou seja, (NMP:MEK) e THF:MEK, sendo a MEK um não solvente para o copolímero não esterificado. As equações utilizadas para extrapolação gráfica foram as de Huggins, Kraemer e Schulz-Blaschke. As equações empregadas em um único ponto foram as de Solomon-Ciuta, Deb-Chanterjee e novamente Schulz-Blaschke. Os resultados obtidos foram comparados e avaliou-se a possibilidade da utilização do método mais rápido, por um único ponto, para os sistemas estudados através dos desvios percentuais tendo como padrão os resultados da equação de Huggins. A equação de Deb-Chanterjee foi a mais adequada aos sistemas em NMP, que foi também o melhor solvente para as amostras. Os resultados obtidos na mistura NMP:MEK sugeriram que a viscosimetria pode ter sido um método sensível às pequenas diferenças estruturais entre os grupos pendentes nas amostras esterificadas. Paralelamente, realizou-se análises de espectroscopia na região do infravermelho (FTIR), análise termogravimétrica (TGA) e ensaios dinâmico-mecânicos (DMA) para a caracterização estrutural e térmica das amostras. Somente os resultados obtidos a partir de DMA indicaram diferenças entre as amostras esterificadas

Responsabilidade socioambiental da empresa de petróleo e gás

O contexto da modernização reflexiva leva ao questionamento sobre o papel das instituições tradicionais, notadamente o Poder Executivo. É possível pressupor que o modelo da sociedade de risco se reflete em alterações no direito e na economia do Brasil. Os riscos tornaram-se uma constante e exigem substituições nas formas de atuação social, o que Ulrich Beck denomina de subpolítica. A mudança é sentida especialmente em setores centrais para economia, tais como o de petróleo e gás, enquanto os riscos globais são sentidos na crise ambiental. Uma vez que as instituições, na Sociedade de Risco, são levadas a rever sua atuação, o empresariado recebe a tarefa de encontrar mecanismos para superar a crise ambiental. A responsabilidade socioambiental da empresa passa a ser exigida como contrapartida pelos lucros obtidos, especialmente para atividades potencialmente poluentes, como a petrolífera. O princípio da precaução, o desenvolvimento sustentável e a equação financeira do contrato podem ser vetores para a adoção da responsabilidade socioambiental pela indústria petrolífera. Mas para que esta possa ser vista como uma nova razão pública deve se demonstrar que ela pode motivar a evolução da sociedade como um todo. A Rio+20 definiu a economia verde como uma nova meta, principalmente para as atividades potencialmente poluentes. O objetivo central do trabalho é investigar a regulação das empresas de petróleo e gás, especialmente sobre a possibilidade de adoção da responsabilidade socioambiental. Ressalte-se que esta visa impor medidas de conservação e ações pró meio ambiente, além daquelas já estipuladas por força legal ou pelo licenciamento ambiental. A pesquisa visa apresentar possíveis soluções para os questionamentos acima, garantindo segurança jurídica para empresas de petróleo e gás, mas ao mesmo tempo visando ampliar a sustentabilidade do setor, propondo novas regras que podem ser adotadas nos editais de licitação e nos contratos de concessão e partilha de produção. Em um momento que se debate a possibilidade de direcionamento dos royalties do petróleo exclusivamente para a educação o estudo de medidas jurídicas para implementar a responsabilidade socioambiental no setor petrolífero torna-se ainda mais necessária.

Processos estocásticos em teoria de campos e aplicação ao universo inflacionário

É conhecido que derivações microscópicas obtidas através de métodos de teoria quântica de campos (TQC) podem conduzir a complicadas equações de movimento (EdM) que possuem um termo dissipativo com memória e um termo de ruído colorido. Um caso particularmente interessante é o modelo que escreve a interação entre um sistema e um banho térmico a temperatura T. Motivado por isso, usamos uma prescrição que nos permite reescrever EdMs não-markovianas semelhantes as obtidas em TQC em termos de um sistema de equações locais, para então confrontarmos a solução desse sistema com a solução aproximada usada correntemente na literatura, a chamada aproximação markoviana. A pergunta chave a qual se pretende responder aqui é: dado um conjunto de parâmetros que descrevem o modelo, a aproximação markoviana é suficientemente boa para descrever a dinâmica do sistema se comparada a dinâmica obtida atravéS da EdM não-markoviana? Além disso, consideramos uma versão linear da ELG de forma que pudéssemos determinar o nível de confiança da nossa metodologia numérica, procedimento este realizado comparando-se a solução analítica com a solução numérica. Como exemplo de aplicação prática do tema discutido aqui, comparamos a evolução não-markoviana do inflaton com a evolução markoviana do mesmo num modelo de universo primordial denominado inflação não-isentrópica (warm inflation).