Neste trabalho apresenta-se um modelo numérico abrangente para análise tridimensional de estruturas de concreto armado submetidas a cargas monótonas e cíclicas, utilizando o método dos elementos finitos. O modelo constitutivo proposto para o concreto é um modelo ortotrópico que segue a teoria de elasticidade não-linear. Neste modelo, aplica-se o conceito de deformação uniaxial equivalente para determinar o comportamento triaxial do concreto, através de três curvas tensão-deformação uniaxiais equivalentes, facilitando a simulação de seu comportamento e permitindo uma boa aproximação com os resultados experimentais. O modelo descreve bem a resposta do concreto solicitado por diversos tipos de carregamentos, sendo capaz de representar o esmagamento e a fissuração do concreto. A implementação em elementos finitos está baseada na consideração de fissuras distribuídas com as fissuras girando segundo as direções de tensões principais. Para a simulação de cargas cíclicas, considera-se a curva de Popovics-Saenz modificada, introduzindo-se algumas alterações nas coordenadas da origem e do pico desta curva tensão-deformação uniaxial equivalente, quando o concreto for submetido a processos de descarga e recarga. O comportamento do aço é descrito por um modelo constitutivo uniaxial. Para modelar este comportamento quando solicitado por cargas monótonas, considera-se uma curva tensãodeformação elástica bilinear com endurecimento. Para simulação da resposta cíclica do aço, utiliza-se um modelo constitutivo não-linear. Este modelo consegue representar bem o comportamento cíclico do aço, simulando a sua resposta não-linear, após alternância de carga. Os modelos constitutivos dos materiais foram implementados em um programa computacional, que utiliza o Método dos Elementos Finitos, para estruturas de concreto armado. Este programa em Elementos Finitos utiliza elementos isoparamétricos hexaédricos lineares e quadráticos para o concreto. Para representar as barras de armadura, utilizou-se um modelo incorporado. Por fim, para comprovar a eficiência do modelo numérico proposto, realizaram-se diversas simulações numéricas com concreto simples e armado. As comparações com resultados de ensaios demonstraram uma ótima concordância das previsões do modelo numérico com os valores experimentais. / In this work, a 3D finite element model for analyses of reinforced concrete structures subjected to monotonic and cyclic loads is presented. The constitutive model proposed for the concrete is an orthotropic model that follows the nonlinear elasticity theory. In this model, the concept of equivalent uniaxial strain is applied for determination of the triaxial behavior of the concrete through three equivalent uniaxial stress-strain curves, which makes it easier to simulate the concrete’s behavior and allows a good approximation with experimental results. The model describes well the response of concrete when subjected to several loading types, being capable of representing crushing and cracking of concrete. The implementation in finite elements is based upon the consideration of smeared cracks with them rotating according to principal stresses’ directions. For simulation of the concrete’s behavior under cyclic loads, the modified Popovics-Saenz curve has been considered, with some changes in the coordinates of the origin and of the peak of this equivalent uniaxial stress-strain curve, when the concrete is subjected to unloading and reloading processes. For simulation of the cyclic response of the steel, a nonlinear constitutive model is used. This model can well represent the cyclic behavior of the steel, simulating its nonlinear response after load reversals. The computational code uses linear and quadratic isoparametric hexahedral finite elements for the concrete. To represent the reinforcing bars, an embedded model is used. Finally, to validate the efficiency of the proposed model, several numerical simulations are carried out in plain and reinforced concrete. The comparisons made with test results show an excellent agreement between the model predictions and experimental data.

