The goal of a condenser is to take a refrigerant stream at a superheated state and cool it at nearly constant pressure so that it condenses to a subcooled refrigerant that can then be throttled to the low-pressure side of the system, exiting the throttling valve at low temperature.
On the refrigerant side, the refrigerant path can be divided into as many as three sections - a superheated section, a two-phase section, and a subcooled section. In practice, operating conditions (particularly low charge), may result in a condenser that does not have a subcooled section, and at extreme conditions, does not have a two-phase section either. The following plot shows the possible configurations for the first two cases described:
(Source code, png, hires.png, pdf)
In both cases, the schematic represents an averaged circuit on the refrigerant side. Overall, the goal of solver for the condenser is to determine how much of the length of an averaged circuit is in the superheated, two-phase, and subcooled phases. The fraction of the circuit length in each of these segments are given by , , and respectively. The sum of these factors must be equal to unity, and if any of the sections do not exist, its respective circuit fraction parameter is set to zero.
Thus if the average length of one circuit is given by , then the lengths of each segment are given by
(1)
and for consistency
(2)
On the air side, the condenser is treated as being pure cross-flow. This is a particularly good assumption for the condition where the condenser is composed of a single bank of tubes, which is a common configuration for condensers in the USA. Thus, since the flow is pure crossflow, the air-side area, as well as the air flow rate for each of the superheated, two-phase, and subcooled sections can be assumed to be proportional to the area for that section (See below).
The analysis here makes a number of assumptions:
On the air side, the air flow rate flow is assumed to be evenly distributed across the face of the coil which has a total air-side area of . Therefore the areas and mass flow rates for each section are given by
(3)
(Source code, png, hires.png, pdf)
Explicitly solve the superheated section to determine
First assume that the refrigerant exits the two-phase portion at a quality of 0, and calculate the required from the heat transfer analysis for the two-phase section. Then either
If the value of is less than 1, there is a subcooled section
- The outlet quality of the two-phase section is therefore known to be 0, which yields the value for
- Use the remainder of the condenser area for the subcooled section (), solve subcooled section for heat transfer rate, pressure drop, etc.
If the value of is greater than 1, there is no subcooled section
- Area fraction for two phase section is known to be , iteratively calculate the outlet quality of the condenser
The analysis for each of these calculations is described in the sections that follow.
In the superheated region, the inlet refrigerant temperature and the outlet refrigerant temperature (dew temperature) are known. Therefore, assuming pure crossflow, the fraction of the area required for the superheated region can be explicitly obtained from
(4)
where
(5)
and where the derivation of this term can be obtained from Condenser Area Derivation. The superheated portion is assumed to exist.
In the two-phase section of the heat exchanger, there are two basic possibilities. Either the outlet of the two-phase section is at some two-phase quality (because there is no subcooled section) or the outlet of the two-phase section is at a quality of 0 (saturated liquid) because there is a subcooled section. The first step is to assume that the outlet of the two-phase section is at a quality of 0 and calculate the required fraction of the circuit length. Before the length fraction can be calculated, the average refrigerant side heat transfer coefficient is required, for which the analysis can be found in section Shah Condensation. The average refrigerant-side heat transfer coefficient is a function of refrigerant outlet quality.
In either configuration of the two-phase section (solving for outlet quality or solving for ), the effectiveness is known, since the parameter cancels out of the solution for , which means that the Ntu are independent of the length fraction . Thus the effectiveness can be obtained directly from
(6)
where the value of is given by
(7)
For a given outlet quality of the two-phase section (), the length fraction can be obtained from
(8)
Otherwise, if the quality is being iterated for, the residual to be driven to zero by altering the outlet quality of the two phase section () is
(9)
where is evaluated from equation (8), and is known from its solution above. Unfortunately, iterative methods are required since the value of the average refrigerant two-phase heat transfer coefficient is dependent on the two-phase section outlet quality.
In the subcooled section, if it exists, the available area is known (since is known), as are the inlet air and refrigerant temperatures. The inlet refrigerant temperature to the subcooled section is the bubble temperature of the refrigerant, and the inlet air temperature is the same inlet air temperature as for all the other sections. Thus the value can be given by
(10)
and the minimum and maximum capacitance rates of air and subcooled refrigerant can be obtained from
(11)
which give the from
(12)
which yields the effectiveness, if the minimum capacitance is on the air side, of
(13)
or, if the minimum capacitance rate is on the refrigerant side, the effectiveness is given by
(14)
The total heat transfer rate in the subcooled section is given by
(15)
which is negative since heat is removed from the refrigerant.
After all the sections of the condenser have been solved and all the requisite terms determined, the results from each section are grouped together, yielding the overall terms
(16)
If there is a subcooled section, the condenser outlet subcooling is evaluated from
(17)
or if there is no subcooled section, an effective subcooling amount is calculated from
(18)
where is the specific heat of saturated liquid. This effective subcooling parameter is primarily needed to allow for continuous behavior during the solving process. Hopefully the subcooled section exists at cycle model convergence.
Minimal Component Test:
from CoolProp.CoolProp import Props
from Condenser import CondenserClass
from FinCorrelations import FinInputs
Fins=FinInputs()
Fins.Tubes.NTubes_per_bank=41 #number of tubes per bank=row
Fins.Tubes.Nbank=1 #number of baks/rows
Fins.Tubes.Ncircuits=5 #number of baks/rows
Fins.Tubes.Ltube=2.286
Fins.Tubes.OD=0.007
Fins.Tubes.ID=0.0063904
Fins.Tubes.Pl=0.0191 #distance between center of tubes in flow direction
Fins.Tubes.Pt=0.0222 #distance between center of tubes orthogonal to flow direction
Fins.Fins.FPI=25 #Number of fins per inch
Fins.Fins.Pd=0.001 #2* amplitude of wavy fin
Fins.Fins.xf=0.001 #1/2 period of fin
Fins.Fins.t=0.00011 #Thickness of fin material
Fins.Fins.k_fin=237 #Thermal conductivity of fin material
Fins.Air.Vdot_ha=1.7934 #rated volumetric flowrate
Fins.Air.Tdb=308.15 #Dry Bulb Temperature
Fins.Air.p=101.325 #Air pressure
Fins.Air.RH=0.51 #Relative Humidity
Fins.Air.FanPower=160
params={
'Ref': 'R410A',
'mdot_r': 0.0708,
'Tin_r': 333.15,
'psat_r': Props('P','T',323.15,'Q',1.0,'R410A'),
'Fins':Fins,
'Verbosity':0
}
Cond=CondenserClass(**params)
Cond.Calculate()
print 'Heat transfer rate in condenser is', Cond.Q,'W'
print 'Fraction of circuit length in subcooled section is',Cond.w_subcool
print 'Fraction of circuit length in twophase section is',Cond.w_2phase
print 'Fraction of circuit length in superheated section is',Cond.w_superheat
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 [3]: execfile('ACHPComponents/ComponentTests/CondenserTest.py')
Heat transfer rate in condenser is -13289.9922 W
Fraction of circuit length in subcooled section is 0.563034708961
Fraction of circuit length in twophase section is 0.380016896545
Fraction of circuit length in superheated section is 0.0569483944945
If not, first stop should be the Frequently Asked Questions
Nomenclature
Variable | Description |
---|---|
Total air side area (fins+tubes) [m^{2}] | |
Air-side area for superheated portion of circuit [m^{2}] | |
Air-side area for two-phase portion of circuit [m^{2}] | |
Air-side area for subcooled portion of circuit [m^{2}] | |
Total refrigerant side area of all tubes [m^{2}] | |
Specific heat of refrigerant [J/kg/K] | |
Specific heat of refrigerant at dewpoint temperature [J/kg/K] | |
Specific heat of humid air per kg of dry air [J/kg_{da}/K] | |
Minimum capacitance rate [W/K] | |
Maximum capacitance rate [W/K] | |
Ratio of capacitance rates [-] | |
Refrigerant latent heat [J/kg] | |
Average circuit length [m] | |
Mass flow rate of dry air [kg/s] | |
Mass flow rate of dry air in superheated section [kg/s] | |
Mass flow rate of dry air in two-phase section [kg/s] | |
Mass flow rate of dry air in subcooled section [kg/s] | |
Mass flow rate of refrigerant [kg/s] | |
Mass of refrigerant charge (total) [kg] | |
Mass of refrigerant charge in superheated section [kg] | |
Mass of refrigerant charge in two-phase section [kg] | |
Mass of refrigerant charge in subcooled section [kg] | |
Number of thermal units [-] | |
Length of superheated portion of circuit [m] | |
Length of two-phase portion of circuit [m] | |
Length of subcooled portion of circuit [m] | |
Refrigerant-side pressure drop (total) [Pa] | |
Refrigerant-side pressure drop in superheated section [Pa] | |
Refrigerant-side pressure drop in two-phase section [Pa] | |
Refrigerant-side pressure drop in subcooled section [Pa] | |
Total heat transfer rate [W] | |
Heat transfer rate in superheated section [W] | |
Heat transfer rate in two-phase section [W] | |
Heat transfer rate in subcooled section [W] | |
Subcooling amount [K] | |
Air inlet temperature [K] | |
Bubble-point temperature of refrigerant [K] | |
Dew-point temperature of refrigerant [K] | |
Refrigerant inlet temperature [K] | |
Refrigerant saturation temperature [K] | |
Overall heat transfer conductance [W/K] | |
Fraction of length of superheated portion of circuit [-] | |
Fraction of length of two-phase portion of circuit [-] | |
Fraction of length of subcooled portion of circuit [-] | |
Quality of refrigerant at outlet of two-phase section [-] | |
Overall surface effectiveness on air side [-] | |
Air side heat transfer coefficient [W/m^{2}/K] | |
Refrigerant side heat transfer coefficient [W/m^{2}/K] | |
Residual term [-] | |
Effectiveness of heat transfer in the two-phase section [-] | |
Effectiveness of heat transfer in the subcooled section [-] | |
Effectiveness term from derivation [-] |