UNIVERSIDADE ESTADUAL DE CAMPINAS

FACULDADE DE ENGENHARIA MECÂNICA

COMISSÃO DE PÓS-GRADUAÇÃO EM ENGENHARIA MECÂNICA
DEPARTAMENTO DE ENGENHARIA TÉRMICA E DE FLUIDOS

Modelagem Numérica de Sistemas

de Refrigeração
por Compressão de Vapor

Autor: Matheus Marcelino Brandão

Orientador: Prof. Dr. Luiz Felipe Mendes de Moura

Curso: Engenharia Mecânica

Área de concentração: Térmica e Fluidos

Dissertação de Mestrado apresentada à comissão de Pós Graduação da Faculdade de Engenharia

Mecânica, como um dos requisitos para obtenção do titulo de Mestre em Engenharia Mecânica.

Campinas, 2005
SP - Brasil

i
Capítulo 1

Introdução

Ao longo das últimas décadas e principalmente dos últimos anos, pudemos notar um
crescente aumento na utilização de sistemas de refrigeração e ar condicionado.

Tais sistemas, cada vez mais, têm sido utilizados, tanto para proporcionar conforto, quanto
para processos industriais.

Os sistemas de condicionamento de ar para conforto, embora ainda considerados um luxo

por algumas pessoas, mostram-se cada vez mais necessários e presentes em nosso dia-a-dia. Tais
sistemas visam o controle da temperatura, umidade, pureza e distribuição do ar, no sentido de
proporcionar conforto aos ocupantes de um determinado recinto. As instalações de ar
condicionado para conforto podem variar em tamanho e complexidade, e podem ter diversas
aplicações, tais como: edifícios comerciais, residências, shopping centers, veículos, etc.

Já os sistemas de condicionamento de ar para finalidades industriais visam o tratamento do

ar, tanto para proporcionar condições mínimas de conforto a trabalhadores em ambientes
insalubres, quanto ao controle das condições do ar em um processo industrial. Algumas das
aplicações mais comuns de sistemas de ar condicionado para finalidades industriais são:
laboratórios com atmosfera controlada, imprensa, indústrias de componentes eletrônicos,
alimentícias, têxteis, salas de computadores, etc.

Diante da existência de tantas aplicações para sistemas de refrigeração e ar condicionado

surge um problema: como prever o comportamento de um sistema quando este é submetido a
condições diferentes das condições ótimas para as quais o mesmo foi projetado?

1
 Uma das formas de responder a essa pergunta é a aplicação de um método para a simulação
numérica que pudesse realizar tal previsão com um certo grau de confiabilidade. Tal método
deverá utilizar-se de um conjunto de equações agrupadas em linguagem computacional e um
conjunto de dados obtidos experimentalmente. Desta maneira, torna-se possível à obtenção de
resultados através da simulação com a utilização de modelos numéricos com uma margem de
erro segura. Tais métodos permitiriam também a simulação do comportamento de equipamentos
de refrigeração e ar condicionado quando estes são submetidos a variações em suas condições de
operação.

1.1 Objetivos

O principal objetivo deste trabalho de revisão bibliográfica foi a pesquisa de trabalhos

disponíveis na literatura a respeito do tema em questão. A pesquisa de tais trabalhos é de vital
importância para a compreensão e o desenvolvimento das etapas seguintes do projeto de
pesquisa. As obras disponíveis na literatura consultada forneceram uma base sólida para a
fundamentação dos conceitos principais do projeto de pesquisa.

Partindo-se do trabalho de revisão, o principal objetivo do projeto de pesquisa é a

elaboração de um modelo numérico capaz de simular com um bom grau de confiabilidade o
comportamento de um sistema de refrigeração operando em regime permanente de serviço.

1.2 Metodologia

A metodologia utilizada para a realização deste trabalho foi a consulta às diversas bases de
dados e livros presentes no acervo das bibliotecas da Unicamp. Os resultados obtidos em tal
pesquisa foram levados ao conhecimento do orientador responsável, Prof. Dr. Luiz Felipe M.
Moura e, posteriormente foi desencadeada uma discussão para a definição das obras mais
relevantes para o posterior desenvolvimento do projeto de pesquisa.

2
Capítulo 2

Revisão da Literatura

Yasuda et al. (1983) descrevem o desenvolvimento de um modelo matemático para prever

o comportamento transiente de sistemas de refrigeração. O sistema de refrigeração utilizado para
o desenvolvimento do trabalho em escala de laboratório é composto por um compressor
alternativo monocilíndrico, condensador tipo carcaça e tubos, válvula de expansão termostática e
um evaporador tipo “seco”. O modelamento de todos os componentes foi baseado em equações
originadas em leis físicas da conservação da massa e da energia. Os modelos de todos os
componentes foram combinados de maneira a obter o modelo global do sistema de refrigeração.
Os resultados da simulação foram comparados a resultados experimentais para três casos
distintos de comportamento transiente, obtendo uma precisão satisfatória e assim, confirmando a
validade do modelo.

Koury et al. (2001) propõem um trabalho onde foram desenvolvidos dois modelos
numéricos para simular os comportamentos transiente e de equilíbrio de sistemas de refrigeração
por compressão de vapor. Em tais modelos, o condensador e o evaporador foram divididos em
um certo número de volumes de controle e as equações utilizadas foram obtidas através de
balanços de energia, massa e momento, para cada um dos volumes de controle em questão. A
válvula de expansão e o compressor foram modelados de acordo com modelos de estado de
equilíbrio, pois tais componentes possuem uma inércia térmica muito pequena. As simulações
foram executadas com o objetivo de verificar a possibilidade de controlar o sistema de
refrigeração e o superaquecimento do gás refrigerante na saída do evaporador pela variação da
velocidade do compressor e pela diminuição da área da secção da válvula. Os gases refrigerantes
utilizados em tais simulações foram o R-12 e o R-134-a. Os resultados da simulação, quando
comparados aos experimentais, apresentaram uma boa fidelidade, validando os modelos. Tais

3
modelos numéricos foram desenvolvidos para que, posteriormente, sirvam como base para
trabalhos de desenvolvimento de algoritmos para controle da capacidade de refrigeração e grau
de superaquecimento do sistema.

Reichler (1999) desenvolve em seu trabalho de mestrado um modelo para a simulação de

sistemas centrais de condicionamento de ar, onde as unidades são instaladas no telhado do local.
Para a determinação dos parâmetros do modelo, foram utilizados dados disponíveis nos catálogos
fornecidos pelos fabricantes dos equipamentos. O modelamento foi baseado na utilização de
modelos matemáticos para a descrição do comportamento de cada um dos componentes do
equipamento (compressor, trocadores de calor, entre outros). O modelo foi desenvolvido no
pacote de softwares EES™ (Engeneering Equation Solver). Os resultados obtidos na simulação
do desempenho do sistema foram satisfatórios, permitindo, inclusive, a extrapolação de alguns
dados e a determinação da performance do equipamento quando submetido a tais condições de
operação.

Murphy e Goldschmidt (1986) desenvolvem um modelo para a simulação do

comportamento de sistemas de condicionamento de ar residenciais nas condições de
desligamento (“shutdown”) do compressor. Tal modelo descreve as transições nos fluxos de
energia e massa resultantes do estado de desequilíbrio causado pela diferença de pressões no
desligamento do compressor. Tal desbalanceamento é resultado das diferenças de temperatura
entre os dois trocadores de calor do sistema (evaporador e condensador) e a massa de gás
refrigerante neles contida. O modelo físico utilizado para a realização de tal modelamento
consiste em dois trocadores de calor ligados por duas linhas de transferência. Uma das linhas de
transferência contém o aparelho utilizado para a expansão (no caso em questão, tubos capilares).
A outra linha, possui a válvula para a equalização da pressão do sistema. Neste modelo, também
foi assumido que não há refluxo de gás refrigerante através das válvulas de descarga do
compressor. Foi verificado que a taxa de mudança das pressões no interior de cada trocador de
calor depende da taxa de fluxo de massa. Foi necessário também, um modelamento preciso do
aparelho de expansão para a obtenção de melhores resultados. Os resultados das simulações
forneceram as taxas de fluxo de massa, concentrações de massa, pressões e temperaturas nos
componentes presentes no equipamento modelado. Tais resultados foram obtidos com boa
precisão comparados a resultados experimentais.

4
 MacArthur e Grald (1989) apresentam um modelo dinâmico para bombas de calor por
compressão de vapor. O evaporador, o condensador e o acumulador recebem um tratamento
matemático mais refinado, devido à importância dos efeitos de expansão térmica e
compressibilidade, notado em tais equipamentos. O modelamento dos demais componentes do
equipamento é descrito de forma breve e a metodologia empregada na solução do sistema
completo é apresentada. Os parâmetros considerados no modelo foram desenvolvidos no
equipamento de expansão, bem como no motor de combustão interna, alimentado por gás natural
e no compressor. É apresentado também um modelo para a determinação da distribuição espacial
da massa de refrigerante, onde foi necessária a utilização da introdução de uma fração de vácuo
adequada. Nos trocadores de calor, a determinação das taxas de variação espaciais temporais de
temperatura, entalpia, fluxo de massa e densidade são determinadas para cada ponto. A
determinação da pressão do refrigerante no condensador e no evaporador e das condições das
fronteiras são necessárias para que o controle do conjunto de equações diferenciais parabólicas
seja satisfeito. As comparações entre os resultados da simulação e os dados experimentais são
favoráveis para a operação transiente de uma bomba de calor água-água hermética e um sistema
ar-ar com compressor aberto, validando assim o modelo.

Joudi et al. (2002) apresentam um modelo computacional com o objetivo de simular o

desempenho de um sistema de condicionamento de ar automotivo na operação com diversos tipos
diferentes de gases refrigerantes. O principal objetivo do modelo foi a determinação do gás
refrigerante alternativo mais recomendável para substituir o R-12. Algumas considerações a
respeito de perdas irreversíveis foram incorporadas no modelo com o objetivo de aproximar mais
os resultados da simulação à realidade. Os efeitos de diversos parâmetros que afetam o
desempenho e a compatibilidade do sistema, tais como a temperatura de evaporação, a
temperatura de condensação e a rotação do compressor, foram investigados. Foi utilizado um
total de cinco gases refrigerantes para o estudo com o auxílio do modelo, incluindo o R-12, R-
134a, R-290, R-600a e uma mistura de propano e isobutano. Os resultados obtidos com a
utilização do modelo apontaram a mistura de propano e isobutano como a mais indicada para
substituir o R-12. Os resultados obtidos nas simulações provaram ter boa fidelidade a resultados
obtidos experimentalmente e a resultados teóricos encontrados na literatura.

5
 Duarte et al. (1999) apresentam um modelamento simplificado para um ciclo de
refrigeração por compressão de vapor que opera somente na região de vapor superaquecido. A
utilização de tal ciclo permite rápidas trocas transientes entre os estados de equilíbrio. A
simulação do teste do ciclo de vapor superaquecido tem como base modelos desenvolvidos para o
compressor, o condensador (“desuperheater”), a válvula de expansão e do tanque regulador de
carga do sistema. O modelamento do ciclo é baseado em modelos matemáticos semi-empíricos
desenvolvidos para cada componente. O principal objetivo das simulações é a obtenção, em
diferentes condições operacionais, dos ajustes necessários para o controle do ciclo, ou seja, sendo
impostas uma determinada carga para o refrigerante, a taxa de fluxo de água e a abertura da
válvula, o modelo deverá ser capaz de determinar as pressões de admissão e exaustão, bem como
a temperatura na entrada do compressor. As propriedades do fluido foram modeladas
considerando um fluido com comportamento real, por meio de um subprograma específico
escrito como uma adaptação do TRNSYS (TRNSYS, 1994) disponível. Isso permite que a partir
da entrada com algumas propriedades, tais como: temperatura e pressão, entre outras, o estado
termodinâmico seja facilmente obtido. Os resultados obtidos nas simulações permitiram a
verificação das influências no ciclo da troca dos três parâmetros de entrada descritos acima e,
quando comparado a resultados experimentais, apresentou uma boa fidelidade.

Fu et al. (2003) propõem um modelo dinâmico para uma bomba de calor ar-água de dupla
modalidade equipada com compressor do tipo parafuso. As linhas de alta e baixa pressão são
divididas em três volumes de controle, os quais incluem o lado do tubo preenchido pelo gás
refrigerante, a parede do tubo e o lado do tubo em contato com o ambiente (água ou ar). As
equações diferenciais ordinárias temporais são obtidas através de balanços de massa e energia
realizados nos volumes de controle. Como o compressor, o corpo da válvula de expansão
termostática e a válvula de reversão possuem inércias térmicas muito baixas. Modelos para
estado de equilíbrio foram aplicados para os processos de compressão, estrangulamento e escape.
As relações entre as temperaturas da mistura saturada líquido-vapor e a temperatura do vapor de
refrigerante na saída do evaporador são descritas por equações diferenciais ordinárias temporais.
A simulação do sistema é finalizada com métodos de integração adaptativa e correção dos
valores. Os resultados obtidos na simulação apresentaram uma boa fidelidade aos resultados
obtidos experimentalmente, o que leva à conclusão que o modelo poderá ser usado como uma
boa ferramenta para o teste e desenvolvimento de produtos.

6
 Dabiri e Rice (1981) apresentam um modelo para compressores com correções no nível de
admissão de vapor superaquecido. O modelo para simulação de bombas de calor (“Heat Pump
Simulation Model” ou HPSM) é utilizado neste trabalho, para descrever uma rotina para a
simulação do compressor. Em lugar de gerar informações para estados de equilíbrio da bomba de
calor pela avaliação de um mínimo de dados experimentais, o HPSM realiza uma rotina na qual
são aplicadas correções para a operação em níveis de superaquecimento do evaporador e dados
fornecidos pelos fabricantes do compressor. A performance de outros componentes da bomba de
calor foi simulada de acordo com o modelo utilizado para descrever a rotina do compressor. Os
métodos utilizados para modelar o compressor são derivados da correção dos valores da taxa de
fluxo de massa e potência, obtidas pelos mapas do compressor e levando-se em conta as
variações na sucção de gás superaquecido. Na maior parte dos casos as correções aplicadas não
foram grandes, entretanto, estas correções são de grande importância para a precisão do modelo.
A magnitude das correções aplicadas foi comparada a dados obtidos experimentalmente pelos
fabricantes e a fidelidade dos resultados da simulação com correção aos resultados experimentais
obtidos em compressores instalados em bombas de calor foram verificadas, validando o modelo.

Tassou e Quereshi (1998) apresentam um trabalho sobre a investigação do desempenho de

compressores de deslocamento positivo para aplicações de controle de capacidade pela variação
da velocidade do compressor. Os equipamentos testados incluem compressores alternativos do
tipo aberto, alternativos do tipo semi-hermético e de aletas giratórias do tipo aberto. Todos os
compressores foram testados para condições de pressão constante e variável. Os resultados
demonstraram que todos os compressores foram projetados para máxima eficiência na velocidade
nominal de operação. Quando testados em condições de pressão constante, somente o compressor
do tipo aberto apresentou uma melhora no COP a velocidades reduzidas. Nos testes com pressões
variáveis, todos os compressores apresentaram um aumento no COP com a redução da
velocidade. O resultado das análises de energia baseadas em resultados experimentais e
condições de clima para dois locais representando climas amenos e temperados, demonstrou que
os três tipos de compressor oferecem uma maior economia de energia operando a velocidades
variáveis. Tais análises também demonstraram que o compressor alternativo do tipo aberto é o
mais eficiente, oferecendo uma economia de energia de até 12% para condições de operação em
climas temperados e de até 24% para condições de operação em climas amenos.

7
 Votsis et al. (1992) apresentam o desenvolvimento de uma investigação do comportamento
das características dinâmicas de uma bomba de calor para uso residencial durante as condições de
partida (“start up”) do equipamento. A simulação do equipamento foi feita através da execução
de testes cíclicos para condições que são muito encontradas em instalações residenciais no Reino
Unido, onde a operação da bomba de calor é restringida somente a ciclos para aquecimento de
ambientes. A bomba de calor utilizada é do tipo “split” ar-água, equipada com compressor do
tipo alternativo totalmente hermético e válvula de expansão, com características de operação para
pressões máximas. Nas simulações, foram admitidas as seguintes condições: não migração de gás
refrigerante do condensador para o evaporador quando o ciclo encontra-se inoperante e a
equalização da pressão através de todo o sistema antes da partida do sistema. Os resultados
sugeriram que as perdas no ciclo, associadas à migração de gás refrigerante durante o tempo em
que o ciclo é desligado e a utilização de um acumulador na sucção, representam respectivamente
5% e 4% de decréscimo no COP do sistema quando comparado a casos em que há controle na
migração do gás refrigerante.