A conversão de energia solar diretamente em energia elétrica a partir do efeito fotovoltaico é uma forma confiável de geração de eletricidade. O crescimento das instalações fotovoltaicas na última década tem provado sua viabilidade quando acompanhada de programas de incentivo apropriados. As células e módulos fotovoltaicos possuem um comportamento elétrico característico e a determinação da curva característica I-V (corrente versus tensão) é essencial para a determinação da qualidade do dispositivo além de fornecer os parâmetros para dimensionamento e estimativa de geração de eletricidade. O presente trabalho apresenta o desenvolvimento de um sistema para traçar a curva característica I-V de módulos fotovoltaicos. O sistema possui uma fonte de quatro quadrantes para polarização do módulo, multímetros para medida simultânea de tensão e corrente formando pares I-V da curva característica, multímetro para medida de irradiância solar por meio de uma célula de referência de silício monocristalino e multímetro para medida de temperatura do módulo fotovoltaico utilizando um sensor Pt100. Todo o sistema é gerenciado por um programa de computador desenvolvido em linguagem Visual Basic®. O programa envia os comandos para os multímetros e fonte através de uma interface USB - GPIB. Os limites do sistema são ± 100 V e ± 10 A, determinados pela característica da fonte utilizada. A tensão de polarização é em forma de rampa com tempo de duração variável e também é ajustada através do programa. O sistema pode medir curvas com até 2000 pontos e o tempo de medida mínimo de cada ponto foi limitado em 50 μs. Durante a medida da curva característica é medida simultaneamente a irradiância solar. O sistema desenvolvido apresentou-se eficiente na medida da curva característica e a densidade de pontos adquirida permite estudos e determinação de parâmetros de maneira precisa. O método de disparo dos multímetros utilizado permite a garantia da simultaneidade da medida do par I-V e a utilização do sistema operacional Windows® permite flexibilidade e facilidade no tratamento dos dados. / The direct conversion of solar energy into electrical energy by means of the photovoltaic effect is a reliable way to generate electricity. The increasing number of photovoltaic installations in the last decade has proved its viability when stimulated by appropriate incentive programs. Photovoltaic cells and modules have peculiar electrical features and the measurement of the characteristic I-V (current - voltage) curve is essential to determine the device's quality. Besides it provides parameters for dimensioning and estimating electricity generation. This work presents the development of a system designed to trace the I-V curve of photovoltaic modules. The system has a four quadrant power supply to polarize the module, two multimeters for measuring I-V pairs of the characteristic curve, one multimeter to measure solar irradiance using a monocristaline silicon reference cell and one multimeter to measure the module temperature using a Pt100. The system is managed by a program developed in Visual Basic® which sends commands to the instruments via a USB - GPIB interface. The system limits are ± 100 V and ± 10 A, limited by the power supply range. The sweeping voltage is a ramp waveform and it is also programmable by the software. The system can measure 2000 points in each curve and the acquisition time of each point is limited to 50 μs. While tracing the characteristic curve, solar irradiance is simultaneously measured. The system is working efficiently and the number of measured points allows studies to obtain precise parameters. The triggering method assures simultaneously measurements of the I-V pair. The use of Windows® operational system allows flexibility and simplicity when dealing with the data in post processing software.

