Evaporator

Overview

In the evaporator, refrigerant enters at some vapor quality between 0% and 100%, and typically exits as a superheated vapor. The analysis for the evaporator shares many features with the analysis employed for the Condenser and the Cooling Coil, but the evaporator’s analysis is significantly more complicated due to the fact that there are moving boundaries on both the refrigerant and air sides. On the refrigerant-side, the moving boundary is between two-phase refrigerant and superheated refrigerant, and on the air-side, there is a moving boundary between wet and dry parts of the coil.

Mathematical Description

The evaporator is treated as being cross-counter-flow. What this means practically is that the mass flow of air to and total air side surface area for each of the superheated and two-phase portions (assuming each exists), are given by

(1)A_{a,superheat}=w_{superheat}A_{a,total}

A_{a,two-phase}=w_{two-phase}A_{a,total}

\dot m_{a,superheat}=w_{superheat}\dot m_{a,total}

\dot m_{a,two-phase}=w_{two-phase}\dot m_{a,total}

which is analogous to the analysis employed for the condenser. On the refrigerant-side, the surface area for superheated and two-phase sections are given by

(2)A_{r,superheat}=w_{superheat}A_{r,total}

A_{r,two-phase}=w_{two-phase}A_{r,total}

Two-Phase Section

In the two-phase section, the target heat transfer rate is given by

(3)\dot Q_{target}=\dot m_r(x_{r,o}-x_{r,i})h_{fg}

and the heat transfer rate from the partially-wet/partially-dry analysis in section Two-phase on working fluid side is given the name \dot Q_{PWPD}. To begin with, the first guess is that all the heat exchanger is in the two-phase region. If using the outlet quality of saturated vapor (x_{r,o}=1) with all the heat exchanger in the two-phase region, and \dot Q_{PWPD} is greater than \dot Q_{target}, too much of the heat exchanger area was given to the two-phase section, and there must be a superheated section as well. Thus, the length fraction w_{two-phase} is iteratively altered using a bounded solver (w_{two-phase} is bounded between 0 and 1). The remaining area is given to the superheated section, so the superheated section circuit length fraction w_{superheat} is given by

(4)w_{superheat}=1-w_{two-phase}

Superheated Section

In the superheated section, for a given w_{superheat}, the analysis is exactly the same as the Cooling Coil, and given by the section Single-phase on working fluid side. As in the cooling coil, the working fluid is single-phase, and the air-side surface may be fully wet, partially wet, or fully dry, and all the inlet conditions for the superheated section are known. The air inlet state to the superheated portion is assumed to be the same as for the two-phase portion because the superheated portion is typically rather small, and both superheated and two-phase portions should see approximately the same inlet air state. The refrigerant properties are calculated using the same single-phase correlations as for the cooling coil.

Minimal Example Code

from CoolProp.CoolProp import Props
from FinCorrelations import FinInputs
from Evaporator import EvaporatorClass

FinsTubes=FinInputs()
    
FinsTubes.Tubes.NTubes_per_bank=32
FinsTubes.Tubes.Ncircuits=5
FinsTubes.Tubes.Nbank=3
FinsTubes.Tubes.Ltube=0.452
FinsTubes.Tubes.OD=0.009525
FinsTubes.Tubes.ID=0.0089154
FinsTubes.Tubes.Pl=0.0254
FinsTubes.Tubes.Pt=0.0219964

FinsTubes.Fins.FPI=14.5
FinsTubes.Fins.Pd=0.001
FinsTubes.Fins.xf=0.001
FinsTubes.Fins.t=0.00011
FinsTubes.Fins.k_fin=237

FinsTubes.Air.Vdot_ha=0.5663
FinsTubes.Air.Tdb=299.8
FinsTubes.Air.p=101.325
FinsTubes.Air.RH=0.51
FinsTubes.Air.FanPower=438

kwargs={'Ref': 'R410A',
    'mdot_r':  0.0708,
    'psat_r':  Props('P','T',282,'Q',1.0,'R410A'),
    'Fins': FinsTubes,
    'hin_r':Props('H','T',282,'Q',0.15,'R410A')*1000, #[J/kg]
    'Verbosity':0
}

Evap=EvaporatorClass(**kwargs)
Evap.Update(**kwargs)
Evap.Calculate()

print 'Evaporator heat trasfer rate is',Evap.Q,'W'
print 'Evaporator capacity (less fan power) is',Evap.Capacity,'W'
print 'Evaporator fraction of length in two-phase section',Evap.w_2phase,'W'
print 'Evaporator sensible heat ratio',Evap.SHR

If you open an IPython(x,y) shell in the root of the documentation (folder Documentation/Web relative to the main trunk), and run the commands below, you should get

In [5]: execfile('ACHPComponents/ComponentTests/EvaporatorTest.py')
Evaporator heat trasfer rate is 13222.5359037 W
Evaporator capacity (less fan power) is 12784.5359037 W
Evaporator fraction of length in two-phase section 0.89179879315 W
Evaporator sensible heat ratio 0.709186116025

If not, first stop should be the Frequently Asked Questions

Nomenclature

Variable Description
A_{a,total} Total air side area (fins+tubes) [m2]
A_{a,superheat} Air-side area for superheated portion of circuit [m2]
A_{a,two-phase} Air-side area for two-phase portion of circuit [m2]
A_{r,superheat} Refrigerant-side area for superheated portion of circuit [m2]
A_{r,two-phase} Refrigerant-side area for two-phase portion of circuit [m2]
h_{fg} Latent heat of refrigerant [J/kg]
\dot m_{a,total} Mass flow rate of dry air [kg/s]
\dot m_{a,superheat} Mass flow rate of dry air in superheated section [kg/s]
\dot m_{a,two-phase} Mass flow rate of dry air in two-phase section [kg/s]
\dot m_r Mass flow rate of refrigerant [kg/s]
\dot Q_{target} Target heat transfer rate [W]
\dot Q_{PWPD} Heat transfer rate from the partially wet/partially-dry analysis [W]
x_{r,i} Refrigerant inlet quality [-]
x_{r,o} Refrigerant outlet quality [-]

Component Class Documentation

class Evaporator.EvaporatorClass(**kwargs)[source]
AirSideCalcs()[source]
Calculate()[source]
Initialize()[source]
OutputList()[source]

Return a list of parameters for this component for further output

It is a list of tuples, and each tuple is formed of items:
[0] Description of value [1] Units of value [2] The value itself
Update(**kwargs)[source]

Table Of Contents

Previous topic

Cooling Coil

Next topic

Multi-Circuited Evaporator

This Page