O presente trabalho foi desenvolvido para avaliar o potencial de utilização de um reator anaeróbico fabricado em PEAD Polietileno de Alta Densidade, produzido pelo processo de rotomoldagem em substituição aos reatores convencionais construídos em concreto e alvenaria, trabalhando em regime de batelada e enterrados no solo. Os estados de tensões e deformações foram avaliados utilizando o programa de Elementos Finitos ABAQUS versão 6.5 e a malha dos nós utilizando o programa MSC PATRAN 2005 formando 7329 nós e 2004 elementos, em uma malha otimizada para as regiões de maior curvatura (pontos concentradores de tensão). O carregamento é formado com uma pressão interna do biogás de 5 kPa acrescido da carga hidrostática de biomassa de 6000 kgf em uma fundação elástica calculada pela razão tensão/recalque a partir do Módulo de Elasticidade equivalente do solo (Esolo). Comparando o estado de tensões avaliado durante o carregamento foi possível constatar que a maior tensão obtida no elemento mais crítico para a utilização mais provável do reator atingiu o valor de 7,46 MPa (não supera 40% do menor valor de resistência à tração e ao cisalhamento do PEAD de 20 MPa) e a maior razão de deformação dR/R foi de 1.0%. O caso mais crítico avaliado foi quando o reator está enterrado, totalmente vazio, em solo com Esolo = 1,55 MPa e o material com EPEAD = 1550 MPa e com uma sobrecarga superficial no terreno de 20kN/m2 gerando uma tensão de 17,80 MPa no elemento 1955 (atingindo 89% do menor valor de resistência à tração e ao cisalhamento do PEAD igual a 20 MPa). Os resultados obtidos comprovam que o reator produzido em PEAD substitui com vantagens os modelos fabricados em concreto ou alvenaria, suportando a pressão interna do biogás e a carga de biomassa. / The present work was developed to evaluate the potential of uses of an anaerobic reactor manufactured in HDPE High Density Polyethylene produced by the rotomolding process in substitution to the conventional reactors built in stonemasonry, working in a batch regime and buried in the soil. The state of tensions and the deformations were assessed using the program of Finite Elements ABAQUS version 6.5 and the mesh of the knots using the program MSC PATRAN 2005 forming 7329 knots and 2004 elements, in an optimized mesh for the areas of larger curvature (tension concentrator points). The loading is formed with an internal pressure of the biogas of 5kPa added of biomass hydrostatic load of 6000 kg in an elastic foundation calculated by the ratio pressure/settling starting from the Module of equivalent Elasticity of the soil (Esolo). Comparing the state of tensions assessed during the loading was possible to verify that the largest tension obtained in the most critical element goes the most probable utilization of the reactor, reached the value of 7, 46 MPa (it doesn't surpass 40% of the smallest resistance value to the traction and to the shearing strain of HDPE of 20 MPa) and the largest ratio of dR/R deformation was of 1.0%. The most critical assessed case was when the reactor is buried in soil with Esolo = 1,55 MPa and material with EPEAD = 1550 MPa, totally empty and with a superficial overload in the land of 20kN/m2 generating a tension of 17,80 MPa in the element 1955 (reaching 89% of the smallest resistance value to the traction and the shearing strain of a 20 MPa HDPE). The obtained results confirmed that the reactor produced in HDPE substitutes with advantages the models manufactured in stonemasonry, supporting the internal biogas pressure and the biomass load.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O objetivo do presente trabalho é apresentar um estudo sobre o equilíbrio de fases sólido-líquido-vapor para hidratos de metano. A análise do equilíbrio trifásico sólido-líquido-vapor tem encontrado diversas aplicações para sistemas hidrocarboneto-água, uma vez que permite, por exemplo, a determinação da região de estabilidade de hidratos de metano e hidratos de gás natural. Inicialmente foi feita uma pesquisa sobre o estado da arte no que diz respeito ao comportamento termodinâmico e equilíbrio de fases de hidratos. Foram implementados os modelos apresentados por Ballard (2002) e Zhang et al. (2005). A proposta de Zhang et al. (2005) é aplicável para equilíbrios de fases a temperaturas abaixo de 300 K. Sua abordagem combinou a teoria de van der Waals e Platteeuw para a fase hidrato com a equação do estado de Peng-Robinson (1976) modificada por Stryjek e Vera (1986) para ambas as fases fluidas (fase vapor e fase aquosa). A abordagem de Ballard (2000) considerou a distorção do hidrato do seu estado padrão, o que fornece uma exata composição do hidrato e melhora a previsão da formação dos hidratos a altas pressões. Ao esclarecer a mudança de volume no hidrato, o raio da gaiola do hidrato é uma função do seu volume. Com isso, Ballard propôs uma nova abordagem considerando tal variação de volume e gerou um equilíbrio de fases em uma rotina de flash multifásico através da minimização da energia livre de Gibbs. Assim, o presente trabalho apresenta as abordagens de Zhang et al. (2005) e Ballard (2002) para o comportamento termodinâmico de hidratos e faz uma análise e comparação entre eles. Para resolver o problema do flash computacionalmente, foi utilizada a ferramenta lsqnonlin (built-in do software MATLAB). O lsqnonlin é um algoritmo baseado no método de Levenberg-Marquadt. / The objective of the present work is to present a study of solid-vapor-liquid three-phase equilibrium for methane hydrates. The analysis of three-phase equilibrium has several applications for water-hydrocarbon systems, since it permits, for example, determination of stability region for methane hydrates and natural gas hydrates. We have started seeking in literature about the state-of-art for thermodynamic behaviour and phase equilibrium for hydrates. And then the models proposed by Ballard (2002) and Zhang et al. (2005) were implemented. Zhang et al. (2005) have proposed a phase equilibrium for single-guest gas hydrates at temperatures below 300 K. Their approach has combined the van der WaalsPlatteeuw theory for the hydrate phase and the PengRobinson equation of state for both fluid phases (vapor and aqueous phase) (1976) modified by Stryjek and Vera (1986). Ballards (2000) approach has allowed the hydrate distortion from its standard state and has gave a more accurate composition of the hydrate and has improved hydrate formation predictions at high pressures. As a direct result of accounting for a changing hydrate volume, the cage radii were functions of the hydrate volume. Thus, Ballard have proposed the hydrate phase equilibrium by Gibbs energy minimization in a multi-phase flash routine. Thus, this work presents the Zhang et al. (2005) and Ballards (2002) approaches for hydrate thermodynamic behavior and makes an analysis and comparison of them. To compute the flash problem, we use the tool lsqnonlin (built-in of MATLAB software). The algorithm lsqnonlin is based on the Levenberg-Marquadt method.

Diversas aplicações industriais relevantes envolvem os processos de adsorção, citando como exemplos a purificação de produtos, separação de substâncias, controle de poluição e umidade entre outros. O interesse crescente pelos processos de purificação de biomoléculas deve-se principalmente ao desenvolvimento da biotecnologia e à demanda das indústrias farmacêutica e química por produtos com alto grau de pureza. O leito móvel simulado (LMS) é um processo cromatográfico contínuo que tem sido aplicado para simular o movimento do leito de adsorvente, de forma contracorrente ao movimento do líquido, através da troca periódica das posições das correntes de entrada e saída, sendo operado de forma contínua, sem prejuízo da pureza das correntes de saída. Esta consiste no extrato, rico no componente mais fortemente adsorvido, e no rafinado, rico no componente mais fracamente adsorvido, sendo o processo particularmente adequado a separações binárias. O objetivo desta tese é estudar e avaliar diferentes abordagens utilizando métodos estocásticos de otimização para o problema inverso dos fenômenos envolvidos no processo de separação em LMS. Foram utilizados modelos discretos com diferentes abordagens de transferência de massa, com a vantagem da utilização de um grande número de pratos teóricos em uma coluna de comprimento moderado, neste processo a separação cresce à medida que os solutos fluem através do leito, isto é, ao maior número de vezes que as moléculas interagem entre a fase móvel e a fase estacionária alcançando assim o equilíbrio. A modelagem e a simulação verificadas nestas abordagens permitiram a avaliação e a identificação das principais características de uma unidade de separação do LMS. A aplicação em estudo refere-se à simulação de processos de separação do Baclofen e da Cetamina. Estes compostos foram escolhidos por estarem bem caracterizados na literatura, estando disponíveis em estudos de cinética e de equilíbrio de adsorção nos resultados experimentais. De posse de resultados experimentais avaliou-se o comportamento do problema direto e inverso de uma unidade de separação LMS visando comparar os resultados obtidos com os experimentais, sempre se baseando em critérios de eficiência de separação entre as fases móvel e estacionária. Os métodos estudados foram o GA (Genetic Algorithm) e o PCA (Particle Collision Algorithm) e também foi feita uma hibridização entre o GA e o PCA. Como resultado desta tese analisouse e comparou-se os métodos de otimização em diferentes aspectos relacionados com o mecanismo cinético de transferência de massa por adsorção e dessorção entre as fases sólidas do adsorvente. / Several important industrial applications involving adsorption processes, citing as an example the product purification, separation of substances, pollution control and moisture among others. The growing interest in processes of purification of biomolecules is mainly due to the development of biotechnology and the demand of pharmaceutical and chemical products with high purity. The simulated moving bed (SMB) chromatography is a continuous process that has been applied to simulate the movement of the adsorbent bed, in a countercurrent to the movement of liquid through the periodic exchange of the positions of input and output currents, being operated so continuous, notwithstanding the purity of the outlet streams. This is the extract, rich in the more strongly adsorbed component, and the raffinate, rich in the more weakly adsorbed component, the method being particularly suited to binary separations. The aim of this thesis is to study and evaluate different approaches using stochastic optimization methods for the inverse problem of the phenomena involved in the separation process in LMS. We used discrete models with different approaches to mass transfer. With the benefit of using a large number of theoretical plates in a column of moderate length, in this process the separation increases as the solute flowing through the bed, i.e. as many times as molecules interact between the mobile phase and stationary phase thus achieving the equilibrium. The modeling and simulation verified in these approaches allowed the assessment and identification of the main characteristics of a separation unit LMS. The application under consideration refers to the simulation of the separation of Ketamine and Baclofen. These compounds were chosen because they are well characterized in the literature and are available in kinetic studies and equilibrium adsorption on experimental results. With the results of experiments evaluated the behavior of the direct and inverse problem of a separation unit LMS in order to compare these results, always based on the criteria of separation efficiency between the mobile and stationary phases. The methods studied were the GA (Genetic Algorithm) and PCA (Particle Collision Algorithm) and we also made a hybridization between the GA and PCA. This thesis, we analyzed and compared the optimization methods in different aspects of the kinetic mechanism for mass transfer between the adsorption and desorption of the adsorbent solid phases.

A fluência é um fenômeno de grande importância no concreto armado. Apesar de ser apontada como a causa de diversas patologias em estruturas, seus efeitos ainda são considerados de forma simplificada pelos projetistas estruturais. Além de estudar o fenômeno em estruturas de concreto armado e sua consideração atual utilizada na análise de estruturas, este trabalho compara as deformações de vigas simplesmente apoiadas de concreto armado, obtidas de forma analítica, e em ensaios experimentais com os resultados de simulações pelo Método dos Elementos Finitos (MEF). Os valores de deformações e deslocamentos calculados de forma analítica foram obtidos a partir de recomendações da Norma Brasileira NBR 6118 (2014) e do Código Modelo Simplificado CEB-FIP (1990), e os resultados experimentais foram extraídos de ensaios disponíveis na literatura. As simulações são realizadas no programa ANSYS Workbench, utilizando o elemento tridimensional SOLID 186 e a simetria das vigas. Posteriormente são realizadas análises de convergência utilizando o elemento bidimensional PLANE 183. Conclui-se ao final que as análises através do MEF são qualitativa e quantitativamente eficazes para estimativa dessa não linearidade e que o método de obtenção dos coeficientes de fluência utilizados é suficientemente preciso. / Reinforced concrete creep is a phenomenon of great importance. Despite being appointed as the main cause of several pathologies, its effects are yet considered in a simplified way by the structural designers. In addition to studying the phenomenon in reinforced concrete structures and its current account used in the structural analysis, this paper compares creep strains at simply supported reinforced concrete beams in analytical and in experimental forms with the finite element method (FEM) simulation results. The strains and deflections obtained through the analytical form were calculated with the Brazilian code NBR 6118 (2014) recommendations and the simplified method from CEB-FIP 90 and the experimental results were extracted from tests available in the literature. Finite element simulations are performed using ANSYS Workbench software, using its 3D SOLID 186 elements and the structure symmetry. Analyzes of convergence using 2D PLANE 183 elements are held as well. At the end, it is concluded that FEM analyses are quantitative and qualitative efficient for the estimation of this non-linearity and that the method utilized to obtain the creep coefficients values is sufficiently accurate.

Redes de sensores sem fio geralmente são compostas por dispositivos que demandam qualidade de transmissão, porém com energia limitada. Desta condição, surge a necessidade do estudo da eficiência energética. Neste documento são apresentadas e comparadas algumas técnicas para melhorar a eficiência energética que englobam o uso da diversidade, tais como o uso de múltiplas antenas no transmissor e receptor bem como o uso de comunicação cooperativa, esta última melhor aproveitada quando um canal de retorno está disponível. São abordados dois modelos de comunicação. No primeiro cenário, apenas a camada física é levada em conta em uma rede MIMO cooperativa com múltiplos relays, em que é proposto um algoritmo distribuído iterativo de alocação de potência baseado no algoritmo de Dinkelbach visando a maximização da eficiência energética. Os resultados demonstram que o uso de múltiplas antenas melhora a eficiência energética, mas que dependendo da técnica MIMO utilizada, o número de relays e o número de antenas podem levar a efeitos diferentes em termos de eficiência energética. Em seguida, é apresentado um cenário que envolve a camada física e a camada de controle de acesso ao meio (MAC), no qual propõe-se um protocolo MAC cooperativo distribuído. Neste caso, além da energia utilizada para transmitir os pacotes de dados, também é analisado o atraso e o consumo para acessar o canal. Pelos resultados obtidos, notase que as técnicas MIMO aumentam a eficiência energética em relação ao caso SISO, das quais pode-se destacar a seleção de antenas. Além disso, outros parâmetros, tais como densidade de nós, eficiência espectral, número de antenas e distância entre nós, são analisados em termos de eficiência energética. / Wireless sensor networks are usually composed of devices that require transmission quality, but with limited energy. From this condition arises the necessity of the energy efficiency study. In this document, we present and compare some techniques to improve the energy efficiency, which include the use of diversity, such as the use of multiple antennas in the transmitter and receiver as well as the use of cooperative communication. The latter has better results when a feedback channel is available. We cover two communication scenarios. In the first scenario, only the physical layer is taken into account in a MIMO cooperative network with multiple relays, where we propose a power allocation distributed iterative algorithm based on the Dinkelbach’s algorithm aiming at the energy efficiency maximization. The results show that the use of multiple antennas improve the energy efficiency, but depending on the used MIMO transmission scheme, the number of relays and the number of antennas can lead to different effects with respect to the energy efficiency. In the following, we present a scenario that contains a cross-layer analysis encompassing physical and medim access control (MAC) layers in a cooperative MIMO network, where we propose a distributed cooperative MAC protocol. In this case, besides the energy used to transmit thte data packets, we also analyze the delay and energy consumption to access the channel. Results show that the MIMO techniques improve the energy efficiency when we compare to SISO, especially when employing antenna selection techniques. In addition, other parameters, such as node density, spectral efficiency, number of antennas and distance between nodes, are analyzed with respect to the energy efficiency.

CAPES / O estudo de um componente estrutural, é mais realístico quando se admite que o componente já possua trincas. A área que estuda este fenômeno é a mecânica da fratura. O componente que possui trinca e é submetido a esforços cíclicos, tende a falhar por fadiga. Este estudo apresenta cotas que “envelopam” a solução numérica aproximada da evolução da trinca. São estimados momentos estatísticos das cotas superior e inferior, afim de se obter resultados mais realísticos com relação a propagação da trinca, visto a existência de incerteza sobre os parâmetros dos modelos de evolução da trinca. As cotas são determinadas via metodologia Fast Crack Bounds, sendo está comparada com a solução numérica aproximada obtida pelo método de Runge-Kutta de quarta ordem. A randomização dos parâmetros do modelo, é executada através de Simulação de Monte Carlo. Para a quantificação da incerteza, da cota superior, inferior e da solução numérica, são considerados exemplos “clássicos” da mecânica da fratura, onde a função de correção do fator de intensidade de tensão é conhecida: placa com largura infinita, placa com largura finita e trinca central e placa com largura finita e trinca na aresta. O trabalho apresenta os desvios relativos do primeiro e segundo momento estatístico, bem como os ganhos computacionais na resolução do problema de valor inicial que descrevem a propagação da trinca. Em todos os casos analisados, a metodologia Fast Crack Bounds apresentou menor tempo computacional, quando comparada à solução numérica do problema, sendo no mínimo 411,23% mais eficaz para o parâmetro a0 , até 8.296,29% para o parâmetro KC . / The study of a structural component is more realistic when it is admitted that the component already has cracks. The area that studies this phenomenon is the fracture mechanics. The component which is cracked and subjected to cyclic stresses tends to fail due to fatigue. This study presents upper and lower bounds that "envelop" the approximate numerical solution of the evolution of the crack. The statistical moments of the upper and lower bounds are estimated, to obtain more realistic results in relation to the crack propagation, considering the existence of uncertainty about the parameters of the evolution models of the crack. Upper and lower bounds are determined using the Fast Crack Bounds methodology, being compared to the approximate numerical solution obtained by the fourth-order RungeKutta method. The randomization of the model parameters and execution through the Monte Carlo Simulation. For the quantification of the uncertainty, the upper and lower bounds and the numerical solution, "classic" examples of fracture mechanics are considered, where the correction function of the tensile strength factor is known: Infinite width plate, finite width plate a centered crack and finite width plate a bordercracked. The work presents the relative deviations of the first and second statistical moments, as well as the computational gains in solving the initial value problem that describe the propagation of the crack. In all cases analyzed, the Fast Crack Bounds methodology presented lower computational time when compared to the numerical solution of the problem, being at least 411.23% more effective for the parameter a0 , up to 8,296.29% for the parameter KC .

A fluência é um fenômeno de grande importância no concreto armado. Apesar de ser apontada como a causa de diversas patologias em estruturas, seus efeitos ainda são considerados de forma simplificada pelos projetistas estruturais. Além de estudar o fenômeno em estruturas de concreto armado e sua consideração atual utilizada na análise de estruturas, este trabalho compara as deformações de vigas simplesmente apoiadas de concreto armado, obtidas de forma analítica, e em ensaios experimentais com os resultados de simulações pelo Método dos Elementos Finitos (MEF). Os valores de deformações e deslocamentos calculados de forma analítica foram obtidos a partir de recomendações da Norma Brasileira NBR 6118 (2014) e do Código Modelo Simplificado CEB-FIP (1990), e os resultados experimentais foram extraídos de ensaios disponíveis na literatura. As simulações são realizadas no programa ANSYS Workbench, utilizando o elemento tridimensional SOLID 186 e a simetria das vigas. Posteriormente são realizadas análises de convergência utilizando o elemento bidimensional PLANE 183. Conclui-se ao final que as análises através do MEF são qualitativa e quantitativamente eficazes para estimativa dessa não linearidade e que o método de obtenção dos coeficientes de fluência utilizados é suficientemente preciso. / Reinforced concrete creep is a phenomenon of great importance. Despite being appointed as the main cause of several pathologies, its effects are yet considered in a simplified way by the structural designers. In addition to studying the phenomenon in reinforced concrete structures and its current account used in the structural analysis, this paper compares creep strains at simply supported reinforced concrete beams in analytical and in experimental forms with the finite element method (FEM) simulation results. The strains and deflections obtained through the analytical form were calculated with the Brazilian code NBR 6118 (2014) recommendations and the simplified method from CEB-FIP 90 and the experimental results were extracted from tests available in the literature. Finite element simulations are performed using ANSYS Workbench software, using its 3D SOLID 186 elements and the structure symmetry. Analyzes of convergence using 2D PLANE 183 elements are held as well. At the end, it is concluded that FEM analyses are quantitative and qualitative efficient for the estimation of this non-linearity and that the method utilized to obtain the creep coefficients values is sufficiently accurate.

Redes de sensores sem fio geralmente são compostas por dispositivos que demandam qualidade de transmissão, porém com energia limitada. Desta condição, surge a necessidade do estudo da eficiência energética. Neste documento são apresentadas e comparadas algumas técnicas para melhorar a eficiência energética que englobam o uso da diversidade, tais como o uso de múltiplas antenas no transmissor e receptor bem como o uso de comunicação cooperativa, esta última melhor aproveitada quando um canal de retorno está disponível. São abordados dois modelos de comunicação. No primeiro cenário, apenas a camada física é levada em conta em uma rede MIMO cooperativa com múltiplos relays, em que é proposto um algoritmo distribuído iterativo de alocação de potência baseado no algoritmo de Dinkelbach visando a maximização da eficiência energética. Os resultados demonstram que o uso de múltiplas antenas melhora a eficiência energética, mas que dependendo da técnica MIMO utilizada, o número de relays e o número de antenas podem levar a efeitos diferentes em termos de eficiência energética. Em seguida, é apresentado um cenário que envolve a camada física e a camada de controle de acesso ao meio (MAC), no qual propõe-se um protocolo MAC cooperativo distribuído. Neste caso, além da energia utilizada para transmitir os pacotes de dados, também é analisado o atraso e o consumo para acessar o canal. Pelos resultados obtidos, notase que as técnicas MIMO aumentam a eficiência energética em relação ao caso SISO, das quais pode-se destacar a seleção de antenas. Além disso, outros parâmetros, tais como densidade de nós, eficiência espectral, número de antenas e distância entre nós, são analisados em termos de eficiência energética. / Wireless sensor networks are usually composed of devices that require transmission quality, but with limited energy. From this condition arises the necessity of the energy efficiency study. In this document, we present and compare some techniques to improve the energy efficiency, which include the use of diversity, such as the use of multiple antennas in the transmitter and receiver as well as the use of cooperative communication. The latter has better results when a feedback channel is available. We cover two communication scenarios. In the first scenario, only the physical layer is taken into account in a MIMO cooperative network with multiple relays, where we propose a power allocation distributed iterative algorithm based on the Dinkelbach’s algorithm aiming at the energy efficiency maximization. The results show that the use of multiple antennas improve the energy efficiency, but depending on the used MIMO transmission scheme, the number of relays and the number of antennas can lead to different effects with respect to the energy efficiency. In the following, we present a scenario that contains a cross-layer analysis encompassing physical and medim access control (MAC) layers in a cooperative MIMO network, where we propose a distributed cooperative MAC protocol. In this case, besides the energy used to transmit thte data packets, we also analyze the delay and energy consumption to access the channel. Results show that the MIMO techniques improve the energy efficiency when we compare to SISO, especially when employing antenna selection techniques. In addition, other parameters, such as node density, spectral efficiency, number of antennas and distance between nodes, are analyzed with respect to the energy efficiency.

O presente trabalho foi desenvolvido para avaliar o potencial de utilização de um reator anaeróbico fabricado em PEAD Polietileno de Alta Densidade, produzido pelo processo de rotomoldagem em substituição aos reatores convencionais construídos em concreto e alvenaria, trabalhando em regime de batelada e enterrados no solo. Os estados de tensões e deformações foram avaliados utilizando o programa de Elementos Finitos ABAQUS versão 6.5 e a malha dos nós utilizando o programa MSC PATRAN 2005 formando 7329 nós e 2004 elementos, em uma malha otimizada para as regiões de maior curvatura (pontos concentradores de tensão). O carregamento é formado com uma pressão interna do biogás de 5 kPa acrescido da carga hidrostática de biomassa de 6000 kgf em uma fundação elástica calculada pela razão tensão/recalque a partir do Módulo de Elasticidade equivalente do solo (Esolo). Comparando o estado de tensões avaliado durante o carregamento foi possível constatar que a maior tensão obtida no elemento mais crítico para a utilização mais provável do reator atingiu o valor de 7,46 MPa (não supera 40% do menor valor de resistência à tração e ao cisalhamento do PEAD de 20 MPa) e a maior razão de deformação dR/R foi de 1.0%. O caso mais crítico avaliado foi quando o reator está enterrado, totalmente vazio, em solo com Esolo = 1,55 MPa e o material com EPEAD = 1550 MPa e com uma sobrecarga superficial no terreno de 20kN/m2 gerando uma tensão de 17,80 MPa no elemento 1955 (atingindo 89% do menor valor de resistência à tração e ao cisalhamento do PEAD igual a 20 MPa). Os resultados obtidos comprovam que o reator produzido em PEAD substitui com vantagens os modelos fabricados em concreto ou alvenaria, suportando a pressão interna do biogás e a carga de biomassa. / The present work was developed to evaluate the potential of uses of an anaerobic reactor manufactured in HDPE High Density Polyethylene produced by the rotomolding process in substitution to the conventional reactors built in stonemasonry, working in a batch regime and buried in the soil. The state of tensions and the deformations were assessed using the program of Finite Elements ABAQUS version 6.5 and the mesh of the knots using the program MSC PATRAN 2005 forming 7329 knots and 2004 elements, in an optimized mesh for the areas of larger curvature (tension concentrator points). The loading is formed with an internal pressure of the biogas of 5kPa added of biomass hydrostatic load of 6000 kg in an elastic foundation calculated by the ratio pressure/settling starting from the Module of equivalent Elasticity of the soil (Esolo). Comparing the state of tensions assessed during the loading was possible to verify that the largest tension obtained in the most critical element goes the most probable utilization of the reactor, reached the value of 7, 46 MPa (it doesn't surpass 40% of the smallest resistance value to the traction and to the shearing strain of HDPE of 20 MPa) and the largest ratio of dR/R deformation was of 1.0%. The most critical assessed case was when the reactor is buried in soil with Esolo = 1,55 MPa and material with EPEAD = 1550 MPa, totally empty and with a superficial overload in the land of 20kN/m2 generating a tension of 17,80 MPa in the element 1955 (reaching 89% of the smallest resistance value to the traction and the shearing strain of a 20 MPa HDPE). The obtained results confirmed that the reactor produced in HDPE substitutes with advantages the models manufactured in stonemasonry, supporting the internal biogas pressure and the biomass load.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O objetivo do presente trabalho é apresentar um estudo sobre o equilíbrio de fases sólido-líquido-vapor para hidratos de metano. A análise do equilíbrio trifásico sólido-líquido-vapor tem encontrado diversas aplicações para sistemas hidrocarboneto-água, uma vez que permite, por exemplo, a determinação da região de estabilidade de hidratos de metano e hidratos de gás natural. Inicialmente foi feita uma pesquisa sobre o estado da arte no que diz respeito ao comportamento termodinâmico e equilíbrio de fases de hidratos. Foram implementados os modelos apresentados por Ballard (2002) e Zhang et al. (2005). A proposta de Zhang et al. (2005) é aplicável para equilíbrios de fases a temperaturas abaixo de 300 K. Sua abordagem combinou a teoria de van der Waals e Platteeuw para a fase hidrato com a equação do estado de Peng-Robinson (1976) modificada por Stryjek e Vera (1986) para ambas as fases fluidas (fase vapor e fase aquosa). A abordagem de Ballard (2000) considerou a distorção do hidrato do seu estado padrão, o que fornece uma exata composição do hidrato e melhora a previsão da formação dos hidratos a altas pressões. Ao esclarecer a mudança de volume no hidrato, o raio da gaiola do hidrato é uma função do seu volume. Com isso, Ballard propôs uma nova abordagem considerando tal variação de volume e gerou um equilíbrio de fases em uma rotina de flash multifásico através da minimização da energia livre de Gibbs. Assim, o presente trabalho apresenta as abordagens de Zhang et al. (2005) e Ballard (2002) para o comportamento termodinâmico de hidratos e faz uma análise e comparação entre eles. Para resolver o problema do flash computacionalmente, foi utilizada a ferramenta lsqnonlin (built-in do software MATLAB). O lsqnonlin é um algoritmo baseado no método de Levenberg-Marquadt. / The objective of the present work is to present a study of solid-vapor-liquid three-phase equilibrium for methane hydrates. The analysis of three-phase equilibrium has several applications for water-hydrocarbon systems, since it permits, for example, determination of stability region for methane hydrates and natural gas hydrates. We have started seeking in literature about the state-of-art for thermodynamic behaviour and phase equilibrium for hydrates. And then the models proposed by Ballard (2002) and Zhang et al. (2005) were implemented. Zhang et al. (2005) have proposed a phase equilibrium for single-guest gas hydrates at temperatures below 300 K. Their approach has combined the van der WaalsPlatteeuw theory for the hydrate phase and the PengRobinson equation of state for both fluid phases (vapor and aqueous phase) (1976) modified by Stryjek and Vera (1986). Ballards (2000) approach has allowed the hydrate distortion from its standard state and has gave a more accurate composition of the hydrate and has improved hydrate formation predictions at high pressures. As a direct result of accounting for a changing hydrate volume, the cage radii were functions of the hydrate volume. Thus, Ballard have proposed the hydrate phase equilibrium by Gibbs energy minimization in a multi-phase flash routine. Thus, this work presents the Zhang et al. (2005) and Ballards (2002) approaches for hydrate thermodynamic behavior and makes an analysis and comparison of them. To compute the flash problem, we use the tool lsqnonlin (built-in of MATLAB software). The algorithm lsqnonlin is based on the Levenberg-Marquadt method.