O potencial dos sistemas de aquecimento de água por energia solar no Brasil é excelente devido às características climáticas favoráveis do país. A performance destes sistemas é bastante influenciada também pelos materiais empregados e dimensões de seus equipamentos e componentes. Diante do aquecimento global, a energia solar ganha maior atenção sendo umas das energias renováveis a ser empregada em grande escala para substituir alguns tipos de energias poluentes. Neste trabalho, apresenta-se o aperfeiçoamento de um programa computacional para simulação de sistemas térmicos de aquecimento de água utilizando energia solar em regime de termossifão ou circulação forçada, programa este que teve seu desenvolvimento iniciado no âmbito do Laboratório de Energia Solar da UFRGS, chamado TermoSim, sendo esta sua versão 3.0. A partir desta versão este software é capaz de realizar simulações com seis diferentes possibilidades de arranjos combinados com o tipo de energia auxiliar utilizada: sistemas com coletores solares com energia auxiliar a gás de passagem, elétrico de passagem, elétrico interno ao tanque, elétrico em série na linha de consumo, a gás em série na linha de consumo, e sem energia auxiliar. O programa é uma ferramenta de auxílio para estudos e análises de sistemas de aquecimento solar, com uma interface amigável e de fácil entendimento e utilização dos resultados. Além disso, esta nova versão permite também realizar simulações levando em conta as perdas de calor à noite, situação em que a circulação reversa pode ocorrer resultando na redução da eficiência, dependendo do sistema simulado. Foram realizadas diversas simulações com os modelos matemáticos empregados e feitas comparações entre os sistemas utilizando dados climáticos da cidade de Caxias do Sul, no Rio Grande do Sul, determinando-se o sistema com a configuração mais eficiente para um dado perfil de consumo de água. Finaliza-se o trabalho com análises econômicas simples com o objetivo de prever o tempo de retorno do investimento, levando-se em conta o preço para tarifa de energia elétrica e as possíveis economias mensais com um sistema utilizando a energia solar. / The potential of solar water heating systems through solar energy in Brazil is excellent due to the climatic features of the country. The performance of these systems is highly influenced also by the materials used to build it and by the dimension of its equipments and components. In face of global warming, solar energy gains more attention, since it´s one of the renewable energy that will be largely used to replace some of the existing polluting types of energy. This paper presents the improvement of a software that conducts simulations of water heating systems using solar energy in thermosiphon regime or forced circulation. TermoSim, as it is called, was initiated at the Solar Labs, and is in its version 3.0. The current version is capable of simulating 6 different arrangements´ possibilities combined with auxiliary energy: systems with solar collectors with auxiliary energy with gas, electric energy, internal electric energy, electric energy in series with the consumption line, and no auxiliary energy. The software is a tool to aid studies and analysis of solar heating systems, it has a friendly interface that is easy to comprehend and results are simple to use. Besides that, this version also allows simulations that consider heat losses at night, situation in which a reverse circulation can occur and mean efficiency loss, depending on the simulated system type. There were many simulations with the mathematical models used and comparisons were maid with the climatic data of the city of Caxias do Sul, in Rio Grande do Sul state, in Brazil. This way we were able to determine the system with the most efficient configuration for the simulated water consume profile. The paper is finalized with simple economic analyses with the intention of foreseeing the time for payback on the investment, taking into account the current prices for electrical energy in the simulated area and the possible monthly economy provided with the use of a solar energy heating system.

O trabalho proposto apresenta um novo método analítico para a resolução da equação de Schrödinger dependente do tempo. Esse método é empregado na simulação de reações químicas envolvendo compostos orgânicos, a fim de estimar o impacto ambiental decorrente do lançamento de substâncias químicas em corpos hídricos. O método é baseado no emprego de transformações de Bäcklund, que fornecem soluções em forma fechada expressas como funções arbitrárias de argumentos conhecidos. Essa formulação proporciona uma redução significativa no tempo de processamento requerido para a obtenção de resultados numéricos, permitindo a simulação de sistemas químicos contendo até 200 átomos. Simulações simbólicas são apresentadas. / This work presents a new analytical method for solving the time dependent Schrödinger equation. This method is employed to simulate chemical reactions involving organic compounds, in order to estimate environmental damages due to the disposal of chemicals in water bodies. The method is based on Bäcklund transformations, which supply closed form solutions as arbitrary functions of known arguments. This formulation provides a significant reduction on the time processing required to obtain numerical results, allowing to carry out the simulation of chemical systems containing up to 200 atoms. Symbolic simulations are reported.