Srinivas e Padmanabhan (2002) propõem um modelo computacional eficiente para

compressores alternativos utilizados tipicamente em sistemas de refrigeração. O forro da
cavidade é modelado através da utilização da teoria de ondas planas acústicas e o processo de
compressão é modelado através de equações unidimensionais para dinâmica de gases. O modelo
dinâmico para as válvulas foi baseado em aproximações de vibrações de um grau de liberdade,
em conjunto com equações para gases dinâmicos e modelos de forro acústico. A solução para
estados de equilíbrio resulta da iteração entre equações não lineares resolvidas por algoritmo de
Warner. O algoritmo de Warner foi aplicado à simulação do compressor devido à sua maior
eficiência computacional quando comparado a métodos convencionais. Tais comparações foram
baseadas no número de iterações necessárias e no tempo levado para atingir a convergência,
provando que o algoritmo de Warner é o mais indicado para este tipo de simulação, pois uma
drástica redução no número de iterações pode der observada, comprovando a eficiência do
modelo.

Popovic e Shapiro (1995) desenvolvem em seu trabalho um modelamento semi-empírico

simples para compressores, para a utilização em simulações de sistemas de refrigeração. Tal
modelo é baseado em princípios termodinâmicos e em uma vasta base de dados experimentais. O

8
principal objetivo do modelo é a redução da quantidade de dados requerida para a caracterização
do desempenho de um compressor com o uso de gases refrigerantes alternativos. Dois parâmetros
principais de desempenho do compressor foram desenvolvidos. A queda efetiva de pressão no
compressor levando em consideração as perdas no fluxo de massa e a dependência funcional
numericamente determinada das perdas do coeficiente de transferência de calor em relação às
perdas de energia do compressor. Estes parâmetros descrevem com razoável precisão a operação
do compressor. Foram utilizados neste projeto dados para um total de 122 pontos, 8 diferentes
refrigerantes e 2 compressores. O modelo é capaz de determinar as taxas de fluxo de massa, o
estado do gás refrigerante na saída do compressor e a potência requerida pelo compressor para
todos os pontos com apenas 10% de erro relativo.

Nyers e Stoyan (1994) apresentam um trabalho no qual é desenvolvido um modelo

dinâmico para um evaporador seco utilizado em bombas de calor adequado para propósitos de
controle do sistema. São levados em conta a válvula de expansão, o compressor e o condensador
para modelar o evaporador e seus arredores. O modelamento é feito para parâmetros do tipo
distribuído e consiste em equações parciais diferenciais e equações algébricas. Para a solução das
equações discretizadas foi utilizado o método de iteração de Newton, juntamente com a solução
de blocos de Gauss para sistemas linearizados. Como o modelamento matemático das variações
temporais e espaciais no evaporador é complicado e, muitas vezes, difícil de ser descrito. Foram
assumidas algumas condições simplificadoras, mas ao mesmo tempo, houve uma preocupação
em manter algumas partes de um modelo mais complexo, consideradas de grande importância
para o controle da bomba de calor. Os métodos numéricos utilizados para modelar o sistema são
de considerável complexidade, tornando impossível uma solução analítica. No entanto, estes
métodos permitem uma solução aproximada. Para a obtenção de tal solução foram aplicados
diversos métodos, destacando-se os métodos da iteração de Newton em conjunto com o método
da iteração dos blocos de Gauss Seidel, pelos quais os resultados mais satisfatórios foram
obtidos. Os resultados obtidos foram apresentados na forma gráfica e correspondem com boa
fidelidade aos resultados obtidos de maneira experimental, tornando o modelo proposto uma
ferramenta versátil no estudo dos fenômenos que ocorrem no evaporador e em seus arredores.

Wang e Touber (1991) propõem um modelo dinâmico distribuído para um tipo de

evaporador comum a muitos sistemas de refrigeração, o resfriador de ar. O modelamento foi

9
realizado por meio da derivação de equações utilizadas para estados de não equilíbrio. Por sua
aplicação puramente distribuída, este modelo pode ser amplamente utilizado em diversos tipos de
arranjo dos circuitos de tubos no interior do evaporador, dando uma maior flexibilidade para a
aplicação do modelo em questão. O modelo foi dividido em duas partes: período em que o
sistema encontra-se em operação e período em que o sistema encontra-se desligado. A razão para
a utilização de tal divisão é que a maioria dos equipamentos de refrigeração encontrados em
operação hoje em dia é equipada com elementos de controle para ligamento e desligamento
automático do sistema. O principal objetivo do modelo foi o estudo e a otimização de
comportamentos de estados de não equilíbrio em sistemas de refrigeração com controle de
capacidade.

Willatzen et al. (1998) apresenta um modelo matemático no qual é descrito o fenômeno

transiente no fluxo bifásico de refrigerante em trocadores de calor. O modelo é baseado em
equações diferenciais parciais unidimensionais, nas leis de conservação da massa e energia e em
equações da energia em paredes de tubulares. Um conjunto relativamente simples de equações
diferenciais ordinárias foi criado, integrando separadamente os três estados do vapor que
geralmente são encontrados no interior do trocador de calor (vapor, líquido-vapor e líquido). Uma
atenção especial foi tomada no estudo de operações em regime transiente causadas pelo controle
de liga-desliga presente em equipamentos de refrigeração por compressão de vapor. O modelo
desenvolvido pode ser facilmente aplicado a refrigeradores domésticos e pode ser uma
ferramenta de grande utilidade na otimização do consumo de energia em sistemas de
refrigeração.

Mithraratne e Wijeysundera (2001) desenvolvem um estudo experimental e numérico do

comportamento dinâmico de evaporadores em contra-fluxo com aquecimento de água. A
distribuição da freqüência das oscilações aleatórias no ponto de transição da mistura de vapor e a
temperatura de superaquecimento na saída do evaporador foram obtidas para a operação em
condições de equilíbrio do sistema. O movimento do ponto de transição foi encontrado para
causar flutuações na temperatura do refrigerante 1 m acima do limite inferior da escala de
movimento. Mudanças de fase nas taxas de fluxo de refrigerante e de aquecimento do fluxo de
água demonstram a não-linearidade das características do evaporador onde as constantes de
tempo para aumentos e decréscimos da mesma magnitude diferem significativamente. O modelo

10
distribuído prevê a variação da temperatura de superaquecimento e a pressão seguindo as
mudanças de fase nas entradas com uma boa precisão.

Zalewski e Gryglaszewski (1997) desenvolvem um modelo visando a discussão dos

processos de transferência de calor e massa em resfriadores de fluido evaporativos que operam
em contracorrente com fluxo de ar. O modelo consiste na utilização de quatro equações
diferenciais ordinárias com suas respectivas condições de contorno em conjunto com algumas
equações algébricas. Foi também desenvolvido um método que permitiu um ajuste dos arranjos
geométricos utilizado em trocadores de calor com tubos desencapados. A fim de aumentar a
precisão dos resultados obtidos pelo modelo, foram introduzidas funções para a correção dos
coeficientes de transferência de massa. Os resultados obtidos pelas simulações utilizando o
modelo foram comparados aos resultados obtidos em testes experimentais, apresentando bons
resultados.

Pfafferrott e Schmitz (2004) apresentam o desenvolvimento de uma biblioteca para

sistemas de refrigeração que operam com CO2 como gás refrigerante baseados em uma biblioteca
para fluidos térmicos existente no programa Modellica ™. O desenvolvimento de tal biblioteca
foi realizado em um projeto de pesquisa na Airbus Deutschland e na TUHH e foi focado
principalmente na obtenção de uma biblioteca para um maior detalhamento de investigações
numéricas de sistemas de refrigeração que operam com CO2. Foram feitas avaliações no conceito
de sistemas de resfriamento a bordo de aviões, bem como o modelamento de trocadores de calor
que operam com CO2 com o auxílio da modelagem no programa Modellica ™ e da biblioteca
criada para o CO2. Foram feitas comparações com resultados obtidos em trocadores de calor em
estado de equilíbrio e os resultados obtidos foram satisfatórios. Os resultados apresentados por
simulações em estados transientes também foram satisfatórios, quando comparados aos
resultados experimentais para tais condições.

11
Capítulo 3

Modelagem Teórica e Matemática

3.1 Modelos de compressores

Compressores podem ser considerados como os principais equipamentos que compõem

sistemas de compressão de vapor. Os principais tipos de compressores utilizados em sistemas
frigoríficos são: alternativo, parafuso, centrífugo e de palhetas. Os compressores alternativos
consistem em um êmbolo, o qual realiza um movimento alternativo no interior de uma cavidade
cilíndrica (cilindro) com válvulas de admissão e descarga abrindo e fechando de maneira
conveniente, de modo a permitir a compressão de gás refrigerante. Já os outros tipos de
compressores, valem-se de elementos rotativos para realizar a compressão do gás refrigerante.
Por ser o mais utilizado em sistemas de refrigeração, o compressor alternativo foi o escolhido
para ser abordado nos modelos a seguir.

Um modelo simples, comumente observado na literatura, é o modelo que supõe a

inexistência na variação no fluxo de massa entre a entrada e a saída do compressor. Tal modelo
permite computar o fluxo de massa e a entalpia na saída do compressor. A vazão mássica de
refrigerante através do compressor pode ser calculada através da equação:

m f = NVρ f η v (3.1)

Onde:
N é a rotação do compressor em rps;
V é o volume deslocado do pistão;
ρ f é a densidade do refrigerante na entrada do compressor;

13
 η v é a eficiência volumétrica.

Para a equação acima, a eficiência volumétrica pode ser obtida através da seguinte equação:

cv
P (c )
η v = 1 + c r − [c r ( c ) p ] (3.2)
Pe

Onde:
c r é o fator de folga do compressor (razão entre volume morto e volume deslocado pelo
pistão);
Pc é a pressão de condensação do refrigerante;

Pe é a pressão de evaporação do refrigerante;

cv é o calor específico para volume constante na entrada do compressor;

c p é o calor específico para pressão constante na entrada do compressor.

Em um primeiro momento, a entalpia na saída do compressor pode ser calculada levando-

se em consideração que o processo de compressão é isoentrópico (s= cte). A diferença entra as
entalpias encontradas na entrada e na saída do compressor (trabalho de compressão) é dividido
pela eficiência de compressão. O valor da eficiência de compressão pode ser adotado baseando-se
em dados experimentais. Então, as temperaturas e entalpias são recalculadas de maneira a
representar o novo valor do trabalho de compressão. Os dados de entrada para este modelo
podem ser obtidos de maneira semi-empírica.

Outro tipo de modelo analisado na literatura é o que leva em consideração o volume morto
no interior do cilindro e algumas perdas. Tal modelo é obtido através de um caminho semi-
empírico e será descrito mais detalhadamente a seguir.

A vazão mássica do refrigerante para este modelo pode ser calculada através da equação
(3.1).

A partir do cálculo da vazão mássica, pode-se calcular a potência através da seguinte

equação:

14
Wcp = Wio + (1 + α )Wis (3.3)

Onde:
Wcp é a potência do compressor;

Wio são as perdas eletromagnéticas do motor do compressor (consideradas constantes);

α é o fator de perdas (experimental);

Wis é a potência de compressão politrópica.

γ −1 γ
P ( ) (
γ −1
)
Wis = m f Pe v 2 [( c ) γ
− 1] (3.4)
Pe

cv
γ = (3.5)
cp

Onde:
v2 é o volume específico do gás refrigerante na saída do compressor;

Wcp = (h2 − h1 )m f (3.6)

Onde:
h1 e h2 são as entalpias do refrigerante na entrada e na saída do compressor, respectivamente.

Para o cálculo da eficiência volumétrica do compressor (ηv), podemos utilizar a equação

(3.2).

T2 = T1 + {[Wio + Wis (1 + α )] / m f * c p12 [ ]} (3.7)

Onde:
T2 é a temperatura do refrigerante na saída do compressor;
T1 é a temperatura na entrada do compressor;
cp12 é o valor médio do calor específico a pressão constante entre a entrada e a saída do
compressor.

15
 Os parâmetros N, c r , V , Wio e α são característicos do compressor utilizado e devem ser

obtidos em catálogos de fabricantes ou por meio de resultados experimentais.

É importante observarmos que, para o modelo acima, a entalpia do refrigerante na saída do

compressor pode ser calculada através da equação (3.6) e que a temperatura na saída do
compressor pode ser calculada através da equação (3.7).

Existem também modelos mais complexos para compressores alternativos, onde as

cavidades de entrada e saída das válvulas são modeladas de acordo com teorias de ondas planas
acústicas. O processo de compressão é modelado por equações para escoamento compressível
unidimensional e as válvulas são modeladas de acordo com equações baseadas em teorias de
vibrações. O modelo opera com um conjunto de equações não-lineares. O princípio de
funcionamento do modelo é a obtenção da oscilação das variáveis envolvidas no estado de
equilíbrio, para uma dada configuração mecânica (tamanho do cilindro, diâmetro das válvulas,
constantes de molas das válvulas de admissão e descarga, etc.). A solução resultante das equações
não lineares para o estado de equilíbrio é colocada como um valor do problema nas fronteiras e
resolvida através de um algoritmo de Warner. Em tal algoritmo, o valor inicial do problema é
convertido em dois pontos nas fronteiras, os quais fornecem condições periódicas dos fenômenos
que ocorrem nas fronteiras. É importante ressaltar que a utilização de tal modelamento depende
de dados obtidos experimentalmente e dados obtidos de fabricantes de compressores.

3.2 Modelos para trocadores de calor

Os trocadores de calor (condensador e evaporador) são os componentes de um sistema de

refrigeração que exercem maior influência no desempenho das trocas térmicas exercidas pelo
sistema. Os modelos para simular o comportamento de evaporadores e condensadores podem ser
classificados de uma maneira geral, em quatro diferentes categorias:
1. Modelos do tipo “black box” (caixa preta);
2. Modelos de zona única;
3. Modelos de duas ou três zonas;
4. Modelos distribuídos.

16
 Os modelos do tipo “black box” são baseados em teorias para automação e controle de
sistemas, sendo mais comumente utilizados para o modelamento de evaporadores. Neste tipo de
modelo, o componente é representado por um conjunto de funções de transferência com diversas
constantes, as quais são determinadas experimentalmente. A principal vantagem na utilização de
tal modelo é a simplicidade do mesmo. No entanto, dados experimentais devem ser obtidos para
garantir o bom funcionamento e a precisão na aplicação desse tipo de modelo.

Os modelos de zona única, duas zonas e distribuídos são baseados em leis físicas da
conservação de massa, quantidade de movimento e energia. Nos modelos de zona única, o
trocador de calor é analisado por meio de variáveis que representam o componente de uma
maneira generalizada. Nestes modelos, somente um valor é utilizado para representar a entalpia
do refrigerante contido no interior do trocador de calor. De maneira semelhante, apenas um valor
é utilizado para representar a densidade e temperatura do fluido secundário (lado externo do
trocador), entre outros. Os coeficientes de transferência de calor entre o refrigerante, o fluido
secundário e a parede dos tubos é suposto como sendo constante ao longo do trocador de calor.

Um modelo de zona única, simples, porém amplamente utilizado na literatura, é o dado

pela equação:

Q = m f (h fi − h fo ) (3.8)

Onde:
Q é o calor trocado,
mf é a vazão mássica de refrigerante no interior do trocador de calor (considerada constante)
hfi e hfo são as entalpias do refrigerante na entrada e na saída do trocador, respectivamente.

Para aplicar tal modelo, deparamo-nos com a dificuldade prática da obtenção das entalpias
na entrada e na saída do trocador de calor. Para obtê-las, é necessário o conhecimento das
temperaturas de entrada e saída do refrigerante no trocador de calor, o que, na prática, torna-se
uma dificuldade para a utilização do modelo. Tal modelo, no entanto, pode ser considerado de
grande utilidade para estimativas teóricas do desempenho de sistemas de ar condicionado.

17
 Outro modelo de zona simples, muito encontrado na literatura, é o modelo da efetividade
