Um guia de política monetária no Brasil durante o regime de metas de inflação

O objetivo deste trabalho é abordar os principais mecanismos utilizados pelo Banco Central do Brasil (BCB) no processo de decisão acerca da condução da política monetária no período durante o regime de metas de inflação. O BCB analisa e reavalia seus modelos econométricos, que buscam captar inter-relações dos principais agregados da economia brasileira, em conjunto com informações externas aos modelos e julgamentos subjetivos de integrantes do conselho. Portanto, inicialmente será discutida a conjuntura histórica que levou o país a adotar o regime de metas de inflação, assim como sua aplicabilidade. Na segunda parte, são feitas as estimações das equações IS e Phillips que mostram que a taxa de juros afeta o hiato do produto com defasagem de um trimestre e que o produto se correlaciona positivamente com a inflação com defasagem de mais um trimestre. Num terceiro momento, são feitas estimações de uma função de reação do BCB, sugerindo que a instituição leva em consideração o desvio da expectativa de inflação do mercado com relação a meta, o comportamento do hiato do produto e a trajetória dos preços livres em seu processo de decisão.

Simulação numérica da transferência simultânea de energia e umidade através do solo em um sistema trocador-armazenador de calor

O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.

Controle de temperatura em reatores de polimerização em fase gasosa

Nos reatores de polimerização em fase gasosa, o controle de temperatura em uma estreita faixa de temperatura é importantíssimo para assegurar que a temperatura na zona de reação seja mantida acima do ponto de orvalho dos reagentes, assim como abaixo do ponto de fusão do polímero. Outro ponto importante é que se a temperatura estiver em malha aberta, estes reatores estão propensos à instabilidade e ao surgimento de ciclos limites, entretanto podem ser facilmente estabilizados com um controlador de temperatura apropriado. Neste trabalho é realizado um estudo detalhado do comportamento dinâmico de um reator industrial de polimerização em leito fluidizado, com a localização de pontos de bifurcação e estabilização do sistema com controlador PID projetado via otimização no domínio da freqüência. O modelo completo do sistema foi implementado no software MATLAB® para ajustar seus parâmetros com os dados da planta industrial, realizar simulações dinâmicas e analisar as estruturas de controle. Comparações entre as simulações e os dados da planta industrial mostram que as predições do modelo proposto são satisfatórias para reproduzir o comportamento da planta industrial e as propriedades do polímero. A adição de um sistema de recirculação de gás e de um sistema de troca térmica ao modelo do reator torna possível o estudo dos comportamentos dinâmicos e a localização dos pontos de bifurcação, os quais não poderiam ser reproduzidos somente com o modelo do reator. O modelo desenvolvido foi implementado no AUTO®, um software para problemas de continuação e bifurcação em equações diferenciais ordinárias.

Três métodos para o cálculo da série zeta(2n) de Riemann

Neste trabalho apresentamos três métodos distintos provando que S(n) = +1 X k=−1 (4k + 1)−n é um múltiplo racional de n para todos os inteiros n = 1, 2, 3, . . . O primeiro utiliza a teoria das função analíticas e funções geradoras. No segundo reduzimos o problema, via mudança de variável devida a E. Calabi, ao cálculo do volume de certos politopos em Rn enquanto que no terceiro usamos a teoria dos operadores integrais compactos. Cada um dos métodos tem um interesse intrínsico e está sujeito a generalizações para aplicações em novas situações.

Solução da equação de difusão-advecção para uma CLP não-homogênea e não-estacionária pelo método GILTT