O estudo do escoamento ao redor de corpos submersos tem sido foco de muitas pesquisas. Escoamentos ao redor de cilindros constituem sistemas de grande interesse prático em vários campos da engenharia e outras ciências. Neste trabalho procurou-se estudar a dinâmica do escoamento ao redor de um cilindro circular submetido a um movimento forçado analisando as alterações ocorridas na formação da esteira de vórtices. Para a realização das simulações foi utilizada a técnica de Simulação Numérica Direta (DNS), com um esquema de diferenças finitas compacto de sexta ordem, em associação com um método de representação virtual das fronteiras do cilindro. Quatro tipos de trajetória foram estudadas: uma trajetória com oscilação transversal, onde o cilindro é forçado a mover-se na direção transversal em relação ao escoamento, uma trajetória elíptica, uma trajetória em arco e uma trajetória em oito. Buscou-se identificar as principais características destes escoamentos, como os modos e freqüências de desprendimento de vórtices na esteira formada, coeficientes hidrodinâmicos e parâmetros estatísticos pertinentes ao fenômeno estudado. As simulações foram realizadas com diferentes amplitudes longitudinais de movimento sendo que a freqüência de movimento do cilindro foi mantida constante para as trajetórias elíptica, em arco e em oito. Os resultados para o caso do cilindro com oscilação transversal representaram bem os fenômenos físicos ocorridos neste tipo de escoamento. Foi possível identificar as fronteiras da região de lock-in primário através das variações bruscas nos coeficientes hidrodinâmicos e do ângulo de fase entre o deslocamento do cilindro e a força de sustentação. No caso das simulações com trajetória elíptica e número de Reynolds Re = 140, alterações significativas foram identificadas na esteira de vórtices à medida que a amplitude longitudinal do movimento era modificada. Nas simulações com Re = 400 e trajetória transversal, notou-se a formação de uma esteira com dois vórtices e duas camadas cisalhantes. Comparando esta simulação com os resultados das simulações em trajetória elíptica, foram observadas alterações no padrão de desprendimento, sendo formado um par de vórtices na região superior da esteira e um vórtice simples na região inferior (padrão assimétrico P+S). Comparando as simulações com trajetória elíptica e sentidos de movimento opostos, notou-se que o escoamento se inverte, evidenciando um efeito de memória deste escoamento que é função do sentido adotado para a trajetória. Tanto a amplitude de deslocamento longitudinal quanto o tipo de trajetória imposta ao cilindro influenciam significativamente os padrões da esteira de vórtices formada. / The study of flow around submerged bodies has been the focus of many researches. Flow around cylinders are systems of great practical interest in many fields of engineering and other sciences. This work aims to study the flow dynamics around a circular cylinder under a forced movement by analysing changes in the patterns of the wake. The simulations were performed through the Direct Numerical Simulation technique (DNS), with a sixth-order compact finite differences scheme associated to a virtual boundary method to represent the cylinder surface. Four trajectories were studied: a trajectory with transversal oscillation; a trajectory on elliptical curve; a trajectory in arc curve and a trajectory in eight curve. The main features of these flows were identified, e.g. modes and vortex shedding frequencies, hydrodynamic coefficients and statistical parameters. The simulations were performed with different streamwise movement amplitude, where the movement frequency was held constant for the elliptical, arc and eight trajectories. The results for the cylinder in transversal oscillation successfully represented the physical phenomena. It was possible to identify the boundaries of the primary lock-in region through jumps in the hydrodynamical coefficients and the phase angle between the cylinder displacement and the lift force. For the simulations with elliptical trajectory and Reynolds number Re = 140, significant changes were identified in the wake as the streamwise displacement of the cylinder was modified. In simulations with Re = 400 and transversal trajectory, two vortex and two shear layers could be observed. Comparing this simulation with simulations in elliptical trajectory, changes in the vortex shedding pattern were observed. A vortex pair was formed in the top of the wake and a single vortex in the lower wake ( asymmetrical mode P+S). Comparing the simulations with elliptical trajectory and opposite directions of movement, it was noted that the flow is reversed showing some effect of memory due to the direction of the trajectory. Both the amplitude of streamwise displacement and the trajectory type may result in different patterns in the wake.