Os sistemas fotovoltaicos autônomos apresentam-se como uma alternativa adequada para a eletrificação rural. Entanto, existem ainda problemas que devem ser resolvidos, fundamentalmente associados ao projeto do sistema e à qualidade técnica dos equipamentos, instalações e os quais influem diretamente sobre a confiabilidade. Para determinar os fatores que incidem sobre a confiabilidade foram estudados os modos de operação e as falhas mais comuns nos sistemas fotovoltaicos autônomos. Pesquisas realizadas no Laboratório de Energia Solar da UFRGS, em conjunto com uma extensa revisão bibliográfica, mostram os problemas técnicos básicos que se produzem em cada um dos elementos da instalação e a dependência entre eles. Estes estudos, combinados com testes sobre baterias automotivas, mostraram que o sistema de acumulação é o elo mais fraco do sistema no aspecto confiabilidade e economia, devido à diminuição de sua capacidade de armazenamento. Isto motivou a considerar os elementos de armazenamento como o foco deste estudo e, através da análise do seu comportamento, apresentar um procedimento que permita dimensionar sistemas fotovoltaicos autônomos com elevada confiabilidade, configuração adequada e baixo custo. Esta Tese tem como objetivo realizar uma análise teórica e prática do impacto das baterias sobre a confiabilidade técnica e viabilidade econômica de sistemas fotovoltaicos autônomos. Para atingir este objetivo, além da realização de testes experimentais, desenvolvimento e ajustes de modelos matemáticos, também foi aperfeiçoada uma ferramenta computacional denominada PVSize. O aperfeiçoamento deste software visa permitir o cálculo de diferentes configurações de sistemas e determinar os riscos de déficit de energia sobre as mesmas, bem como as figuras de mérito econômicofinanceira associados ao projeto escolhido. Como parte deste estudo foi instalado um sistema fotovoltaico autônomo e foram montadas bancadas de ensaio de componentes do sistema. Medições realizadas nestes sistemas permitiram verificar os modelos matemáticos que descrevem o comportamento individual de cada dispositivo e caracterizar os elementos que compõem o sistema. Através de análises experimentais e o seguimento do comportamento do banco de baterias do sistema experimental ao longo de um ano, determinou-se que a conexão de baterias em paralelo acelera o processo de degradação das baterias, e esta degradação tem impacto diferenciado sobre a perda de capacidade de cada bateria. Como resultado deste estudo foi verificado o impacto que a perda de capacidade das baterias tem sobre a confiabilidade de um sistema fotovoltaico autônomo. Além disso, foi aperfeiçoada uma ferramenta de análise e modelagem computacional que determina o nível de confiabilidade de um sistema fotovoltaico autônomo associado à aparição de falhas e ao dimensionamento com elementos reais, à característica estocástica da radiação solar, ao perfil de consumo de energia elétrica e ao efeito de diminuição da capacidade do sistema de acumulação com o tempo. / Stand-alone photovoltaic systems are a suitable alternative for rural electrification. However, there are still problems to be solved, mainly related to the system design and the technical quality of the equipment and facilities, which have impact on the systems reliability. To determine the factors that affect the reliability of these systems were studied the most common configurations and associated failures. The Laboratory experimental research, together with an extensive literature review, show the basic technical problems that occur to each of the elements of the installation and the dependence between them. These studies have shown that the storage system, considering system reliability and economy, is the weakest element due to the decrease of their storage capacity. This fact led to consider the storage systems as the focus of this study and, through the analysis of their behavior, to develop a procedure to size systems with high reliability, lower cost and appropriate configuration. The goal of this Thesis is determinate the impact of batteries on the technical reliability and economic viability of photovoltaic systems. This goal was achieved through experimental testing and the development and adjustment of mathematical models. These models were implemented to preexisting software called PVSize. The improved software allows the calculation of different configurations of systems and to determine the loss of load probability and the figures of merit associated to the chosen economic-financial project. On this work was installed a photovoltaic system and was developed a battery testing system. The values measured in these systems allow to verify the mathematical models that describe the behavior of each device and characterize the components of the system. Experimental analysis of the behavior of a bank of batteries along a year showed that the connection of batteries in parallel accelerates the batteries degradation process, and this degradation has differentiated impact on the loss of capacity of each battery. This study allowed to evaluate the impact that the loss of capacity of batteries has on the reliability of an stand-alone photovoltaic system. It was also improved a software that allow to determine the level of technical reliability due to the design, the stochastic characteristic of solar radiation, the profile of electric power consumption and the effect of the accumulation loss capacity over the time.