(NUT). A principal vantagem de tal modelo é a possibilidade da análise de trocadores de calor
quando as temperaturas de saída dos fluidos quente e frio (fluido secundário e refrigerante, no
caso) são desconhecidas, ou difíceis de serem determinadas.

Para definir a efetividade de um trocador de calor, primeiramente devemos determinar a

máxima taxa de transferência de calor possível para o trocador (qmax). Tal taxa, a princípio,
poderia ser obtida através de um trocador de calor em contracorrente com comprimento tendendo
a infinito.

A partir daí, podemos formular a seguinte expressão geral:

q max = c p min (Tq ,e − T f ,e ) (3.9)

Onde:
qmax é a taxa máxima de transferência de calor admitida no trocador
cpmin é o calor específico à pressão constante que apresentar o menor valor entre o fluido quente
para a temperatura da entrada e o fluido frio para a temperatura da entrada.
Tq.e é a temperatura do fluido quente na entrada do trocador de calor;
Tf,e é a temperatura do fluido frio na entrada do trocador de calor.

A efetividade do trocador de calor, então, nada mais é do que a razão entre a taxa real de
transferência de calor e a taxa máxima de transferência de calor possível para um determinado
trocador de calor.

q
ε= (3.10)
q max

Onde:
ε é a efetividade do trocador de calor

Então, podemos também definir a efetividade através das seguintes equações:

18
 c pq (Tq ,e − Tq , s )
ε= (3.11)
c p min (Tq ,e − T f ,e )

ou

c pf (T f ,e − T f , s )
ε= (3.12)
c p min (Tq ,e − T f ,e )

Onde:
cpq e cpf são os calores específicos à pressão constante dos fluidos quente e frio, respectivamente
(média dos calores específicos entre a entrada e a saída do trocador).
Tq,s e Tf,s são as temperaturas do fluido quente e do fluido frio na saída do trocador de calor,
respectivamente.

Por definição, a efetividade, que é um parâmetro adimensional, deve estar no intervalo

compreendido entre 0 ≤ ε ≤ 1. A efetividade pode vir a ser muito útil, uma vez que, se os valores
de ε , Tq,e e Tf,e forem conhecidos, a taxa real de transferência de calor pode ser determinada
através da seguinte expressão:

q = εc p min (Tq ,e − T f ,e ) (3.13)

Para qualquer trocador de calor, pode ser demonstrado que:

 C 
ε = f  NUT , min  (3.14)
 C max 

Onde Cmin/Cmax pode ser igual a cpf / cpq ou cpq / cpf , dependendo das magnitudes das taxas
de capacidades caloríficas dos fluidos quente e frio.

Já o Número de Unidades de Transferência (NUT) é um parâmetro adimensional

amplamente utilizado na análise de trocadores de calor, sendo definido através da seguinte
expressão:

19
 UA
NUT = (3.15)
C min

Onde:
U é o coeficiente global de transferência de calor
A é a área do trocador de calor

Relações Efetividade-NUT:

Para a determinação da forma específica da relação efetividade-NUT (Eq. 3.14) podemos

considerar um trocador de calor com escoamento em paralelo, no qual cpmin=cpq. Desta maneira,
podemos obter a seguinte expressão:

Tq ,e − Tq , s
ε= (3.16)
Tq ,e − T f ,e

Sabemos, também que:

C min m q c p , q T f , s − T f ,e
= = (3.17)
C max m f c p , f Tq ,e − Tq , s

Onde:
mq e mf são as vazões mássicas do fluido quente e do fluido frio, respectivamente.

  C  
1 − exp− NUT 1 +  min  
   C max  
ε= (3.18)
C 
1 +  min 
 C max 

A expressão 3.18 foi desenvolvida através de um processo de rearranjo e substituição de

equações, o qual é apresentado de maneira mais clara na literatura (Incropera F.P. & DeWitt D.P.,
Fundamentos de transferência de calor e massa, 4ª Edição, Editora LTC, Capítulo 11, pp.327).

20
 De maneira semelhante, foram desenvolvidas expressões aplicáveis a outros tipos de
trocadores de calor. Tais expressões podem, facilmente, ser encontradas em sua forma final em
tabelas apresentadas na literatura (Incropera F.P. & DeWitt D.P., Fundamentos de transferência
de calor e massa, 4ª Edição, Editora LTC, Capítulo 11, pp.328,329).

Os modelos de duas ou três zonas para trocadores de calor são os mais utilizados para
simular o comportamento deste tipo de equipamento devido à relativa facilidade do modelamento
matemático e da maior fidelidade ao comportamento real do equipamento, quando comparado
aos modelos de zona única. Em tais modelos, o evaporador é dividido basicamente em duas zonas
(evaporação e superaquecimento) e o condensador é subdividido basicamente em três zonas
(superaquecimento, condensação e sub-resfriamento). Para este tipo de modelamento, as
equações da continuidade, energia e quantidade de movimento podem ser estabelecidas e
aplicadas com maior precisão e, considerando a utilização de um coeficiente de transferência de
calor para cada uma das zonas, o problema físico pode ser melhor descrito e matematicamente
representado.

Os modelos distribuídos têm sido cada vez mais utilizados ao longo dos últimos anos,
devido à maior fidelidade ao comportamento real do equipamento em relação aos modelos
apresentados anteriormente. Neste tipo de modelamento, o trocador de calor é dividido em vários
pequenos volumes de controle. A aplicação das equações para balanços de massa, energia e
quantidade de movimento em cada um dos volumes de controle, juntamente com a utilização de
coeficientes de transferência de calor locais permitem uma maior fidelidade aos fenômenos
físicos que realmente ocorrem nesses equipamentos.

A escolha do tipo de modelo a ser utilizado deve pode ser feita de acordo com a geometria
do trocador de calor. Por exemplo: para trocadores de calor do tipo carcaça e tubos, com o
processo de evaporação ocorrendo na superfície externa do tubo, os modelos de zona única e de
duas zonas podem ser utilizados para representar os fenômenos físicos adequadamente.
Entretanto, para evaporadores de tubos concêntricos, por exemplo, tais modelos não fornecem
uma representação adequada dos fenômenos físicos que ocorrem nos mesmos, sendo necessária à
aplicação de modelos mais complexos, tais como os modelos distribuídos.

21
3.3 Modelos para dispositivos de expansão

Os dispositivos de expansão são importantes componentes para o funcionamento de

sistemas de refrigeração. Tais dispositivos têm como principal função reduzir a pressão na
entrada do evaporador e regular a vazão de refrigerante. Os dispositivos de expansão mais
comumente utilizados são os tubos capilares e as válvulas de expansão termostática. Os tubos
capilares consistem em tubos de pequeno diâmetro que atuam como uma restrição ao fluxo,
diminuindo a pressão. Já as válvulas de expansão termostática controlam a vazão de refrigerante,
por meio de um sensor localizado na linha de descarga do evaporador mantendo constante o grau
de superaquecimento na saída do evaporador.

A maioria dos modelos para dispositivos de expansão é muito simples. Considera-se que a
variação de entalpia em tais dispositivos é nula, ou seja, h entrada = h saída. O fluxo de massa
através de tais dispositivos pode ser modelado através da equação:

mf = K ∆P * ρ (3.19)

Onde:
K é uma constante característica do dispositivo de expansão;
∆P é a diferença entra as pressões de entrada e saída da válvula;
ρ é a densidade do refrigerante na entrada da válvula.

22
Capítulo 4

Modelagem Numérica

O modelo matemático utilizado para a simulação de sistemas de refrigeração como um

todo, na verdade, nada mais é do que um conjunto de equações formado pelas equações presentes
nos modelos de cada componente do sistema (compressor, condensador, evaporador e válvula de
expansão). Tais modelos interagem ente si, de maneira que as variáveis de saída do modelo de
um determinado componente servem como variáveis de entrada para a simulação de outros
componentes.

Neste capítulo serão inicialmente apresentadas as equações representativas de cada

componente. Em seguida será apresentado o procedimento numérico utilizado para a solução
iterativa do conjunto de equações.

4.1 Modelo para a simulação do compressor

O modelo escolhido para a simulação do compressor é relativamente simples. Tal modelo é

regido pelas seguintes equações:

 cv 
 
P   c p 
η v = 1 + c r − c r  c  (4.1)
 Pe 

Onde:
η v é o rendimento volumétrico do compressor;
c r é o fator de folga do compressor;

23
 Pc é a pressão na saída do compressor (pressão de condensação);
Pe é a pressão na entrada do compressor (pressão de evaporação);
cv é o calor específico do refrigerante para volume constante na entrada do compressor;
cp é o calor específico do refrigerante para pressão constante na entrada do compressor.

(η vVd )
mf = N (4.2)
v1

Onde:
mf é a vazão mássica de refrigerante.
Vd é o volume de refrigerante deslocado pelo compressor.
v1 é o volume específico do refrigerante na entrada do compressor.
N é a rotação do compressor.

γ −1 γ
Pc ( ) ( )
Wcomp = m f Pe v 2 [( ) γ
− 1] γ −1 (4.3)
Pe

Onde:
Wcomp é a potência de compressão politrópica;
v2 é o volume específico do refrigerante na saída do compressor;
cp
ϒ é a razão:
cv

Wcp = Wi + (1 + α )Wcomp (4.4)

Onde:
Wcp é a potência real do compressor;
Wi é a potência dissipada por perdas eletromagnéticas no motor do compressor;
α é o fator de perdas do compressor.

Wcp
h2 = h1 + (4.5)
mf

24
 Onde:
h2 e h1 são as entalpias do refrigerante na saída e na entrada do compressor, respectivamente.

Wcp
T2 = T1 + (4.6)
m f c p12

Onde:
T2 é a temperatura do refrigerante na saída do compressor;
T1 é a temperatura na entrada do compressor;
cp12 é o valor médio do calor específico a pressão constante entre a entrada e a saída do
compressor.

É interessante observarmos que em tal modelo temos como variáveis de entrada as pressões
de evaporação e condensação (entrada e saída do compressor, respectivamente), a entalpia e a
temperatura na entrada do compressor e o volume específico do refrigerante na entrada e na saída
do compressor. Para que possamos utilizar o modelo de maneira efetiva, necessitamos ainda
definir alguns parâmetros, tais como os calores específicos para pressão e volume constante, o
fator de folga do compressor, o volume de refrigerante deslocado, a rotação, o fator de perdas e
as perdas eletromagnéticas do compressor.

A partir das variáveis de entrada e dos parâmetros geométricos e operacionais, o modelo

calcula algumas variáveis intermediárias, tais como: o rendimento volumétrico e a potência de
compressão politrópica, necessários para o cálculo das variáveis de saída do modelo (potência
real do compressor entalpia e temperatura do refrigerante na saída do compressor, e vazão
mássica de refrigerante).

4.2 Modelo para a simulação dos trocadores de calor (condensador e evaporador)

O tratamento de ambos os trocadores de calor foi feito através do modelo da efetividade (ε-
NUT). Tal modelo, além de ser relativamente simples e facilmente encontrado na literatura,
oferece como principal vantagem a possibilidade de estimar a taxa real de transferência de calor

25
do trocador sem que haja um prévio conhecimento das temperaturas do fluido secundário e do
refrigerante na saída do trocador de calor.

Evaporador

As equações que regem o modelo para o evaporador são as seguintes:

C fe = m f c pre (4.7)

Onde:
Cfe é a capacidade térmica do refrigerante no evaporador;
mf é a vazão mássica de refrigerante;
cpre é o calor específico à pressão constante do refrigerante na entrada do evaporador.

C qe = m fse c par (4.8)

Onde:
Cqe é a capacidade térmica do fluido secundário no evaporador;
mfse é a vazão mássica de fluido secundário que atravessa o evaporador;
cpfs é o calor específico à pressão constante do fluido secundário na entrada do evaporador.

C min e
Ce = (4.9)
C max e

Onde:
Ce é a razão de capacidades térmicas;
Cmine é o menor valor entre Cfe e Cqe;
Cmaxe é o maior valor entre Cfe e Cqe.

U e Ae
NUTe = (4.10)
C min

Onde:
NUTe é o número de unidades de transferência térmica;

26
 Ue é o coeficiente global de transferência de calor;
Ae é a área do trocador de calor.

 1  
ε e = 1 − exp { [ ] }
( NUT )(0, 22 ) exp − C e ( NUT )(0, 78 ) − 1  (4.11)
 C e  

Onde:
εe é a efetividade do trocador do evaporador.

A expressão (4.11) foi obtida através de um procedimento semelhante ao procedimento de

obtenção da equação (3.18). Tal equação permite-nos calcular a efetividade para trocadores de
calor com fluxo cruzado, sem mistura de fluidos.

1 − exp[− NUTe (1 − C e )]
εe = (4.12)
1 − C e exp[− NUTe (1 − C e )]

A expressão (4.11) também pode ser substituída por outras correlações de efetividade,
como, por exemplo, a correlação para trocadores de calor de tubos concêntricos com escoamento
em contra corrente, descrita pela equação acima.

q max e = C min (t fs − t 4 ) (4.13)

Onde:
qmaxe é a máxima taxa de transferência de calor no evaporador;
tfs é a temperatura do fluido secundário na entrada do evaporador;
t4 é a temperatura do refrigerante na entrada do evaporador

q e = ε e q max e (4.14)

Onde:
qe é a taxa real de transferência de calor no evaporador

27
 qe
t fsse = t fs − (4.15)
m fse c pfs

Onde:
tsse é a temperatura do ar na saída do evaporador

q 
h1 =  e  + h4 (4.16)
m 
 f 

Onde:
h1 é a entalpia do refrigerante na saída do evaporador;
h4 é a entalpia do refrigerante na entrada do evaporador.

A temperatura t1 pode, então, ser obtida através de tabelas termodinâmicas, conhecendo-se

os valores da entalpia (h1) e da pressão de evaporação.

Cálculo da temperatura de saturação através do grau de superaquecimento do refrigerante

na saída do evaporador:

t1sat = t1 − ∆t (4.17)

Onde:
t1sat é a temperatura de saturação do refrigerante no evaporador;
∆t é o grau de superaquecimento do refrigerante na saída do evaporador

A partir do cálculo da temperatura de saturação, através da especificação de um grau de

superaquecimento na saída do evaporador, podemos determinar a pressão de saturação do
refrigerante, com o auxílio de tabelas termodinâmicas. A pressão de saturação, por sua vez, será
admitida como sendo igual à pressão de evaporação (ou pressão na entrada do compressor, Pe).

É importante observarmos que o modelo possui como variáveis de entrada, as temperaturas

do fluido secundário (tfs) e do refrigerante (t4) na entrada do evaporador, a vazão mássica de
refrigerante (mf) e a vazão mássica de fluido secundário que atravessa o evaporador (mfse). Para

28
que o modelo seja capaz de efetuar os cálculos, também são necessários alguns parâmetros de
entrada, tais como: o calor específico do fluido secundário (cpfs), o calor específico do
refrigerante (cpre), o coeficiente global de transferência de calor (Ue), a área do trocador de calor
(Ae) e o grau de superaquecimento do refrigerante na saída do evaporador (∆t). A partir das
variáveis de entrada e da definição dos parâmetros acima apresentados, torna-se possível a
realização dos cálculos das variáveis intermediárias do modelo ( Cqe,Cfe, NUT, qmaxe, εe). Estes
cálculos, quando realizados, permitem o cálculo das variáveis de saída do modelo, tais como: a
taxa real de transferência de calor (qe), a temperatura do refrigerante na saída do evaporador (t1),
a temperatura do fluido secundário na saída do evaporador (tfsse), a temperatura do refrigerante no
estado de vapor saturado no interior do evaporador (t1sat), a pressão de evaporação (Pe) e a
entalpia na saída do evaporador (h1). É importante notar, que a variável de saída Pe é obtida
através de tabelas termodinâmicas através da temperatura de saturação.

Condensador

As equações que regem a simulação da operação do condensador são, basicamente, as

mesmas utilizadas para simular a operação do evaporador, conforme descrevemos mais
detalhadamente a seguir:

C fc = m fsc c pfs (4.18)

Onde:
Cfc é a capacidade térmica do fluido secundário no condensador;
mfsc é a vazão mássica do fluido secundário no condensador;
cpfs é o calor específico à pressão constante do fluido secundário.

C qc = m f c prc (4.19)

Onde:
Cqc é a capacidade térmica do refrigerante no condensador;
mf é a vazão mássica do refrigerante;
cprc é o calor específico a pressão constante na entrada do condensador.