Este trabalho apresenta uma metodologia para a simulação numérica de eventos de capotamento de veículos automotores. Pois apesar de o rolamento veicular não ser o mais comum tipo de acidente, ele é o de maior importância com respeito a lesões e traumas provocados nos ocupantes. Assim, reduzir a incidência de mortes e lesões graves tem aumentado a importância de simulações computacionais e testes físicos. Inicia-se com o desenvolvimento de um modelo computacional de multicorpos para entendimento e avaliação dos parâmetros da suspensão que influenciam no fenômeno de capotamento. Os resultados deste modelo de multicorpos são comparados com testes experimentais. Na segunda etapa são desenvolvidos modelos numéricos computacionais de veículos automotivos para simulação e comparação de diferentes normas e procedimentos que avaliam a resistência da estrutura do teto em suportar as cargas decorrentes do fenômeno de capotamento. Este estudo também apresenta modelos computacionais de elementos finitos para simular eventos de capotamento e predizer possíveis lesões ocasionadas à cabeça, pescoço, tórax e coluna. Modelos numéricos de um veículo de passageiros e de um veículo utilitário esportivo (SUV) são simulados com manequim antropomórfico para representar o condutor. São verificados os riscos de lesões e traumas ao condutor se o mesmo estiver usando ou não cinto de segurança. O risco de lesões ao condutor de um veículo de passageiros é comparado ao risco de um condutor de um SUV. A metodologia computacional desenvolvida mostrou-se bastante eficiente na avaliação da estrutura dos veículos em eventos de capotamento. Resultados computacionais demonstram que algumas normas e procedimentos, que atualmente são amplamente utilizados na indústria automotiva para avaliar a estrutura do teto dos veículos, podem não prevenir a ocorrência de lesões e traumas aos ocupantes. Além disso, os resultados das simulações sugerem que a carga utilizada para avaliar a resistência do teto dos veículos seja aumentada. Demonstrando que métodos numéricos podem ser utilizados na avaliação da estrutura dos veículos e verificação de riscos de lesões ao condutor, reduzindo a quantidade de testes experimentais a serem realizados. / This work presents a methodology for the numerical simulation of vehicle’s rollover events. Even though the rollover is not the most frequent type of accident, it is of the greatest significance with respect to injury and trauma caused to the vehicle occupants. The need to reduce death incidence and serious injuries has increased the importance of computational simulations and prototype testing. Firstly, a multi-body model is developed for understanding and evaluating the parameters of the suspension that most influence in the rollover phenomenon. The results of this numerical model are compared with experimental tests. In the second part, numerical models of automotive vehicles are developed for simulation and comparison of different standards and procedures that evaluate the roof resistance with respect to withstand the rollover loads. This study also presents finite element models to simulate rollover events and to predict possible injuries caused in the head, neck, thorax and cervical spine. Numerical models of passenger's vehicle and a sport utility vehicle (SUV) are simulated including anthropomorphic dummy to represent the driver. The injury risks and traumas are verified to the driver considering belted and unbelted dummies. The driver's risk of injury when driving a passenger car is also compared with the risks when driving a SUV. The computational methodology developed proved to be efficient for the evaluation of the vehicle's roof structure in rollover events. Regarding the computational results, some standards and procedures which are broadly used by the automotive industry to evaluate the vehicle's roof resistance, cannot prevent the occurrence of injuries and traumas to the occupants. Moreover, according to the results of numerical simulations, it is suggested that the load required evaluating the car's roof resistance should be increased. It is shown that numerical methods can be used to assess the structure of the vehicle and verifying the driver's injury risk, reducing the amount of experimental tests to be performed.