Neste trabalho é apresentada a solução da equação de difusão-advecção transiente para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método analítico GILTT (Generalized Integral Laplace Transform Technique) e da técnica de inversão numérica da quadratura de Gauss. A validação da solução é comprovada utilizando as concentraçãos obtidas a partir do modelo com as obtidas experimentalmente pelo Experimento de Copenhagen. Nesta comparação foram utilizados os perfis de vento potencial e logaritmo e os parâmetros de turbulência propostos por Degrazia et al (1997) [19] e (2002) [17]. Os melhores resultados foram obtidos utilizando o perfil de vento potencial e o coeficiente de difusão propostos por Degrazia et al (1997). A influência da velocidade vertical é mostrada através do comportamento das concentrações de poluentes na pluma. Além disso, as velocidades verticais e longitudinais geradas pelo Large Eddy Simulation (LES) foram colocadas no modelo para poder simular uma camada limite turbulenta mais realística, a qual apresentou resultados satisfatórios quando comparados com os disponíveis na literatura.

Modelo esférico quântico de vidro de spin na aproximação de recozimento

Apresentamos aqui o modelo esférico quântico de vidro de spin usando a aproximação de recozimento. São calculadas a energia livre, bem como a temperatura crítica em função do momentum de inércia e a entropia. São consideradas interações aleatórias de longo alcance (campo médio) com distribuição normal de média zero, e a energia cinética de cada spin. O cálculo é feito utilizando o formalismo funcional de Feynman de integrais de caminhos. O limite clássico é apresentado e coincide com o limite conhecido de teorias anteriores.

Doença infecciosa bursal: estudo sobre amostras vacinais e de campo, imunidade materna e desafio com amostra muito virulenta do vírus.

As aves foram separadas em quatro grupos, dois vacinados e dois não vacinados no 1º dia de vida. As aves foram desafiadas no 1º,4º,7º,10º,13º,16º,19º e 22º dias de vida. A cada dia de desafio eram coletados e medidos, antes e depois, os seguintes itens do grupo: peso relativo da bolsa de Fabrício, diâmetro da bolsa, bolsa de Fabrício para exame histológico e soro para medir os anticorpos contra a DIB, através de Elisa e avaliação clínica da DIB. Os resultados mostraram uma queda de anticorpos sem diferença significativa entre aves vacinadas e não vacinadas. Analisando os diferentes resultados, foi estabelecido que um título de Elisa (log10) de 3,4, foi o ponto de corte entre aves saudáveis e aves doentes. Foram construídas equações de regressão, para o estabelecimento do melhor dia para a vacinação ou do título de Elisa (log10), que as aves podem ter numa determinada idade. Assim sendo, os pintos da Empresa A deveriam receber vacina contra a DIB a partir do 6º - 7º dia, e os da Empresa B deveriam receber a vacina entre o 11º-12º dia de idade. Com os resultados obtidos neste experimento fica claro que as aves não devem ser vacinadas no 1º dia de vida, que o esquema de vacinação das reprodutoras ocasiona diferenças na proteção da progênie e que não se deve aplicar o mesmo esquema de vacinação indiscriminadamente para todos os lotes de frangos.

Características do escoamento ao longo de uma escada de peixes do tipo ranhura vertical

Sucessivos barramentos ao longo de um rio impedem a migração, fenômeno característico de algumas espécies de peixes. Essa interrupção pode provocar a extinção local de espécies migratórias de peixes e acentuada queda da produção pesqueira. Mecanismos de Transposição para Peixes (MTP) são estruturas capazes de mitigar os efeitos negativos desses barramentos, possibilitando a transposição segura dessas espécies através dos barramentos. Esta pesquisa visou a compreensão do funcionamento de um MTP conhecido como escada para peixes do tipo ranhura vertical. Para tanto, foram realizados experimentos em uma estrutura de laboratório, geometricamente semelhante à Escada de Peixes do tipo Ranhura Vertical do reservatório da UHE de Igarapava/MG. Foram realizados experimentos em diversas vazões para a verificação do regime de escoamento ao longo da estrutura e para a determinação dos parâmetros hidráulicos de vazão adimensional, de coeficiente de descarga e de coeficiente de cisalhamento, que foram comparados aos encontrados na bibliografia. Por meio desses ensaios foi possível sugerir equações simplificadas para esses parâmetros Também foram executados ensaios a vazão constante para gerar mapas de distribuição de velocidades médias e de pressões dentro de um tanque da estrutura. A vazão constante também foram medidos valores de altura de lâmina d’água ao longo de dois eixos de um tanque e realizadas visualizações do escoamento por meio do uso de traçadores. Os resultados puderam demonstrar a existência de um jato e de duas zonas de recirculação de água, à esquerda e à direita do tanque, assim como a alta variação de valores de pressões no jato e a existência de velocidades no sentido vertical, principalmente na zona de recirculação à esquerda do tanque.