Apesar do crescente número de publicações estudando as diferentes propriedades e o comportamento peculiar inerente dos materiais hiperelásticos submetidos a deformações, as quais muitas vezes propõem modelos matemáticos para caracterizá-los, muito pouco esforço é observado no sentido de conferir a este modelo matemático propriedades básicas como garantia de unicidade das soluções, estabilidade numérica ou mesmo comportamento físico condizente com a realidade. Este trabalho apresenta uma revisão das diferentes restrições matemáticas que devem ser impostas às equações da função energia de deformação dos modelos hiperelásticos para garantir-lhes estas propriedades básicas desejadas aplicadas à família de modelos de Hoss- Marczak [Hoss, 2009]. Como resultado, obtém um conjunto de desigualdades que devem ser satisfeitas pelas constantes constitutivas destes modelos. Posteriormente estes mesmos modelos são implementados em um programa comercial de elementos finitos e casos de deformações homogêneas e não-homogêneas são simulados para comparar os seus resultados com aqueles obtidos utilizando-se modelos consagrados pela literatura. Por fim, este trabalho apresenta uma metodologia de otimização multi-critério das constantes hiperelásticas baseada na minimização do erro entre dados teóricos e experimentais não somente de um ensaio, como usualmente é feito, e sim da soma dos erros de dois ou mais ensaios. As curvas dos ajustes teóricos assim obtidas são confrontadas com as experimentais e, através da comparação direta entre os coeficientes de correlação não-linear das curvas, concluise que a metodologia proposta fornece resultados superiores a outras metodologias. / In spite of the increasing in the number of publications proposing hyperelastic models to better characterize different types of rubber-like materials, the number of them concerned with aspects related to unicity of solution, numerical stability or plausibility of the physical behavior is much smaller. This work presents a literature review of the various types of mathematical restrictions that should be verified in strain energy density functions for hyperelastic constitutive models. These restrictions are verified for the HMI constitutive models [Hoss, 2009]. As a result, a set of inequalities that must be respected by the constitutive parameters of those models is obtained. In order to further assess the performance of the HMI models, their strain energy functions are implemented in a commercial finite element software. Different cases of homogeneous and non-homogeneous deformations are analyzed and the results are compared to those obtained when using some classical constitutive models. Finally, a multi-objective optimization procedure is implemented to obtain the best possible set of constitutive parameters for a general hyperelastic model. The method allows one to fit the chosen model to more than one set of experimental data. Theoretical predictions obtained using the proposed method are compared with experimental curves through non-linear fitting coefficients. The results show that the methodology generates better predictions than those obtained by adjusting the parameters against a single set of experimental data or combining multiple single-objective fittings.

O presente trabalho apresenta uma modelagem matemática simplificada para projeção do módulo de elasticidade de materiais de fricção. A modelagem foi desenvolvida a partir de equações apresentadas por outros autores. Estas equações foram complementadas para que fossem adequadas aos materiais de fricção. A partir da equação base, foi inserida a Lei da Mistura simples na equação a fim de possibilitar a consideração de todas as fases presentes no material de fricção e sua contribuição para a rigidez do mesmo. Foi proposto um fator de forma α para as matérias-primas particuladas empregadas na fabricação dos materiais de fricção. Este fator visou corrigir a eficiência da contribuição das fases particuladas para a rigidez do material de fricção em virtude da sua irregularidade geométrica. Da mesma forma, foi incluído, ainda, no modelo matemático, um fator exponencial que quantifica a contribuição de fases viscoelásticas de baixo módulo de elasticidade. Conforme observado, a introdução deste fator exponencial é imprescindível devido à magnitude peculiar da influência das fases de baixo módulo de elasticidade sobre a rigidez do material de fricção. Para correta alimentação do modelo matemático proposto, foi realizada a caracterização do fator de forma das fases particuladas utilizadas na confecção dos materiais de fricção que foram utilizados para teste do modelo. Da mesma forma, foi feita a caracterização do módulo de elasticidade da resina fenólica, fabricada sob os mesmos parâmetros de processo dos materiais de fricção. Foram produzidas oito composições representativas de materiais de fricção semelhantes aos utilizados para fabricação de lonas de freio, aplicadas em veículos pesados com sistema de freio a tambor. As mesmas foram testadas em flexão a 3 pontos, de acordo com a norma ASTM D790-92, 1992, método de teste I, segundo procedimento A. Os resultados experimentais obtidos para o módulo de elasticidade dos materiais de fricção mostraram bom grau de confiabilidade, evidenciada pelo desvio padrão obtido nos resultados. A comparação entre os resultados teóricos, calculados através do modelo proposto e a partir das propriedades das matérias-primas, e os resultados experimentais mostrou boa correlação, encorajando a utilização do modelo desenvolvido neste trabalho para o projeto de materiais de fricção. / This work shows a simplified mathematical modeling for friction materials elasticity modulus design. The modeling was developed from other authors equations. These equations were modified to be adequate to friction materials. In the base equation, was introduced a simple Rule-of-Mixtures in order to consider all components in the friction material composition, and their contribution to the rigidity of it. A shape factor (α) was proposed to be applied to the particulate raw materials used in the fabrication of friction material. The shape factor function is to correct the particulate phases contribution efficiency to the rigidity of the friction material, due to its geometrical irregularities. Also, was introduced in the mathematical modeling an exponential factor to quantify the contribution of viscoelastic phases that show low modulus of elasticity. It was observed that the introduction of the exponential factor is very important due to the peculiar magnitude of the influence of these phases with low modulus of elasticity on the friction material rigidity. To obtain the necessary support to the mathematical modeling, the shape factor characterization was proceeded for the particulate phases used in the manufacture of the friction materials tested to check modeling validation. Also, the modulus of elasticity characterization of the phenolic resin was made with pure phenolic resin pieces manufactured with the same process parameters used for the friction materials manufacture. Eight compositions, representative of friction materials used for heavy vehicles drum brake linings, were produced. They were tested under 3-points bending test procedure, according to ASTM D790-92, 1992 standard, using test method I, procedure A. Experimental results obtained for the modulus of elasticity of the friction materials show reliable confidence, evidenced by the standard deviation of the results. Comparison between theoretical results, calculated with the proposed mathematical modeling and the raw materials properties, and the experimental results show good correlation, encouraging the use of the proposed modeling for friction materials design.