Calamidades naturais são cada dia mais frequentes e a maior parte delas está relacionada aos danos ocasionados por tempestades severas. Desta forma, são necessários estudos que auxiliam tanto na compreensão desses fenômenos quanto na previsão dos mesmos para que sejam desenvolvidas estratégias mitigatórias e edificações cada vez mais preparadas para as adversidades meteorológicas. Intensas rajadas de vento provocam um dos principais danos provenientes de tempestades severas e podem ser originadas a partir de fenômenos tais como downbursts ou tornados. Downbursts são fenômenos meteorológicos que apresentam distinto perfil de velocidade de vento e turbulência que causam diversos acidentes no mundo e, por meio deste trabalho, verificou-se que são a causa de diversos acidentes também no Brasil, em especial na região da Bacia Amazônica, na Região Sul e Sudeste brasileira e na região próxima ao litoral do nordeste brasileiro. A previsibilidade de downbursts também é abordada nessa pesquisa; buscou-se verificar a relação existente entre altos valores do índice CAPE e baixos valores de CINE na ocorrência de ventos provenientes de downbursts, tendo sido constatado que esses índices não são eficientes se analisados de forma isolada na previsão de downbursts. Por fim, verificou-se que estudos mais aprofundados a respeito das ações de downbursts sobre edificações fazem-se necessários para o Brasil. Por fim desenvolveu-se uma revisão de literatura a respeito da simulação física e numérica do fenômeno, buscando gerar um embasamento teórico para futuros trabalhos nesta área. / Natural disasters are becoming more frequent every day and most of them are related to damages caused by severe storms, thus studies that help both in understanding these phenomena and in predicting are necessary to the development of mitigation strategies and to the development of building increasingly prepared to the weather adversity. Intense wind gusts are responsible for the major damage from severe thunderstorms and can arise from phenomena such as downbursts or tornadoes. Downbursts are meteorological phenomena that exhibit distinctive profiles of wind speed and turbulence causing many accidents in the world and through this work are found to be the cause of many accidents in Brazil too, especially in the Amazon Basin region in the South and Southeast Brazilian’s Regions and in the region next to the coast of Brazilian’s Northeastern. Predictability of downbursts is also addressed in this study, it was valued the relationship between high values of CAPE index and low values of CINE index for the occurrence of winds from downbursts, it was found that these index are not efficient when analyzed in isolation in predicting downbursts. Finally, it was concluded that more intensive studies on the effects of downbursts on structures are up needed in Brazil. Finally it was developed an investigation on literature about the physical and numerical simulation of this phenomena, seeking to generate a theoretical basis for future work in this area.

Desde a exploração até o planejamento de lavra, metodologias para quantificação do risco associado ao empreendimento mineiro podem melhorar drasticamente o processo de tomada de decisão. Pequenas variações nas condições de contorno de um projeto podem ter impacto significativo sobre o retorno final do mesmo. Dentro deste contexto, na mineração, a incerteza geológica é vista como um dos principais fatores que contribuem para o insucesso do projeto. A medição e gerenciamento do risco geológico na avaliação do projeto e tomada de decisão pode ser realizada em varias etapas do empreendimento mineiro, desde a quantificação dos recursos, sequenciamento de lavra até a exaustão das reservas. Além disso, a incerteza geológica pode ser traduzida no risco financeiro associado, proporcionando uma análise sobre a viabilidade econômica do empreendimento. A necessidade da quantificação da incerteza geológica já era reconhecida pela indústria de mineração desde a década de 1970, onde estimativas globais e locais eram consideradas insuficientes para a otimização do planejamento de produção, sequenciamento de lavra ou estratégias de homogeneização de pilhas. Desta forma, a simulação estocástica vem desempenhando um papel crucial na construção de modelos de incerteza de teores em depósitos minerais, proporcionando uma ferramenta para obter uma análise do risco associado. Este trabalho tem como objetivo o desenvolvimento de uma metodologia que contribua para a construção de cenários de planejamento de lavra e avaliação do impacto provocado pelas incertezas geológicas fornecidas a partir de modelos simulados. As respostas buscadas, usando a abordagem sugerida, compreendem a definição de uma metodologia para controle e planejamento de lavra incorporando a incerteza local associada ao minério, avaliando a sensibilidade do retorno financeiro e o impacto no planejamento de lavra, considerando variações e incertezas relacionadas ao modelo geológico e estimativa de teores. / From exploration to mine planning, methodologies for quantifying the risk associated with mining enterprise can dramatically improve decision making process. Small variations in the boundary conditions of a project may have significant impact on the final return of the same. Within this context, in mining, geological uncertainty is seen as a major factor contributing to the failure of the project. The need to quantify and manage the geological risk for project evaluation and decision making can be translated into the need to assess the geological risk of any parameter at all stages, from the resources quantification, mining sequencing to reserves exhaustion. In addition, the geological uncertainty can be translated into financial risk, providing an analysis of the economic feasibility of the project. However, the need for quantification of geological uncertainty was already recognized by the mining industry since the 1970, where local and global estimates were considered insufficient for the optimization of production planning, sequencing mining, strategies of homogenization. Thus, the stochastic simulation has played a crucial role in the construction of grade uncertainty models in mineral deposits, providing a tool to perform risk analysis. This study aims to develop a methodology able to contribute in building scenarios of mine planning and evaluation of the impact caused by geological uncertainties supplied from simulated models. Answers searched using the suggested approach include understand the definition of a methodology to control and mine planning incorporating the local uncertainty associated with the ore, assessing the sensitivity of financial return and mine planning considering variantions and uncertainty related to the geological model and grade estimation.