Um estudo da otimização da geometria de um pára-quedas simplificado

Neste trabalho desenvolve-se um estudo numérico do fluxo de ar em torno da geometria de um pára-quedas tradicional simplificado, para alguns valores de Reynolds. O método baseia-se na solução das equações incompressíveis de Navier- Stokes discretizadas pelo método de diferenças finitas e integradas pelo método de Runge-Kutta. Utiliza-se o método dos contornos virtuais para representar a geometria numa malha cartesiana e o método de otimização não-linear dos poliedros flexíveis para otimização do coeficiente de arraste calculado através do código de dinâmica de fluidos computacional; esteé um método de busca multivariável, onde o pior vértice de um poliedro com n + 1 vérticesé substituído por um novo.

Análise experimental e numérica do comportamento térmico de um coletor solar acumulador

Esta tese apresenta um estudo do comportamento térmico de um coletor solar acumulador e desenvolve uma metodologia para medir a sua eficiência diária. O coletor solar acumulador está instalado na face norte do prédio de Energia Solar da UFRGS e possui cerca de 26 m2. É constituído de uma massa espessa de concreto com uma superfície absorvente feita de tijolos, possuindo uma cobertura dupla de vidros colocada de modo a deixar um espaço para a circulação de ar. Os raios solares atravessam a cobertura de vidro e aquecem a massa absorvente de tijolo, a qual aquece o ar que é introduzido no interior da construção por efeito de termossifão. Uma das principais características do coletor solar acumulador consiste no fato de que a resposta do coletor é defasada no tempo. Este fenômeno permite que o coletor entregue calor ao ambiente mesmo após o término da radiação solar. Essa defasagem dos picos de energia térmica ocorre devido ao baixo valor da difusividade térmica do concreto. Este trabalho foi dividido em duas etapas. A primeira etapa consistiu na montagem de um calorímetro para controle e monitoração das variáveis envolvidas. No interior do calorímetro foram instaladas 36 garrafas com água. As temperaturas dos conteúdos das garrafas, do coletor solar e as radiações envolvidas foram monitoradas através de 26 sensores de temperatura de CI, 8 sensores resistivos PT100 e dois sensores de radiação fotovoltaicos. Para obter as medidas dos sensores instalados foi feita a montagem de um sistema de aquisição de dados interfaceado a um microcomputador A segunda etapa consistiu na produção de um programa computacional, escrito em linguagem Fortran 90, para simular o comportamento térmico dos diversos elementos constituintes do coletor, determinar a potência térmica do coletor solar e sua eficiência diária. Para a simulação numérica do coletor solar acumulador, adotou-se um modelo simplificado bidimensional do mesmo. Foi integrada, através do Método dos Volumes Finitos, a equação de difusão de calor transiente em 2 dimensões. Na formulação das equações lineares optou-se pelo emprego das diferenças centrais no espaço e formulação explícita no tempo. Ao todo foram produzidas 4 malhas computacionais, com distintos refinamentos e foi realizado o estudo da estabilidade numérica das diversas malhas. Através da montagem experimental obtiveram-se várias características térmicas do comportamento do sistema, entre as quais, a transmitância da cobertura, curvas de temperatura do ar fornecido ao calorímetro e curva da eficiência diária do coletor solar . Através da simulação numérica foi possível determinar a potência térmica que o coletor entrega para o laboratório, a eficiência do coletor, os campos de temperatura e a vazão mássica nos diversos canais interiores do coletor solar.

Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas

No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.

Ambiente de alto desempenho com alta exatidão para a resolução de problemas