Dado a larga aplicação dos modelos matemáticos, a procura pelo aprimoramento dos métodos de identificação de modelos sempre é bem vinda. Foi em cima disto que se desenvolveu o presente trabalho. Considerando a problemática de que os distúrbios desconhecidos nas entradas comprometem a qualidade da estimação dos parâmetros, se estudou a estimação combinada com o estimador de estados CEKF (Constrained Extended Kalmam Filter). Isto porque se supôs que com o emprego do CEKF o efeito do distúrbio desconhecido seria reduzido. Para tanto, se tomou um caso de estudo de seis tanques esféricos interligados e formam simulados dados experimentais computacionalmente. Foi escolhido um ponto de operação para a estimação com base na estimabilidade dos parâmetros. Antes da estimação os dados experimentais foram filtrados. Para a estimação dos parâmetros foi elaborada uma metodologia iterativa, caracterizada pelos três procedimentos: (1) estimação dos parâmetros, (2) ajuste de função objetivo: possibilidade de incorporar as derivadas medidas, bem como os estados e derivadas estimados pelo CEKF e (3) melhoria dos dados com a passagem do CEFK. Em outras palavras, a metodologia desenvolvida consiste em estimações de parâmetros seqüenciais acompanhadas com a melhoria progressiva dos dados devido à aplicação do CEKF. Para o estudo de caso tomado foi perceptível a melhoria na estimação dos parâmetros, tornando interessante a continuidade do estudo que foi iniciado no presente trabalho. / Given the wide application of mathematical models, the demand for improvement of methods of identification of models is always welcome. It was upon this that has developed in this work. Considering the problem of the unknown disturbances in the inputs affect the quality of the estimation of parameters, we studied the estimation combined with the estimator states CEKF (Extended Constrained Filter Kalmar). This is because it was assumed that with the use of CEKF the effect of unknown disturbance would be reduced. As such, it took a case study of six spherical tanks interconnected to form simulated experimental data computationally. Was chosen an operating point for the estimation based on the estimated parameters. Before estimating the experimental data were filtered. For the estimation has been elaborated iterative methodology, characterized by three procedures: (1) parameter estimation, (2) adjustment of the objective function: the possibility of incorporating the derived measures, and states and derive the estimated CEKF and (3 ) improvement of data with the passage of CEFK. In other words, the methodology consists of sequential estimation of parameters together with the gradual improvement of data due to application of CEKF. For the case study was taken a noticeable improvement in the estimation of parameters, makes it interesting to continue the study that was initiated in this work.