29
 C min c
Cc = (4.20)
C max c

Onde:
Cc é a razão das capacidades térmicas no condensador;
Cminc é o menor valor entre Cfc e Cqc;
Cmaxc é o maior valor entre Cfc e Cqc.

U c Ac
NUTc = (4.21)
C min c

Onde:
NUTc é o número de unidades de transferência térmica do condensador;
Uc é o coeficiente global de transferência de calor do condensador;
Ac é a área do condensador.

 1  
ε c = 1 − exp { [ ] }
( NUTc )(0, 22 ) exp − C c ( NUTc )(0, 78 ) − 1  (4.22)
 C c  

Onde:
εc é a efetividade do condensador.

A expressão (4.21) foi obtida através de um procedimento semelhante ao procedimento de

obtenção da equação (3.18). Tal equação permite-nos calcular a efetividade para trocadores de
calor com fluxo cruzado, sem mistura de fluidos.

A expressão 4.21 também pode ser substituída por outras correlações de efetividade
encontradas na literatura (Incropera F.P. & DeWitt D.P., Fundamentos de transferência de calor e
massa, 4ª Edição, Editora LTC, Capítulo 11, pp.328,329), de maneira a representar outros tipos
de configuração geométrica para o condensador.

q max c = C min c (t 2 − t fs ) (4.23)

30
 Onde:
qmaxc é a máxima taxa de transferência do calor no condensador
t2 é a temperatura do refrigerante na entrada do condensador (aproximadamente igual á
temperatura do refrigerante na saída do compressor).
tfs é a temperatura do fluido secundário na entrada do condensador.

q c = ε c q max c (4.24)

Onde:
qc é a taxa real de transferência de calor no condensador.

q 
h3 = h2 −  c  (4.25)
m 
 f 

Onde:
h3 é a entalpia do refrigerante na saída do condensador.

A temperatura t3 pode ser obtida com o auxílio de tabelas termodinâmicas, através do

conhecimento prévio da entalpia (h3) e da pressão de condensação.

qc
t fssc = t fs + (4.26)
m fsc c fsr

Onde:
tfssc é a temperatura do fluido secundário na saída do condensador.

Podemos observar que as variáveis de entrada do modelo do condensador, quando

comparadas com as do evaporador, são colocadas de maneira mito semelhante, ou seja, o modelo
possui como variáveis de entrada as temperaturas do refrigerante (t2) e do fluido secundário (tfs)
na entrada do condensador, bem como a vazão mássica de refrigerante (mf) e a vazão mássica do
fluido secundário (mfsc). Os parâmetros utilizados para o condensador também são semelhantes
aos usados para o evaporador, sendo cpfs o calor específico do fluido secundário, cpref o calor

31
específico do refrigerante, Uc o coeficiente global de transferência de calor do condensador e Ac a
área do condensador. A partir das variáveis de entrada e dos parâmetros acima descritos,
podemos efetuar os cálculos das variáveis intermediárias do modelo ( Cqc,Cfc, NUTc, qmaxc, εc).
Estes cálculos, quando realizados, permitem o cálculo das variáveis de saída do modelo, tais
como: a taxa real de transferência de calor no condensador (qc), a temperatura do refrigerante no
estado de líquido sub-resfriado na saída do condensador (t3), a temperatura do fluido secundário
na saída do condensador (tfssc) e a entalpia na saída do condensador (h3).

4.3 Modelo para a simulação da válvula de expansão

O modelo utilizado para a simulação do dispositivo de expansão admite que o processo de

expansão é isentálpico.

h3 = h4 (4.27)

Onde:
h3 é a entalpia do refrigerante na entrada da válvula de expansão
h4 é a entalpia do refrigerante na saída da válvula de expansão

A determinação da pressão de condensação é feita através do modelamento da válvula de

expansão descrito pela equação (3.19). Através da especificação de um valor para a constante
característica da válvula (K), podemos determinar o valor da pressão de condensação.

A vazão mássica do refrigerante é relacionada com a diferença de pressão na válvula

através da equação:

mf = K ∆P * ρ (4.28)

Onde:
K é uma constante característica do dispositivo de expansão;
∆P é a diferença entra as pressões de entrada e saída da válvula;
ρ é a densidade do refrigerante na entrada da válvula.

32
 Através da equação acima, pode-se determinar a pressão de condensação (Pc), uma vez que
os valores da pressão de evaporação (Pe), constante da válvula (K), vazão mássica de refrigerante
e densidade do refrigerante na entrada da válvula de expansão (ρ) são previamente calculados
pelo programa.

O modelo numérico apresentado tem como variáveis de entrada as temperaturas iniciais do

refrigerante na entrada do compressor (t1i) e na saída do condensador (t3i). A partir da estimativa
inicial de tais variáveis e da parametrização de cada componente do sistema, o programa inicia o
processo iterativo de solução.

As principais variáveis de saída do programa são as temperaturas do refrigerante na sucção

do compressor (t1), descarga do compressor (t2), saída do condensador (t3) e descarga da válvula
de expansão (t4), as pressões de condensação (Pc) e de evaporação (Pe), as potências térmicas
trocadas no condensador (qc) e no evaporador (qe), a potência consumida pelo compressor (Wcp) e
as temperaturas do fluido secundário na saída do condensador (tfsc) e do evaporador (tfse).

4.4 Algoritmo para a resolução do modelo matemático no programa EES

Para a construção do algoritmo, foi definido um conjunto de equações representando os

modelos de cada componente do sistema, conforme descrito nos itens anteriores.

Uma vez definido o conjunto de equações, deve-se especificar o fluido de trabalho (R-134a,
por exemplo) e o conjunto de parâmetros operacionais e geométricos do sistema de refrigeração.
Tais parâmetros devem ser especificados para cada sistema de refrigeração diferente que possa
vir a ser simulado. Nas tabelas a seguir, podemos observar os parâmetros operacionais e
geométricos a serem especificados para cada um dos componentes do sistema de refrigeração.

33
Tabela 4.1 - Parâmetros operacionais e geométricos do compressor

PARÂMETROS DO COMPRESSOR

N Rotação
Wi Potência dissipada
Vd Volume deslocado pelo pistão
cr Fator de folga
α Fator de perdas

Tabela 4.2 - Parâmetros operacionais e geométricos do condensador

PARÂMETROS DO CONDENSADOR

Ac Área
Uc Coeficiente Global de Transferência de Calor
Vazão mássica do fluido secundário que
mfsc
atravessa o condensador
tfs Temperatura do fluido secundário

Tabela 4.3 - Parâmetros operacionais e geométricos do evaporador

PARÂMETROS DO EVAPORADOR

Ae Área
Ue Coeficiente Global de transferência de calor
Vazão mássica do fluido secundário que
mfse
atravessa o evaporador
Grau de superaquecimento na saída do
∆te
evaporador
tfs Temperatura do fluido secundário

Tabela 4.4 - Parâmetros operacionais e geométricos da válvula de expansão

PARÂMETROS DA VÁLVULA DE EXPANSÃO

K Constante da válvula de expansão

Uma vez definido o conjunto de parâmetros operacionais e geométricos a serem utilizados

nos modelos de cada um dos componentes do sistema de refrigeração, existe a necessidade da

34
estimativa inicial das variáveis de entrada para que o programa possa realizar os primeiros
cálculos. O conjunto de variáveis de entrada do programa a serem especificadas pode ser
visualizado na tabela a seguir.

Tabela 4.5 - Variáveis de entrada

VARIÁVEIS DE ENTRADA

Temperatura inicial do refrigerante na saída do

t1,i
evaporador
Temperatura inicial do refrigerante na saída do
t3,i
condensador

A partir da estimativa inicial, o programa realiza uma seqüência de cálculos preliminares

das entalpias para cada ponto do sistema, das pressões de condensação e de evaporação e de mais
alguns outros cálculos de propriedades do refrigerante como, por exemplo, o volume específico e
o título na entrada do compressor. Tais cálculos preliminares foram feitos baseados em funções
internas do próprio programa (E.E.S.) e no comportamento do ciclo padrão de funcionamento de
um ciclo frigorífico (compressão isentrópica, expansão isentálpica). Vale a pena, neste ponto,
ressaltar, que os valores obtidos através destes cálculos preliminares são apenas estimativas para
iniciar os cálculos.

A partir da segunda iteração e nas outras iterações seguintes, os valores das pressões e
entalpias, entre outros, são atualizados a partir dos valores obtidos como resultados da iteração
anterior.

O critério de parada do programa para a obtenção dos resultados finais da simulação foi
feito através de uma análise comparativa da variação dos valores obtidos para as pressões de
condensação e evaporação a cada iteração. Tal critério compara a variação dos valores das
pressões calculadas entre duas iterações consecutivas com um valor pré-estabelecido. Quando a
variação dos valores das pressões de evaporação e de condensação é menor que o valor pré-
estabelecido, o programa é finalizado e são apresentados os resultados finais.

O algoritmo com a seqüência de resolução das equações, bem como a ordem com a qual
elas foram resolvidas, será demonstrada a seguir.

35
Seqüência de equações para a simulação do modelo

Compressor

Inicialmente, o programa calcula o rendimento volumétrico do compressor através da

equação:

 cv 
 
P   c p 
η v = 1 + c r − c r  c  (4.29)
 Pe 

Para o cálculo do rendimento volumétrico do compressor, são necessários o fator de folga,

as pressões de condensação e de evaporação e os calores específicos para volume e para pressão
constante na entrada do evaporador. Tanto as pressões como os calores específicos são estimados
através de cálculos preliminares para a primeira iteração e são substituídos a cada iteração pelos
valores obtidos na iteração anterior.

Após o cálculo do rendimento volumétrico, podemos calcular a vazão mássica de

refrigerante, partindo-se do rendimento volumétrico (equação 4.29), volume deslocado, rotação
do compressor e do volume específico do refrigerante na entrada do compressor.

(η vVd )
mf = N (4.30)
v1

Na equação acima, o volume deslocado e a rotação são parâmetros do compressor e o

volume específico na entrada do compressor é estimado através do cálculo preliminar para a
primeira iteração e substituído a cada iteração pelo valor obtido na iteração anterior.

A potência de compressão é calculada supondo-se um processo politrópico partindo-se da

vazão mássica de refrigerante (equação 4.30), pressão de evaporação, pressão de condensação,
volume específico do refrigerante na saída do compressor e de γ que nada mais é do que a razão
entre os calores específicos à pressão constante e ao volume constante na entrada do compressor.

γ −1 γ
P ( ) (
γ −1
)
Wcomp = m f Pe v 2 [( c ) γ
− 1] (4.31)
Pe

36
 O cálculo da potência real do compressor, leva em consideração as perdas eletromagnéticas
no motor do compressor e o fator de perdas mecânicas (folgas, perdas nas válvulas, entre outras).

Wcp = Wi + (1 + α )Wcomp (4.32)

Na equação acima, tanto as perdas eletromagnéticas como o fator de perdas, são parâmetros
do compressor.

O cálculo da entalpia na saída do compressor é feito a partir da entalpia na entrada do

compressor, da potência do compressor (equação 4.32) e da vazão mássica de refrigerante
calculada (equação 4.30).

Wcp
h2 = h1 + (4.33)
mf

A entalpia na entrada do compressor é estimada através dos cálculos preliminares para a

primeira iteração e substituída pelo valor obtido na iteração anterior a cada nova iteração.

A temperatura na saída do compressor (t2), estimada através de cálculos preliminares para a

primeira iteração, é calculada com base na entalpia na saída do compressor (equação 4.33) e na
pressão de condensação, através de funções internas do próprio E.E.S, para as iterações seguintes.

Condensador

O cálculo da capacidade térmica do condensador (Cfc ) é dado através da vazão mássica do

fluido secundário que atravessa o condensador (parâmetro do condensador) e do calor específico
à pressão constante do fluido secundário.

C fc = m fsc c pfs (4.34)

37
 O calor específico à pressão constante do fluido secundário é calculado pelo programa
através da temperatura do fluido secundário (parâmetro) através de funções internas.

Para o cálculo da capacidade térmica do refrigerante no condensador (Cqc) são necessários:

a vazão mássica do refrigerante (equação 4.30) e o calor específico à pressão constante do
refrigerante na entrada do condensador.

C qc = m f c prc (4.35)

O calor específico do refrigerante à pressão constante na entrada do condensador é

calculado através de funções internas do programa baseando-se na pressão de condensação e na
temperatura do refrigerante na entrada do condensador.

Os valores de Cmin e Cmax são determinados através de funções de mínimo e máximo do

programa. Para tal, são comparados os valores das constantes Cfc e Cqc, sendo que o menor valor
obtido é atribuído à constante Cminc e o maior valor é atribuído à constante Cmaxc.

Partindo da determinação de tais valores, calcula-se a razão de capacidades térmicas Cc,

conforme a equação abaixo.

C min c
Cc = (4.36)
C max c

O cálculo do número de unidades de transferência térmica do condensador (NUTc) é feito

baseado no coeficiente global de transferência de calor, na área do condensador (parâmetros
operacionais e geométricos) e no valor de Cminc.

U c Ac
NUTc = (4.37)
C min c
O cálculo da efetividade do condensador é realizado a partir do valor de Cc (equação 4.36) e
do número de unidades de transferência térmica do condensador, NUTc (equação 4.37).

38
  1  
ε c = 1 − exp { [ ] }
( NUTc )(0, 22 ) exp − C c ( NUTc )(0, 78 ) − 1  (4.38)
 C c  

É importante ressaltar, a este ponto do trabalho, que a equação 4.38 descreve o cálculo da
efetividade para um trocador de calor com fluxo cruzado entre o refrigerante e o ar e convecção
forçada.

O cálculo da efetividade pode ser realizado para diferentes tipos de trocadores de calor,
como, por exemplo, trocadores de tubos concêntricos com fluxo cruzado; bastando, apenas, a
substituição da equação 4.38 no corpo do programa pela equação adequada ao tipo de trocador a
ser simulado.

O cálculo da máxima taxa de transferência de calor no condensador é realizado partindo-se

da constante Cminc e das temperaturas do refrigerante na saída do compressor e do fluido
secundário (parâmetro operacional).

q max c = C min c (t 2 − t fs ) (4.39)

A taxa real de transferência de calor no condensador é calculada partindo-se da efetividade

(equação 4.38) e da máxima taxa de transferência de calor no condensador (equação 4.39).

q c = ε c q max c (4.40)

O cálculo da entalpia do refrigerante na saída do condensador é realizado a partir da

entalpia do refrigerante na saída do compressor (equação 4.33), da taxa real de transferência de
calor no condensador (equação 4.40) e da vazão mássica de refrigerante (equação 4.30).

q 
h3 = h2 −  c  (4.41)
m 
 f 

A temperatura de saída do fluido secundário do condensador é calculada partindo-se da

temperatura do fluido secundário (parâmetro), da taxa real de transferência de calor no

39
condensador (equação 4.40) da vazão mássica do fluido secundário que atravessa o condensador
(parâmetro) e do calor específico do fluido secundário.

qc
t fssc = t fs + (4.42)
m fsc c pfs

A temperatura do refrigerante na saída do condensador é estimada através de cálculos

preliminares para a primeira iteração e calculada através de funções internas do programa E.E.S.
com base na pressão de condensação e na entalpia do refrigerante na saída do condensador para
as iterações seguintes.

Evaporador:

O cálculo da capacidade térmica do refrigerante no evaporador (Cfe ) é realizado partindo-se

da vazão mássica de refrigerante (equação 4.30) e do calor específico do refrigerante à pressão
constante na entrada do evaporador.

C fe = m f c pre (4.43)

O calor específico do refrigerante à pressão constante na entrada do evaporador é calculado

através de funções internas do programa baseando-se na pressão de evaporação e considerando
que o refrigerante entra no estado líquido saturado na entrada do trocador de calor.

A capacidade térmica do fluido secundário no evaporador (Cqe) é calculada a partir da

vazão mássica do fluido secundário que atravessa o evaporador (parâmetro) e do calor específico
à pressão constante do fluido secundário.

C qe = m are c par (4.44)

A razão de capacidades térmicas Ce é calculada a partir dos valores de Cmine e Cmexe, as

quais são determinadas de maneira semelhante à já descrita para o condensador.

40
 C min e
Ce = (4.45)
C max e

O cálculo do número de unidades de transferência térmica para o evaporador é realizado