O desenvolvimento de protocolos distribuídos é uma tarefa complexa. Em sistemas tolerantes a falhas, a elaboração de mecanismos para detectar e mascarar defeitos representam grande parte do esforço de desenvolvimento. A técnica de simulação pode auxiliar significativamente nessa tarefa. Entretanto, existe uma carência de ferramentas de simulação para investigação de protocolos distribuídos em cenários com defeitos, particularmente com suporte a experimentos em configurações “típicas” da Internet. O objetivo deste trabalho é investigar o uso do simulador de redes NS (Network Simulator) como ambiente para simulação de sistemas distribuídos, particularmente em cenários sujeitos à ocorrência de defeitos. O NS é um simulador de redes multi-protocolos, que tem código aberto e pode ser estendido. Embora seja uma ferramenta destinada ao estudo de redes de computadores, o ajuste adequado de parâmetros e exploração de características permitiu utilizá-lo para simular defeitos em um sistema distribuído. Para isso, desenvolveu-se dois modelos de sistemas distribuídos que podem ser implementados no NS, dependendo do protocolo de transporte utilizado: um baseado em TCP e o outro baseado em UDP. Também, foram estudadas formas de modelar defeitos através do simulador. Para a simulação de defeito de colapso em um nodo, foi proposta a implementação de um método na classe de cada aplicação na qual se deseja simular defeitos. Para ilustrar como os modelos de sistemas distribuídos e de defeitos propostos podem ser utilizados, foram implementados diversos algoritmos distribuídos em sistemas síncronos e assíncronos. Algoritmos de eleição e o protocolo Primário-Backup são exemplos dessas implementações. A partir desses algoritmos, principalmente do Primário-Backup, no qual a simulação de defeitos foi realizada, foi possível constatar que o NS pode ser uma ferramenta de grande auxílio no desenvolvimento de novas técnicas de Tolerância a Falhas. Portanto, o NS pode ser estendido possibilitando que, com a utilização dos modelos apresentados nesse trabalho, simule-se defeitos em um sistema distribuído.

Este trabalho tem por objetivo simular e analisar uma usina termelétrica a carvão em várias condições de funcionamento. A usina simulada neste trabalho é a AVV 1 localizada em Copenhague, Dinamarca. A AVV 1 é uma usina de geração de potência e aquecimento distrital que pode funcionar em várias condições de carga. A simulação da usina supracitada foi tema de um concurso de simuladores proposto no congresso ECOS 2003 realizado em Copenhague, Dinamarca. Para a realização deste trabalho foi construído um programa na linguagem FORTRAN 90. Cada componente da usina é modelado através de equações de balanço de massa e energia, e o sistema completo tem sua solução obtida pelo método de substituição sucessiva. Para viabilizar essa solução é necessário também implementar uma rotina de cálculo de propriedades do fluido de trabalho. No caso estudado, o fluido de trabalho da usina é a água e a formulação utilizada para o cálculo de suas propriedades nos diversos estados é a IAPWS IF-97. A usina é simulada em dois modos de operação: modo de condensação, onde ocorre apenas geração de eletricidade, e em modo de contrapressão, onde há geração de eletricidade e aquecimento distrital, conforme nomenclatura sugerida pela organização do concurso No modo de condensação, são feitas quatro séries de simulações variando a carga de 100% a 40%. Cada série contém um conjunto de hipóteses quanto à variação das eficiências isoentrópicas e pressões das turbinas em função da vazão mássica. No modo de contrapressão, a usina é simulada funcionando com 100% da carga. O programa desenvolvido calcula as propriedades para qualquer ponto de trabalho ao longo da planta, assim como a eficiência da mesma, a potência gerada, e todas as vazões mássicas pertinentes. Além disso, é feita também uma análise exergética da planta. A simulação demonstrou que a planta possui uma eficiência global de 42,02% com uma geração de 250,2 MW em 100% de carga no modo de condensação. Nessas mesmas condições, do ponto de vista exergético, a eficiência encontrada é de 37,21%. No modo de contrapressão, a usina apresenta uma eficiência exergética de 40,19% com um aproveitamento energético de 90,55%. Por fim, é possível também verificar a comportamento da eficiência da planta e a variação de água de resfriamento do condensador com a carga. Os resultados gerados são próximos àqueles encontrados pelos diversos pesquisadores que abordaram o problema.

