Aproximação numérica – analítica para a modelagem da conversão termoquímica de combustíveis sólidos

Um algoritmo numérico foi criado para apresentar a solução da conversão termoquímica de um combustível sólido. O mesmo foi criado de forma a ser flexível e dependente do mecanismo de reação a ser representado. Para tanto, um sistema das equações características desse tipo de problema foi resolvido através de um método iterativo unido a matemática simbólica. Em função de não linearidades nas equações e por se tratar de pequenas partículas, será aplicado o método de Newton para reduzir o sistema de equações diferenciais parciais (EDP’s) para um sistema de equações diferenciais ordinárias (EDO’s). Tal processo redução é baseado na união desse método iterativo à diferenciação numérica, pois consegue incorporar nas EDO’s resultantes funções analíticas. O modelo reduzido será solucionado numericamente usando-se a técnica do gradiente bi-conjugado (BCG). Tal modelo promete ter taxa de convergência alta, se utilizando de um número baixo de iterações, além de apresentar alta velocidade na apresentação das soluções do novo sistema linear gerado. Além disso, o algoritmo se mostra independente do tamanho da malha constituidora. Para a validação, a massa normalizada será calculada e comparada com valores experimentais de termogravimetria encontrados na literatura, , e um teste com um mecanismo simplificado de reação será realizado.

Variação radial da retratibilidade e densidade básica da madeira de Eucalyptus saligna Sm.

A ultra-estrutura e a composição química da madeira, bem como suas propriedades físicas e mecânicas, variam significativamente entre espécies, entre árvores de uma mesma espécie e, mesmo, entre diferentes partes de uma mesma árvore. Com este trabalho objetivou-se o estudo dos parâmetros de retratibilidade e de densidade básica da madeira Eucalyptus saligna, com idade de 16 anos, proveniente de talhões experimentais da EMBRAPA Florestas, de Colombo, Paraná. As amostras foram retiradas à altura do DAP de quatro posições eqüidistantes a partir da medula em direção à periferia, correspondendo a 0, 33, 66 e 100%, com dimensões nominais de 1,0 x 2,0 x 3,0 cm, sendo a última dimensão no sentido longitudinal. Elas foram mantidas em câmara fechada com ventilação, próximo de soluções salinas supersaturadas, com o objetivo de proporcionar diferentes condições de umidade relativa. Uma vez atingidas as distintas condições de umidade de equilíbrio, as amostras foram secas em estufa a 105 ºC e obtidos os dados de retratibilidade e densidade básica da madeira nas posições mencionadas. Constataram-se valores de contração volumétrica mais baixos na região medular, apresentando um acréscimo para as demais posições. Comportamento semelhante foi observado para os coeficientes das contrações lineares nas direções tangencial e radial. O fator anisotrópico foi consideravelmente mais elevado na região medular, decrescendo substancialmente em direção ao alburno. A densidade básica não mostrou sinais efetivos de estabilidade, apesar de mostrar tendência de aumento em direção à periferia do tronco.

Estudo da corrosão dos aços AISI 1010 e AISI 316L em soluções com elevado teor de cloreto e em ácido naftênico