partindo-se do coeficiente global de transferência de calor, da área do evaporador (parâmetros) e
da constante Cmine.

U e Ae
NUTe = (4.46)
C min

A partir do cálculo do número de unidades de transferência térmica do evaporador e da

constante Ce podemos calcular a efetividade.

 1  
ε e = 1 − exp { [ ] }
( NUT )(0, 22 ) exp − C e ( NUT )(0, 78 ) − 1  (4.47)
 C e  

Assim como ocorre para o cálculo da efetividade condensador, a equação acima é utilizada
para a simulação de um evaporador com fluxo cruzado e convecção forçada.

Para a simulação de outros tipos de evaporador, é necessária a substituição no corpo do

programa da equação (4.47) pela equação da efetividade pertinente ao tipo de evaporador a ser
simulado.

O cálculo da máxima taxa de transferência de calor no condensador é dado a partir do valor

de Cmine e das temperaturas do refrigerante na entrada do evaporador e do fluido secundário
(parâmetro).

q max e = C min (t fs − t 4 ) (4.48)

A temperatura do refrigerante na entrada do evaporador é estimada através de cálculos

preliminares para a primeira iteração e calculada com base na entalpia na entrada do evaporador e
na pressão de evaporação, através de funções internas do programa E.E.S., para as iterações
seguintes.

41
 O cálculo da taxa real de transferência de calor no evaporador é dado de maneira
semelhante à descrita para o condensador, partindo-se da efetividade (equação 4.47) e da máxima
taxa de transferência de calor no evaporador (equação 4.48).

q e = ε e q max e (4.49)

A temperatura do fluido secundário na saída do evaporador é calculada com base na vazão

mássica do fluido secundário que atravessa o evaporador (parâmetro), no calor específico do
fluido secundário à pressão constante, na temperatura do fluido secundário (parâmetro) e na taxa
real de transferência de calor no evaporador (equação 4.49).

qe
t fsse = t fs − (4.50)
m fse c pfs

A entalpia do refrigerante na entrada do compressor é calculada partindo-se da taxa real de

transferência de calor no evaporador (equação 4.49), da vazão mássica de refrigerante (equação
4.30) e da entalpia do refrigerante na saída do condensador (equação 4.41).

q 
h1 =  e  + h4 (4.51)
m 
 f 

A entalpia na entrada do evaporador é estimada através de cálculos preliminares para a

primeira iteração e substituída pelo valor obtido na iteração anterior para as iterações seguintes.

A temperatura de saturação do refrigerante é calculada baseando-se na temperatura na

entrada do compressor e no grau de superaquecimento do refrigerante na entrada do compressor
(parâmetro).

t1sat = t1 − ∆t (4.52)

42
 A temperatura do refrigerante na entrada do compressor é calculada através de funções
internas do programa E.E.S., baseando-se na entalpia do refrigerante na entrada do compressor e
na pressão de evaporação.

Válvula de expansão

A válvula de expansão é admitida isoentálpica, portanto, o cálculo da entalpia na saída da

válvula depende apenas da entalpia do refrigerante na saída do condensador (equação 4.41).

h3 = h4 (4.53)

O cálculo da pressão de condensação é feito através da pressão de evaporação (determinada

através do grau de superaquecimento na saída do evaporador), do volume específico do
refrigerante na entrada da válvula (determinado através de funções internas do programa E.E.S.),
na vazão mássica de refrigerante (equação 4.30) e na constante da válvula de expansão
(parâmetro).

mf = K ∆P * ρ (4.54)

Partindo-se do algoritmo inicial descrito acima, foram possíveis as primeiras simulações

com o programa, feitas ainda de modo manual, calculado-se apenas uma iteração por vez, a fim
da obtenção dos primeiros resultados, correção de possíveis erros e verificação da convergência
do programa.

A partir da determinação dos parâmetros iniciais (temperatura na sucção do compressor,

grau de superaquecimento na saída do evaporador, temperatura do fluido secundário, áreas dos
trocadores de calor, coeficientes globais de transferência de calor para os trocadores de calor,
rotação do compressor e fator de folga do compressor) foi possível a obtenção dos primeiros
resultados (pressão de condensação, pressão de evaporação, calor trocado no condensador, calor
trocado no evaporador, temperatura de saída do fluido secundário do condensador, temperatura
de saída do fluido secundário do evaporador, efetividade do condensador e do evaporador, entre
outros). O primeiro conjunto de resultados foi novamente inserido como dados de entrada no

43
programa (ainda de maneira manual, executando-se apenas uma iteração por vez), de maneira a
obter-se um segundo conjunto de resultados e assim sucessivamente. Através de tal
procedimento, tornou-se possível um monitoramento do comportamento do programa. Através do
monitoramento de alguns resultados, tais como: pressões e temperaturas nos diversos pontos do
sistema, foi possível observar se o programa tendia a convergir ou divergir dos valores
apresentados a cada iteração. Tal verificação de convergência foi feita através da comparação dos
resultados obtidos a cada nova iteração com os resultados obtidos nas iterações anteriores.

Após a verificação manual da convergência do programa, o mesmo foi automatizado,

através de rotinas de programação, que fazem com que o mesmo realize a mesma seqüência de
cálculos diversas vezes sempre se baseando nos valores obtidos como resultados na iteração
anterior à que está sendo calculada, até que a convergência seja obtida, através de um critério de
parada previamente estabelecido.

O fluxograma descrito abaixo permite uma visualização geral do funcionamento do

programa.

44
 Parâmetros de Entrada

Cálculos Preliminares

Procedimentos Iterativos

Verificação da
Convergência

Variação dos resultados Variação dos resultados

entre duas iterações entre duas iterações
consecutivas menor que o consecutivas maior que o
critério de convergência critério de convergência

Parada dos
Procedimentos
Iterativos

Apresentação
dos resultados
Finais

Para uma análise mais detalhada do programa segue, em anexo, no final deste trabalho, uma
listagem completa das rotinas de programação utilizadas.

45
Capítulo 5

Resultados Numéricos

Após a verificação das equações e do algoritmo para a resolução do modelo numérico, foi
realizada a validação do modelo através da comparação dos resultados obtidos numericamente
com resultados experimentais . Tal validação é de vital importância para a verificação da
confiabilidade do modelamento utilizado.

Dada a validação do programa, foram feitos dois estudos de casos, para a verificação da
influência de cada parâmetro sobre o comportamento do sistema e otimização dos sistemas de
refrigeração propostos para cada caso.

Tanto a validação do modelo numérico quanto os estudos de caso propostos serão descritos
mais detalhadamente ao longo deste capítulo.

5.1 Validação do programa através da comparação com dados experimentais

A validação do programa foi feita através da comparação dos resultados obtidos em testes
na bancada frigorífica experimental da Faculdade de Engenharia Mecânica da Universidade
Estadual de Campinas (Unicamp), com os resultados obtidos no programa.

Para tal comparação, o programa foi parametrizado de maneira a simular os equipamentos

que compõe a bancada experimental. Tal bancada é composta por um compressor alternativo de
um cilindro, com potência de ¼ de hp, condensador composto por tubos aletados, de fluxo
cruzado com ar em convecção forçada, válvula de expansão do tipo eletrônica e evaporador
composto por tubos concêntricos de fluxo em contracorrente com água (fluido secundário),

47
operando com fluido R-134a. Para tal sistema, o grau de superaquecimento do refrigerante na
sucção do compressor é mantido constante através do controle eletrônico da válvula de expansão.

5.1.1 Parâmetros operacionais e geométricos da bancada experimental

Os parâmetros operacionais e geométricos da bancada experimental foram definidos com

base em catálogos técnicos dos equipamentos que a compõe e medições experimentais realizadas
previamente a este trabalho.

A tabela abaixo mostra os parâmetros geométricos e operacionais de cada componente do

sistema e os valores médios obtidos de um conjunto de dados experimentais para os coeficientes
globais de transferência de calor do evaporador e do condensador. A simplificação da adoção de
valores médios para os coeficientes globais de transferência de calor no evaporador e no
condensador pode ser explicada pela pouca variação dos valores destes para as faixas de
temperaturas nas quais os ensaios experimentais foram realizados.

É importante, ainda, ressaltar que os parâmetros descritos na tabela a seguir foram mantidos
constantes para todas as simulações numéricas feitas para a validação do programa.

Tabela 5.1 - Parâmetros do sistema de refrigeração

Parâmetros do compressor

Wi (W) Vd (m3) cr α
0,0968 7,95 x 10-6 0,05 0,178
Parâmetros do evaporador

Ae (m2) Ue ( kW/m2°C) ∆te (ºC) Ae (m2)

0,025 1,52 6,5 0,025
Parâmetros do condensador

Ac (m2) Uc ( kW/m2°C) mar (kg/s) Ac (m2)

0,1214 0,0470 0,188 0,1214
Parâmetro da válvula de expansão

K
0,14435 x 10-6

48
 A tabela a seguir, mostra o conjunto de parâmetros do sistema que sofreram variações em
cada ensaio experimental. Os fluidos secundários do evaporador e do condensador são água e ar,
respectivamente.

Tabela 5.2 - Parâmetros do sistema de refrigeração

Compressor Evaporador Condensador

Ensaio N (rpm) magua (kg/s) tagua (ºC) tar (ºC)

1 3415 0,0257 25,0 27,0

2 3923 0,0515 26,3 28,4
3 3635 0,0515 29,7 32,5
4 2073 0,0223 31,6 34,1
5 2349 0,0223 28,8 30,5

5.1.2 Estimativas iniciais para o sistema de refrigeração

Para dar início à simulação numérica, para cada ponto experimental foram escolhidos os
valores iniciais das temperaturas do refrigerante na sucção do compressor (t1) e na saída do
condensador (t3).

Tabela 5.3 - Temperaturas iniciais na saída do evaporador e na saída do condensador

Ensaio t1 (ºC) t3 (ºC)

1 3,10 49,10
2 2,80 57,90
3 3,80 58,80
4 3,50 56,30
5 1,70 45,90

5.1.3 Comparação entre os resultados numéricos e experimentais

As tabelas e gráficos a seguir fornecem dados para a análise comparativa dos dados de
temperatura do refrigerante (R-134a) calculados pelo programa em relação aos dados obtidos

49
experimentalmente. As temperaturas analisadas são as dos quatro principais pontos do sistema:
saída do evaporador (t1), descarga do compressor (t2), saída do condensador (t3) e descarga da
válvula de expansão (t4).

Tabela 5.4 - Temperaturas na entrada do compressor

Bancada Frigorífica Programa

Ensaio
t1 (ºC) t1 (ºC)
1 3,10 5,20
2 2,80 4,02
3 3,80 5,77
4 3,50 4,92
5 1,70 2,88

Temperaturas na Entrada do Compressor

7,00

6,00

5,00

4,00
Experimental
ºC

Programa
3,00

2,00

1,00

0,00
0 1 2 3 4 5 6
Ensaio

Figura 5.1 - Temperaturas na entrada do compressor

50
 Tabela 5.5 - Temperaturas do refrigerante na descarga do compressor
Bancada Frigorífica Programa
Ensaio
t2 (ºC) t2 (ºC)
1 49,10 57,39
2 57,90 65,01
3 58,80 66,50
4 56,30 59,70
5 45,90 50,10

Temperaturas na descarga do compressor

70,0

60,0

50,0

40,0
Experimental
ºC

Programa
30,0

20,0

10,0

0,0
0,0 1,0 2,0 3,0 4,0 5,0 6,0
Ensaio

Figura 5.2 - Temperaturas do refrigerante na descarga do compressor

51
 Tabela 5.6 - Temperaturas do refrigerante na saída do evaporador
Bancada Frigorífica Programa
Ensaio
t3 (ºC) t3 (ºC)
1 39,90 48,87
2 43,30 56,32
3 44,30 58,37
4 42,60 54,27
5 37,20 46,87

Temperaturas na saída do condensador

70,0

60,0

50,0

40,0
Experimental
ºC

Programa
30,0

20,0

10,0

0,0
0 1 2 3 4 5 6
Ensaio

Figura 5.3 - Temperaturas na saída do condensador

52
 Tabela 5.7 - Temperaturas do refrigerante na descarga da válvula de expansão
Bancada Frigorífica Programa
Ensaio
t4 (ºC) t4 (ºC)
1 -3,50 -3,18
2 -4,10 -3,13
3 -3,10 -3,47
4 -3,30 -2,49
5 -4,50 -3,28

Temperaturas na descarga da válvula

0,00

-0,50

-1,00

-1,50

-2,00
Experimental
ºC

-2,50
Programa
-3,00

-3,50

-4,00

-4,50

-5,00
0 1 2 3 4 5 6
Ensaio

Figura 5.4 - Temperaturas do refrigerante na descarga da válvula de expansão

A análise dos gráficos e tabelas fornecidas acima permite uma comparação dos dados de
temperatura para os quatro principais pontos do sistema. Podemos notar, que os dados calculados
pelo programa tendem a seguir o comportamento dos dados obtidos experimentalmente. As
variações de valores das temperaturas calculadas em relação às obtidas experimentalmente
devem-se a simplificações presentes nos modelos numéricos utilizados para cada componente do
ciclo e a pequenas imprecisões dos equipamentos de medição presentes na bancada frigorífica
experimental.

53
 As tabelas a seguir fornecem dados de pressões de condensação e evaporação obtidos
experimentalmente e calculados pelo programa para cada um dos seis ensaios realizados na
bancada frigorífica.

Tabela 5.8 - Pressões de condensação

Bancada Frigorífica Programa
Ensaio
Pcb (Mpa) Pc (Mpa)
1 1,04 1,28
2 1,14 1,54
3 1,17 1,65
4 1,17 1,60
5 0,96 1,18

Pressões de Condensação

1,8

1,6

1,4

1,2

1,0 Experimental
MPa

0,8 Programa

0,6

0,4

0,2

0,0
0 1 2 3 4 5 6
Ensaio

Figura 5.5 - Pressões de condensação

54
 Tabela 5.9 - Pressões de evaporação
Bancada Frigorífica Programa
Ensaio
Peb (Mpa) Pe (Mpa)
1 0,27 0,26
2 0,26 0,26
3 0,27 0,28
4 0,27 0,27
5 0,26 0,26

Pressões de Evaporação

0,29

0,29

0,28

0,28
Experimental
MPa

0,27
Programa
0,27

0,26

0,26

0,25
0 1 2 3 4 5 6
Ensaio

Figura 5.6 - Pressões de evaporação

A análise das tabelas e gráficos acima mostra que a concordância das pressões calculadas
pelo programa em relação às obtidas experimentalmente é semelhante à concordância da
temperaturas. As pressões calculadas tendem a variar em cada ponto nas mesmas proporções que
as obtidas experimentalmente. Nota-se ainda, que a variação entre as pressões de evaporação
calculadas e experimentais é bem reduzida, quando comparada à variação das pressões de
condensação em cada ensaio.

55
 A tabela e o gráfico a seguir, fornecem dados experimentais e calculados pelo programa
para a comparação dos calores trocados no evaporador.

Tabela 5.10 - Calores trocados no Evaporador

Bancada Frigorífica Programa
Ensaio
qe (kW) qe (kW)
1 0,52 0,58
2 0,66 0,77
3 0,63 0,67
4 0,60 0,63
5 0,21 0,23

Calor trocado no evaporador

0,90

0,80

0,70

0,60

0,50 Experimental
kW

0,40 Programa

0,30

0,20

0,10

0,00
0 1 2 3 4 5 6
Ensaio

Figura 5.7 - Calor trocado no evaporador

Analisando a tabela e o gráfico apresentados acima, podemos notar que os valores obtidos
para os calores trocados no evaporador nos ensaios experimentais são inferiores aos calculados
pelo programa. Tal variação pode ser explicada por simplificações que os modelos de efetividade

56
para trocadores de calor adotam. No entanto, nota-se que a variação entre os resultados obtidos
através do modelo e os resultados experimentais não é grande, e o modelo tende fortemente a
seguir o comportamento experimental.

A tabela e o gráfico abaixo fornecem dados para a análise comparativa dos valores
calculados pelo programa com os obtidos experimentalmente.

Tabela 5.11 - Potência consumida pelo compressor

Bancada Frigorífica Programa
Ensaio
Wcp (W) Wcp (W)
1 0,26 0,31
2 0,32 0,37
3 0,29 0,35
4 0,28 0,32
5 0,15 0,17