Este trabalho visa a disponibilização de um ambiente de alto desempenho, do tipo cluster de computadores, com alta exatidão, obtida através da utilização da biblioteca C–XSC. A alta exatidão na solução de um problema é obtida através da realização de cálculos intermediários sem arredondamentos como se fossem em precisão infinita. Ao final do cálculo, o resultado deve ser representado na máquina. O resultado exato real e o resultado representado diferem apenas por um único arredondamento. Esses cálculos em alta exatidão devem estar disponíveis para algumas operações aritméticas básicas, em especial as que possibilitam a realização de somatório e de produto escalar. Com isso, deseja-se utilizar o alto desempenho através de um ambiente de cluster onde se tem vários nodos executando tarefas ou cálculos. A comunicação será realizada por troca de mensagens usando a biblioteca de comunicação MPI. Para se obter a alta exatidão neste tipo de ambiente, extensões ou adaptações nos programas paralelos tiveram que ser disponibilizadas para garantir que a qualidade do resultado final realizado em um cluster, onde vários nodos colaboram para o resultado final do cálculo, mantivesse a mesma qualidade do resultado que é obtido em uma única máquina (ou nodo) de um ambiente de alta exatidão. Para validar o ambiente proposto foram realizados testes básicos abordando o cálculo do produto escalar, a multiplicação entre matrizes, a implementação de solvers intervalares para matrizes densas e bandas e a implementação de alguns métodos numéricos para a resolução de sistemas de equações lineares com a característica da alta exatidão. Destes testes foram realizadas análises e comparações a respeito do desempenho e da exatidão obtidos com e sem o uso da biblioteca C–XSC, tanto em programas seqüenciais como em programas paralelos. Com a conseqüente implementação dessas rotinas e métodos será aberto um vasto campo de pesquisa no que se refere ao estudo de aplicações reais de grande porte que necessitem durante a sua resolução (ou em parte dela) da realização de operações aritméticas com uma exatidão melhor do que a obtida usualmente pelas ferramentas computacionais tradicionais.

Aplicação do princípio da correspondência elasto-viscoelástica para previsão de deformabilidade de misturas asfálticas

Para a pesquisa, utilizou-se um concreto asfáltico com ligante modificado com borracha ensaiado nas dependências do LAPAV-UFRGS, sob compressão diametral de cargas cíclicas, cujo procedimento é brevemente descrito, associando-se os seus procedimentos com as hipóteses assumidas na modelagem, o que é importante na concepção do modelo empregado. Apresenta-se e discute-se o manuseio dos dados experimentais para o modelo, bem como a escolha do modelo constitutivo adotado, no qual se optou por um modelo viscoelástico linear anisotrópico simplificado, cujas equações base e métodos para calibração são apresentados, discutidos e comparados, bem como os resultados das calibrações dos três primeiros ciclos de carga-descanso, no qual se obteve uma excelente descrição do comportamento do material, o que foi comprovado nas previsões dos dois ciclos seguintes. As previsões de resultados a ensaios sob diferentes condições aos das calibrações mostraram resultados qualitativamente bons às freqüências de ciclo de 1 Hz, comprovados na comparação das curvas de deslocamentos normalizadas. Para a pesquisa, utilizou-se um concreto asfáltico com ligante modificado com borracha ensaiado nas dependências do LAPAV-UFRGS, sob compressão diametral de cargas cíclicas, cujo procedimento é brevemente descrito, associando-se os seus procedimentos com as hipóteses assumidas na modelagem, o que é importante na concepção do modelo empregado. Apresenta-se e discute-se o manuseio dos dados experimentais para o modelo, bem como a escolha do modelo constitutivo adotado, no qual se optou por um modelo viscoelástico linear anisotrópico simplificado, cujas equações base e métodos para calibração são apresentados, discutidos e comparados, bem como os resultados das calibrações dos três primeiros ciclos de carga-descanso, no qual se obteve uma excelente descrição do comportamento do material, o que foi comprovado nas previsões dos dois ciclos seguintes. Para a pesquisa, utilizou-se um concreto asfáltico com ligante modificado com borracha ensaiado nas dependências do LAPAV-UFRGS, sob compressão diametral de cargas cíclicas, cujo procedimento é brevemente descrito, associando-se os seus procedimentos com as hipóteses assumidas na modelagem, o que é importante na concepção do modelo empregado. Apresenta-se e discute-se o manuseio dos dados experimentais para o modelo, bem como a escolha do modelo constitutivo adotado, no qual se optou por um modelo viscoelástico linear anisotrópico simplificado, cujas equações base e métodos para calibração são apresentados, discutidos e comparados, bem como os resultados das calibrações dos três primeiros ciclos de carga-descanso, no qual se obteve uma excelente descrição do comportamento do material, o que foi comprovado nas previsões dos dois ciclos seguintes.