A mudança no mercado global do petróleo nos últimos anos, com o declínio das reservas de óleo leve, têm forçado a busca por novos campos petrolíferos em ambientes mais remotos, como nos campos localizados na camada pré-sal, e a exploração de óleos pesados que possuem elevado teor de ácidos naftênicos. Isso acarreta em grandes desafios para a previsão do desempenho de materiais frente às novas condições ambientais em que estão inseridos. No presente trabalho, o comportamento da corrosão do aço carbono AISI 1010 e do aço inoxidável AISI 316L foi estudado em soluções aquosas com elevado teor de cloreto e em solução de ácido naftênico ciclopentanóico a fim de ter melhor entendimento da ação dessas espécies no processo de corrosão e simular a corrosão pela água de produção na indústria petrolífera. Foram aplicadas as técnicas de potencial de circuito aberto, polarização potenciodinâmica, voltametria cíclica, espectroscopia de impedância eletroquímica, espectroscopia Raman, microscopia eletrônica de varredura e microscopia de força atômica, usadas, em cada caso, de acordo com a conveniência. O aumento da [Cl-] na faixa de 1,2–2,8 mol.L-1 não altera os processos catódicos e anódicos perto do Ecorr para os aços AISI 1010 e AISI 316L. Em condições de sobrepotenciais afastados do Ecorr, o aumento de [Cl-] aumenta os processos oxidativos de corrosão, o que é expresso pelas maiores densidades de corrente e carga anódica e aumento da perda de massa sofrida pelos eletrodos de ambos os aços. Portanto, os danos da corrosão são mais intensos quando se aumenta a [Cl-]. O aço AISI 1010 é ativo nas soluções de NaCl e a corrosão se propaga livremente de forma uniforme. Para o aço AISI 316L, uma ampla faixa de passividade pode ser vista nas soluções de NaCl; no Epit ocorre a ruptura do filme passivo e o crescimento de pites estáveis. Após 24 h de imersão em soluções de sulfato de sódio (branco) e de ácido naftênico ciclopentanóico ocorre crescimento de filme de óxido e as fases α-Fe2O3, Fe3O4 e δ-FeO(OH) foram identificadas nos espécimes de aço AISI 1010 e Fe3O4 foi identificado nos defeitos do filme prévio presente na superfície do aço AISI 316L. Os filmes formados em solução de ácido ciclopentanóico possuem menor resistência à polarização, maior rugosidade e maior taxa de corrosão quando comparado aos filmes crescidos na solução branco, para ambos os aços. A presença do ácido naftênico muda a forma como a reação de corrosão se procede e contribui para o aumento da corrosão. A corrosão naftênica foi mais pronunciada no aço carbono porque a presença dos elementos de liga no aço inox reduzem o número de sítios ativos ricos em Fe e tornam menos oportuna a ligação do Fe com o naftenato.

Comparação de desempenho entre a formulação direta do Método dos Elementos de Contorno com Funções Radiais e o Método dos Elementos Finitos em problemas de Poisson e Helmholtz

O presente trabalho objetiva avaliar o desempenho do MECID (Método dos Elementos de Contorno com Interpolação Direta) para resolver o termo integral referente à inércia na Equação de Helmholtz e, deste modo, permitir a modelagem do Problema de Autovalor assim como calcular as frequências naturais, comparando-o com os resultados obtidos pelo MEF (Método dos Elementos Finitos), gerado pela Formulação Clássica de Galerkin. Em primeira instância, serão abordados alguns problemas governados pela equação de Poisson, possibilitando iniciar a comparação de desempenho entre os métodos numéricos aqui abordados. Os problemas resolvidos se aplicam em diferentes e importantes áreas da engenharia, como na transmissão de calor, no eletromagnetismo e em problemas elásticos particulares. Em termos numéricos, sabe-se das dificuldades existentes na aproximação precisa de distribuições mais complexas de cargas, fontes ou sorvedouros no interior do domínio para qualquer técnica de contorno. No entanto, este trabalho mostra que, apesar de tais dificuldades, o desempenho do Método dos Elementos de Contorno é superior, tanto no cálculo da variável básica, quanto na sua derivada. Para tanto, são resolvidos problemas bidimensionais referentes a membranas elásticas, esforços em barras devido ao peso próprio e problemas de determinação de frequências naturais em problemas acústicos em domínios fechados, dentre outros apresentados, utilizando malhas com diferentes graus de refinamento, além de elementos lineares com funções de bases radiais para o MECID e funções base de interpolação polinomial de grau (um) para o MEF. São geradas curvas de desempenho através do cálculo do erro médio percentual para cada malha, demonstrando a convergência e a precisão de cada método. Os resultados também são comparados com as soluções analíticas, quando disponíveis, para cada exemplo resolvido neste trabalho.