Potência Consumida pelo Compressor

0,40

0,35

0,30

0,25
Experimental
kW

0,20
Programa
0,15

0,10

0,05

0,00
0 1 2 3 4 5 6
Ensaio

Figura 5.8 - Potência consumida pelo compressor

57
 A análise da tabela e do gráfico acima, permite observar que a potência calculada pelo
programa sofre pouca variação em relação aos dados obtidos experimentalmente para cada ponto
ensaiado. Para as simulações de cada um dos cinco ensaios, foi adotada a hipótese simplificadora
da média dos valores obtidos experimentalmente nos ensaios.

A tabela abaixo fornece dados experimentais e calculados pelo programa para a comparação
do COP do ciclo

Tabela 5.12 - COP do sistema de refrigeração

Bancada Frigorífica Programa E.E.S.
Ensaio
COP COP
1 2,10 1,87
2 2,10 2,05
3 2,10 1,92
4 2,20 2,00
5 1,40 1,34

COP do sistema

2,50

2,00

1,50
Seqüência1
COP

Seqüência2
1,00

0,50

0,00
0 1 2 3 4 5 6
Ensaio

Figura 5.9 - COP do sistema

58
 Podemos notar que o COP calculado pelo programa segue a mesma tendência dos obtidos
experimentalmente. As variações dos valores calculados para os obtidos experimentalmente pode
ser explicada pela variação dos valores calculados para o calor trocado no evaporador e para a
potência consumida pelo compressor em relação aos obtidos experimentalmente.

Após a validação do programa através da comparação dos dados obtidos com os dados
experimentais, foram definidos dois estudos de caso, para a verificação da influência da variação
dos principais parâmetros (tais como: área do condensador, área do evaporador, vazão mássica de
ar que atravessa o condensador, vazão mássica de ar que atravessa o evaporador, grau de
superaquecimento na saída do evaporador, grau de sub-resfriamento na saída do condensador e
temperatura do ar ambiente).

Após a verificação da influência dos principais parâmetros para o primeiro estudo de caso
(otimização e projeto de uma câmara frigorífica) foi definido um novo conjunto de parâmetros,
mas desta vez para a análise da simulação de um sistema de ar condicionado. Uma vez definido o
novo conjunto de parâmetros para o sistema de ar condicionado, foram verificadas as influências
dos novos principais parâmetros sobre o comportamento do sistema na simulação.

5.2 Estudo de caso: projeto e otimização de uma câmara frigorífica industrial

O sistema de refrigeração idealizado para ser simulado neste estudo de caso é uma câmara
frigorífica industrial, com capacidade de refrigeração de 18,5 KW, composta por um compressor
alternativo, condensador de fluxo cruzado com ar em convecção forçada, evaporador de fluxo
cruzado com ar e convecção forçada, e válvula de expansão. É interessante observarmos que o
modelamento adotado para a válvula de expansão independe do tipo de equipamento utilizado,
pois leva em consideração apenas as condições do refrigerante na entrada e na saída do
equipamento. O sistema opera com o refrigerante R-134a. A escolha do refrigerante para a
operação do sistema (R-134a) foi feita a partir da observação da grande utilização deste
refrigerante atualmente. O programa também nos permite a execução de simulações com outros
fluidos de trabalho, bastando a alteração da nomenclatura do fluido a ser utilizado no corpo do
programa.

59
 A princípio, existe a necessidade da definição de alguns parâmetros e estimativas iniciais
para a simulação. Tais parâmetros e estimativas foram definidos para um caso-teste, a partir do
qual foram obtidos os primeiros resultados numéricos do modelo.

5.2.1 Parâmetros geométricos e operacionais da câmara frigorífica industrial

Parâmetros iniciais para o compressor:

Wi = 0.
N = 3000 rpm.
Vd = 0,0005 m3.
cr = 0,05.
α = 0.

Parâmetros iniciais para o condensador:

Ac= 40 m2
Uc = 0,05 KW/m2°C
marc = 2,0 Kg/s
tar = 20 ºC

Parâmetros iniciais para o evaporador:

Ae = 20 m2
Uc = 0,05 KW/m2°C
mare = 1,0 Kg/s
∆te = 5 ºC
tar = 20 ºC

Parâmetro inicial para a válvula de expansão:

K=0,0000042

60
5.2.2 Estimativas iniciais para a câmara frigorífica industrial

Temperatura inicial na saída do evaporador:

t1i = -5 °C

Temperatura inicial na saída do condensador:

t3i= 50 ºC

5.2.3 Resultados obtidos na simulação do estudo de caso da câmara frigorífica industrial

Os principais resultados obtidos na simulação deste estudo de caso, com a utilização dos
parâmetros descritos acima serão demonstrados na tabela a seguir:

Tabela 5.13 - Resultados numéricos

Pc COP
mf (Kg/s) Pe (Mpa) qc (KW) qe (KW) Wcp (KW)
(Mpa)
0,15 1,33 0,16 25,78 18,42 7,29 2,50

5.2.4 Análise da influência dos principais parâmetros nos resultados obtidos na câmara
frigorífica industrial

A partir dos resultados obtidos para o caso-teste, pôde-se fazer a análise comparativa dos
efeitos da variação dos principais parâmetros iniciais nos resultados.

É importante ressaltar, que as análises foram feitas para a variação de apenas um parâmetro
em questão, sendo que nos demais parâmetros foram mantidos os valores utilizados para a
simulação do caso teste.

Os efeitos da variação da área do condensador sobre os resultados podem ser visualizados

através da tabela abaixo:

61
Tabela 5.14 - Análise do efeito da variação da área do condensador
Ac mf Pc Pe qc qe Wcp
COP
(m2) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
60,0 0,14 1,14 0,15 25,58 18,83 6,68 2,57
55,0 0,14 1,18 0,15 25,48 18,77 6,79 2,63
50,0 0,14 1,22 0,15 25,73 18,66 6,94 2,64
45,0 0,15 1,27 0,16 25,95 18,56 7,11 2,62
40,0 0,15 1,33 0,16 25,78 18,42 7,29 2,50
35,0 0,15 1,42 0,16 25,73 18,24 7,57 2,43

Analisando os dados fornecidos pela tabela acima, pode-se observar que há uma queda na
pressão de condensação com o aumento da área do condensador. Tal queda deve-se ao fato de
que com o aumento da área para a troca de calor, torna-se possível a rejeição de calor em uma
temperatura mais baixa e conseqüentemente a uma menor pressão de condensação. A queda de
pressão diminui a diferença entre as pressões de condensação e evaporação, diminuindo a
potência necessária para a compressão do refrigerante (R-134a), e conseqüentemente,
melhorando o COP do sistema.

Observa-se ainda, através dos dados fornecidos acima, que o aumento da área do
condensador não exerce grande influência sobre o calor trocado no condensador e no evaporador,
e que a vazão mássica de refrigerante deslocada pelo compressor permanece praticamente
inalterada.

Os efeitos da variação da área do evaporador sobre os resultados obtidos no caso-teste

padrão podem ser analisados através da tabela abaixo:

62
Tabela 5.15 - Análise da influência da variação da área do evaporador
Ae mf Pc Pe qc qe Wcp
COP
(m2) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
40,0 0,19 1,58 0,21 32,26 23,05 9,33 2,50
35,0 0,19 1,54 0,20 31,14 22,29 8,98 2,51
30,0 0,18 1,49 0,19 29,70 21,26 8,54 2,52
25,0 0,16 1,42 0,18 27,88 19,97 7,97 2,52
20,0 0,15 1,33 0,16 25,78 18,42 7,29 2,50
15,0 0,13 1,23 0,14 22,85 16,35 6,46 2,51

Uma análise dos dados fornecidos pela tabela acima, permite constatar que a variação da
área do evaporador não exerce uma influência significativa sobre o COP do sistema. No entanto,
pode-se notar um aumento nas pressões de condensação e evaporação, nos calores trocados pelo
condensador e pelo evaporador, na potência de compressão e na vazão mássica de refrigerante.
Como a área do evaporador é aumentada gradativamente até chegar ao valor de 100% da área
original, a efetividade do trocador aumenta, e conseqüentemente o calor trocado aumenta. O
aumento da pressão de evaporação provoca um aumento no volume específico do refrigerante.
Tal aumento de pressão, também favorece a sucção do compressor, que passa a deslocar uma
maior vazão mássica de refrigerante, porém, o trabalho de compressão também aumenta, devido
ao aumento na pressão de condensação e o conseqüente aumento da resistência à descarga do
refrigerante na linha de condensação, exigindo uma maior potência do compressor.

A análise deste comportamento permite concluir que o aumento da área do evaporador

implica em um aumento de todo o sistema de refrigeração, mantendo as mesmas características
iniciais.

Os efeitos da variação da vazão mássica de ar que atravessa o condensador sobre os

resultados obtidos no caso-teste padrão podem ser analisados através da tabela abaixo.

63
Tabela 5.16 - Análise da influência da variação da vazão mássica de ar que atravessa o
condensador
marc mf Pc Pe qc qe Wcp
COP
(Kg/s) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
4,0 0,14 1,21 0,15 25,52 18,71 6,88 2,75
3,5 0,15 1,22 0,16 25,60 18,67 6,92 2,69
3,0 0,15 1,25 0,16 25,66 18,61 7,03 2,64
2,5 0,15 1,28 0,16 25,67 18,53 7,14 2,60
2,0 0,15 1,33 0,16 25,78 18,42 7,29 2,50
1,5 0,15 1,41 0,17 25,85 18,22 7,63 2,39

Analisando os dados da tabela acima podemos notar que a variação da vazão mássica de ar
que atravessa o condensador não influi significativamente na vazão mássica de refrigerante
deslocada pelo compressor e na pressão de evaporação. A pressão de condensação e o calor
trocado no condensador sofrem uma pequena queda com o aumento da vazão mássica de ar que
atravessa o condensador. A queda na pressão de condensação causa uma queda na temperatura de
condensação, diminuindo a diferença entre a temperatura do refrigerante no interior do trocador e
a temperatura do ar ambiente, conseqüentemente provocando uma queda no calor trocado. O
COP do sistema, por sua vez, sofre um aumento, devido à diminuição da diferença entre as
pressões de condensação e de evaporação. Como a pressão de evaporação permanece
praticamente inalterada e a pressão de condensação diminui, o trabalho de compressão também
diminui, reduzindo a potência consumida pelo compressor.

Os efeitos da variação da vazão mássica de ar que atravessa o evaporador sobre os

resultados obtidos no caso-teste padrão podem ser analisados através da tabela abaixo:

64
Tabela 5.17 - Análise da influência da vazão mássica de ar que atravessa o evaporador
mare mf Pc Pe qc qe Wcp
COP
(Kg/s) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
2,0 0,17 1,46 0,18 29,12 20,76 8,33 2,49
1,4 0,16 1,40 0,18 27,47 19,70 7,85 2,52
1,2 0,15 1,36 0,17 26,74 19,15 7,57 2,52
1,0 0,15 1,33 0,16 25,78 18,42 7,29 2,50
0,8 0,14 1,29 0,15 24,27 17,44 6,89 2,55
0,6 0,13 1,22 0,14 22,54 16,14 6,35 2,52

Na tabela acima, podemos observar que a variação da vazão mássica de ar que atravessa o
evaporador não influi significativamente na vazão mássica de refrigerante deslocada pelo
compressor.

O aumento da vazão mássica de ar que atravessa o evaporador acarreta em um aumento no

calor trocado no evaporador e por conseqüência, um aumento no trabalho de compressão, na
pressão de condensação e no calor rejeitado no condensador. Tal efeito é semelhante ao aumento
da área do evaporador. O COP do sistema permanece praticamente inalterado, pois a diferença
entre as pressões de condensação e evaporação não varia de maneira significativa.

Os efeitos causados pela variação do grau de superaquecimento do refrigerante na saída do

evaporador podem ser observados na tabela a seguir:

Tabela 5.18 - Análise da influência da variação do grau de superaquecimento do refrigerante na

saída do evaporador
∆te mf Pc Pe qc qe Wcp
COP
(ºC) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
10,0 0,14 1,24 0,16 25,34 18,51 6,89 2,71
8,0 0,14 1,28 0,16 25,59 18,47 7,06 2,60
6,0 0,15 1,31 0,16 25,78 18,43 7,22 2,51
5,0 0,15 1,33 0,16 25,78 18,42 7,29 2,50
4,0 0,15 1,35 0,16 25,88 18,39 7,38 2,46
2,0 0,15 1,39 0,17 25,91 18,64 7,58 2,42

65
 Podemos notar, pela análise da tabela acima, que o aumento do grau de superaquecimento
na saída do evaporador provoca uma queda na pressão de condensação e na potência do
compressor. Como a pressão de evaporação permanece quase constante, e a temperatura do
refrigerante na saída do evaporador torna-se maior com o aumento do grau de superaquecimento,
o volume específico do refrigerante na sucção do compressor também é maior. Como a vazão
mássica de refrigerante deslocada pelo compressor permanece quase inalterada, e o volume
específico aumenta, podemos concluir que o volume de refrigerante deslocado pelo compressor
diminui, causando a diminuição da pressão de condensação, diminuindo o trabalho de
compressão. O COP do ciclo sofre uma melhora, pois a diferença entre as pressões de
condensação e evaporação diminui com o aumento do grau de superaquecimento do refrigerante
na saída do evaporador.

Os efeitos causados pela variação da temperatura do ar ambiente podem ser observados na

tabela a seguir.

Tabela 5.19 - Análise da influência da variação da temperatura do ar ambiente

K mf Pc Pe qc qe Wcp
COP
(x 10-6) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
27,5 0,16 1,45 0,17 26,95 18,96 7,90 2,37
25,0 0,15 1,33 0,16 25,78 18,42 7,29 2,50
22,5 0,14 1,23 0,15 24,62 17,85 6,78 2,64
20,0 0,13 1,13 0,14 23,74 17,30 6,31 2,69
17,5 0,12 1,02 0,13 22,51 16,73 5,80 2,90
15,0 0,11 0,94 0,12 21,63 16,15 5,38 2,94

A partir da análise da tabela acima, podemos notar que a variação da temperatura do ar

ambiente tem grande influência sobre o comportamento do sistema de refrigeração. Quando a
temperatura ambiente é diminuída, a temperatura e a pressão de evaporação tendem a cair. Tal
queda ocorre devido à maior facilidade de rejeição de calor do condensador para o ambiente,
causando uma diminuição na temperatura e na pressão de condensação. Como podemos observar,
a vazão mássica de refrigerante deslocada pelo compressor, também tende a diminuir.A queda na
pressão de condensação dá se de maneira muito mais acentuada que a da pressão de evaporação,
ocasionando uma diferença menor entre as pressões de condensação e de evaporação,

66
melhorando o COP do sistema. A diferença entre as pressões, também ocasiona uma diminuição
na potência de compressão, o que acarreta em um menor consumo de potência pelo compressor.

A tabela abaixo demonstra a influência da variação da constante da válvula de expansão (K)

sobre o sistema.

Tabela 5.20 - Efeitos da variação da constante característica da válvula de expansão

K mf Pc Pe qc qe Wcp
COP
(x 10-6) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
4,2 0,15 1,32 0,16 25,80 18,41 7,34 2,49
4,0 0,14 1,37 0,16 25,65 18,64 7,29 2,53
3,8 0,14 1,42 0,16 25,57 18,35 7,25 2,54
3,6 0,13 1,65 0,15 25,14 18,00 7,14 2,52
3,4 0,13 1,57 0,15 25,29 18,15 7,17 2,54
3,2 0,12 1,65 0,15 25,14 18,00 7,14 2,52

A partir da análise da tabela acima podemos notar que a pressão de condensação aumenta
com a diminuição da constante da válvula. Isso se deve ao fato de que a diminuição de tal
constante representa uma restrição física (diminuição do diâmetro de um tubo capilar, por
exemplo). Tal restrição acarreta em uma diminuição na vazão de refrigerante que atravessa a
válvula, aumentando a pressão de condensação e diminuindo a vazão mássica de refrigerante
deslocada pelo compressor. Nota-se, também, que a variação da constante da válvula tem pouca
influência sobre os calores trocados pelo condensador e pelo evaporador.

A pressão de evaporação permanece praticamente inalterada, mas nota-se uma diminuição