Controle preditivo de sistemas mistos lógicos-dinâmicos

O comportamento dinâmico de um sistema é tradicionalmente descrito por um modelo, em geral associado a um conjunto de equações diferenciais ou de equações de diferenças, onde as variáveis representam grandezas físicas. Sistemas complexos, principalmente na indústria de processos, incorporam elementos com comportamento dinâmico lógico, tais como atuadores e sensores ON-OFF (estados aberto/fechado) ou proposições lógicas (estados verdadeiro/falso). Estes sistemas são denominados “Sistemas Híbridos”, “Sistemas Mistos Lógicos-dinâmicos” ou, simplesmente, “Sistemas Mistos”. Neste trabalho, são apresentadas técnicas que, associando variáveis lógicas a estes elementos, e mediante a incorporação de restrições sobre as variáveis, permitem obter um modelo matemático do sistema misto. Neste caso, técnicas clássicas de controle não permitem a incorporação destas novas variáves e restrições. Como opção de controle de sistemas mistos, é então proposta e estudada uma técnica de controle preditivo baseado em modelo. São apresentados tanto a formulação teórica do problema de controle, quanto exemplos e simulações bem como um estudo de caso de sua aplicação sobre um sistema de equalização de uma planta de tratamento de efluentes.

Análise energética e exergética de plantas de co-geração

O presente trabalho apresenta uma comparação das eficiências energética e exergética de três alternativas diferentes de montagem de sistemas de co-geração de energia que empregam o gás natural como combustível. O caso estudado é uma situação real de um curtume que utiliza lenha como combustível para geração de sua energia térmica e compra energia elétrica da distribuidora. A primeira alternativa é composta por um motor de combustão interna a ciclo Otto, a segunda emprega uma turbina a gás a ciclo Brayton e a terceira usa um gerador de vapor superaquecido juntamente com uma turbina a vapor em contrapressão. Para simulação das plantas é utilizado um software comercial que foi alterado, ficando capacitado para calcular exergias físicas e químicas de plantas de utilidades, cálculos exergoeconômicos, tais como taxas de custo das correntes e eixos, e também de determinar para cada equipamento os valores das variáveis usadas para a avaliação e o melhoramento do desempenho de plantas. Após as modificações o programa foi testado resolvendo um problema conhecido e cuja solução encontrada foi comparada com o da bibliografia. Foram obtidos resultados muito próximos ao padrão usado para comparação, os desvios encontrados são desprezíveis para os objetivos deste trabalho. A planta com melhor eficiência energética e exergética é aquela com motor ciclo Otto, seguida pela montagem que emprega a turbina a gás e finalmente, o conjunto gerador de vapor e turbina a vapor. São calculados os custos de insumos: investimentos necessários para construção das plantas, os custos de operação e manutenção e os custos com combustível. A partir destes valores são resolvidos os sistemas de equações para determinar as taxas de custos das correntes e eixos. Para as três plantas as variáveis de decisão são modificadas buscando aumentar a eficiência exergética e diminuir os custos dos produtos de co-geração. A metodologia para avaliação e melhoramento do desempenho das plantas é empregada e observa-se que é na terceira configuração é onde as modificações tiveram maior impacto, seguida pela primeira planta e por último a segunda alternativa.