Este trabalho apresenta uma sistemática para realizar a otimização numérica de pré-formas e de matrizes em problemas de forjamento axissimétricos e em estado plano de deformações. Para este fim, desenvolveu-se um código computacional composto basicamente de três módulos: módulo de pré-processamento, módulo de análise e módulo de otimização. Cada um destes foi elaborado acrescentando rotinas em programas comerciais ou acadêmicos disponíveis no GMAp e no CEMACOM. Um programa gerenciador foi desenvolvido para controlar os módulos citados no processo de otimização. A abordagem proposta apresenta uma nova função objetivo a minimizar, a qual está baseada em uma operação booleana XOR (exclusive or) sobre os dois polígonos planos que representam a geometria desejada para o componente e a obtida na simulação, respectivamente. Esta abordagem visa eliminar possíveis problemas geométricos associados com as funções objetivo comumente utilizadas em pesquisas correlatas. O trabalho emprega análise de sensibilidade numérica, via método das diferenças finitas. As dificuldades associadas a esta técnica são estudadas e dois pontos são identificados como limitadores da abordagem para problemas de conformação mecânica (grandes deformações elastoplásticas com contato friccional): baixa eficiência e contaminação dos gradientes na presença de remalhamentos. Um novo procedimento de diferenças finitas é desenvolvido, o qual elimina as dificuldades citadas, possibilitando a sua aplicação em problemas quaisquer, com características competitivas com as da abordagem analítica Malhas não estruturadas são tratadas mediante suavizações Laplacianas, mantendo as suas topologias. No caso de otimização de pré-formas, o contorno do componente a otimizar é parametrizado por B-Splines cujos pontos de controle são adotados como variáveis de projeto. Por outro lado, no caso de otimização de matrizes, a parametrização é realizada em termos de segmentos de reta e arcos de circunferências. As variáveis de projeto adotadas são, então, as coordenadas das extremidades das retas, os raios e centros dos arcos, etc. A sistemática é fechada pela aplicação dos algoritmos de programação matemática de Krister Svanberg (Método das Assíntotas Móveis Globalmente Convergente) e de Klaus Schittkowski (Programação Quadrática Sequencial – NLPQLP). Resultados numéricos são apresentados mostrando a evolução das implementações adotadas e o ganho de eficiência obtido.

O presente trabalho é parte integrante de um projeto multidisciplinar cujo objetivo é construir uma ferramenta computacional capaz de simular a operação de uma Usina Termoelétrica tratando, especificamente, da câmara de reação (local da caldeira onde se desenvolve a queima do combustível). Foi desenvolvido um modelo matemático dinâmico, pseudo-homogêneo, unidimensional e, com propriedades físicas variáveis, com o objetivo de descrever a conversão do calor gerado pela combust em energia térmica, que está baseado em balanços de massa e energia. Este modelo foi implementado no simulador gPROMS e seus parâmetros foram estimados utilizando o módulo gEST deste software (método da máxima verossimilhança), com base em dados de operação da Usina e dados cinéticos de literatura. Para a simulação dinâmica da câmara de reação foi utilizado um método de integração implícita de passos múltiplos com ordem e passo variáveis. A câmara de combustão foi dividida em três zonas, sendo realizados, em cada zona, os balanços de massa e energia, adotando-se a hipótese de mistura perfeita. O modelo desenvolvido prevê o perfil de temperatura ao longo da câmara, as correntes de saída e a composição dos gases de combustão em cada zona; incluindo uma previsão aproximada da concentração dos poluentes formados, principalmente monóxido de nitrogênio e dióxido de enxofre. Os resultados apresentados mostraram que o modelo prediz satisfatoriamente a composição do gás de combustão na saída da câmara de reação e que a divisão da câmara em diferentes zonas de combustão é apropriada para avaliar a geração e o aproveitamento energético na caldeira.