da vazão mássica de refrigerante. Tal diminuição na vazão pode ser explicada pelo aumento da
restrição na vazão pela válvula de expansão, o que acarreta em um menor volume de refrigerante
deslocado para o evaporador. Como a vazão mássica a ser deslocada é menor, a potência de
compressão diminui. Tal diminuição de potência, somada ao calor trocado no evaporador (o qual
não sofre grande influência da variação da constante da válvula) provoca uma melhora no COP
do sistema.

67
5.2.5 Otimização da câmara frigorífica industrial

Em um primeiro momento, a otimização da câmara frigorífica foi feita a partir da

modificação dos parâmetros do sistema. Os valores escolhidos para cada parâmetro foram os que
apresentaram o melhor COP para cada uma das análises feitas acima, e podem ser visualizados na
tabela abaixo. A temperatura ambiente foi fixada em 20 ºC para a otimização.

Tabela 5.21 - Parâmetros para a otimização da câmara frigorífica

Ac Ae marc mare ∆te Vd
K (x 10-6)
(m2) (m2) (kg/s) (kg/s) (ºC) (m3)
50,0 30,0 4,0 0,8 10 0,5 x 10-3 3,80

Uma vez que os parâmetros acima foram inseridos no programa, obteve-se o primeiro
conjunto de resultados numéricos otimizados, conforme o demonstrado na tabela a seguir.

Tabela 5.22 - Resultados numéricos otimizados