Este trabalho visa o uso da função de transferência, a qual relaciona distribuição de temperatura e fluxo de calor, na comparação no domínio freqüência, entre o modelo de difusão usual (parabólico) e um modelo ondulatório (hiperbólico) que inclue o efeito de propagação do calor não instantâneo, sendo avaliados os casos de meio semi-infinito e finito. Para o caso de meio semi-infinito, são determinadas as expressões para as características de amplitude e de fase, considerando tanto a abordagem parabólica quanto hiperbólica. É observada a relação entre estas duas abordagens, mostrando que a abordagem parabólica é uma caso particular da abordagem hiperbólica, podendo ser obtida através de um processo de limite envolvendo o tempo de relaxação r . Para o caso de meio finito, são determinadas as expressões para as caracteríısticas de amplitude de ambas as faces da placa unidimensional, considerando tanto a abordagem parabólica quanto hiperbólica. Estas expressões são transformadas para a forma adimensional quando então são deduzidas as expressões correspondentes das características de amplitude. Mais uma vez, todos os resultados para o caso parabólico podem ser determinados a partir dos resultados do caso hiperbólico, através de um processo de limite envolvendo o tempo de relaxação r São apresentados resultados numéricos referentes às características de amplitude, onde é apontada a existência de uma freqüência limite, acima da qual a diferença entre os dois modelos, do tipo parabólico ou hiperbólico, aumenta rapidamente. Também é apresentada uma forma alternativa de cálculo da distribuição de temperatura transiente que faz uso da função de transferência do sistema.

A modelagem matemática é uma ferramenta que apresenta diversas vantagens no estudo de processos industriais como, por exemplo, entender e desenvolver tecnologias, avaliar o impacto de variáveis nas propriedades do produto, estudos de redução de custos de produção e de impacto ambiental. O objetivo deste trabalho é desenvolver e validar um modelo matemático para um reator de bancada de polimerização em lama de eteno, enfatizando o desenvolvimento do modelo cinético. O modelo do reator contemplou: 1) modelo de troca térmica; 2) modelo termodinâmico; 3) modelo de transferência de massa gás-líquido e 4)modelo cinético. Para uma melhor predição do modelo foi realizada a estimação de alguns dos seus parâmetros conforme uma metodologia que compreendeu a análise de sensibilidade paramétrica e das variáveis de entrada do modelo e os efeitos de um planejamento de experimentos para a geração dos dados experimentais. A metodologia utilizada mostrou-se eficiente na avaliação do modelo quanto as suas características de predição, ajudando na identificação de possíveis falhas e evidenciando as vantagens de uma adequada metodologia experimental. Uma etapa determinante para o processo de estimação dos parâmetros é a escolha dos dados de saída a utilizar para a estimativa de determinado parâmetro. A modelagem do reator experimental mostrou-se satisfatória em todos os aspectos do processo (troca térmica, cinética, termodinâmica e transferência de massa), pois o modelo prediz com precisão as características encontradas nos experimentos. Assim, este pode ser utilizado para avaliar o mecanismo reacional envolvido de forma a aproximar-se das reais características do processo.

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.

Nas indústrias de processo, algumas variáveis como as composições de topo de colunas de destilação são quantificadas através de análises laboratoriais e ou de analisadores em linha. Este tipo de controle apresenta algumas limitações tais como custos de manutenção e o tempo de amostragem (em geral análises laboratoriais ocorrem de 1 a 3 vezes ao dia). Para melhoria destes métodos, as variáveis podem ser continuamente estimadas a partir de sua relação com as variáveis medidas diretamente no processo. Através de um algoritmo matemático, um analisador virtual de propriedades pode ser construído, sendo necessário para o seu desenvolvimento um modelo simplificado do processo. Este trabalho teve como objetivo a construção de dois analisadores virtuais para estimação das composições de topo de duas colunas fracionadoras de tolueno presentes em uma indústria petroquímica. Para tal, estudou-se as metodologias existentes para construção de analisadores voltados para aplicações em colunas de destilação. O desenvolvimento de analisadores virtuais tem como base três etapas: seleção de variáveis secundárias, construção de modelos e correção/adaptação de modelos. Tais etapas, baseadas principalmente em métodos estatísticos, foram estudadas e as técnicas que melhor se adaptaram ao caso em questão foram empregadas com resultados satisfatórios. Realizaram-se também estudos comparativos entre modelos estacionários e dinâmicos e modelos construídos a partir de dados de simulação e de processo. As simulações foram conduzidas nos softwares Aspen PlusÒ e Aspen DynamicsÒ e o software usado para implementação dos analisadores na indústria foi o Aspen IQÒ.