mf Pc Pe qc qe Wcp
COP
(Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,12 0,91 0,12 23,03 17,46 5,49 3,18

Após a primeira otimização, a influência de cada um dos parâmetros sobre o COP do

sistema foi novamente estudada, permitindo fazer uma nova parametrização do sistema, de
maneira a obter-se um COP um pouco mais elevado.

A tabela a seguir mostra o novo conjunto de parâmetros obtidos para a otimização do

sistema.

Tabela 5.23 - Parâmetros para a otimização do sistema

Ac Ae marc mare ∆te Vd
K (x 10-6)
(m2) (m2) (kg/s) (kg/s) (ºC) (m3)
40 38 4,2 1,3 10 0,5 x 10-3 3,90

68
 Utilizando o conjunto de parâmetros descritos na tabela acima, obtivemos um novo
conjunto de resultados numéricos otimizados do sistema de refrigeração, conforme o descrito na
tabela a seguir.

Tabela 5.24 - Resultados numéricos otimizados

mf Pc Pe qc qe Wcp
COP
(Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,11 0,89 0,12 22,51 18,62 5,37 3,46

A tabela a seguir fornece os resultados numéricos obtidos para a comparação do

comportamento sistema antes da otimização e após a primeira e a segunda otimizações.

Tabela 5.25 - Análise comparativa dos resultados numéricos da otimização

mf Pc Pe qc qe Wcp COP
(Kg/s) (Mpa) (Mpa) (KW) (KW) (KW)
0,15 1,33 0,16 25,78 18,42 7,29 2,50
0,12 0,91 0,12 23,03 17,46 5,49 3,18
0,11 0,89 0,12 22,51 18,62 5,37 3,46

A análise da tabela acima permite concluir que a primeira fase da otimização foi a que
refletiu efeitos mais significativos na melhora do desempenho do sistema, ficando a segunda fase
como apenas um ajuste mais refinado dos parâmetros afim da obtenção do melhor desempenho
possível.

5.3 Estudo de caso: projeto e otimização de um sistema de ar condicionado

O sistema de refrigeração idealizado para a simulação do segundo caso-teste é um sistema

de ar condicionado de 0,4 kW, composto por um compressor alternativo, condensador de fluxo
cruzado com ar em convecção forçada, evaporador de fluxo cruzado com ar em convecção
forçada e válvula de expansão do tipo tubo capilar, operando com refrigerante R-134a. O sistema
acima descrito, foi escolhido por se tratar de um tipo de equipamento muito comumente utilizado
na prática para a climatização de ambientes (salas, quartos, consultórios, etc).

69
5.3.1 Parâmetros geométricos e operacionais do sistema de ar condicionado

Parâmetros iniciais para o compressor:

Wi = 0,025 KW
N = 3000 rpm.
Vd = 7,95 x 10-6 m3.
cr = 0,05.
α = 0,05.

Parâmetros iniciais para o condensador:

Ac= 0,025 m2
Uc = 1,75 KW/m2°C
marc = 0,2 Kg/s
tar = 25 ºC

Parâmetros iniciais para o evaporador:

Ae = 0,025 m2
Uc = 1,10 KW/m2°C
mare = 0,035 Kg/s
∆te = 4 ºC
tar = 25 ºC

Parâmetro inicial para a válvula de expansão:

K=1,5 x 10-7

70
5.3.2 Estimativas iniciais para o sistema de ar condicionado

Temperatura inicial na saída do evaporador:

t1i = 2 °C

Temperatura inicial na saída do condensador:

t3i= 40 ºC

5.3.3 - Resultados obtidos na simulação do caso-teste de um sistema de ar condicionado

Tabela 5.26 - Resultados do caso-teste para o sistema menor

mf Pc Pe qc qe Wcp COP
(Kg/s) (Mpa) (Mpa) (KW) (KW) (KW)
4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41

A partir dos resultados obtidos para o caso-teste do sistema de ar condicionado, pôde-se

fazer a análise comparativa dos efeitos da variação dos principais parâmetros iniciais nos
resultados.

É importante ressaltar, que as análises foram feitas para a variação de apenas um parâmetro
em questão, sendo que nos demais parâmetros foram mantidos os valores utilizados para a
simulação do caso teste do sistema de ar condicionado, já descritos acima.

Os efeitos da variação da área do condensador sobre os resultados podem ser visualizados

através da tabela abaixo:

71
Tabela 5.27 - Análise do efeito da variação da área do condensador
Ac mf Pc Pe qc qe Wcp
COP
(m2) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,040 4,16 x 10-3 0,90 0,24 0,42 0,43 0,15 2,86
-3
0,035 4,23 x 10 0,94 0,24 0,45 0,43 0,16 2,68
0,030 4,48 x 10-3 1,04 0,26 0,48 0,41 0,17 2,41
-3
0,025 4,48 x 10 1,04 0,26 0,43 0,41 0,17 2,41
-3
0,020 4,72 x 10 1,14 0,27 0,44 0,40 0,18 2,22
0,015 5,15 x 10-3 1,39 0,30 0,50 0,37 0,21 1,76

Analisando a tabela acima, podemos notar que o aumento da área do condensador causa
uma diminuição na pressão de condensação. A diminuição da pressão de condensação pode ser
explicada pelo aumento no volume do trocador, e conseqüente diminuição na resistência à
descarga do fluido na saída do compressor. A pressão de evaporação sofre uma ligeira
diminuição devido à diminuição da pressão na entrada da válvula de expansão. A diminuição da
pressão de condensação causa um aumento no volume específico do refrigerante na sucção do
compressor, causando a diminuição da vazão mássica de refrigerante deslocada. Como a pressão
de condensação diminui de maneira mais acentuada que a de evaporação, a diferença entre
pressões tende a cair, diminuindo o trabalho de compressão (e conseqüentemente a potência
consumida pelo compressor) e melhorando o COP do ciclo.

Os efeitos da variação da área do evaporador sobre os resultados obtidos no caso-teste

padrão podem ser analisados através da tabela abaixo.

Tabela 5.28 - Análise da influência da variação da área do evaporador

Ae mf Pc Pe qc qe Wcp
COP
(m2) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,040 4,90 x 10-3 1,23 0,28 0,63 0,50 0,19 2,63
-3
0,035 4,77 x 10 1,17 0,28 0,59 0,47 0,18 2,61
-3
0,030 4,61 x 10 1,09 0,26 0,56 0,45 0,18 2,50
0,025 4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41
-3
0,020 4,45 x 10 1,04 0,26 0,46 0,36 0,17 2,11
-3
0,015 4,21 x 10 0,93 0,24 0,34 0,32 0,16 2,00

72
 Na tabela acima, podemos observar que o aumento da área do evaporador produz uma
melhora na quantidade de calor trocado pelo mesmo. O aumento da quantidade de calor trocada
no evaporador, produz um aumento na temperatura do refrigerante em seu interior, acarretando
na diminuição do volume específico do refrigerante na sucção do compressor, causando um
aumento na vazão mássica deslocada para o condensador. Como o a área do condensador
permanece inalterada, a pressão de condensação tende a aumentar causando maior resistência à
descarga do compressor e aumentando o trabalho de compressão e a potência consumida pelo
compressor. O aumento da pressão de condensação contribui para o aumento da pressão de
evaporação, pois, com o aumento da pressão na entrada da válvula de expansão, esta tende a
permitir a passagem de uma maior vazão de refrigerante. O COP do ciclo aumenta, pois o
aumento da potência consumida pelo compressor é menor do que o aumento do calor trocado no
evaporador.

Os efeitos da variação da vazão mássica de ar que atravessa o condensador sobre os

resultados obtidos no caso-teste padrão podem ser analisados através da tabela abaixo:

Tabela 5.29 - Análise da influência da variação da vazão mássica de ar que atravessa o

condensador
marc mf Pc Pe qc qe Wcp
COP
(Kg/s) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,6 4,47 x 10-3 1,05 0,25 0,54 0,41 0,17 2,41
0,5 4,48 x 10-3 1,05 0,26 0,52 0,41 0,17 2,41
-3
0,4 4,47 x 10 1,05 0,26 0,51 0,41 0,17 2,41
-3
0,3 4,48 x 10 1,06 0,26 0,52 0,41 0,17 2,41
0,2 4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41
-3
0,1 4,47 x 10 1,04 0,26 0,42 0,41 0,17 2,41

Na tabela acima, podemos observar que o aumento da vazão mássica de ar que atravessa o
condensador tem pouca influência sobre o calor trocado no evaporador, potência consumida pelo
compressor, COP e pressão de evaporação.

A pressão de condensação sofre um sensível aumento devido ao aumento no calor trocado

no condensador, e conseqüente diminuição da temperatura do refrigerante no interior do trocador,

73
o que acarreta em um aumento no volume específico. A vazão mássica de refrigerante deslocada
pelo compressor permanece praticamente inalterada.

Os efeitos da variação da vazão mássica de ar que atravessa o evaporador sobre os

resultados obtidos no caso-teste padrão podem ser analisados através da tabela abaixo.

Tabela 5.30 - Análise da influência da vazão mássica de ar que atravessa o evaporador

mare mf Pc Pe qc qe Wcp
COP
(Kg/s) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,065 4,68 x 10-3 1,13 0,27 0,50 0,46 0,18 2,55
-3
0,055 4,59 x 10 1,10 0,26 0,47 0,45 0,18 2,50
-3
0,045 4,57 x 10 1,08 0,26 0,47 0,43 0,18 2,39
0,035 4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41
-3
0,025 4,44 x 10 1,03 0,26 0,43 0,38 0,17 2,23
-3
0,015 4,31 x 10 0,97 0,25 0,41 0,32 0,16 2,00

Na tabela acima, podemos notar que o calor trocado no evaporador aumenta com o aumento
da vazão mássica de ar que atravessa o evaporador. Isso acarreta um aumento na temperatura do
refrigerante na sucção do compressor, causando uma diminuição do volume específico do mesmo
e aumentando a vazão mássica deslocada pelo compressor. Tal aumento ocasiona um aumento na
pressão de condensação e conseqüentemente, um pequeno aumento na potência consumida pelo
compressor, devido à maior resistência ao deslocamento do refrigerante causada pelo aumento
da pressão de condensação. O COP do ciclo aumenta, pois o aumento na quantidade de calor
trocado no evaporador é mais significativo que a potência consumida pelo compressor.

Os efeitos causados pela variação do grau de superaquecimento do refrigerante na saída do

evaporador podem ser observados na tabela a seguir.

74
Tabela 5.31 - Análise da influência da variação do grau de superaquecimento do refrigerante na
saída do evaporador
∆te mf Pc Pe qc qe Wcp
COP
(ºC) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
6,0 4,45 x 10-3 1,02 0,25 0,42 0,41 0,17 2,41
-3
5,0 4,45 x 10 1,04 0,25 0,43 0,41 0,17 2,41
4,0 4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41
-3
3,0 4,64 x 10 1,10 0,27 0,48 0,40 0,18 2,22
-3
2,0 4,75 x 10 1,16 0,27 0,58 0,40 0,19 2,10
1,0 4,95 x 10-3 1,26 0,29 0,57 0,39 0,20 1,95

Podemos notar, observando os dados fornecidos pela tabela acima, que a diminuição do
grau de superaquecimento do refrigerante na saída do evaporador acarreta em um aumento na
pressão de condensação e uma conseqüente diminuição no volume específico do refrigerante no
condensador. A diminuição do volume específico acarreta em um aumento da vazão mássica de
refrigerante que atravessa a válvula de expansão, provocando um aumento na pressão de
evaporação. Tal aumento de pressão ocasiona um aumento no volume específico do refrigerante
no evaporador, aumentando a vazão mássica deslocada pelo compressor. No entanto, como a
pressão de condensação sofre um aumento significativamente maior que a pressão de evaporação,
o trabalho de compressão e a potência consumida pelo compressor tendem a aumentar, enquanto
que o calor trocado no evaporador tende a diminuir levemente, causando uma diminuição no
COP do sistema.

Os efeitos causados pela variação da temperatura do ar ambiente podem ser observados na

tabela a seguir:

75
Tabela 5.32 - Análise da influência da variação da temperatura do ar ambiente
tamb mf Pc Pe qc qe Wcp
COP
(ºC) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
30,0 4,94 x 10-3 1,27 0,29 0,49 0,45 0,20 2,25
-3
27,5 4,79 x 10 1,18 0,28 0,48 0,43 0,19 2,26
25,0 4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41

22,5 4,27 x 10-3 0,93 0,24 0,39 0,40 0,16 2,50

20,0 4,04 x 10-3 0,86 0,22 0,39 0,38 0,15 2,53
-3
17,5 4,04 x 10 0,83 0,22 0,40 0,35 0,14 2,50

Podemos notar, na tabela acima, que com o aumento da temperatura do ar ambiente, a

pressão de evaporação tende a aumentar, acarretando na diminuição do volume específico do
refrigerante na entrada do compressor e no conseqüente aumento da vazão mássica deslocada
para o condensador. O aumento da vazão mássica deslocada causa o aumento da pressão de
condensação, o que gera um aumento no trabalho de compressão e na potência consumida pelo
compressor.

A pressão de evaporação tende a aumentar, pois com o aumento na pressão de condensação,

uma maior vazão de refrigerante tende a deslocar-se pela válvula de expansão. O COP do ciclo
tende a diminuir com o aumento da temperatura ambiente, devido ao fato de que o crescimento
da potência consumida pelo compressor é maior do que o aumento do calor trocado no
evaporador.

A tabela abaixo demonstra a influência da variação da constante da válvula de expansão (K)

sobre o sistema.

76
Tabela 5.33 - Efeitos da variação da constante característica da válvula de expansão
K mf Pc Pe qc qe Wcp
COP
(x 10-6) (Kg/s) (MPa) (MPa) (KW) (KW) (KW)
0,17 5,02 x 10-3 1,03 0,28 0,50 0,40 0,17 2,35
-3
0,16 4,69 x 10 1,00 0,26 0,40 0,40 0,17 2,35
0,15 4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41
-3
0,14 4,31 x 10 1,08 0,25 0,46 0,42 0,17 2,47
-3
0,13 4,25 x 10 1,19 0,25 0,51 0,42 0,18 2,33
0,12 3,92 x 10-3 1,20 0,23 0,53 0,43 0,18 2,39

A análise da tabela acima permite observar que a pressão de condensação tende a aumentar
com a diminuição da constante característica da válvula de expansão. Tal aumento de pressão
pode ser explicado pelo fato de que a diminuição da constante representa uma redução no
diâmetro físico do tubo da válvula de expansão, causando uma maior resistência ao escoamento
do refrigerante em seu interior. Isso também explica a sensível redução na pressão de evaporação,
com um conseqüente aumento do volume específico do refrigerante na sucção do compressor,
que passa a deslocar uma menor vazão mássica e a consumir uma maior potência devido ao
aumento da diferença entre as pressões de evaporação e condensação, reduzindo o COP do
sistema.

5.3.4 Otimização do sistema de ar condicionado

Em um primeiro momento, a otimização do sistema de ar condicionado foi feita a partir da

modificação dos parâmetros do sistema e da fixação da temperatura ambiente a 20 ºC. Os valores
escolhidos para cada parâmetro foram os que apresentaram o melhor COP para cada uma das
análises feitas acima, e podem ser visualizados na tabela a seguir.

Tabela 5.34 - Parâmetros para a otimização do sistema de ar condicionado

Ac Ae marc mare ∆te Vd
K (x 10-6)
(m2) (m2) (kg/s) (kg/s) (ºC) (m3)
0,04 0,04 0,2 0,06 4,0 7,95 x 10-6 0,14

77
 Utilizando o conjunto de parâmetros descritos na tabela acima, obtivemos um novo
conjunto de resultados numéricos otimizados do sistema de ar condicionado, conforme o descrito
na tabela a seguir.

Tabela 5.35 - Resultados numéricos otimizados

mf Pc Pe qc qe Wcp
COP
(Kg/s) (MPa) (MPa) (KW) (KW) (KW)
4,12 x 10-3 0,98 0,24 0,59 0,54 0,16 3,37

Após a primeira otimização, a influência de cada um dos parâmetros sobre o COP do

sistema foi novamente estudada, permitindo fazer uma nova parametrização do sistema de ar
condicionado, de maneira a obter-se um COP um pouco mais elevado.

A tabela a seguir mostra o novo conjunto de parâmetros obtidos para a otimização do

sistema.

Tabela 5.36 - Parâmetros para a otimização do sistema

Ac Ae marc mare ∆te Vd
K (x 10-6)
(m2) (m2) (kg/s) (kg/s) (ºC) (m3)
0,03 0,04 0,20 0,03 4,0 0,5 x 10-6 0,014

Utilizando o conjunto de parâmetros descritos na tabela acima, obtivemos um novo

conjunto de resultados numéricos otimizados do sistema de ar condicionado para a potência no
evaporador requerida no projeto (0,4 kW), conforme o descrito na tabela a seguir.

Tabela 5.37 - Resultados numéricos otimizados

mf Pc Pe qc qe Wcp
COP
(kg/s) (MPa) (MPa) (kW) (kW) (kW)
3,79 x 10-3 0,87 0,33 0,35 0,40 0,11 3,64

A tabela a seguir fornece os resultados numéricos obtidos para a comparação do

comportamento sistema antes da otimização e após a primeira e a segunda otimizações.

78
Tabela 5.38 - Análise comparativa dos resultados numéricos da otimização
mf Pc Pe qc qe Wcp COP
(kg/s) (MPa) (MPa) (kW) (kW) (kW)
4,48 x 10-3 1,04 0,26 0,43 0,41 0,17 2,41
-3
4,12 x 10 0,98 0,24 0,59 0,54 0,16 3,37
3,79 x 10-3 0,87 0,33 0,35 0,40 0,11 3,64

A análise da tabela acima permite concluir que para este estudo de caso, a primeira fase foi
fundamental para uma pré otimização do sistema, ficando a segunda fase apenas como um ajuste
fino dos parâmetros afim da otimização do COP do sistema para a capacidade de refrigeração
requerida no projeto.

79
Anexo
"Programa de Simulação de Sistemas de Ar Condicionado"

"Cálculo das estimativas iniciais"

PROCEDURE
INICIAL(t_1_i;t_3_i;DELTAt_e:t_1_i_sat;P_e_i;v_1_i;h_1_i;s_1_i;c_p_1_i;i;X_1_i;P_c_i;v_3_i;h_3_i;s_3_i;
c_p_3_i;c_v_3_i;X_3_i;t_2_i;v_2_i;h_2_i;s_2_i;c_p_2_i;c_v_2_i;X_2_i;t_4_i;h_4_i;X_4_i)

"Temperatura de saturação no evaporador:" t_1_i_sat=t_1_i-

DELTAt_e
"Pressão de evaporação:"
P_e_i=P_SAT(R134a;T=t_1_i_sat)
"Volume específico no ponto 1 (saída do evaporador):"
v_1_i=VOLUME(R134a;T=t_1_i;P=P_e_i)
"Entalpia específica no ponto 1 (saída do evaporador):"
h_1_i=ENTHALPY(R134a;T=t_1_i;P=P_e_i)
"Entropia específica no ponto 1 (saída do evaporador):"
s_1_i=ENTROPY(R134a;T=t_1_i;P=P_e_i)
"Calor específico a pressão constante no ponto 1 (saída do evaporador):"
c_p_1_i=CP(R134a;T=t_1_i;P=P_e_i)
"Calor específico a volume constante no ponto 1 (saída do evaporador):"
c_v_1_i=CV(R134a;T=t_1_i;P=P_e_i)
"Título no ponto 1 (saída do evaporador):"
X_1_i=QUALITY(R134a;T=t_1_i;P=P_e_i)

"Pressão de condensação:"
P_c_i=P_SAT(R134a;T=t_3_i)
"Volume específico no ponto 3 (saída do condensador):"
v_3_i=VOLUME(R134a;x=0;P=P_c_i)
"Entalpia específica no ponto 3 (saída do condensador):"
h_3_i=ENTHALPY(R134a;x=0;P=P_c_i)
"Entropia específica no ponto 3 (saída do condensador):"
s_3_i=ENTROPY(R134a;x=0;P=P_c_i)
"Calor específico a pressão constante no ponto 3 (saída do condensador):"
c_p_3_i=CP(R134a;x=0;P=P_c_i)
"Calor específico a volume constante no ponto 3 (saída do condensador):"
c_v_3_i=CV(R134a;x=0;P=P_c_i)
"Título no ponto 3 (saída do condensador):" X_3_i=0

"Entropia específica no ponto 2 (saída do compressor):" s_2_i=s_1_i

"Temperatura no ponto 2 (saída do compressor):"
t_2_i=TEMPERATURE(R134a;S=s_2_i;P=P_c_i)

87
"Volume específico no ponto 2 (saída do compressor):"
v_2_i=VOLUME(R134a;T=t_2_i;P=P_c_i)
"Entalpia específica no ponto 2 (saída do compressor):"
h_2_i=ENTHALPY(R134a;T=t_2_i;P=P_c_i)
"Calor específico a pressão constante no ponto 2 (saída do compressor):"
c_p_2_i=CP(R134a;T=t_2_i;P=P_c_i)
"Calor específico a volume constante no ponto 2 (saída do compressor):"
c_v_2_i=CV(R134a;T=t_2_i;P=P_c_i)
"Título no ponto 2 (saída do compressor):"
X_2_i=QUALITY(R134a;T=t_2_i;P=P_c_i)

"Entalpia específica no ponto 4 (saída da válvula):" h_4_i=h_3_i

"Temperatura no ponto 4 (saída da válvula):"
t_4_i=TEMPERATURE(R134a;H=h_4_i;P=P_e_i)
"Título no ponto 4 (saída da válvula):"
X_4_i=QUALITY(R134a;H=h_4_i;P=P_e_i)

END

"Cálculos iterativos"

PROCEDURE
ITERAÇÕES(t_1_i_sat;P_e_i;v_1_i;h_1_i;s_1_i;c_p_1_i;c_v_1_i;;P_c_i;v_3_i;h_3_i;s_3_i;c_p_3_i;c_v_3_
i;X_3_i;t_2_i;v_2_i;h_2_i;s_2_i;c_p_2_i;c_v_2_i;X_2_i;t_4_i;h_4_i;X_4_i;W_i;N;V_d;c_r;alpha;c_val;A_c;U
_c;m_ar_c;t_3_i;A_e;U_e;m_ar_e;t_1_i;DELTAt_e;t_ar;c_p_ar;t_agua;c_p_agua;E:P_c_atual;P_e_atual)

"Inicialização das variáveis"

t_1_anterior=t_1_i
t_1_sat_anterior=t_1_i_sat
P_e_anterior=P_e_i
v_1_anterior=v_1_i
h_1_anterior=h_1_i
s_1_anterior=s_1_i
c_p_1_anterior=c_p_1_i
c_v_1_anterior=c_v_1_i
X_1_anterior=X_1_i
t_3_anterior=t_3_i
P_c_anterior=P_c_i
v_3_anterior=v_3_i
h_3_anterior=h_3_i
s_3_anterior=s_3_i
c_p_3_anterior=c_p_3_i
c_v_3_anterior=c_v_3_i
X_3_anterior=X_3_i
t_2_anterior=t_2_i
v_2_anterior=v_2_i
h_2_anterior=h_2_i
s_2_anterior=s_2_i
c_p_2_anterior=c_p_2_i
c_v_2_anterior=c_v_2_i
X_2_anterior=X_2_i
t_4_anterior=t_4_i
h_4_anterior=h_4_i
X_4_anterior=X_4_i

88
"Contador de iterações" contador=0

Repeat

"Compressor"
gama=c_p_1_anterior/c_v_1_anterior
eta_v=1+c_r-(c_r*((P_c_anterior/P_e_anterior)^(1/gama)))
m_dot=((eta_v*V_d)/v_1_anterior)*N
W_comp=m_dot*P_e_anterior*v_1_anterior*((P_c_anterior/P_e_anterior)^((gama-1)/gama)-
1)*(gama/(gama-1))*1000
W_cp=W_i+(1+alpha)*W_comp
h_2_atual=h_1_anterior+W_cp/m_dot
t_2_atual=TEMPERATURE(R134a;h=h_2_atual;P=P_c_anterior)

"Condensador"
C_f_c=m_ar_c*c_p_ar
c_p_r_c=SPECHEAT(R134a;x=1;P=P_c_anterior)
hlv_c=ENTHALPY(R134a;x=1;P=P_c_anterior)-ENTHALPY(R134a;x=0;P=P_c_anterior)
C_q_c=m_dot*c_p_r_c+m_dot*hlv_c/(t_2_atual-t_ar)
C_min_c=min(C_f_c;C_q_c)
C_max_c=max(C_f_c;C_q_c)
C_c=C_min_c/C_max_c
NUT_c=(U_c*A_c)/C_min_c
epsilon_c=1-exp((1/C_c)*(NUT_c^(0,22))*(exp((-C_c)*(NUT_c^(0,78)))-1))
q_max_c=C_min_c*(t_2_atual-t_ar)
q_c=epsilon_c*q_max_c
h_3_atual=h_2_atual-(q_c/m_dot)
t_3_atual=TEMPERATURE(R134a;h=h_3_atual;P=P_c_anterior)
t_ar_s_c=t_ar+(q_c/(m_ar_c*c_p_ar))
x_3_atual=QUALITY(R134a;h=h_3_atual;P=P_c_anterior)
t_3_sat_atual=TEMPERATURE(R134a;x=1;P=P_c_anterior)
v_3_atual=VOLUME(R134a;x=0;P=P_c_anterior)

"Válvula de expansão"
h_4_atual=h_3_atual
t_4_atual=TEMPERATURE(R134a;h=h_4_atual;P=P_e_anterior)
x_4_atual=QUALITY(R134a;h=h_4_atual;P=P_e_anterior)

"Evaporador"
C_q_e=m_ar_e*c_p_ar
c_p_r_e=CP(R134a;x=1;P=P_e_anterior)
hlv_e=ENTHALPY(R134a;x=1;P=P_e_anterior)-ENTHALPY(R134a;x=0;P=P_e_anterior)
C_f_e=m_dot*c_p_r_e+m_dot*hlv_e/(t_ar-t_4_atual)
C_min_e=min(C_f_e;C_q_e)
C_max_e=max(C_f_e;C_q_e)
C_e=C_min_e/C_max_e
NUT_e=(U_e*A_e)/C_min_e
epsilon_e=1-exp((1/C_e)*(NUT_e^(0,22))*(exp((-C_e)*(NUT_e^(0,78)))-1))
"epsilon_e=(1-exp(-NUT_e*(1-C_e)))/(1-C_e*exp(-NUT_e*(1-C_e)))"
q_max_e=C_min_e*(t_ar-t_4_atual)
q_e=epsilon_e*q_max_e
t_ar_s_e=t_ar-((q_e)/(m_ar_e*c_p_ar))
h_1_atual=(q_e/m_dot)+h_4_atual
t_1_atual=TEMPERATURE(R134a;h=h_1_atual;P=P_e_anterior)
t_1_sat_atual=t_1_atual-DELTAt_e
x_1_atual=QUALITY(R134a;h=h_1_atual;P=P_e_anterior)

89
"Cálculos das pressões - modelo da válvula de expansão"
P_e_atual=0,2*(P_SAT(R134a;T=t_1_sat_atual))+(1-0,2)*P_e_anterior
P_c=P_e_anterior+0,000001*v_3_atual*(m_dot/c_val)^2
v_1_atual=VOLUME(R134a;T=t_1_atual;P=P_e_atual)
c_p_1_atual=CP(R134a;T=t_1_atual;P=P_e_atual)
c_v_1_atual=CV(R134a;T=t_1_atual;P=P_e_atual)

"Variáveis para a verificação de parada"

Verif_Pc=P_c_atual-P_c_anterior
Verif_Pe=P_e_atual-P_e_anterior

"Atribuição de novos valores"

t_1_anterior=t_1_atual
t_1_sat_anterior=t_1_sat_atual
P_e_anterior=P_e_atual
v_1_anterior=v_1_atual
h_1_anterior=h_1_atual
c_p_1_anterior=c_p_1_atual
c_v_1_anterior=c_v_1_atual
X_1_anterior=X_1_atual
t_3_anterior=t_3_atual
t_3_sat_anterior=t_3_sat_atual
P_c_anterior=P_c_atual
h_3_anterior=h_3_atual
X_3_anterior=X_3_atual
t_2_anterior=t_2_atual
h_2_anterior=h_2_atual
t_4_anterior=t_4_atual
h_4_anterior=h_4_atual
X_4_anterior=X_4_atual

"Contador de iterações"
contador=contador+1

"Verficação de parada"
Until ((abs(Verif_Pc)<E) and (abs(Verif_Pe)<E))

END

"Parâmetros iniciais"

"Potência dissipada no compressor:" W_i=0,025

"Número de rotações por segundo do compressor:" N=3000/60
"Volume deslocado pelo compressor:" V_d=0,0000050
"Fator de folga do compressor:" c_r=0,05
"Fator de perdas do compressor:" alpha=0,05

"Coeficiente da válvula de expansão" c_val=0,00000015

"Área do condensador:" A_c=0,035

"Coeficiente global de transferência de calor do condensador:" U_c=1,75
"Vazão mássica de ar que atravessa o condensador:" m_ar_c=0,2
"Temperatura no ponto 3 (saída do condensador):" t_3_i=40

90
"Área do evaporador:" A_e=0,04
"Coeficiente global de transferência de calor do evaporador:" U_e=1,1
"Vazão mássica de ar que atravessa o evaporador:" m_ar_e=0,035
"Temperatura no ponto 1 (saída do evaporador)" t_1_i=2,0
"Grau de superaquecimento na saída do evaporador:" DELTAt_e=4,0

"Temperatura do ar ambiente:" t_ar=25,0

"Calor específico do ar ambiente:"
c_p_ar=CP(Air;T=t_ar)
"temperatura da água"
t_agua=22,6
"calor específico da água"
c_p_agua=CP(H2O;T=t_agua)

"Critério de parada:" E=0,002

"Chamada dos procedimentos de cálculo"

CALL
INICIAL(t_1_i;t_3_i;DELTAt_e:t_1_i_sat;P_e_i;v_1_i;h_1_i;s_1_i;c_p_1_i;c_v_1_i;X_1_i;P_c_i;v_3_i;h_3_
i;s_3_i;c_p_3_i;c_v_3_i;X_3_i;t_2_i;v_2_i;h_2_i;s_2_i;c_p_2_i;c_v_2_i;X_2_i;t_4_i;h_4_i;X_4_i)

CALL
ITERAÇÕES(t_1_i_sat;P_e_i;v_1_i;h_1_i;s_1_i;c_p_1_i;c_v_1_i;X_1_i;P_c_i;v_3_i;h_3_i;s_3_i;c_p_3_i;c
_v_3_i;X_3_i;t_2_i;v_2_i;h_2_i;s_2_i;c_p_2_i;c_v_2_i;X_2_i;t_4_i;h_4_i;X_4_i;W_i;N;V_d;c_r;alpha;c_val;
A_c;U_c;m_ar_c;t_3_i;A_e;U_e;m_ar_e;t_1_i;DELTAt_e;t_ar;c_p_ar;t_agua;c_p_agua;E:P_c_atual;P_e_
atual)

91