to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, ... I. A Steady-State. Computer. Design Model for Air-to-Air. Heat Pumps. S. K. Fischer. C. K. Rice ..... consulting engineers, research ...
Printed
in the United States of America. Available from National Technical Information Service U.S. Department of Commerce 5285 Port Royal Road, Springfield, Virginia 22161 NTIS price codes-Printed Copy: ,410 Microfiche ,401
This report was prepared as an account of work sponsored by an agency of the UnitedStatesGovernment. NeithertheUnitedStatesGovernment noranyagency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its usewould not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United StatesGovernment or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United StatesGovernment or any agency thereof.
ORNL/CON-80/Rl
Contract
No. W-7405-eng-26
Energy
Division
The Oak Ridge Heat Pump Models: I. Design Model for Air-to-Air
'/ i
A Steady-State Heat Pumps
S. K. Fischer C. K. Rice
:
DEPARTMENTOF ENERGY Division
,
of Building
Date Published:
Equipment
Auqust
1983
OAK RIDGE NATIONAL LABORATORY Oak Ridge, Tennessee 37830 operated by UNION CARBIDE CORPORATION for the DEPARTMENTOF ENERGY
Computer
.
-
TABLE OF CONTENTS
r
Page LIST OF FIGURES ................................................
vii
LIST OF TABLES .................................................
ix
ABSTRACT .......................................................
1
1.
INTRODUCTION AND BACKGROUND...............................
2.
UTILIZATION
3.
CALCULATIONAL PROCEDUREAND ORGANIZATION .................. 3.1 Modeling Procedure for the Vapor Compression Cycle 3.2 Organization of the Computer Program ................. 3.3 Input and Output Description .........................
4.
5.
6.
3
...............................................
COMPRESSORMODELS ......................................... 4.1 Introduction ......................................... 4.2 General Calculational Scheme ......................... 4.3 Map-Based Compressor Model ........................... 4.4 Loss and Efficiency-Based Compressor Model
7 *I 9 ... 1: 16
...........
FLOW CONTROL DEVICES ...................................... 5.1 Introduction ......................................... 5.2 Capillary Tube Model ................................. 5.3 Thermostatic Expansion Valve Model ................... 5.4 Short Tube Orifice Model .............................
37 37
CONDENSERAND EVAPORATORMODELS ........................... 6.1 Introduction ......................................... 6.2 Single-Phase Heat Transfer Coefficients in Heat Exchanger Tubes ...................................... 6.3 Two-Phase Heat Transfer Coefficients in Heat Exchanger Tubes ................................................ .................. 6.4 Air-Side Heat Transfer Coefficients 6.5 Heat Exchanger Performance ........................... 6.6 Fan Motor and Compressor Shell Heat Losses ...........
47 47
PRESSURE DROPS AND FAN POWERS ....................
65
7.
AIR-SIDE
8.
PRESSURE AND ENTHALPY CHANGES IN REFRIGERANT LINES ........
iii
z; 44
50 51 55 2:
Page .
9.
10.
MODEL VALIDATION .......................................... 9.1 Introduction ........................................ ................................. 9.2 Compressor Modeling .......................... 9.3 Heat Exchanger Calibration 9.4 Heat Pump Validation at 5.4"C (41.7OF) Ambient 9.5 Heat Pump Validation at 10.6"C (51.0°F) Ambient'::::: ................................ 9.6 Cooling Mode Results
73 73 73 74
RECOMMENDATIONS . . . . . . . . . . . . . ..*..........................
83
REFERENCES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
:: 81
85
APPENDIX A.
DEFINITIONS
OF INPUT DATA . . . . . . . . . . . . . . . . . . . . . . . .
89
APPENDIX B.
SAMPLE INPUT AND OUTPUT DATA . . . . . . . . . . . . . . . . . . . . .
95
APPENDIX C.
DEFINITIONS OF CONSTANTS ASSIGNED IN BLOCK DATA SUBROUTINE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
101
APPENDIX D.
DEFINITIONS
OF VARIABLES IN COMMONBLOCKS . . . . . . . .
105
APPENDIX E.
CROSS-REFERENCE OF COMMONBLOCKS . . . . . . . . . . . . . . . . .
131
APPENDIX F.
SUBROUTINE DESCRIPTIONS AND CROSS-REFERENCE OF SUBROUTINE CALLS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
133
ALGEBRAIC NOTATION . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
139
APPENDIX H. DERIVATION OF DEHUMIDIFICATION SOLUTION . . . . . . . . . . H.l Introduction ........................................ H.2 Heat Transfer Equations - Air to Wall . . . . . . . . . . . . . . . H.3 Heat Transfer Equations - Air to Mean Surface Conditions .......................................... - Mean Surface to Refrigerant H.4 Heat Transfer Equations Conditions .......................................... H.5 Derivation of Coil Characteristic and Total Heat Flow Equation ............................................ H.6 Derivation of the Solution for the Exit Effective Surface Temperature, Tf, ,(x0) ...................... H.7 Derivation of the Soluti6n for the Exit Air Dry Bulb Temperature T (x ) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . References for App%aixOH ...............................
147 147 147
APPENDIX I. DESCRIPTION OF THE DEHUMIDIFICATION ALGORITHM . . . . References for Appendix I ...............................
163 170
APPENDIX G.
APPENDIX J.
RATIOS OF GEOMETRIC PARAMETERSUSED IN HEAT EXCHANGERCALCULATIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . .
iv
149 150 151 153 156 161
171
. .
Page .
APPENDIX K. .
AVAILABILITY AND DISTRIBUTION OF THE HEAT PUMP MODEL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . OF COMPUTERMODEL ........................
179 185
APPENDIX L.
LISTING
APPENDIX M.
SAMPLE INPUT DATA ................................
185
APPENDIX N.
LISTING
PROGRAM. . . . . . . . . . . . . . .
185
APPENDIX 0.
SAMPLE COMPRESSORMAP DATA . . . . . . . . . . . . . . . . . . . . . . .
185
APPENDIX P.
LISTING
OF THE INTERACTIVE SUBROUTINES . . . . . . . . . . .
185
APPENDIX Q.
LISTING
OF THE INTERACTIVE PREPROCESSOR. . . . . . . . . .
185
OF THE MAP FITTING
. .
V
.
LIST OF FIGURES Page c
Fig. Fig.
3.1 3.2
Pressure vs. Enthalpy Diagram for the Heat Pump Cycle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10
Block Diagram of Iterative Loops in the Main ....................................... Program
15
Models
21
........
27
for the Model . . . . .
31
Fig.
4.1
Computational
Sequence
in the
Fig.
4.2
Components
Fig.
4.3
Iteration on Internal Energy Balance Loss and Efficiency-Based Compressor
of Compressor
Energy
Compressor Balance
Fig.
6.1
General
Structure
of the
Condenser
Fig.
6.2
General
Structure
of the
Evaporator
Fig.
6.3
Computational
..
57
Fig.
6.4
Computational Sequence in the Evaporator Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
58
Sequence
Fig.
I.1
Block
Fig.
J-1
Sample Tube-and-Fin
Diagram
of the
in the
Model
Condenser
Dehumidification Heat Exchanger
vii
Model
.......
48
......
49
Model
Algorithm
164
.............
172
.
*
. .
.
LIST OF
TABLES
P c
Page .
c
Table Table Table
9.1 9.2 9.3
Comparison Predictions
of Experimental of Compressor
Data with Map Model
............
75
Heat Pump Model Validation at 5.4'C (41.7"F) ........................................ Ambient
77
Heat Pump Model Ambient Without
79
Validation Corrections
at 5.4"C (41.7OF) to Compressor Map . .
Table
9.4
Heat Pump Validation
Table
E-1
Cross-Reference
of Common Blocks
Table
F-1
Cross-Reference
of Subroutine
Table
K-1
Sample Pages of Input Data for Interactive Version of Heat Pump Model . . . . . . . . . . . . . . . . . . ...183-184
at 10.6OC (51°F)
Calls
Ambient
..
80
...............
132
............
137
I. A Steady-State The Oak Ridge Heat Pump Models: Design Model for Air-to-Air Heat Pumps* S. K. Fischer C. K. Rice Energy Division Oak Ridge National Laboratory Oak Ridge, Tennessee 37830
Computer
ABSTRACT The ORNL Heat Pump Design Model+ is a FORTRAN-IV computer program to predict the steady-state performance of conventional, vapor compression, air-to-air heat pumps in both heating and cooling electrically-driven, This model is intended to serve as an analytical design tool for modes. use by heat pump manufacturers, consulting engineers, research instituand universities in studies directed toward the improvement of tions, The Heat Pump Design Model allows-the user to heat pump performance. specify: 0
system
l
compressor
0
refrigerant
0
fin-and-tube
l
fan
l
any of ten
The model will
operating
conditions,
characteristics, flow
control
heat
and indoor
exchanger
duct
parameters,
characteristics,
and
refrigerants.
compute:
system
capacity
compressor
and COP (or
and fan motor
coil
outlet
air
dry-
air-
and refrigerant-side
a summary of the the cycle, and overall compressor effectiveness.
.
devices,
EER),
power consumptions,
and wet-bulb pressure
refrigerant-side efficiencies
This report provides thorough documentation This is a revision of an earlier the model. corrections and information on availability including an interactive version. *Research sponsored by the Office U.S. Department of Energy under Carbide Corporation. 'The ORNL Heat Pump Design Laboratory.
temperatures, drops, states and heat
throughout exchanger
of how to use and/or modify report containing miscellaneous and distribution of the model -
of Buildings and Community contract W-7405-eng-26 with
Model was developed 1
at Oak Ridge
Systems, the Union National
.
3
INTRODUCTION AND
1.
The ORNL Heat Pump Design developed
to predict
compression,
an analytical
consulting
design
engineers
directed
toward
the
The current
[6].
routines
in Ref.
in
tool
for
of heat
has evolved
served
1978,
improved
version
informal
documentation.
model.
indicated
Many of these
This
report
rections
is
identified
the
greatest
extent
correlations
derived
and algorithms critically appropriate
correlations
The model earlier
has thus
versions
been reorganized be able
by vertical
bars
principles
to readily
[7,8]
and documented
as well
An limited, to a
as the
capabilities with
and generalized the
experiences in the
1979 and 1981. miscellaneous
cor-
or additions
have
limitations
and efficient
standpoint.
into
format specific
The correlations
to include
from an engineering
to
of empirical
of the model
improved
pro-
correlations
literature.
or replaced
to their
pump modeling
has been to base the
been significantly
it
and Kusuda
in the left margin.
and more accurate
adapt
heat
corrections
releases
a more modular
[5],
between
so as to avoid
and improved
Flower
makes use
in 1979 with
version
from manufacturers'
reviewed
and also
for
additional
where
used in the previous
[3]
have been distributed
programs
The places
possible,
to
institutions.
of an earlier
physical
is
at ORNL [l,Z]
ORNL Heat Pump Model.
of the model development
gram on underlying
model
in studies
was modified
needs were incorporated
a revision
The philosophy
point
of the
the need for
and additions.
been made are
[4],
program
,and research
heating
pump manufacturers,
written
code was made available
Our own use of the previous users
programs
starting
version
vapor
pump efficiency.
These two versions
number of manufacturing
of this
and universities
and Erth
as the
program
pumps in both
of Technology
the original
of this
heat development
from
Institute
computer
of conventional,
use by heat
improvement
1 as a preliminary
of other
the
by Kartsounes
The MIT model
at ORNL, where
for
institutions,
and at the Massachusettes of selected
performance
, research
model
Fortran-IV
a
air-to-air
The purpose
modes.
is
steady-state
electrically-driven,
and cooling provide
the
Model
BACKGROUND
have been
more recent algorithms.
and expanded Furthermore,
so that needs.
or
reseachers
from
the
it will
has
4
The ORNL Heat Pump Design e
.
System Operating 1. 2. .
0
0
0
0
l
Model
allows
the
user
to specify: II
Conditions
the desired indoor relative humidities,
and outdoor
air
temperatures
*
the arrangement of the compressor and fans up or downstream of the flow stream, i.e.,
Compressor
in the air heat exchangers;
Characteristics
1.
either a map-based equipment,
2.
or an efficiency compressors;
Refrigerant
model
for
and loss
Flow Control
designs
model
a capillary a short-tube
2.
a specified value of refrigerant at the condenser exit (in this the equivalent capillary tube, parameters);
Fin-and-Tube
with
for
tube, thermostatic orifice, or
Heat Exchanger
advanced
expansion
reciprocating
valve
(TXV),
or
subcooling (or quality) case the program calculates TXV, and short-tube orifice
Parameters
1.
tube size, spacing, parallel circuits,
and number of rows,
2.
fin pitch, thickness, fins (smooth, wavy,
and thermal or louvered),
3.
air
flow
available
Devices
1.
Fan and
and
and number of
conductivity;
type
of
rates;
Indoor
Duct Characteristics
1.
overall fan fans, or
efficiency
2.
a specified
fan
3.
the
diameter
values
efficiency
of one of six
for
curve
indoor for
the
equivalent
and outdoor outdoor
fan,
ducts;
Refrigerants 1.
either
2.
one of six additional refrigerants (Rll, R113, and C318) by adding the appropriate constants given by Downing [9];
R12, R22, R114, or R502 as a standard
option,
R13, R21, R23, thermodynamic
5 Refrigerant
0
1.
lengths
2.
pipe
specifications
3.
heat
losses
The user is
Lines
cannot
implicity
specify
satisfactory
the
specified
model
for
heat
a low level
charge
over
unless
the
known.
a range level
of ambient
and liquid
with
the
having
instead,
proper This
a suction
excess
is
line
the
performance
superheat
however,
model
and output
has evolved
system
from
a accumulator
refrigerant with
cannot
and
a given
be accurately
at each ambient
single
of measurement
of the
computer and it
design
would
The specific
parameter
in Appendices
the
user
output
are given to control
includes
the
refrigerant modeled
temperature studies
can still
English
units.
(The
Sensible Outlet Air-
dry-
effort
to convert
each input
Four
output
output.
British
and output
options
The basic
allow program
heat
power consumption transfer
and wet-bulb
isentropic
of program
ratio
temperatures pressure
and volumetric effectiveness
capillary
Summary of the
A.
used for
A and B.
and refrigerant-side
Equivalent parameters
Appendix
units
and fan motor
Heat exchanger
levels
based on the
have been a major
amount of printed
to total
Overall
Other
were
and COP
Compressor
l
are-in
tube,
output
drops compressor
efficiencies
and UA values TXV, and short-tube
refrigerant-side are discussed
states
is
point
(as appropriate):
Capacity
0
program codes that
them to SI units).
it
amount
be made. The input
mode,
lines.
charge;
conditions.
temperatures
of evaporator
or cooling
superheat.
systems,
For such systems,
charged
can store
of evaporator
For charge-sensitive
is
operating
which
of heating
refrigerant
pump systems
systems)
pipes,
discharge,
system
system
for
(charge-insensitive
suction,
the the
of interconnecting independent
from
assumed that
of refrigerant
maintain
and diameters
orifice
throughout
in Section
3.2
the and
cycle
6
This from
its
second
generation
predecessors
heat
pump model contains
and many additional
significant
capabilities.
This
changes report
is
intended: e
to serve as a rigorous documentation (where correlations, and algorithms the assumptions, the program, to provide equations
e
A basic
a complete list of references were obtained, and
to serve as a user's manual for intact and to aid those who will understanding
been presumed
for
of vapor-compression this
report.
appropriate) of that are used by
from which
data
and
those who will use the code need to add new capabilities. air-to-air
heat
pumps has
7 2. UTILIZATION The physical outdoor
operating
described
in
less
complete
listing
Although program,
software. numerical optimization
than
are read
A; tape
The program
executes
of the
conditions
in Appendix
required.
alone
specification
drives
requires
5
of the
and other
seconds
on the
program
is
provided
is generally
and graphical
display [7,8].
file,
and as
are not
computer
in Appendix
of the model studies
data
peripherals
IBM 3033
Model
indoor
of memory and typically
version
and parametric
pump and the a single
has been used in conjunction
An earlier analysis
from
56 K words
the Heat Pump Design it
heat
with
at ORNL.
A
K. run as a stand-
other
computer
has been used at ORNL with
programs
in
heat
pump design
P CALCULATIONAL PROCEDUREAND
3. 3.1
Modeling
The heat sections. flow
pump model device
routines
into
low-side
unit,
Calculations
the proceed
overall
independent
the
the
into
compressor,
an interrelated contains between
is obtained.
of whether
Compression
functionally
combines
iteratively
balance
the Vapor
is organized
section
section,
desired
for
The first
control
second
Procedure
ORGANIZATION
two major condenser,
high-side
the these
is
operating
the
basic
and
unit.
evaporator
model.
two sections
until
The calculational
unit
Cycle
in
the
The the
scheme is heating
or cooling
mode. Figure
3.1
an exaggerated The user
is
pressure-enthalpy required
vapor
compression that
diagram,
cycle,
is modeled
shown on
by the
program.
to specify:
l
the
l
design parameters for a flow-control device condenser exit subcooling (or quality),
l
condenser
l
dimensions
l
heat
The user
level
of evaporator
exit
and evaporator of components
losses
must also
from
provide
superheat
inlet
air
(or
quality), or the
estimates
of
temperatures,
and interconnecting
interconnecting
level
pipes,
and
pipes.
for:
e
the
refrigerant
saturation
temperature
at the
evaporator
exit,
l
the refrigerant and
saturation
temperature
at the
condenser
inlet,
e
the
mass flow
These
three
final
results
values
represents
refrigerant
estimates
used
are
of the model
are chosen;
as starting
do not
however,
rate.
better
points
for
iterations.
depend on how well guesses
will
The
or how poorly
yield-shorter
these
running
times.
High-Side
Computations
The computations refrigerant which
state
is defined
for
at the by the
the exit
high-side of the
specified
system evaporator
superheat
balance
begin
(point
a in
and the
estimate
with Fig.
the 3.1),
of the
10
I(
.
ORNL- DWG 81-4749
ENTHALPY a
EVAPORATOR EXIT b COMPRESSOR SHELL INLET C COMPRESSOR SHELL OUTLET d CONDENSER INLET e CONDENSER OUTLET
f g
INLET TO EXPANSION DEVICE EVAPORATOR INLET Fig.
3.1
Pressure vs. Enthalpy Heat Pump Cycle
Diagram
for
the
11
saturation the
This
temperature.
compressor,
condenser,
The compressor
model
point
remains
expansion
uses
state
fixed
device,
one iteration
and evaporator
a along
point
for
the estimates of the refrigerant mass fl ow rate temperature at the condenser inlet, and
l
the specification of the dimensions suction and discharge lines the
condenser
Fixed
state
at the
d, as well
inlet,
shell
calculations.
with:
l
to calculate
and heat
and saturation
losses
b, she1 1 outlet,
inlet,
of
as a new estimate
for c,
the
and the
of the mass flow
rate.
Condenser Subcooling The remainder
flow
control
exit
subcooling
described
of the
device (or
quality)
since
l
the
physical
l
the
calculated
0
the
condenser
l
the
refrigerant
the
of condenser compared a fixed
states
inlet
compressor
b, c,
condenser
device,
f,
is
value.
computed
submodel
and relative inlet,
of state
the
two values
do not
the
saturation
calculations
are
achieved,
is
e, and agree
within at point
pressure),
changed,
mass flow updates the
The level
repeated.
temperature model
e.
entrance
and
d,
temperature
a new condenser
condenser
uses:
point
knowledge
C),
is
humidity,
from
a new refrigerant is
case
rate,
condenser
saturation
subcooling
The latter
outlet,
specifying
calculates
of condenser
condenser
If
a
value
at the
is computed
and d and the
desired
at the
on whether
exchanger,
temperature
state
condenser
heat
mass flow
air
and condenser
the
model
of the
state
effect
a desired
The condenser
refrigerant
(see Appendix
(in
Each time
simpler.
description
specified
tolerance
compressor
is
subcooling
d, is changed and the
it
or if
depend
has been specified.
refrigerant
to the
calculations
has been chosen
first
to evaluate
high-side
state
the
rate
and new
e.
Once the
state
at the
flow
control
using:
l
the
state
at the
l
the
dimensions
l
the most recent
condenser
and heat calculation
loss
exit, for
e, the
of the
liquid
refrigerant
line,
and mass flow
rate.
Refrigerant
state
the
equivalent
that
would
f and the mass flow
capillary produce
Specified
When condenser has a flow the the
the
exit
heat
pressure
is
expansion
rate
This
balance
characteristic
refrigerant turn,
is
flow
rate
found
obtained
flow is
control fixed.
inlet the
at the
at the
condenser
liquid
0
the
line
dimensions
computed
mass flow entrance
required,
the
rate
outlet
expansion
at the
evaporator
pressure
at point
in mass
exit,
subcooling
for
the
Since
estimate the
proceed
of
condenser
to the
using:
e,
the
loss,
and
rate.
flow
calculated
control
device
using
state
valve
or short-tube
entrance,
point
g is
a specified
and condenser
inlet.
mass flow
is
which,
with
f is evaluated
state
which
conditions
case of a thermostatic pressure
shell
and heat
compressor
on the
condenser
the calculations
state
0
low-side is
been specified,
the
the mass flow
refrigerant
balance
condenser
where
the
the
with
high-side
device
case where
compressor
entrance
of the
model.
to the
device
e
using
a high-side
at the
high-side
dependent
e, d, and e, are computed
states temperature
control
control
model
compressor
The states
The refrigerant the
the
is
by
of the
the
since
device
to flow
to achieve
has not
based on its
entrance
are similar
subcooling flow
flow-control
balance
capability
proper
iteratively
is,
(constrained
model must be coupled the
it
that
by the
characteristics
to achieve
as if
opening,
specified,
drop
control
device
saturation
is
acts
condenser
pumping
device
condenser
from
and outlet,
of the
must be found
the
The calculations
determined
routines
at the
from
is completely
pressure
flow
of the
state
parameters
system
variable
refrigerant
by the
so the
the
an infinitely
control
and condenser
balance.
specified,
capability
and the
controlled
device,
compressor
with
rejection
When a flow
compressor.
is
pressure
subcooling),
orifice
Device
device
system
used to calculate
subcooling.
subcooling
control
high-side
between
specified
Flow Control
are then
TXV, and short-tube
tube,
the
rate
initially
g.
will
f and, orifice,
If
pass
the
in the
latter
assumed to be equal
13
to that then the
at point
later
calculations If
estimate.
and flow
in the
the
control
device
a new estimate
balance
do not
of the
computing
Low-Side
agree
two mass flow
rates
of the
calculations
temperature are
compressor tolerance
are repeated
at the
condenser
in agreement,
analysis
the
a high-side
can now be completed
by
performance.
or low-side,
computations
l
the
refrigerant
condition
l
the
refrigerant
enthalpy
l
the
refrigerant
mass flow
These values
have all
temperature
and specified
the
inlet
the
assumed exit is
air
temperature
executed
one execution evaporator
exit
would
quality)
pressure
drop
at the
at
evaporator yield
a,
inlet,
point
g,
the
the
and
inlet
saturation evaporator point
specified
inlet
air heat
since
the
temperature
from
at the
exchanger exit
at
The evaporator
of superheat
the
g, and
superheat
unknown.
a value
across
the
inlet,
are still
to calculate
and a pressure
hence a saturation
(or
varying
next,
state
based on the estimated
temperature
iteratively,
evaporator
exit,
rate.
at the
which
saturation
to the
at the
superheat
pressure
are based on:
at the evaporator
been computed
The saturation
model
are used to improve
device rates
and
Computations
The evaporator,
exit.
model
evaporator),
a specified
One cycle
the evaporator
the
within
saturation
has been obtained.
in
mass flow
and expansion
Once the
entrance.
drop
evaporator
refrigerant
condenser,
compressor, with
no pressure
a (i.e.,
(and
conditions
are
fixed). A system air
temperature,
calculated given
tolerance. flow
has been completed
though
exit
thermostatic the
solution
not
superheat
agrees
(The highexpansion
control
device
pressure
is used during
solution
is found
changing
the
for
saturation
necessarily
valve
the
with
the
or a short-tube that
high-side
desired
temperature
some evaporator
the desired
and low-side
to ensure the
for
one, value
loops
repeated
orifice
an accurate
is
at the
inlet evaporator
within
a
once if
specified
evaporator
calculations.)
evaporator
when the
specified are
inlet
a
as
inlet
The system air
temperature exit,
point
by a,
14 and repeating state
point
agrees
with
sequence inlet
entire
calculational until
the
the
value
within
is
summarized
desired
air
temperature
is
temperature
outermost
loop
as the
distinct into
distinct
for
example,
the
of these
other.
This
subroutines,
is
shown in
components
one for
the
the
accomplished
flow
in the
the
compressor
l
modeling
the
condenser
l
modeling
the
flow
*
iterating
iterations
properties,
pressure
(Appendix
F contains
control
on condenser
each of these
inner
The evaporator of the
exit
refrigerant over
Program
3.2
are divided
these
computations
of the
system
performance control
into are divided
high-side
of the
compressor,
device
and then
interconnecting
pipes
heat
pump model
balance,
using
balancing with
each
individual
each task:
modeling
perform
Fig.
and the
l
In addition
The
one or two iterations
to perform
computing
and (optionally)
temperature
3.2.
Computer
The calculation
condenser,
on
tolerance.
function
of the
subroutines
requires
air
in Fig. only
iteration
are required.
Organization
modules.
output
3.2
inlet
a specified
usually
calculations
sections,
computed
a linear
so that
in Fig.
3.2 Just
nearly
This
process.
a continues
of calculations
saturation the
the
device saturation
subroutines and calculate
drops , and heat a cross-reference
temperature
calls
other
thermodynamic transfer of all
service
routines
to
and thermophysical
coefficients, subroutines
etc. called
from
each subroutine.) The advantages
to this
program
and identifiable
structure tasks
are that:
0
separate
are performed
by each subroutine,
0
modules can be easily modified to make improvements, replaced to use models for a different component type, or added to increase system analysis capability
l
familiarization with the code can be done at different levels, depending upon the needs of a specific user.
‘
ORNL-DWG
81-4741
READ INDOOR AND OUTDOOR AMBIENT TEMPERATURES AND ESTIMATES OF THE CONDENSING AND EVAPORATING TEMPERATURES
TXV, CAP-TUBE ,OR ORFICE
SPECIFIED
FIND THE CONDENSING TEMPERATURE T sti,o, WHICH BALANCES THE COMPRESSOR AND EXPANSION DEVICE MASS FLOW RATES
SUBCOOLING
FIND THE CONDENSING TEMPERATURE T so+,e, WHICH GIVES THE SPECIFIED SUBCOOLING
IENT TEMPERATURE
Fig.
3.2
Block
Diagram
of
Iterative
Loops in the
Main Program
r
x
I
16
One consequence large
of this
program
variables E are
passed
between
them through to identify
which
common blocks
As part data
with
that
a single
a large
common blocks. variables
number of
Appendices
D and
used in common blocks
are used to pass data
input
and output of the There
output. not
read
output
than
data,
between
amount
of data
the
output
results
The program
Appendix
A.
The first
computed
results energy
summary.
in cases out
where
the
proper
been made.
Intermediate
with
of
either
Appendix
A, is
standard
iteration
sequence
to another
other
operation
B is a system in
are for
3.2,
of computations.
variable,
to print
the
heat
Further,
in
pump and
diagnostic
pur-
or for
use in
changes
have
are printed
out
as defined abbreviated
B, a summary after
or continuous
has
one page of
of the
after
3.2
as defined
to converge
An input
summary as in Appendix
Fig.
to only
in the computations whether
summary and
options,
two options
two options.
in Fig.
loop
output
fails
the
them as necessary.
of the program
results
to reduce
can change
COP and capacity loop
are assigned These parameters
the output
The other
to be
each run.
for
outermost
limits
used to specify
output
of each major entire
these
has three
an iterative
and
parameters,
are unlikely
program
the user
the
the
input
of printed
tolerance
which
must specify
C so that
gives
sample
levels
(BLOCK DATA) in order
option
which
(e.g.,
pressure)
after
A contains
used by the model which
shown in Appendix
computed
of the
B.
constants
user
all
Appendix
are other
one use of the
in Appendix
The printed
almost
B contains
constants
and there
subroutine
that
are described
checking
other
atmospheric from
in a single
converged.
and Appendix
and physical
number,
Description
of the program,
shown in Appendix
frequently
a brief
data
as input
values
poses
input
however,
that
refrigerant
contains
and Output
are done by two subroutines.
are,
Parameters changed
Input
of the modularization
definitions
the
is
subroutines. 3.3
the
however,
subprograms
as aides
and to identify
are
approach,
becomes many short
included
which
modular
in output,
completion
output
throughout
the output
statements
17 in most of the from print
a particular control
subroutines
are arranged
subroutine statements
so that
can be obtained
in that
routine.
detailed,
by removing
diagnostic one or two
output
19 4.
COMPRESSORMODELS 4.1
Since primary
the
user
compressor
of electrical
to good system be tempered
of the
models
design
parameters.
models
depending The first
Model
manufacturers' type;
for
Accurate
available
to
the
has built-in
superheat
those
for
for
use with
rotary,
simulation
between
two simpler
of existing
compressor
corrections
to adjust compressors
generated. it
or for
which
Although
should
or centrifugal
compressors
the
use of compressor
compressors,
screw,
heat
which
hardware
a specific
the maps were
reciprocating
reasons
detailed
in reciprocating
which
of
a subroutine
based on the
maps) for
types
For these
using
needs.
the model
different
incorporate
specific
model was written modify
information
upon their
of refrigerant from
of
can choose
(compressor
important must
users
data
is
however,
may be used.
is
modeling
criterion,
performance
model
and the
This
Instead,
compressor
pump system
compressor
and of
does not
compressor
a heat
type
program
the program
ORNL Heat Pump Design
of
accurate
of the
in which
different
power,
by consideration
pump studies
levels
heart
prediction.
users
compressor
the
performance
most potential
rigorously
is
Introduction
are this
be easy
to
compressors.
is
possible
with
this
model. The second is
intended
changes It ant.
for
model,
a loss
use in heat
in compressor
can also This
more general
efficiency
values.
This
over
the
based model without
performance
routine
values.
improvements about
using
same range
is well
suited,
and interactions
a particular
the predict
design
energy
heat
conditions
for
and their point.
loss
compressor
of some of the however,
system with
internal
model,
to predict
affect
performance
of operating
adjustment
e.g.,
terms
user-supplied
model cannot
local It
models
compressor
studies,
and efficiency
compressor
as accurately
compressor
pump design
be used to model compressor
a reciprocating
and loss
loss
and efficiency-based
how
performance.
a new refrigerbalances
in
and internal performance as the map-
input
efficiency
studying
internal
effect
on system
20
4.2 Both
compressor
models
with
rest
and interact general
the
calculational
tions
of those
temperature refrigerant
the
sequence
frigerant
degree
at the at
pressure
volume
refrigerant
procedures
inlet
exit
Figure
the evaporator
exit
is
the
estimated
shows simulation.
outlet
the
condenser
entrance.
using
and the
refrigerant
re-
corresponding
identified
or quality
the refrigerant
4.1
and evaporator and the
superheat
and on
of compressor
are used to calculate
in the
suction
8 using
the
the
calculated
saturation
temperature,
enthalpy,
suction
point,
gain
and discharge current
refrigerant
at the‘compressor
of heat
From this
exit.
exit
and average
calculated
value
evaporator
saturation
entropy,
are computed.
drops
state
refrigerant
evaporator
(from
in Section rate
condenser
by descrip-
are different.
on the
method
The
same way.
followed
and the
either
in common
in the
which
evaporator
for
from which
as described
the
at the
temperatures
The pressure mass flow
entrance
of evaporator
exit
and specific
two routines dependent
of the
state
temperature),
estimates
temperatures shell
line
inlet
pressure
in the
each of the
lines
is
drop
suction
for in
are computed the
the
refrigerant
lines.
then
identified
and the
specified
The
line.
compressor
models
use different
to calculate:
0
the
refrigerant
0
the
compressor
motor
input
0
the
compressor
shell
heat
0
the
refrigerant
the
previously
(input)
of the
pressures
evaporator
the
first,
estimates
The refrigerant
(input)
scheme is discussed
condenser
saturation
specified
of the
of calculations
refrigerant
using
model
superheat
The current
Scheme
have a number of calculations
are functionally
at the
the
Calculational
system
sections
Both models
Finally
General
value
refrigerant calculated of discharge
mass flow
state state
power, loss
at the
(optional), compressor
at condenser
discharge line
rate,
heat
line loss.
entrance pressure
and shell is drop
exit. computed
using
and the
specified
21 ORNL-DWG
(
COMP OR CMPMAP
Sb4745R
J
CALCULATE THE PRESSURES AT THE CONDENSER INLET AND THE EVAPORATOR OUTLET /
c 1
r
3
CALCULATE REFRIGERANT PROPERTIES AT THE EVAPORATOR EXIT _I
\
r .
1 , CALCULATE b P IN THE DISCHARGE LINE AND TqcmpAND Pcut,cmp J t 7 CALCULATE AP IN THE SUCTION LINE AND Tin,cmp AND Pin,cmp + CALCULATE REFRIGERANT CONDITIONS AT THE SHELL INLET
LOSS AND
I
EFFICIENCY-BASED
6 SUCTION GAS ENTHALPY ITERATION FROM FIG. 4.3
,
MAP -BASED t CALCULATE POWER AND RATE FROM AND SHELL
ti CALCULATE THE REFRIGERANT STATE AT THE SHELL OUTLET
.
AND CORRECT MASS-FLOW CURVE FITS HEAT LOSS I
CALCULATE THE REFRIGERANT STATE AT THE SHELL OUTLET .
i 1
f
I
I
CALCULATE THE REFRIGERANT STATE AT THE CONDENSER INLET
r CALCULATE NEW AVERAG TEMPERATURES AND t- SPECIFIC VOLUMES IN THE COMPRESSOR LINES
CALCULATE OVERALL VOLUMETRIC, ISENTROPIC,AND SUCTION GAS HEATING EFFICIENCIES \ i END
Fig.
4.1
Computational
Sequence
in the
Compressor
Models
22 Upon completion and specific
of these
volumes
in the
These are used with flow
rate
suction
latest
to recalculate The entire
drops. the
the
drops
agree
After
completion
new average
and discharge
calculation
the
process
pressure
calculations,
suction
is within
lines
of the
are computed.
refrigerant
and discharge
repeated,
temperatures
line
pressure
as shown in Fig.
tolerances
from
mass
4.1,
until
one iteration
to the
next. ficiency
indices
of the
pressure
are computed.
In both
efficiency
indices
are calculated
efficiency
and the
overall
is
used to refer
appropriate)
The overall the
outlet
compressor
compressor
models,
two basic
isentropic
compression
The term
efficiency.
based on compressor
shell
ef-
"overall"
inlet
and (when
conditions.
isentropic
compression
efficiency,
nCm is
given
by
equation
n cm "k
flow
Qcm,ideal
rate
represents (ideal)
the
outlet
isentropic
compressor specific
specific inlet
and the
a given
mass per unit
outlet, completely
ratio
be emphasized houtlets
is
insulated
The overall
not
)
inlet
(4.1)
of
that
refrigerant
would
time)
used
in Eq. 4.1; would
pressure
i r,actual %ol = fir,ideal
refrigerant
nvol,
_
is
tilr,actual
given
'inlet DS
case, of
outlet.
(to
required
the
an
entropy
at shell
enthalpy
were this
irom
refrigerant
have an ncm value
efficiency,
isen
be obtained
to the actual
actual
mass
houtlet
the minimum power required
the
volumetric
The term
based on the
that
compressor
power and refrigerant
enthalpy.
process
the
refrigerant
motor
enthalpy
compression
shell
-h
R cm,actual
Thus T-Icm represents should
(houtlet,isen
cm,actual
and h represents
at compressor
lilr,actual
=
i cm and mr represent
where
It
iteration,
- the overall
volumetric
to a value
shell
drop
compress power.
at shell any 1.0.
by
(4.2)
23 where
v
inlet, rated
inlet D is
is
the
the
total
compressor
with
actual
motor
and efficiency conditions
compressor as discussed
The next compressor methods,
for
measurements provide and/or
performed
saturation condenser maps,
"evaporator"
1 in Appendix
loss port
the
two
calculational
is
of
for
the
superheat
to six
shell
outlet).
fixed
with
Given
in
gas temperature the
or fixed
of
For some
changes
suction
"condenser"
values
superheat.
to vary
constant.
saturation
four
for
curves
rate
"evaporator"
inlet
while
curves
calorimeter
mass flow
compressor
allowed
held
performance
compressor
inlet)
at the
inlet
rate
condenser
suction
exit
gas temperature,
values for
can be derived from refrigerating ,, performing these calculations is given
N. routine
(kW) and the
of compressor
shell
the
published
performance
Cl'
c,,
**a,
the
Model
map is generated
is
A routine
The map-based input
shell
temperature inlet
mass flow data.
for
between
refrigerant
and compressor
and compressor
refrigerant
power
(i.e.,
saturation
subcooling
be confused
These performance
as functions
superheat
shell
from
input,
at compressor
the
at compressor
capacity
power
subcooling
however,
model.
and discharge
uses empirical
obtained
each performance exit
not
requirements,
by the manufacturers.
temperatures
Usually
should
efficiencies
Compressor
model
capacity
(i.e.,
and efficiency
differences
input
Map-Based
motor
refrigerating
temperature
3450 or
output.
compressors
compressor
the
in
compressor
reciprocating
(e.g.,
the
4.4.
describe
4.3 The map-based
loss
based on suction
- differences
and resulting
the
ncm and nVol
in Section
two sections
models
model
compressor
model
shell
the map-based
speed for
internal
at compressor and S is either
indices,
defined
volume
displacement,
speed for
two efficiency
similarly
specific
compressor
motor
3500 rpm) or the These
refrigerant
for
uses curve refrigerant
inlet
bi-quadratic
to the
mass flow
and outlet data.
fits
saturation
The user functions
compressor
rate
(lbm/h)
temperatures
must provide
of the
form
sets
given
motor as functions to model of coefficients;
by Eq. 4.3
24 for
the
outlet
power
input
saturation
and mass flow
rate
as functions
of the
inlet
and
temperatures, .
f(T outlet'
T inlet)
=
'1 Tiutlet '4
I
A short
computer
compute
these
below,
these
rated
must also
or temperature
at the
based)
compressor
for
the
displacement
is
also
map-based
model
represent
a compressor
as the
original
fits model the
specify
the
total
actual
and the
fixed
speed,
an input
to scale
with
Dabiri compressor rate,
the but
and Rice motor
which
factors
at actual power
presented
input,
they
displacement,
refrigerant
superheat the map is
The desired
value
is
compressor
used by the
curves
performance
for
the
linearly
to
characteristics
biquadiatic superheat
were generated. empirical
level
or
The map-based
curve
fits
to model
conditions. a technique
i
for
and the cm,map' of superheat or suction
than give
their
factors
compressor.
capacity.
m for values r,map those for which the maps were generated. correction
as noted
compressor
maps and the only
to
Except
performance
to the
operating
[lo]
this
of a different
applicable for
N.
(upon which
same general
(4.3)
'6
algorithm
modeled.
compressor
'
+
one particular
inlet
being
the compressor
correction
compressor
is
only
parameter;
the
gas temperature applies
shell
which
to them are strictly
suction
in Appendix for
compressor
earlier,
squares
apply
compressor
As noted
listed
Tinlet
will motor
+ '3 TLet
+ '5 Toutlet
uses a least
is
coefficients
compressor
TinIet
which
coefficients
The user the
program
+ '2 Toutlet
to account
for
correcting
refrigerant
the mass flow
gas temperature Equations
non-standard
4.4
other and 4.5
superheat
values,
rilr,actual
= [l
t F (Vmap v v actual
- 1)-J lil
r ,maP
(4.4)
25 and
(4.5) where of
and Ah represent
v
the
"actual",
conditions,
port
suction
1 respectively, a value
"map",
map superheat
estimated
and "isen"
conditCons
in the
the suction
shell
inlet
conditions
using
gain
to the
gas between
Ah
from
of a constant inlet
(assigned
C).
conditions
shell
pressure, factor
see Appendix
assumption compressor
outlet
from
the
enthalpy and suction
, of 21 kJ/kg (9 Btu/lbm), due to motor and port by the suction gas [lo]. This assumption has been
inlet,suction compressor cooling improved
the
superheat
correction
port
The
process
shell
efficiency
respectively,
conditions. actual
to compressor
Data subroutine,
estimated
port,
represent
Block
and Rice
port
change,
and an isentropic
a volumetric
Dabiri
suction
and enthalpy
suction
conditions,
and Fv is
of 0.75
volume
based on estimated
refrigerant
subscripts
specific
upon in our map-based
model
by assuming
instead
that
.
Ah
where
Fs,
is
inlet,suction
an appropriate
a value
of 0.33
in the
Eq. 4.6
assumes
that
inlet
to suction
temperatures compressor
input
has been observed Once the to the shell
Block will
and higher work
values
suction
gas heating
condensing per unit
corrections
at lower saturation
mass flow reciprocating
for
actual
houtlei, h
outlet
is = (A
compuied
cm,actual
Fsh is
assigned
(see Appendix
C).
Use of
change (from evaporating
rate
from
- 'can"'r
as the
increases.
compressors level
enthalpy
shell
saturation
temperatures
superheat
of icm map and mr map, the
outlet,
factor.
gas temperature
increase
in typical
r4w
Data subroutine
the suction
port)
= Fs, +
port
Such a trend [ll].
have been applied
at the
compressor
Eq. 4.7
, actual
+ hinlet
(4.7)
26 where
is
i,,
the
specified
by the
fraction
of actual
are
not
usually
in
at the
the
loss
user
at the
power,
by the
calculations
then
proceed
as described
in Section
outer
balances
supplied
design , internal of the 4-2.
of (j
can
manufacturers.) at this
thermodynamic
are computed pressure
Loss and Efficiency-Based
energy
Fig.
is
or as a specified (Values
relevant
entrance
to the
value
i,,
properties
next.
drop
The
convergence
loop
4.2.
internal
shown in
input
shell.
has been identified
the
and efficiency-based
representation
compressor
icrn actual.
exit
and all
and condenser
The loss
the
compress&
compressor
exit
4.4
from a fixed
compressor
calculations
shell
rate
as either
provided
The state point
heat
Compressor
compressor
in a reciprocating efficiency,
loss
routine
models
the
compressor
from
user-
and heat-loss
and efficiency-based
The user
is
Model
required
values.
A conceptual
compressor
to specify
the
model
is
following
values: D
-
total
c
-
actual
S input
-
synchronous appropriate),
i
s,fl
-
4.8
motor
maximum compressor
n.lsen
isentropic discharge
low-side
compressor
are also
full
mechanical
loss
gas,
specified
from
compressor
motor
speed
load
power
(optional),
shaft
efficiency compressor
and the user
(as
efficiency,
efficiency,
the
i hilo, by the
ratio,
motor
compression port.
heat
refrigerant
volume
or actual
nmot,maxnmech -
displacement,
clearance
compressor
The interna 6 can,
compressor
and from
discharge
compressor
according
suction
shell
to the
to
tube
to the
heat
loss,
choices
in Eqs.
and 4.9. specified . 4 hilo
=
Qhilo 0.03
directly,
'cm,actual' i
cm,actual
or
(see Appendix
A)
(4.8)
27
ORNL-DWG
1
IIII I
77-19443R
-DISCHARGE TUBE
hr
SUCTION TUBE
houtlet
. my
hsuction
mr
port
Fig.
4.2
Components
of Compressor
hdischorge
Energ.y Balance
port
28
specified
I
=
' can
'can
'cm,actual'
Oeg (lso
The explicit the
program
directly,
losses,
- rlmotor i,
as described
The unknowns
which
ih
suction
h h
port
discharge
port
outlet
(see Appendix .
Or 'mech)
(4.9)
'cm , actual ~1, are part
or fractions, in Appendix
A.
of the
input
data
to
'
need to be calculated
are:
refrigerant
rate
mass flow
enthalpy
at the
suction
port
enthalpy
at the
discharge
enthalpy
at the
shell
port
outlet
cl r
work done on the
refrigerant
ii
work
done on the
shaft
work
input
ii
A)
S
cm
they
compressor
rate of heat loss due to cooling compressor and motor
4 cooling and if
to the
are not
specified
explicitly,
the:
6 can
rate
of compressor
'hilo
rate of heat transfer from gas to the suction gas
These ten unknowns are calculated in Eqs.
from
of
the
shell
three
heat
loss,
and
the discharge
energy bahace
equations
4.10 - 4.12,
(. (hsuction
port
ir (hdischarge
-h
inlet)
port
-h
'= 'hilo
suction
+ 'cooling
port )=c$
-'can
3
'
(4'10)
(4.11)
29
and
'r
(h outlet
- hdischarge
and the seven through
port)
equations
defining
= -'hilo
(4.12)
'
given
in Eqs.
suction
port
4.8
and 4.9
and Eqs.
4.13
4.17.
fir = I& (h isen,discharge
port
-h
s = irhmech
Bcm=
nisenl
(4.13)
3
(4.14)
,
's'Tlmotor
=(1---l
il cooling
)/(3413
(4.15)
'
mech'motor"cm
(4.16)
'
and Ii
r
The combination compressor
i
cm
which
port
S oper
D'u suction
- 4.16
yields
Eq. 4.18
vol,suction
of Eqs.
4.10
(4.17)
port for
the
overall
balance,
(4.18)
= ~~r(houtlet - hinlet) + EIcan1/34’3 is
map-based c
energy
=n
refrigerant dependent
identical model.
to Eq. 4.7, Equations
properties), on hsuction
the
corresponding
4.8 - 4.15 coupled
port-
set
form
a non-linear
of equations
The enthalpy
expression
at the
where suction
for
the
(through ir
the
and mr are port
is
I
30
of W r’ calculational
The iterative is
equations
shown in Fig. is
h
suction port is calculated
h
from
suction
values
4.3.
the
= h
port
conditions
using
Eq. 4.20
(which
form
and the
discharge
where
port
al 1 of
to the
shell
the
outlet
4.20
compression
refrigerant
gas.
actual
modeled the
the
process
result
This
is,
process
compressor
body to the
Davis
and Scott
heat [12],
efficiency
term
frictional
sources.
Once h conditions at
port
by Eq. 4.19. (4.19)
are computed
pressure discharge
using
at the port
shell
is
computed
and 4.13),
isen,discharge
loss
discharge Dart at the discharge the discharge
is
through
in Eq. 4.16
assumption
that
in additional is
of polytropic
extra
hsuction
- 'can)"r
refrigerant
of course,
as sequence This
calculations,
given
port
at the
of
-h
port n.men
suction
port(4.20)
pressure).
contains
compression
shell.
pressure
port
+
system
is the enthalpy after an isentropic compression port conditions to discharge port pressure (assumed
isen,discharge the suction port Equation
from
suction
h
equal
the
=
h
suction
Eq. 4.11
this
Thereafter
+ 'cooling
enthalpy
h
in
and tcan.
the
of Eq. 4.10,
at the
combines
to begin
to hinlet.
+ (0 hilo
inlet
suction port the refrigerant
Next,
In order
rearranged
of h
inlet.
from
, thilo,
scheme used to solve
assumed to be equal
The refrigerant
h
tcooling
can be dependent
which
the
1 osses,
on mr and on theeheat
dependent
not
(and
processes).
There
in a manner might
has been calculated,
the
port
are computed is equal
to that
for
the since
is more properly is
heat the
similar
increase
which
port
that
gain correct
gas inside
an artifical
inefficiencies
enthalpy
nonadiabatic
handled
the
theoretically
refrigerant
above
all
in
transfer
compressor to that
the mechanical
be expected remaining
assuming at the
of
the shell
due to refrigerant
refrigerant outlet.
. b.
31 ORNL- DWG
v
a
CALCULATE THE REFRIGERANT CONDITIONS AT THE SUCTION PORT
.
81- 4743
4 ,
I c
COMPUTE THE ENTHALPY AT THE DISCHARGE PORT AFTER ISENTROPIC COMPRESSION
J
1
+ COMPUTE
I
CALCULATED
SPECIFIED
c
+ CALCULATE THE COMPRESSOR MOTOR SPEED c
4
COMPUTE REFRIGERANT MASS Fu)W RATE AND WORK DONE ON THE < REFRIGERANT r
\
4
COMPUTE MASS FLOW RATE OF THE REFRIGERANT 4
j
/
L
,
4
\
0
\
4.3
/
I
1,
I
&r Fig.
,
CALCULATE THE FULL- LOAD MOTOR SHAFT POWER
DETERMINE THE FRACTION OF FULLLOAD MOTOR POWER
CALCULATE THE COMPRESSOR MOTOR EFFICIENCY
\
DETERMINE WORK DONE ON THE REFRIGERANT
. CONVERGENCE
Iteration on Internal Energy Balance Efficiency-Based Compressor Model
for
the
Loss and
32
The volumetric from
efficiency
Eq. 4.21
given
nvol,suction
based on suction
by Davis
and Scott
=
port,actual
port
conditions
is computed
[12].
'vol,suction
port,theoretical
-b
(4.21)
where
11 vol,suction
=
port,theoretical
1-c
-
refrigerant
pressure
at discharge
-
refrigerant
pressure
at suction
Y
-
specific evaluated
heat ratio at suction
C
-
actual
P discharge P suction
and "a"
port
and "b"
compressors Data).
valid
for
clearance
Davis
clearance
empirical
capacities
port
derived
refrigerant
ratio, from
test
data
but of the same series
[12]
similar
port
for the specific port conditions
volume
constants
and Scott
"basically
note
units
that
which
their
differ
(defined
in
should
be
equation
only
for
in stroke
or
volume."
Once the is
power,
are
of different
Block
point
port
is
reached fl,
of compreisor
operating the
refrigerant,
are
then
on the efficiency.
motor
4.3.
If
for
compressor
the
, and the from is motkr
Eqs. fl
is
a decision
compressor-motor as input
speed and efficiency
S oper The refrigerant i
has been evaluated,
been specified
i.e.,
computed value
efficiency
in Fig.
has not
values,
respectively.
computed
volumetric
data,
full the
are taken
input
load
shaft
values
to be the actua2
= S and n input motpr = 'motor max' mass flow rate, m , the work input required
4.17, the
full
4.13, rated
to achieve
load
iotor
and 4.14, (full)
load
power,
is
to fl
respectively.' that
The
must be imposed
the assumed maximum rated
motor
33 Alternatively,
if
is
speed, Soper'
and efficie;cy,
of fractional
rated
fl
of ffl,
motor
specified
as input
nmot , are
(full)
load,
f
The operating
is
fl
speed,
f,,,
computed
functions
in Eq. 4.23.
(4.23)
's,fl)
is assumed to be a linear
function
i.e.,
Soper = (l-0
+ afl ffl)
is the synchronous motor speed (e.g. S sync specified in the input data, i.e., S sync = S. The value of afl is determined by solvi?::. conditions
full-load
(f,,
motor
= 1.0)
speeds
using
known values
afl =
(i.e.,
'fl
(4.24)
Ssync
where
load
the motor
as built-in
as defined
= 'r'(?nech Soper,
data,
3600 or 1800 rpm) as . 4.24
under
of synchronous
fulland
- 'sync). s
Since 4.13), not
and s%~
directly
used with results
S
is
a function
G
is dependent
solva;rble Eqs.
4.13
for
of i
S
in Eq. 4.25, S oper
onrS oper .
and 4.1?':0'
(tE%Ggh
Eqs.
(from
However, Eqs. solve for S oper
4.24,
4.23,
Eq. 4.17), 4.23
Eq. 4.24
and 4.24
algebraically.
=
and is
can be This
(4.25) G&-
afl
($)(g5J
where
60~1 vol,suction
discharge
port
port
D(l728v
-h
suction
suction
port
port >
3
(4.27)
34
and
H f
fl
=
T/(3413
(4.28)
rlmech ws,fJ
iJ r
and the
constants
Once 'oper are calculated After
has been evaluated from ffl
calculated
1728 and 3413 are used to maintain Eqs.
from
the,fit
=
The coefficients a fit
slot
size
Eq. 4.25,
4.27,
the compressor efficiency
nmotor,max
5 i=O
published
i
1;r'
C).
efficiency
motor
'Imotor,i
and ffl
r'
efficiency
given
(4.29)
(ffl.'i
are assigned
[13]
can be
by Eq. 4.29.
values
These coefficients
by Johnson
and maximum motor
for
to f,,,
i in Eq. 4.29
(see Aipendix
to data
values
units.
and 4.28.
of motor
nmotor
Data subroutine
from
4.26,
has been computed,
rlmotor
Block
from
consistent
for
the
were obtained
a motor
of 86% at a load
in
with
"medium"
of 6.8
N-m (80
oz-ft). The next
step
and efficiency
in the
calculations,
have been computed
value
has been calculated
input
from
for
i r,
value
of i
(optionally)
ihilo
and can
A new refrigerant
next
using the is
compressor two methods
is
the
to calculate
motor
speed
and a
compressor
power
Eq. 4.30.
The computed
until
the
by one of these
GJ cm = fir/h
culated
after
Eq. 4.19
change
smaller
in the than
is
then
mech'motor
used to evaluate
using
Eqs.
enthalpy
at the
and the
sequence
value
a specified
of h
4.16,
4.8,
compressor
suction tolerance.
(4.30)
>
4
suction
shown in Fig. port
and
cooling and 4.9. 4.3
port is
is
repeated
from one iteration The refrigerant
calto the
enthalpy
at
-
35
compressor
shell
rearranged
form
exit,
outlet
outlet
= hdischarge
the
associated
and the
refrigerant
been computed, vergence
the
state
Once all
of these
iterative
efficiency-based
model
to the
values
in Section
defined
Eq. 4.31
(a
by Shaffer
in Section loops
of overall
isentropic
The suction
and Lee [14]
Ah
=
Ah
compressor
at condenser outer
shell
inlet
pressure
have
drop
con-
4.2.
have converged,
a suction
n sh
(4.31)
at the
to the
calculates
4.2.
d hilo hr
--
properties
proceed
earlier
is
using
and properties
calculations
described
described
computed
port
refrigerant
loop
addition
then
of Eq. 4.12), h
After
is
houtlet,
the
gas heating
loss
efficiency
and volumetric
gas heating
and in
efficiency
efficiency,
nsh,
as
isen,inlet
(4.32)
isen,suction
port
where
Ah
specific enthalpy change for an isentropic compression from she22 inZet conditions to shell outlet pressure
F
isen,inlet
and
Ah
isen,suction
The overall
isentropic
rl sh by Eq. 4.33
=
port
specific enthalpy change for an isentropic compression from suction port conditions to shell outlet pressure
compression
efficiency,
no,
is
related
to
[7,14-j. 11
cm
=
rl
sh 'isen
'mech
'motor
(4.33)
36
The terms five
on each side
efficiency
terms
magnitude input
values
in the
lossthan
note, for
n,,,h,
of the heating
with
be checked suction
four
between
for
by the
gas enthalpy
Eq. 4.16
this
dividing
through
model.
the
caution
the
loss
check
relative
This
the compressor
+ 'cooling
in suction
gas cooling
inlet
and suction
in termination
and the
requirement
port.
of the
Such a situation
or zero,
- 'can
(oracan) certain
Data Set."
be positive
selecting
is because
shell
and results
in and iC,
can result
use of Eq. 4.10 change
for
as an independent
(orahilo),
"Inconsistent
a message
i hilo
or,
must exercise
situation
output
All
inefficiencies.
parameters
dhilo Using
to identify
nmOtoTP ihilo
is an unrealistic
program
and also
user
independently.
in the
serve
compressor the
computed
are printed
and efficiency-based
combinations rather
various
thus
These values
calculations
of the
As a final
This
model.
compressor
are
in Eq. 4.33
and efficiency-based on the
of Eq. 4.33
that
can the
i.e.,
(4.34)
2 '
gives
+ ('
- n,,,hn,oto,)
w,
by i cm, an alternative
- km
(4.35)
2 '
expression
is
(4.36)
37
5.
FLOW CONTROL DEVICES 5.1
Introduction
The ORNL Heat Pump Design fixed
level
expansion high
of condenser device
which
employ
valves devices held
of the
each of these
to compute
the
design
TXV, and orifice
that
involves
using
to obtain
rate
from
the
described
three
would
the
between
compressor
One of the
at the empirical for
requires
inlet
pressure through
an equivalent
or more, the
subcooling
(this
from
BLOCK
and condenser
models
Capillary
corresponding in an iterative
refrigerant expansion
mass flow device
as
Tube Model
capillary
refrigerant tube
given
capillary
and subcooling the capillary
tube,
the
the
is
capillary
is specified, the calculated
three
in such a way as
characteristics
through
expansion
subcooling
called
in the Heat Pump Model
to curves
a standardized
is
the
subroutines
any of these
condenser
a particular
between
thermostatic
the specified
and that
of the capillary fits
flow
a
3.1.
options
of one,
model
produce
compressor
5.2
tube
for
device
a balance
in Section
operation
that
submodels
parameters
control
is used with
loop
event
for
contains
so that
some assumed operating
When a flow
DATA). model
In the
to specify
refrigerant tubes,
orifices
user
parameters
The program
capillary
and short-tube
can be modeled.
fixed,
the
system. for
the
or design
to control
correlations
(TXV),
allows
subcooling
in order
and low sides
Model
for tube
tubes
is
to simulate
in parallel.
pressure
and degree
or tubes.
The model
in the ASHRAE Equipment tube
flow
rate
R-12 and R-22. is given
as a function The actual
the
The capillary of subcooling consists
of
Handbook
[15]
of inlet mass flow
rate
by
in = +N r cap "k
(5.1)
38 where Ii;
-
standardized
N cap
-
number of capillary
@
-
flow
The capillary Equipment inside
flow
factor,
Handbook
[15]
for
tubes
each particular
mass flow
rate
from
the ASHRAE
combination is
of N
also
mi isciyven
P, and two functions
pressure,
in parallel,
4, must be obtained
The value
diameter.
The standardized
rate,
factor.
tube
tube
inlet
mass flow
of length
required
by Eq. 5.2
of the degree
and
as input. in terms
of
of subcooling,
m
0
and k.
(>
P k = m. 1500
ii;
mo, and exponent,
The coefficient, subcooling, within
AT, as given
1.5% of the
(5.2)
k, have been fit
by Eq. 5.3 and 5.4,
curves
across
the
range
as functions
such that of Fig.
of
mC agrees
40 in Ref.
to
15
3.56 m
0
= 356 + 0.641
k = 0.4035
In the level
situation
of condenser the
capillary
compute
the
flow
pressure fit
tube
factor,
tube
additional
curves
be added by the
user
critical
conditions.
in
Design lieu
Model
is used with
of a specific
would
give
pressure. capillary
control
the desired
subcooling.
assumes that However,
in the ASHRAE Equipment to predict
flow
a specified
are used in such a way as to
used in the program critical
(5.4)
e-0*04AT
correlations 0, which
model
is at or below
+ 0.4175
the
subcooling
device,
The capillary
where
(5.3)
tube
Handbook
the downstream
Smith
[16]
[15]
which
performance
under
has could non-
39
The current capillary the
tube
final the
to zero will
not
incomplete the
For all
flow
cases
zero
is
0 is
5.3
state
set
Expansion
l
specific empirical nozzle and tubes,
l
additional empirical equations tube liquid line temperatures, loadings.
used by Smith
ii
where
lil
r,TXV
in
r,TXV
=c TXIJtAToper
on the
keeping
with
the
rate
AT
- actual
operating
AP TXV
entrance. value
Model
one size
expansion of distributor
to correct for nonstandard lengths, and nozzle and tube TXV is a form
- available
ATstatic)
of the
equation
that flow
liquid
line
density,
of the
TXV will area
p,,
standard
in Eq. 5.5 are defined
has orifice
as follows:
pass, coefficient,
superheat,
superheat (superheat open), and pressure
(5.5)
(p??AplXV)1’2
genera 1 form
parameters
- mass flow
static barely
-
refrigerant
orifice
static-
for
of a cross-charged
- general
AT
required
thermostatic
correlations and
C TXV oper
the
Valve
model of a cross-charged
The remaining
equation.
tube
[16]:
a dependence
been added,
the
,'
a general valve,
model
that
capillary
contains:
0
The general
liquid
to zero.
Thermostatic
The TXV subroutine
in the
as input,
that
subcooling.
to make sure
at the
specified
assuming
of condenser
"flashing
at the
condensation,
continue
degrees
be checked
in a two-phase
factor
condensation
of incomplete
concerning
should
condensation
incomplete
and calculations
have
results
is
handle
are printed
final
refrigerant for
set
messages
line," If
is
conditions
error
cannot
entrance.
subcooling
1 the If
model
drop
at which
from TXVinlet
the
valve
to TXVexit.
is just
40
The assumption
AT
oper' that is the
p, and AP,,,,
independent
of evaporator
efficiency
such valves
tend
heat
the
effective
actual
*Tmax eff,
= ATstatic
be used to compute
AT AT
defined
static'
b fat,
the
The given factor
valve
in the
mass flow
greater
(ATrated
mass flow the
rating
bleed
over
a
the maximum
(5.6)
- ATstatic)
input
before
rate. data
pressure
drop
Eq. 5.5 can
The value
of CT,,
[17]: (tons),
of the
valve
(at
75% of
and
factor.
rating
TXV rating Block
than
of the TXV valve
earlier,
capacity
and standard
(specified
used in above since
superheat
have to be calculated
from
capacity
is
the rated operating superheat the maximum valve opening),
rated'
described
of AT in Eq. 5.5. oper C,,, and the available
a refrigerant
variables
the
6 rated'
type
opening
where
+ 1*33
AP,,,,
the TXV valve, on four
valve
Most TXV valves constant
superheat
is used in place ATmax eff The flow area coefficient,
depends
a given
a certain
pressure.
a relatively
operating
superheat,
across
provide
that
conditions.
ATmax , eff
then
will
implies
pumps are of the
to maintain
of operating
If
TXV valve
at a given
higher
range
of a cross-charged
is used in conjunction conditions
Data routine)
for
with
a specific
to calculate
the
the
bleed
refrigerant rated
refrigerant
rate:
f2,OOO 0 r-i r,rated
=
h
rated
evap out,rated
bfac -h
liquid
(5.7)
line,rated
where the
constant
Btu/hr.
12,000
is
used to convert
from
tons
of capacity
to
41 h
is the evaporator outlet refrigerant enthalpy rated evaporator pressure and superheat,
evap out,rated
at
and h
liquid
is
line,rated
CTXV is calculated
'TX,
the refrigerant line temperature.
from
= --)(hT
Eq. 5.8
enthalpy
using
at rated
liquid
mr rated: ,
Iir,rated rated
,p
static
----TPr,rated
rated
at rated
liquid
(5.8)
)V
where =
Pr,rated
refrigerant
density
line
temperature
and =
AP rated
rated
The available remaining
pressure
unknown
the TXV inlet
pressure drop
nozzle
The pressures pump model tube for
pressure
drops
one set
additiona
empirical
data
A,11 the
equations
pressure
are computed
da~~'~~~pthe
the drops
the
ptfessures
at
in the
to correct
and distributor
non-standard
fits
heat
and
R-22 fZow and
for
loading from
the
empirica7Y fits
sizes for
in
The nozzle
specific
and tube
were obtained
(5.9)
elsewhere
TXV routine. using
nozzle
equations lengths,
using
is
for
liquid
R-22,
to published
line
22, and performance
[18]. Distributor
the
tube
AP,,,,
- APnoz - APtube
in,evap
are calculated
of distributor
temperature,
502.
input
computed
and the
Pin TXV and P. and ari
device.
i.e.,
-P
aP P TXV = in,TXV
expansion
the TXV valve,
is
inlet
and tubes,
across
across
This
in Eq. 5.5.
and evaporator
distributor
drop
refrigerant
nozzles flow
and tubes
are
in each evaporator
often circuit
used with
TXV's
and to assure
to equalize proper
42 TXV operation.
At standard
have a pressure
drop
(15 psi)
for
R-121
69 kPa (10 psi) the
user
nozzle
culation
of the
in the
actual
capacity
tube do not
=
i no2
r;lr(h
",
is
-h
evap out,rated
0tube
~,,,,
for
components. is
the
cal-
and of each tube
the mass flow
liquid
across
the
cj,,
rated
,
nozzles
calculated
and tubes,
acor
(5.10)
(5.11)
= qnoz'Ncircuits
rate
= (3)12000
>
line,actual
by the
compressor
is the number of circuits in the evaporator. N circuits Next, the capacities are computed which yield the drops
allows
computations
have these
nozzle,
kPa
which
of nPnoz and nPtube
of the
will
each have
has been provided
which
nozzle
i.e.,
tube'
where
will
and distributor evaluation
distributor
R-22 and R502 [103
tubes
option
to systems
step
the
for
distributor
An input
applications
The first tj
and the
the
conditions
172 kPa (25 psi)
drop.
to omit
specific
of
rating
model
rated
and
pressure
i.e.,
10
(T
- 40)/201.0
(5.12)
sat,evap
and
4tube,rated
= (1.1)12000
aCOr .sor
10
(T
sat*evaP
- 40)/177.64
(5.13)
where the
constants
specifications acor
=
3 and 1.1 are tonnage
ratings
from
the
equipment
[18], correction factor than lOO"F, i.e.,
for
liquid
line
temperature
other
43
= lo(‘OO a
car I
- Tin
T
,)/155.18
T
(100 -Tin:,,,),140.19 10
in,TXV in,TXV
'.'O"
(5 ,4J
,100
.
and =
Bcar
correction i.e.,
factor
tube lengths
for
other
than
30 inches,
I3car = (30/&be)1’3 Equations
5.12
and 5.13,
a Sporlan
[ll]
nozzle
respectively,
No. 2.5
represent
and l/4
other nozzZe or tube, or refrigerant, replaced (the correction terms acor do not
the
5.11
and the
tube
loadings
drops,
inch
the
OD copper
R-22 capacity tubes.
of
FOP any
Eqs. 5.22 and 5.13 must be and BcOr, Eq. 5.14
and 5.15,
need to be replaced).
Given
0 tube"
(5.15)
actual
rated
nozzle
capacities
Lnoz and Ltube,
and tube from
capacities
Eqs.
defined
5.12
from
Eqs.
and 5.13,
the
as ~,,,/~,,,
respectively, and the tube,rated APnoz and aPtube, can be found.
nozzle
rated and'tube
5.10 nozzle
and and
and pressure
For R-22 or R-502
I
25.0 (L
AP = no2
29.4
)1.838
AP tube
.
Il.2 no2 > 1.2 L no2
(5.16)
5 1.1 noz > 1.1 L noz
(5.18)
L
noZ 0.9547 (Lnoz)
= 10.0(Ltube)'.8'22
For R-22
AP = noz
I
15.0(Lno2)'*8'7 )I.265 15.8l(L noz
L
44
A!'
=
tube
I
Equation
5.9
Eq. 5.5
lo.o(Ltube)'.772
L tube
'
12.265(Ltube)1'3377
L tube
>lS6
is
then
to evaluated
subroutine
from
the
which,
(5.19)
in turn,
is
used in
mr TXV.
In a case where the
AP,
used to obtain
l-6
a +ixed
calculates
level
the
of condenser
required
rated
subcooling
capacity
is
of the
specified, TXV valve
equations ii C TXV
rir,rated
=
6
AT
oper Or by ATrated tube pressure cooling
is
is
(AT oper
sta:ic)(","PTXV'1'2
CTXV(ATrated
replaced
- ATstatic)(Pr,ratedAPrated'1'2
(5*21)
out,rated - hliquid 12000 bfac
(5 . 22)
by ATmax eff
if AToper 5 aT,t,ti~* drop computations is model
is
not
In such cases,
calculations
continue
as described
incomplete
condensation
values
is
The short-tube Heat Pump Design
ratios
also
AToper > ATmax,eff
if
to omit
available
the
nozzle
when condenser
and sub-
data
incomplete
the
subcooling
in the
capillary
the
program
specified,
five
from 7.5
Short
orifice Model
for
diameters from
to handle
is
set
tube will
condensation to zero
If
case.
set
the
and
TXV
to zero. 5.4
with
Eq. 5.6
The option
intended
at the TXV entrance.
obtained
from
line,rated)
specified.
The current
rating
(5.20)
- AT
il r,rated(hevap
=
rated
where
=
subroutine
which
uses correlations 0.127
1.067
to 11.9)
Tube Orifice
m (0.5
to 1.694 using
R-22.
Model is
included
developed in.)
long
mm (0.0420
in
by Mei [19].
Carrier
ORNL He
Accurators
to 0.0667
He observed
the
pressure
inches) drops
(L/D
45 between
620 and 1515 kPa (90 and 220 psi)
levels
of subcooling
observed
from
refrigerant
470 lbm/h).
The standard
was solved with
both
level
for
r
orifice
orifice
of subcooling
rates
ranged
at the
across
0.63683 - 0.019337
x
=
C which short
tube
Equation
(P
-P
2Prgc(Pin
- 0.00325AT
4 i
to correlate
difference
occuring
the at
data
orifice,
AP, and the
+ 0.006aT I
,
AT < 40
1 l/2
2Prgc(Pin
- 'in,sat)
to within
low subcooling
7.5% (with
levels), inlet
the maximum
where
P in' to the orifice,
pressure
(5.24)
AT > 40 ,
at the
P in,evap' the inlet
inlet
to the
respectively.
If
the
rather
than
solved
for
As in
correlated
5.24,
q/2
- 'in
(5.23)
orifice,
in,evap
I[
are the pressures at the and P in,sat to the evaporator, and the saturation
(150 to
,
for
was found
The
to 213 kg/h
for
l-ID2
inlet.
was then
AT.
in
and
in Eq. 5.23,
C
the
inlet,
68.0
c(2g PAP)"~
=
coefficient
drop
restrictor
at the
from
relation
$/Aorifice
the pressure
IID2 4
r;l
the
=
the
0 to 28 Co (0 to 50 F")
mass flow
G r
across
user
has specified
selecting the
orifice
the case of the
specified,
the
a fixed
a particular diameter other
corresponding
flow that
control
would
two models, orifice
of condenser
level
device
provide
Eq. 5.24
that
if
incomplete
diameter
is
set
subcooling is
subcooling. condensation
to zero.
is
46
The user at the the
flow
three
results extended
should
always
control flow
should beyond
device
control not their
check
that
entrance
device
models.
be used because intended
is the
range.
the
computed
refrigerant
equal
to zero when using
If
is
flow
it
not,
control
quality any of
the calculated models
have been
47
CONDENSERAND EVAPORATORMODELS
6.
6.1
Introduction
The ORNL Heat Pump Design air-to-refrigerant
and evaporators
performance
a modified version of the effective surface approach when there is dehumidification,
0
the Thorn correlation the Moody friction phase refrigerant
0
friction partially
vs.
principal
NtU correlations
for
the
differences
there
heat
exchangers
with
unmixed
and the
flow
superheated
differences
on both
the
, and subcooled
6.1
is
a general
organization
and iterative
6.2
is
diagram
a similar
The air-side density
air
calculated the
Universal
the
air
temperature.
inlet
several being
equivalent specified
estimated
divided
by the
flow
for
the
refrigerant
Since
the
for the
data),
mass flow
number of circuits
the
The for
the
the
condenser.
chart,
outlining
model.
Figure
calculations. in air
heat
exchanger
is calculated
the
input
data and the
using
atmospheric
(given
in BLOCK DATA),
exchangers
are modeled
circuits the rate
[3].
condenser
gas equation
refrigerant input
the
or flow
specified ideal
they circuits
and for
for
each heat
gas constant
parallel with
rate
are
computations
evaporator
diagram,
for
from the
pressure,
the
into
the evaporator rate
out where
sides
dry,
algorithm.
refrigerant
regions
loops
for
mass flow
volumetric
block
for
identical,
and evaporator
parallel
for
drop
are not
and refrigerant
regions
the
the
air
drops and the single-
have been used assume that
are separated
and two-phase
Figure
which
a
temperature
pressure
are pointed
of equivalent,
two-phase
for
the dehumidification
the condenser
calculations
superheated,
transfer
and evaporator
around
methods
consist
refrigerant-side
from
condenser
center
The calculational
occur.
heat
for the air-side wetted coils.
are many similarities, together
for
for two-phase refrigerant pressure factor chart plus momentum terms for pressure drop, and
factor equations wetted, or fully
the models
of
by using:
0
discussed
the
the
effectiveness dry coil,
Since
air
calculates
0
Although the
condensers
Model
(the
air-side from
to obtain
the
values
actual
mass flow compressor
for
and
as
number rate model
each circuit.
and are
48 ORNL-DWQ
Sb4750R
.
CALCULATE AIR-AND REFRIGERANTSIOE MASS FLOW k4TE PER CIRCUIT AND QEOMttRlC PARAMETERS DETERMINE SUBCOOLED AND SUPERHEAlEO REFRlGERANt PROPERTIES AND REFRIQERANTSIDE HEAT TRANSFER COEFFICIENTS
*
4 DETEf3MI~E REFRIGERANT PROPERTIES At@ HEAT TRANSFER ,3%F&IENT FOR TWO-PHASE 4 ’
k AOJUST THE TWO-PHASE HEAT TRANSER COEFFlCkNT FOR THE REQION WHERE THE BULK TEMPERATURE EXCEEDS THE SATURATION TEMPERATURE
COMPUTE AIR-SIDE PROPERTIES AND PRESSURE DROPS, OPTIONALLY ADD FAN AND/OR COMPRESSOR HEAT LOSS TO INCOMING AIR .AND COMPUTE AIRSIDE HEAT TRANSFER COEFFICIENT
RETURN CODE
t 1'
v
1
COMPUTE REFRIGERANT- SIDE PRESSURE DROP FOR SUPERHEATED ,TWO-PHASE, AND SUGCOOLEO REGION I
QPTIONALLY ADD FAN AND/OR COMPRESSOR HEAT LOSS TO OUTLET AIR AN0 COMPUTE UA VALUES
.
Fig.
6.1
General
Structure
of the Condenser
Model
ORNL-
DWG
864751
. f INITIALIZE ITERATION I PARAMETERS 1 c CALCULATE AIR- AND REFRIGERANTSIDE MASS FLOW RATE PER CIRCUIT AND GEOMETRIC PARAMETERS
.
I
f
DETERMINE SUPERHEATED REFRIGERANT PROPERTIES AND SUPERHEATED VAPOR HEAT TRANSFER COEFFICIENT
*
t DETERMINE ENTERING QUALITY AND AVAILABLE HEAT TRANSFER IN THE TWO-PHASE REGION * DETERMINE REFRIGERANT PROPERTIES AND HEAT TRANSFER COEFFICIENT FOR THE TWO-PHASE REGION L CALCULATE AIR-SIDE PROPERTIES AND PRESSURE DROPS, OPTIONALLY ADD FAN AND/OR COMPRESSOR HEAT LOSS TO INCOMING AIR,AND COMPUTE AIR-SIDE HEAT TRANSFER COEFFICIENT
COMPUTE REFRIGERANT-SIDE PRESSURE DROP FOR SINGLE AND TWO-PHASE REGION A I
+
COMPUTE NEW AVERAGE EVAPORATING TEMPERATURE
OPTIONALLY ADD FAN AND/OR COMPRESSOR HEAT LOSS TO OUTLET AIR AND COMPUTE U A VALUES
Fig.
6.2
General
Structure
of
the
Evaporator
Model
50
The average
densities
two-phase
region
calculated
from
peratures
in the
the
of each coil the
current
heat
single-phase
for
conductivity,
the
6.2
computes
the
Single-Phase
region
in the
condenser
by Eq. 6.1.for throughout
gas flow
this
report
average
i.e.,
of ah at
p
an abrupt
heat
=
for
the
condenser
thermal according
for
the
by Kartsounes
to a
vapor
region,
and Erth
[4]
cp,vm in Heat Exchanger
coefficient
is calculated for
viscosity,
heat
tem-
properties for
the are
saturation
are calculated
Coefficients
transfer
from
the
in
of vaporization
and superheated
developed
Heat Transfer heat
of the
The specific
by a routine
The refrigerant
heat
heat,
[5]. value
and vapor
latent
(subcooled specific
local
liquid
The thermophysical
evaporator),
by Flower
Cp V, is determined which
estimates
exchangers.
and liquid
written
refrigerant
and the
refrigerant
and superheated routine
of the
for
the
Tubes
superheated
(as reported
in the
MIT Model
contraction
entrance
[ZO]
transfer
coefficients
and h for
[3])
(h is used enthalpy
variables). h
where
N,
=
GrD/p
=
c,
=
=
-2/3 *I+
N '2 Re
(6.1)
and,
I
for
NRe < 3500
0.01080
for for
3500 5 NRe < 6000 NRe 2 6000
-0.78992
for
NRe < 3500
1.03804
for
3500 5 NRe < 6000
-0.13750
for
NRe 2 6000
1.10647
c,
C G c 1 r p,v
3.5194
x 10 -7
(6.2)
(6.3)
51 The constants Eqs.
6.2
Appendix
used to define
and 6.3
the
are assigned
step
in the
functions
C1 and C2 defined
BLOCK DATA subroutine
in
(see
C).
The heat
transfer
condenser
and the
using
Dittus-Boelter
the
coefficients
superheated
for
region
correlation
the
subcooled
in the
evaporator
for
fully
region
of the
are computed
developed
flow
[Zl]:
(6.4)
where
"C" is 0.3
when being
when the
Two-Phase
The two-phase the
evaporator
local
values
where
tube
wall.
to assuming coefficient,
being
cooling
Refrigerant Heat Transfer in Heat Exchanger Tubes
heat
transfer
are average over
the
the
of tube
change
Assuming
over
the
and 0.4
range
related
heat
flux)
is given
Z
ave
= J Z
the
practically
[3,22],
the
s
=
' dz
and
integrating The actual quality
in quality
[3] of dx by
refrigerant
and the
constant
(as opposed
average
X
0
condenser
refrigerant
by
lz i
of
between
h(x)dz h
the
region.
to a change
AT remains
Coefficients
by effectively
two-phase
in temperature that
for
obtained
of the
dz is
a constant have,
coefficients
values
length
is performed
a length
and AT is
is
heated. 6.3
integration
refrigerant
x Odx i ,
heat
transfer
which
reduced
The local calculated Travis
to
heat
transfer
as a function
coefficient,
of the
quality
h(x),
for
according
the
condenser
is
to a correlation
by
[22].
I
0 9 F(Xtt) kl
h(x)
*Pr *Re ’
F(Xtt)
DF2
< 1.0
= 4
(6.5) kL *Pr
o g *Re ’
F(Xtt)1*15
1.0 < F(Xtt)
DF2
< 15
where (6.6)
Gr D (1 -xJ N = Re
F(xtt)
= 0.15
Yl
(6.7)
,
+ 2.85Xtt-o'476
(6.8)
, and
>
N
0.707 Npr NReoe5 5.0 Npr + 5.0 F2
lntl.0
Re
+ Npr(0.09636
=
NReo'585
l.O)l --l.O)]
50 < N
Re
LnCl.0 5.0 Npr + 5.0 Ln[l.O
+'5.0
c: 50
NPr] + 2.5 1nC0.00313
< 1125
NRe0s8121 NRe0S812]
NRe > 1125 i
(6.9)
53
The integration
is
carried
Ax and has been changed
MIT Model [3]. to liquid adjusted where bulk
slightly
heat
from
condensation
given
temperature
in
steps
the
to the
heat
and Glicksman
manner
at
the
transfer
of extrapolating
end points
coefficient
is occurring
exceeds
of a constant
Hiller
by Rohsenow [23]
(desuperheating)
refrigerant
over
coefficients
two-phase
a correlation
that
pertain
transfer
The average using
equation
by a summation
The modifications
and vapor
integration.
out
for
is
the
region
at the wall
the saturation
of the
but
the
temperature.
The
>
(6.10)
is
T h
The bulk begins
=
tP
refrigerant at the
v,ds
h
wall
h
P9v
is calculated
v’T T v,ds
tp,in r-l
=
sat,in
fg
TV ds, where
temperature,
tube
l/4
-T
refrigerant
from Eqs. -T
6.11
condensation
and 6.12
where
a,in
(6.11)
and
(6.12)
T
tp,in phase
and Ta in are the region
And the
The average the
evaporator
is for
derived
[24]
and from
the general
region.
Equation of the
inlet
air
local
heat
for
from a correlation
expresses transfer
12 in
of Sthapak the
heat
coefficients
the
the
two-phase
developed
of refrigerant
observations
entering
two-
respectively.
coefficient
evaporation 6.13
temperature
temperature,
heat-transfer
Noerager
terms
refrigerant
[25]
transfer
region
of
by Chaddock
and
horizontal for
the
coefficient
as functions
tubes dry-out in of refrigerant
54 quality
from
and from
the
dry-out
inlet
to the
to exit
quality, 'do - s xi
1 h
dry-out h2(x)
the
for
Xtt
liquid
is
between
from
xi
to xdo, .
xdo and x . 0
..
o dx s + Xdo h2(x)
(6.13)
tP
two integrals
xxodx i
separately,
h&x) where
$(x)
X
dx hlTx) s
Examining
point,
as defined
= 3*OhQ(Xtt
in Eq. 6.6
refrigerant
obtained
hi(X)
hi(x)
given
Eq. 6.4.
'R l/3
by Eq. 6.14
1 )2/3
(6.14)
and hL is the
from
= 3.0$(;;-)
is
'v $1
heat
transfer
coefficient
Thus
0.0667(;
x x)O.6
V
and
J X
'do k i
Xdo (1 - x)06 xO.6
J X
dx
(6.15)
i
where
0 0 Pt
C
Equation yield
6.15
=
3*0ht
can be expanded,
l/3
-7
truncated,
lJ
V
0.0667
5
and integrated
term
by term
to
55
The second nonlinear phase pure
integral
in
interpolation
value
at the h
vapor,
of the dry-out
from
V'
the
numerator heat
point,
Eq. 6.4,
of Eq. 6.13
transfer
hl(xdo),
from
of the
The second Eq. 6.17,
evaporator
integral
from
Eq. 6.14
the
to value
twoof
i.e.,
x do and x0 are refrigerant
outlet
based on a
coefficient
h’ltxdo) where
is
qualities
or the
in Eq. 6.13
at the
two-phase
region
can be evaluated
- “VI dry-out
(6.17) point
(whichever in closed
and the
comes first). form
using
i.e.,
(h, (Xdo)Ch(Xdo)- hvlF
X
0
+j-=
/
X
2(xo
- 'do)
'do
A value
of 0.65
is
assumed for
6.4
Air-Side
The air-side McQuiston the
heat
[26,27],
correlation
h
given
Ch,fXdo)
- h,])
{ h&xdo)
b'l(xdo)
- hvl}
xdo in the model Heat Transfer
transfer
Yoshii
{ hl(Xdo)
[28],
coefficients and Senshu
[25].
Coefficients are based on the work of [29]
and are calculated
a
by
by Eq. 6.19.
1 - 1280 N, NRe-lm2 =
(6.18)
C G c N _ -2'3 0 a pa Pr
j 1 - 5120 NRe-ls2
1
(6.19)
where j = 0.0014 N Re=
+ 0.2618
(,
1 FJ-o*15
(y
)
-OB4 ,
GawT I lJ
and C
0
= 1.0, 1.45, or 1.75 depending on whether the fins are smooth, wavy, or louvered (see Appendix G for parameter definitions).
/
Equation
6.19
over
Reynolds
the
transfer
was obtained for
approximately
extensive of 3,000
constant
C .
a number of fin-and-tube
surfaces
are assumed to be referenced
increases smooth
for
fin
heat
values.
Equation
number of geometric spacing,
and the
calculated which
using
is
Heating,
a modified
version
based on correlations
evaporator dry
heat
which
is wetted
coefficient,
The fin and for
on the work a tube
the
of Schmidt
surrounded fin
arrangement).
This
tube
changes
spacings
[32]
shape
calculate
the
surface
a the
fin
properties
are
[5]
Society
of
[30]. portion
of the
is calculated
from
[31]. (6.20)
effectiveness
evaporator
as reported fin
segment
surrounding
that
occur
as the
fin
longitudinal
[33]
for
(the
in a staggered to properly
conbased
and Tree
of a sheet
each tube
has been generalized
for.the
are calculated
by McQuiston
work
tube
handle
and transverse
are varied. 6.5
Figures
the
for
($)“.lolha
of the
by a hexagonal
representative orientation
region
for
rows,
Engineers
from Myers
= 0.626
a,w
thus
from
by Flower
due to dehumidification
and overall
dry
to adjust
by the American for
area;
area
The air-side
coefficient
fin
to account
of a subroutine
h a, by Eq. 6.20
efficiency
surface
intended
Conditioning
transfer
available.
wavy and louvered
number of tube
published
increased
assumption
are not
terms
spacing.
and Air
h
denser
such as the tube
equation
and surface
includes
transverse
Refrigeration,
are
The heat
this
for
coefficient 6.19
effects
The air-side the
transfer
fin
to smooth-fin fins
fins
are assumed to be
however,
so calculated
wavy and louvered
in both
on smooth
Data to verify
transfer
Co values
fins smooth
geoietries,
coefficients
data
2 NRe 5 15,000.
wavy and louvered
The heat the
test
by the use of the
by the multiplicative for
from
number range
coefficients
predicted
56
6.3
and 6.4
performance
Heat Exchanger are block of the
Performance
diagrams condenser
for
the
subroutines
and evaporator,
which
respectively.
57 ORNL-DWG
El-4742R
CALCULATE THE REFRIGERANT TEMPERATURE AT THE STARTOF THE DESUPERHEATING REGION
f
1
COMPUTE THE EFFECTIVE DRIVING ENTHALPY DIFFERENCE FOR TWO-PHASE REGION
ITERATE TO FIND FRACTION OF COIL IN VAPOR REGION f” \\
f SET
ftp= 1 - f,
AND
CALCULATE
THE
EXIT QUALITY
I
>
COMPUTER FRACTION OF COIL IN TWO-PHASE REGION f
tp
NO
1
SETf
=f-(f
AND
C&PUTE?,q
u&
tf,, AND
FOR SUBCOOLED
REGION
1
1 ‘COMPUTE C, FROM ESTIMATED cp FOR THE SUBCOOLED REGION t
I
,
CALCULATE EFFECTIVENESS AND NTU FOR SUBCOOLED REGION AND A NEW EXIT STATE , r
\ CALCULATE ESTIMATE
A NEW OF cP
FOR THE REGION
SUECOOLED
CONVERGENCE
ON
/ , 2. \
Fig.
6.3
Computational
CALCULATE HEAT TRANSFER RATES AND EXIT AIR TEMPERATURE FOR EACH FLOW REGION
Sequence
in the
Condenser
Model
ORNL-DWG
CALCULATE
81-4744R
PROPERTIES
OF MOIST AIR AND SUkFACE EFFECTIVENESS AND THERMAL RESISTANCE OF DRY COIL
COMPUTE AIR TEMPERATURE WHEN DEHUMIDIFICATION BEGINS,
Toid
I
\ CALCULATE
NTU
AND
EFFECTIVENESS TWO-PHASE ASSUMING REMOVAL
FOR REGION
NO
MOISTURE
OCCURS
COMPUTE EXIT AIR TEMPERATURE FROM PHASE REGION NO MOISTURE
TWO-
ASSUMING REMOVED,
/ CALCULATE HEAT TRANSFER RATE AND THE FRACTION OF THE COIL IN TWO-PHASE \
YES
.
FLOW,f+P
MOISTURE REMOVAL OCCURS SEE FIG.5.1
I
AND
FOR THE
ESTIMATE
COMPUTE
VAPOR
REGION
OF
YES \ CALCULATE CONDITIONS c
6.4
Computational
. ORIFD Diameter of the orifice GERANT LINES
Card 5
1
1.0
3 0.0544
BeTB
FORMAT(5F10.5) Inside diameter of liquid line from indoor coil to outdoor coil (in). XLEQLL Equivalent length of liquid line from indoor coil to outdoor coil (ft).
DLL
5.44;
0.1900 30.4
91 QDISLN Heat loss rate in the compressor discharge line (Btu/h) QSUCLN Rate of heat gain in the compressor suction line (Btu/h) from the indoor coil QLIQLN Heat loss rate in the liquid-line to the outdoor coil (Btu/h). c
Card 6
Card 7
FORMAT(4F10.5) DLRVOC Inside diameter of refrigerant line from the reversing valve to the outdoor coil (in). length of the refrigerant line from the XLRVOC Equivalent reversing valve to the outdoor coil (ft). DLRVIC Inside diameter of refrigerant line from the reversing valve to the indoor coil (in). XLRVIC Equivalent length of the refrigerant line from the reversing valve to the indoor coil (ft). FORMAT(4F10.5) Inside diameter of the refrigerant line from the compressor shell inlet to the reversing valve (in). length of the refrigerant line from the XLEQLP Equivalent compressor suction port to the reversing valve (ft). DDLRV Inside diameter of the refrigerant line from the compressor shell outlet to the reversing valve (in). XLEQHP Equivalent length of the refrigerant line from the compressor discharge port to the reversing valve (ft).
DSLRV
2000.0 300.0
200.0 0.68 6.00 .5556 30.40
0.680 2.00 0.556 2.00
S INITIAL, ESTIMATES m BEGIN SITERATION Card 8
FORMAT(3F10.2) TSATEO Estimate of the of the evaporator TSATCI Estimate of the of the condenser XMR Estimate of the
saturation (F). saturation (F). refrigerant
temperature
at the outlet
temperature
at the inlet
mass flow rate
(lbm/h).
30.0 115.0 400.0
COMPRESSOR Q&E,A Card 9 ICOMP
DISPL SYNC
QCAN
FORMAT(I10,7F10.0) switch to specify which compressor submodel is to be used, =I, for the efficiency and loss model, ~2, for the map-based compressor model. Total compressor piston displacement (cu in). when ICOMP = 1, synchronous compressor motor speed when FLMOT (on card 11) is specified (rpm); actual compressor motor speed when FLMOT (on card 11) is to be calculated (rpm). when ICOMP = 2, rated compressor motor speed (rpm); Heat rejection rate from the compressor can, used if CANFAC= 0.0 (Btu/h).
2
4.520 3600.0 3450.0 3450.0
0.0
92 CANFAC Switch to control the method of specifying QCAN ~0.0, to specify QCANexplicitly, C1.0, to set QCAN= CANFAC* compressor power, ~1.0, to set QCAN= 0.9*(1.0 - actual motor efficiency mechanical efficiency)*compressor power (not used with the map-based model). The variables is chosen. Efficiency Card 10
Card 11
0.350 *
on the next two cards depend on which compressor model
and Loss Model: FORMAT(8FlO.O) VR Compressor actual clearance volume ratio. EFFMMX Maximum compressor motor efficiency. ETAISN Isentropic compressor efficiency. ETAMEC compressor mechanical efficiency. FORMAT(I10,7F10.0) MTRCLC Switch to determine whether compressor full load motor power is specified as input data or if it is calculated, =0, compressor full load motor power is specified as input, =I, compressor full load motor power is calculated. FLMOT Compressor motor output at full load, i.e. (0.7457 kW/hp)x(brake hp at rated load), used only if MTRCLC=O(kW). rate from the compressor discharge line QHILQ Heat transfer to the inlet gas, used if HILOFC=O.O (Btu/h). HILOFC Switch to control method of specifying QHILO ~0.0, to specify QHILO explicitly, .Thickness of the fins in the indoor coil (in). Fin pitch for the indoor coil (fins/in). Thermal conductivity of the fins in the indoor coil (Btu/h-ft-F). Frontal area of the indoor coil (sq ft).
FORMAT(2FlO.O) Number of row of tubes in the indoor coil in the direction of the air flow. refrigerant circuits in the NSECT Number of parallel indoor coil.
NT
2.0
0.400 0.336 0.00636 14.00 128.00 3.1667
3.0 3.0
94 Card 14
FORMAT(3F10.5) HCONT Contact conductance between fins and tubes for the indoor coil (Btu/h-sq ft-F). ST The vertical spacing of tube passes (in). WT The spacing of tube rows in the direction of air flow (in).
Card 15 QA TAI RH RTB FANEF
Card 16 FINTY
DEA DER DELTA FP XKF AAF Card 17 NT NSECT Card 18
30000.0
1.00 0.875
FORMAT(5F10.4) Air flow rate for the indoor coil (cfm). Inlet air temperature for the indoor coil (F). Relative humidity of the inlet air for the indoor coil. The number of return bends in the indoor coil. Combined fan and fan motor efficiency for the indoor coil.
FORMAT(7F10.6) Switch to specify the type of fins, ~1.0, for smooth fins, ~2.0, for wavy fins, ~3.0, for louvered fins. Outside diameter of tubes in the outdoor coil (in). Inside diameter of the tubes in the outdoor coil (in). The thickness of the fins in the outdoor coil (in). The fin pitch for the outdoor coil (fins/in). The thermal conductivity of the fins in the outdoor coil (Btu/h-ft-F). The frontal area of the outdoor coil (sq ft). FORMAT(2FlO.O) The number of rows of tubes in the outdoor direction of air flow. The number of parallel refrigerant circuits outdoor coil.
coil
in the
in the
FORMAT(3F10.5) HCONT The contact conductance between the fins and tubes for the outdoor coil (Btu/h-sq ft-F). The vertical spacing of tube passes for the outdoor ST coil (in). The spacing of tube rows in the direction of air flow WT for the outdoor coil (in).
Card 19 QA TAI RH RTB FANEF
FORMAT(5F10.5) The air flow rate for the outdoor coil (cfm). The inlet air temperature for the outdoor coil (F). The relative humidity of the air entering the outdoor coil. The number of return bends in the outdoor coil. The combined fan and fan motor efficiency for the outdoor coil, not used if MFANFT=l on Card 3.
1200.0 70.0 0.50 72.00 0.200
2.0
0.400 0.336 0.00636 14.00 128.00 5.040
3.0 4.0
30000.0 1 .oo 0.875
2300.0 47.0 0.70 64.00 0.160
95 Appendix B: Sample Input and Output Data SAMPLEDATA SET FOR ORNLHEAT PUMPDESIGNMODEL-- @ 47.OF AMBIENT 15000.0
8.00 2
* o.,,o: 0.6800 0.6800 30.00
f .5
2
10.00 1
30.40 45.0 6.00 2.00 130.00 4.520
1 2000.00 0.5550 0.5550 400.0 3450.00
2 300.00 30.40 2.00
1
0
200.00
0.350
-1.5093-04 4.0893-02-1.3383-04 5.86OE% 3.6383-04 9.7593-05 -2.6753-02 4.6333+00 4.703E-02 9.6403+00-1.8683-02 1.207E-04 3.z
l
4.523+00
0.4000 3.000
0.3360
0.00636
14.0
128.0
3.1667
30000.00 1200.0 3.;;:
1.00 70.0 0.4000 4.000
0.875 0.500 0.3360
72.00 0.00636
0.200 14.0
128.0
5.0400
30000.00 2300.0
1.0 47.00
0.875 0.700
64.00
0.160
****
2O.OE+OO
INPUT DATA l **** SAMPLEDATA SET FOR ORNLHEAT PUMPDESIGNMODEL-- @ 47.OF AMBIENT HEATINGMODEOF OPERATION(NCORHz2) THE HOUSELOADIS 15000. BTU/H DIAMETEROF 6 EQUIVALENTDUCTS- 8.000 IN COMPRESSOR CAN HEAT LOSS ADDEDTO AIR BEFORECROSSINGTHE OUTDOOR COIL. POWER TO THE INDOORFAN ADDEDTO AIR AFTERCROSSINGTHE COIL. POWER TO THE OUTDOOR FAN ADDEDTO AIR BEFORECROSSINGTHE COIL. CONDENSER SUBCOOLING IS HELDFIXED AT 45.00 F EVAPORATOR SUPERHEAT IS HELDFIXED AT 10.00 F
L
DESCRIPTIONOF CONNECTING TUBING: LIQUID LINE FROMINDOORTO OUTDOOR HEAT EXCHANGER ID 0.19000 IN EQUIVALENTLENGTH 30.40 FT FROMINDOORCOIL TO REVERSINGVALVE FROMOUTDOOR COIL TO REVERSINGVALVE ID ID 0.55500 IN 0.68000 IN EQUIVALENTLENGTH 30.40 FT EQUIVALENTLENGTH 6.00 FT FROMREVERSINGVALVE TO COMPRESSOR INLET FROMREVERSINGVALVE TO COMPRESSOR OUTLET ID ID 0.68000 IN 0.55500 IN EQUIVALENTLENGTH 2.00 FT EQUIVALENTLENGTH 2.00 FT HEAT LOSS IN DISCHARGELINE HEAT GAIN IN SUCTIONLINE HEAT LOSS IN LIQUID LINE
2000.0 BTU/H 300.0 BTU/H 200.0 BTU/H
ESTIMATEOF: 400.000 LBM/H REFRIGERANT MASSFLOWRATE 130.000 F SATURATIONTEMPERATURE INTO CONDENSER SATURATIONTEMPERATURE OUT OF EVAPORATOR 30.000 F COMPRESSOR CHARACTERISTICS: 4.520 CUBIC INCHES TOTALDISPLACEMENT SYNCHRONOUS MOTORSPEED3450.000 RPMS MAP-BASEDCOMPRESSOR INPUT: POWER CONSUMPTION. -1.509%04*CONDENSING TEMPERATURE**2+ 4.089E-02'CONDENSING TEMPERATURE + -1.338E-04'EVAPORATINGTEMPERATURE**2+ 5.860E-04*EVAPORATINGTEMPERATURE + 3.638E-04'CONDENSING TEMPERATURE*EVAPORATING TEMPERATURE,+9;i59E:05 MASSFLOWRATE-
-2.675E-02*CONDENSING TEMPERATURE**2+ 4.633E+OO*CONDENSINGTEMPERATURE + 4.703E-02*EVAPORATINGTEMPERATURE**2+ 9.640E+OO*EVAPORATING TEMPERATURE + -1.868E-02*CONDENSING TEMPERATURE*EVAPORATING TEMPERATURE + 1.2073-04
CORRECTION FACTORFOR SUCTIONGAS CORRECTION FACTORFOR VOLUMETRIC EFFICIENCY BASE SUPERHEAT FORCOMPRESSOR MAP BASE DISPLACEMENT FOR COMPRESSOR MAP
0.330 0.750 20.000 4.520
F CU IN
HEAT REJECTEDFROMCOMPRESSOR SHELL IS 0.35 TIMES THE COMPRESSOR POWER
INDOOR UNIT: CONDENSER 0.40000 OD OF TUBES IN HX ID OF TUBES IN HX 0.33600 FRONTAL AREA OF HX 3.167 3.00 NUMBEROF PARALLEL CIRCUITS 3.00 NUMBEROF TUBES IN DIRECTION OF AIR FLOU 72.00 NUMBEROF RETURN BKNDS AIR PLOWRATE 1200.00 70.000 INLET AIR TEMPERATURE WAVY FINS OUTDOORUNIT: EVAPORATOR 0.40000 OD OF TUBES IN HX 0.33600 ID OF TUBES IN HX FRONTAL AREA OF HX 5.040 4.00 NUMBEROF PARALLEL CIRCUITS NUMBEROF TUBES IN DIRECTION OF AIR FLOW 3.00 64.00 NUMBEROF RETURN BENDS 2300.00 AIR FLOWRATE 47.000 ISLET AIR TKMPERATURE WAVY FINS
IN IN SQ FT
CFM F
IN IN SQ
FT
CFM F
FIN THICKNESS FIN PITCH THERMAL CONDUCTIVITY OF FINS CONTACT CONDUCTANCE HORIZONTAL TUBE SPACING VERTICAL TUBE SPACING FAN BFFICIENCY RELATIVE HUMIDITY
FIN THICKNESS FIN PITCH THERMAL CONDUCTIVITY OF FINS CONTACT CONDUCTANCE HORIZONTAL TUBE SPACING VERTICAL TUBE SPACING FAN EFFICIENCY RELATIVE HUMIDITY
0.00636 14.00
IN
1.000 0.20000 0.50000
IN
0.00636
IN
FINS/IN 128.000 BTU/H-FT-F 30000.0 BTU/H-SQ FT-F 0.875 IK
14.00 FINS/IN 128.000 BTU/H-FT-F 30000.0 BTU/H-SQ FT-F 0.875 IN 1.000 IN 0.16000 0.70000
i
’ I
;
l
****
CALCULATED HEAT PUMP PERFORMANCE*****
COMPRESSOR OPERATING CONDITIONS: COMPRESSORPGWER MOTORSPEED
4.011 3450.000
REFRIGERANT MASS FLOW RATE COMPRESSOR SHELL HEAT LOSS SYSTEM SUMMARY COMPRESSOR SUCTION LINE INLET SHELL INLET SHELL OUTLET CONDENSERINLET OUTLET
EFFICIENCY VOLUMETRIC OVERALL
Km
RPH
LBWH BTU/H
413.828
4791.414
REFRIGERANT TEMPERATURE 39.470 F 43.687 206.498
0.6311 0.5053
POWERPER UNIT MASS FLOW POWERCORRECTIONFACTOR MASS FLOW RATE CORRECTIONFACTOR
183.889 F 81.217
F
126.558
126.196
EXPANSION DEVICE
79.607
F
120.765
F
EVAPORATORINLET OUTLET
34.916 F 39.470
34.916
F
TOTAL SEAT EXCHANGEREFFECTIVENESS
NTU HEAT EXCHANGEREFFECTIVENESS CR/CA FRACTION OF HEAT EXCHANGER HEAT TRANSFER RATE OUTLET AIR TEHPERATURE AIR SIDE: MASS FLOW RATE PRESSUREDROP HBAT TRAESFER COEFFICIENT UA VALUES: VAPOR REGION (BTU/H-F) REFRIGERANT SIDE AIR SIDE conB1NED
0.0 0.0 0.0
29.577
BTlJ/LBM
1:0089
SATURATION REFRIGERANT REFRIGERANT TEMPERATURE ENTHALPY QUALITY 108.958 BTU/LBM 1 .oooo 29.577 F 29.263 109.683 1.0000 126.899
PERFORMANCEOF EACH CIRCUIT IN THE CONDENSER INLET AIR TEMPERATURE 70.000 OUTLET AIR TEMPERATURE 100.147 HEAT LOSS FROM FAN 1257.1 AIR TEMPERATURECROSSINGCOIL 99.193
33;3;9";;
REFRIGERANT PRESSURE 69.054 PSlA
131.185
1.0000
68.660 299.682
126.352 BTU/LB!4
1 .oooo 0.0
298.411 297.060
PSIA
33.477
32.994 BTU/LBM
0.0
277.316
PSIA
32.994 BTU/LBM
0.1483 1.0000
'V$E~ .
PSlA
108.958
70.000 100.147
F
47.000 38.908
F
F F BTU/H F
W U
0.6949
SUPERHBATER REGION 0.0 1.0000 0.0 0.0 0.0 BTU/H 70.000
AIR TEMPERATURE
9
1798.1 0.5224
LBH/H IN H20
13.136
BTU/H-SQ FT-F
TWO-PHASE REGION 1.1351 0.6786 0.6453 10838.4 108.270
BTU/H F
suBcooLED RRGION 1.9987 0.8005 0.2816 0.3547 1973.0 82.676
BTU/H F
REFRIGERANT SIDE: MASS FLOW RATE 137.9 PRESSURE,DROP 1.351 HEAT TRANSFER CORFFICIENT VAPOR REGION 102.905 453.291 TWO PHASE REGION SUSCOOLBD REGION 110.787
TV0 PBASE REGION (BTU/B-P) REFRIGERANT SIDE 2933.419 AIR SIDE 1436.706 COKBINED 964.380
LBWH
PSI BTU/H-SQ FT-F BTU/B-Sp
FT-F
BTU/H-SQ PT-F
SUB~OGLEDREGION (BTU/H-F) REFRIGERANT SIDE 394.035 AIR SIDE 789.615 COMBINED 262.862
n .ow
CONTROLDRVICE - CONDENSERSUBCOOLINGIS SPECIFIED AS 45.000 F CORRESPONDINGTXV RATING PARAMRTERS: CORRESPONDINGCAPILLARY TUBE PARAMETERS: RATRD OPERATING SUPERHEAT 11.000 F NUMBEROF CAPILLARY TUBES 1 STATIC SUPERHEATRATING 6.000 F CAPILLARY TUBE FLOW FACTOR 2.626 PERMANENTBLEED FACTOR 1.150 1.880 TONS TXV CAPACITY RATING INCLUDING NOZZLE AND TUBES
PERFORMANCEOF EACH CIRCUIT IN THE EVAPORATOR INLET AIR TEMPERATURE 47.000 HEAT LOSS FROM COMPRESSOR 4791.4 HEAT LOSS FROM FAN 2274.5 AIR TEMPERATURECROSSING COIL 49.696 OUTLET AIR TEMPERATURE 38.908
CORRESPONDINGORIFICE PARAMRTER: ORIFICE DIAHETER 0.0544 IN
F BTU/H BTU/H F F
MOISTURE REMOVAL OCCURS SUNMARYOF DEHUMYDIFICATION PERFORMANCE(TWO-PHASE REGION) LEADING EDGE OF COIL AIR DRY BULB TEMPERATURE 49.696 F 0.00475 HUMIDITY RATIO 17.090 BTU/LBH ENTHALPY
POINT WHEREMOISTURE REMOVAL BEGINS AIR WALL 37.740 F 44.912 F 0.00475 0.00475 14.191 BTU/LB!+! 15.930 BTU/LBM
RATE OF MOISTURE REMOVAL FRACTION OF NAPORATOR THAT IS WET LATENT HEAT TRANSFER RATE IN TWO-PHASEREGION SENSIBLE HEAT TRANSFER RATE IN TWO-PHASEREGION SENSIBLE TO TOTAL HEAT TRANSFER RATIO FOR TWO-PHASE REGION
0.7345 LBM/H 0.6941 789. BTU/H 6896. BTU/H 0.8974
OVERALL SENSIBLE TO TOTAL HEAT TRANSFER RATIO
0.8996
OVERALL CONDITIONS ACROSS COIL ENTERING AIR DRY BULB TEMPERATURE 49.696 F 43.924 F WET BULB TMPERATURE 0.633 RELATIVE HUMIDITY 0.00475 HUMIDITY RATIO
AIR SIDE MASS FLOW RATE PRESSUREDROP HEAT TRANSFER COEFFICIENT DRY COIL WET COIL DRY FIN EFFICIENCY WET FIN EFFICIENCY
0.6476
SUPERHEATED REGION 0.7608 0.4834 0.4025 0.0682 172.2 ,-.- BTU/H ~~_~~ 184.36 LBM/H 45.845 F 2702.1 LBM/H 0.395 IN HZ0 14.462 BTU/H-SQ FT-F 16.803 BTU/H-SQ FT-F 0.791 0.763
COIL WALL 35.434 F 0.00433 13.183 BTU/LBM
EXITING AIR 38.908 F 37.659 F 0.901 0.00448
TOTAL HEAT EXCHANGEREFFECTIVENESS (SENSIBLE)
NTU HEAT EXCHANGEREFFECTIVENESS CR/CA FRACTION OF HEAT EXCHANGER HEAT TRANSFER RATE AIR MASS FLOW RATE OUTLET AIR TEMPERATURE
LEAVING EDGE OF AIR 38.400 F 0.00446 14.036 BTU/LBM
TWO-PHASE REGION 1.0590 0.6532 0.9318 7684.9 BTU/H 2517.70 LBM/H 38.400 F REFRIGERANT SIDE MASS FLOW RATE 103.5 PRESSUREDROP 7.005 HEAT TRANSFER COEFFICIENT VAPOR REGION 63.625 TN0 PHASE REGION 631.979
LBM/H PSI ,- BTU/H&Q Fr-r BTU/H-S9 FT-F
UA VALUES: REFRIGERANT SIDE AIR SIDE DRY COIL WET COIL COMBINED DRY COIL WET COIL
VAPOR REGION 69.286 261.250 54.763
TWO PHASE REGION 9398.352 BTU/H-F 1091.526 2772.596
BTU/H-F BTU/H-F
791.186 BTU/H-F 1945.610 BTU/H-F
SUMMARYOF ENERGYINPUT AND OUTPUT: SAMPLE DATA SET FOR ORNL HEAT PUMP DESIGN MODEL -- @ 47.0F AMBIENT AIR TBMPERATUREINTO EVAPORATOR 47.00 F HEAT FROM CONDENSERTO AIR 38434. BTU/H HEAT TO EVAPORATORFROM AIR 71428. BTU/H POWERTO INDOOR FAN -1257. BTU/H POWERTO OUTDOORFAN BTU/H 2275. POWERTO COMPRESSORMOTOR 13690. BTU/H COMPRESSOR SHELL HEAT LOSS 4791. BTU/H TOTAL HEAT TO/FROM INDOOR AIR 39691. BTU/H SYSTEM EFFICIENCY: cop (HEATING) COP (WITH RESISTANCE HEAT)
2.305 2.305
HEAT OUTPUT: HEAT FROM HEAT PUMP RESISTANCE HEAT HOUSELOAD
39691.29 BTU/H 0.0 BTU/H 15000.00 BTU/H
101
C - Definitions
Appendix
of Constants
A number of variables Design
Model
that
quently
they
are
simply
of input
data.
each set
DATA subroutine routines into
several
they
order.
Physical
Air-Side
values
values
which
are used via
common blocks.
determined
by-use;
and consespecified in the
passed
are defined
with
BLOCK
to the
These data
they
sub-
are divided in that
Parameters of filter
AHEATR
cross-sectional 1.28 ft*
area
of resistance
PA
atmospheric
RACKS
number of resistance
RAU
universial
pressure,
14.7 heater
gas constant,
Tubing
on indoor heaters
unit,
2.78
in
indoor
ft* unit,
lbf/in* racks,
53.34
3.0
ft-lbf/lbm-OR
Parameter
roughness Motor
Heat Pump
being
together
are in turn
area
Compressor
often,
than
They have been brought
cross-sectional
E
very
rather
AFILTR
Refrigerant
in BLOCK DATA Subroutine
are used by the
to be changed
assigned
categories
"natural"
and constants
are unlikely
and given
where
Assigned
of interior
tube
walls,
5
x
low6 ft
Efficiency
CETAM
coefficients for the 0th through 5th order terms of the fit of the compressor motor efficiency as a function of the fractional motor load (Eq. 4.29), 0.4088, 2.5138, -4.6289, 4.5884, -2.3666, and 0.48324
RPMSLR
slope of linear fit speed as a function (see Eq. 4.24)
for fraction of fraction
of no-load compressor motor of full load power, -0.042
102 --Data for COFAN
Outdoor
Fan Efficiency
constant term for the fit of outdoor fan specific speed, -3.993
fan
static
efficiency
to
. .
ClFAN
coefficient for static efficiency
the linear term of the fit of outdoor to the fan specific speed, 4.266
C2FAN
coefficient for static efficiency
the quadratic term of the fit of outdoor to fan specific speed, -1.024
EFFMOT
constant
motor
efficiency
RPMFAN
constant
motor
speed for
Tolerance
Parameters
AMBCON
convergence 0.20 F".
parameter
CMPCON
convergence calculations,
parameter for the Ah iteration 0.05 Btu/lbm
CNDCON
convergence parameter for iterations and refrigerant temperature at the 0.20 F"
CONMST
tolerance parameter temperatures, 0.003
for F"
the
EVPCON
tolerance 0.5O"F.
for
iterations
FLOCON
convergence rate through
TOLH
tolerance parameter used in calculating properties heated vapor when enthalpy is known, 10e3
of super-
TOLS
tolerance parameter used in calculating properties heated vapor when entropy is known, 10 4
of super-
Refrigerant NR
the outdoor
the outdoor
for
fan,
the outdoor
iterations
fan,
fan
0.55
825 rpm
ambient
air
interation,
in the
compressor
on condenser subcooling flow control device, 0.20 on evaporator on evaporator
criteria for iterations the expansion device,
fan
tube
wall
superheat,
on refrigerant 0.2 lbm/h
mass flow
Specification refrigerant
Refrigerant-Side ClR,
parameter
for
.
number,
Heat Transfer
C3R, and C5R (Eq. 6.2),
coefficients 1.10647,
22 Correlation for 3.5194
for
Single-Phase
Regions
single-phase heat transfer x 10m7, and 0.01080
coefficient
.^
103
I
exponents -0.78992,
CZR, C4R, and C6R (Eq. 6.3),
for single-phase heat 1.03804, and -0.13750
transfer
coefficient
-.
. I
ULR
lower limit refrigerant,
for the 6,000
XLLR
upper limit ant, 3,500
on the
Reynolds
Valve
Constants
Thermostatic BLEEDF
Expansion
Reynolds
bypass or bleed meters
factor when condenser
number for number for
coefficient subcooling
is
turbulent laminar
flow flow
of
of refriger-
used to compute TXV paraheld fixed, 1.15
DPRAT
rated pressure drop across the TXV at the design 100 psi for R-22 and R-502, 50 psi for R-12
conditions,
NZTBOP
switch to bypass nozzle and distributor tube calculations when calculating TXV parameters subcooling is held fixed, 0
pressure drop if the condenser
STATIC
static superheat setting used to compute the when the condenser subcooling is held fixed,
TXV parameters
SUPRAT
rated operating superheat used to compute the TXV parameters when the condenser subcooling is held fixed, 1l.O.F"
TERAT
rated
TLQRAT
rated liquid lOO.O"F
Capillary
Tube Parameter
NCAP
number of capillary tubes factor, $, when condenser
6.0 F”
f
Compressor
evaporating
Efficiency
temperature
refrigerant
for
the TXV, 4O.O"F
temperature
at the
inlet
to the
used to compute capillary subcooling is held fixed,
tube 1
TXV,
flow
Parameters
ETAVLA
intercept volumetric parameter
of the fit of theoretical minus actual compressor efficiencies as a linear function of the correlating of Davis and Scott, -0.0933
ETAVLB
slope of the fit of theoretical minus actual volumetric efficiencies as a linear function parameter of Davis and Scott, 0.733
SUCFAC
suction
VOLFAC
volumetric Eq. 4.4,
gas heating efficiency 0.75
factor
Fsh used
correction
compressor of the correlating
in Eq. 4.6,
factor
0.33
Fv used in
105 Appendix
D:
Definitions
of Variables
in Common Blocks
Variable A
Definition constant used in a first order approximation of the saturation temperature in the logarithm of the refrigerant pressure. assigned a value in TABLES and used in TSAT.
Common SUPER
AAFC
the frontal area of the condenser, by CALC and COND (sq ft).
CONDEN
AAFE
the frontal area of the evaporator, by CALC and EVAPR (sq ft).
read by DATAIN, used read by DATAIN, used
EVAPOR
ACV,BCV,..., FCV constants used in calculating refrigerant heat capacity with temperature and specific volume, assigned values in TABLES and used in SATPRP.
OTHER
AFILTR
AIR
frontal area of filter for the indoor unit, in BLOCK DATA and used in PDAIR (sq ft).
assigned
AHEATR cross-sectional area of the resistance heaters, value in BLOCK DATA and used in PDAIR (sq ft).
a value
assigned
a
AIR
AL.BL,...,GL coefficients for expressing density of the liquid refrigerant as a function of absolute temperature (R), assigned values in TABLES and used in SATPRP (lbm/cu ft).
DENSIT
ALFAAC ration of the air-side heat transfer area for the condenser to the total volume of the condenser, calculated in CALC, used by COND, EFFCT, and EXCH (/ft).
CONDEN
ratio of the air-side heat transfer area for the evaporator to the total volume of the evaporator, calculated in CALC, used by EVAP and EVAPR (/ft).
EVAPOR
ALFARC ratio of the total refrigerant-side heat transfer area for the condenser to the volume of the condenser, calculated in CALC, used by COND, EFFCT, and EXCH (/ft).
CONDEN
ALFARE ratio of the total refrigerant side heat transfer area for the evaporator to the volume of the evaporator, calculated in CALC, used by EVAP and EVAPR (/ft).
EVAPOR
ALPHA
STATEQ
AALFE
constant used in refrigerant equation of state correlating pressure with temperature and specific volume, assigned a value in TABLES, used by SATPRP, SPVOL, and VAPOR.
AMBCON tolerance parameter for the ambient air temperature iteration in the main program, assigned a value in BLOCK DATA and used in the main program (F).
MPASS
ARFTC
CONDEN
inside cross-sectional area of the tubes in the condenser calculated in CALC and used in COND (ft).
,
106
ARFTE
inside cross-sectional area of the tubes in the evaporator calculated in CALC and used by EVAPR (sq ft).
EVAPOR
ARHTC
total refrigerant-side heat transfer area for each circuit in the condenser, calculated and used in CALC, used by COND, EFFCT, and EXCH (sq ft).
CONDEN
ARHTC
total refrigerant side heat transfer area for each circuit in the evaporator, calculated and used in CALC, used by EVAP and EVAPR (sq ft).
EVAPOR
used in correlating refrigerant AVP,BVP,..., FVP constants pressure with absolute temperature (R), assigned in TABLES, used in SATPRP and TSAT. A2,B2,C2
vapor values
constants used in refrigerant equations of state correlating pressure with temperature and specific volume, assigned values in TABLES, used by SATPRP, SPVOL, and VAPOR.
A3,B3,C3
SAT
STATEQ
constants used in refrigerant equations of state correlating pressure with temperature and specific volume, assigned values in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
constants used in refrigerant equations of state correlating pressure with temperature and specific volume, assigned values in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
A5,B5,C5
constants used in refrigerant equations of state correlating pressure with temperature and specific volume, assigned values in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
A6,B6,C6
constants used in refrigerant equations of state correlating pressure with temperature and specific volume, assigned values in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
B
constant used in a first order approximation to the refrigerant saturation temperature as a function of the logarithm of the saturation pressure, assigned value in TABLES and used in TSAT.
SUPER
A4,B4,C4
BLEEDF by-pass or bleed factor coefficient expansion valve, assigned a value DATAIN, used by TXV. Bl
for the thermostatic by BLOCK DATA or read by
constant used in refrigerant equation of state correlating pressure with temperature and specific volume, assigned a value in TABLES, used by SATPRP, SPVOL, and VAPOR.
* c
1
TXVDAT
STATEQ
_
107 CA c
.
air-side heat capacity rate of the condenser, calculated
for the superheated portion of the in EFFCT, used by EXCH (Btu/h-F).
capacity of the evaporator at the operating conditions with an assumed 10 F of superheat at the evaporator exit, calculated in and used in TXV, printed by OUTPUT (Btu/hr),
FLOWBA
CAPNOZ heating capacity valve, calculated
of the nozzle of the thermostatic expansion and used in TXV, printed by OUTPUT (Btu/hr).
TXVDAT
CAPTUB heating capacity of each tube from the thermostatic expansion valve to the evaporator, calculated and used in TXV, printed in OUTPUT (Btu/hr).
TXVDAT
CARC
ratio of the fin heat transfer area to the contact area between the fins and the tubes in the condenser, calculated in CALC and used in COND.
CONDEN
CARE
ratio of the fin heat transfer. between the fins and the tubes in CALC and used by EVAP.
area to the contact area in the evaporator, calculated
EVAPOR
CETAM
coefficients for a fit to the motor efficiency as a function of the fraction of full load power output, assigned values in BLOCK DATA and used by COMP.
CMPMOT
CLCXMR refrigerant calculated (lbm/hr).
mass flow rate through each capillary in CAPTUB, used by FLOBAL, and printed
CMPCON tolerance parameter for entering the compressor in BLOCK DATA (Btu/h-F).
the iteration suction port
tube, in OUTPUT
on the enthalpy in COMP, assigned
COIL
evaporator coil resistance to heat transfer, used in EVAP, used by TWISOL.
calculated
PRNT2
MPASS a value
parameter for the condenser subcooling iterations CNDCON tolerance in the main program and in COND, and the refrigerant temperature at the flow control device in FLOBAL, assigned value in BLOCK DATA, used in the main program, COND, and FLOBAL (F). II
TXVDAT
CAPFLO capillary tube flow factor, default value assigned in BLOCK DATA, read by DATAIN if capillary tube submodel is used, calculated or used by CAPTUB, printed by OUTPUT.
2.
‘
COMPR
as CANFAC factor used to express compressor shell heat rejection a fraction of compressor power consumption, read by DATAIN end used by COMP. CAP
EFFARG
and
MPASS a
TWI
108 parameter for the iterations CONMST tolerance temperature in EVAP, assigned a value in EVAP (F).
on the tube wall in BLOCK DATA and used the heat pump with in the main program
MPASS
COP
operating coefficient of performance for any necessary resistance heat, calculated and printed by OUTPUT.
COPHP
operating coefficient of performance for the heat pump, does PRNT8 not include any resistance heat, calculated by the main program and printed in-OUTPUT.
CPA
specific heat of dry air, calculated COND and EVAPR, used by COND, EVAPR, XMOIST, and WBF (Btu/lbm-F).
AIR
CPM
specific heat of moist air, calculated by CONDand EVAPR and used by COND, EFFCT, EVAP, and EXCH (Btu/lbm-F).
AIR
CPOW
coefficients for compressor power consumption as a function MAPFIT of the condensing and evaporating temperatures, read by DATAIN if the map-based model is used and used by CMPMAP.
CPR
constant used in refrigerant equation of state correlating pressure with temperature and specific volume, assigned a value in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
CPSCC
specific heat for the subcooled refrigerant in the condenser, calculated by MUKCPand EXCH, used by COND, EFFCT, EXCH, and SPHTC (Btu/lbm-F).
CONDS
CPSP
specific heat for the superheated refrigerant in the condenser calculated in COND, modified in EXCH, and used by COND, EFFCT, and EXCH (Btu/lbm-F).
CONDS
CPSPC
specific heat for the condensing superheated refrigerant in the condenser, calculated in COND and used by COND and EXCH (Btu/lbm-F).
CONDS
CPSPE
specific heat for the superheated refrigerant ator calculated and used by EVAPR (Btu/lbm-F).
EVAPS
CR
refrigerant-side heat capacity rate for the superheated region of the condenser, calculated in EFFCT and used in EXCH (Btu/hr-F).
in the evapor-
PRNT8
EFFARG
CSCOOL ratio of refrigerant to air-side heat capacity rates for the subcooled region of the condenser, calculated in EXCH and printed in OUTPUT.
PRNT4
CSUPER ratio of the refrigerant to air-side heat capacity rates for the superheated region of the evaporator, calculated in EVAP and printed in OUTPUT.
PRNT5
e *
109 CUASC
.
overall thermal conductance times the heat transfer area (UA) for the subcooled region of the condenser, calculated in COND and printed in OUTPUT (Btu/hr-F).
thermal conductance times the heat transfer CUASCA air-side for the subcooled region of the condenser, calculated and printed in OUTPUT (Btu/hr-F).
area (UA) in COND
thermal conductance times the heat transfer CUASCR refrigerant-side area (UA) for the subcooled region of the condenser, calculated in COND and printed in OUTPUT (Btu/hr-F). CUATP
overall thermal conductance times heat transfer area (UA) for the two-phase region of the condenser, calculated in COND and printed in OUTPUT (Btu/hr-F).
thermal conductance times the heat transfer CUATPA air-side for the two-phase region of the condenser, calculated and printed in OUTPUT (Btu/hr-F).
area (UA) in COND
UAS
UAS
UAS
UAS
UAS
CUATPR refrigerant-side thermal conductance times the heat transfer area (UA) for the two-phase region of the condenser, calcuated in COND and printed in OUTPUT (Btu/hr-F).
UAS
CUAV
overall thermal conductance times the heat transfer area for the vapor region of the condenser, calculated in COND and printed in OUTPUT (Btu/hr-F).
UAS
CUAVA
air-side thermal conductance times the heat transfer area (UA) for the vapor region of the condenser, calculated in COND and printed in OUTPUT (Btu/hr-F).
UAS
CUAVR
refrigerant-side thermal conductance times heat transfer for the superheated region of the condenser, calculated COND and printed in OUTPUT (Btu/hr-F).
UAS
CXMR
coefficients for the refrigerant mass flow rate as a function of the condensing and evaporating temperatures, read by DATAIN, if the map-based compressor model is chosen, and used by CMPMAP.
area in
MAPFIT
COFAN,ClFAN,C2FAN coefficients for the fan static efficiency for the outdoor heat exchanger as a function of the static speed of the fan, assigned values in BLOCK DATA and used by FANFIT.
FANMOT
coefficients and exponents for ClR,C2R ,...,C6R heat transfer coefficient as a function Reynolds and Prandtl numbers, assigned and used in SPHTC.
VAPHTC
fit to single-phase of refrigerant values in BLOCK DATA
/
110
DDL
inside diameter of the compressor reversing valve to the condenser, and used by COMP (ft).
discharge line from the assigned a value by DATAIN
LINES .
DDLRV
inside diameter of the line from the compressor discharge port to the reversing valve, read by DATAIN and used by COMP (ft).
LINES
DDUCT
diameter of each of 6 identical air ducts with equivalent lengths of 100 ft, read by DATAIN and used by COND and EVAPR (ft).
BLANK
DEAC
outside diameter of the tubes in the condenser, DATAIN, used by CALC and COND (ft).
CONDEN
DEAE
outside diameter of the tubes in the evaporator, DATAIN, used by CALC, EVAP and EVAPR (ft).
of the fins DELTAC thickness by CALC and COND (ft).
in the condenser,
read by read by
read by DATAIN, used
of the fins in the evaporator, DELTAE thickness by CALC and EVAPR (ft).
read by DATAIN, used
-
EVAPOR CONDEN EVAPOR
DERC
inside diameter of the refrigerant tubes in the condenser, read by DATAIN, used by CALC and COND (ft).
CONDEN .
DERE
inside diameter of the refrigerant tubes in the evaporator, read by DATAIN, used by CALC and EVAPR (ft).
EVAPOR *
DISPL
total compressor by COMP (cu in).
COMPR
DLL
inside diameter of the liquid line connecting the indoor and outdoor heat exchangers, read by DATAIN and used by FLOBAL (ft).
LINES
DP
air-side pressure drop across the outdoor coil, calculated in FANFIT and printed by OUTPUT (inches of water).
PRNTS
DPDL
pressure drop in the compressor discharge line from the discharge port to the condenser, calculated and used by COMP (psi).
LINES
DPLL
pressure drop in the liquid line connecting the heat exchangers, calculated and used in FLOBAL (psi).
LINES
DPNOZ
pressure drop in the nozzle for the thermostatic expansion valve, calculated and used in TXV, print in OUTPUT (psi).
TXVDAT -
DPRAT
rated pressure drop across the thermostatic at the design conditions, assigned a value used by TXV (psi).
TXVDAT *
piston
displacement,
read by DATAIN and used
expansion valve in BLOCK DATA and
111 DPSL
LINES
pressure drop in the suction line connecting the evaporator and compressor, calculated and used in COMP (psi).
expanDPTUBE pressure drop across the tubes from the thermostatic sion valve to the evaporator, calculated and used in TXV, printed in OUTPUT (psi).
h s
TXVDAT
DPTXV
pressure drop across the valve of the thermostatic valve (does not include nozzle and tube losses), and used in TXV, printed in OUTPUT (psi).
expansion calculated
DSL
inside diameter to the reversing
DSLRV
inside diameter of the line from the compressor suction port to the reversing valve, read by DATAIN and used by COMP (ft).
LINES
DTROC
specified condenser subcooling, read by DATAIN used by CNDNSR, FLOBAL, and FLODEV, and printed by OUTPUT (F).
FLOWBA
DZC
equivalent calculated
length of the refrigerant tubing in CALC and used by COND.
in the condenser,
CCNDEN
DZE
equivalent calculated
length of the refrigerant tubing in CALC and used by EVAPR.
in the evaporator,
EVAPOR
E
roughness factor for all of the refrigerant lines, assigned a value in BLOCK DATA, used by COMP, COND, EVAPR, and FLOBAL (ft).
TXVDAT
of the suction line from the evaporator valve, read by DATAIN and used by COMP (ft).
LINES
LINES
and the mechanical EFFMMX product of the maximum motor efficiency efficiency for the compressor, read by DATAIN and used by COMP. for the outdoor EFFMOT fan motor efficiency BLOCK DATA and used in FANFIT.
fan,
assigned
COMPR
a value
fan, calculated printed by
in
EINDF
rate of power consumption for the indoor COND or EVAPR, used by the main program, OUTPUT (Btu/hr).
in
EOUTF
rate of power consumption for the outdoor fan, calculated COND or EVAPR, used by the main program, printed by OUTPUT (Btu/hr).
ETAISN
isentropic efficiency from the suction port to the discharge port of the compressor, read by DATAIN and used by COMP, printed and used by OUTPUT.
in
FANMOT PRNT8
PRNT8
COMPR
112
ETAMEC mechanical efficiency used by COMP.
of the compressor,
read by DATAIN and
ETAMOT operating motor efficiency for used in COMP, used and printed
the PRNTlessor, by OUTPUT.
ETASUP compressor and printed
efficiency,
suction gas heating by OUTPUT.
calculated
calculated
and
by COMP
COMPR
COMPR * PRNTI
isentropic efficiency from the shell inlet ETATOT compressor overall to the shell outlet, calculated by COMP, used and printed by OUTPUT.
PRNTl
of the fit of theoretical minus actual compressor ETAVLA intercept volumetric efficiency as a linear function of the correlating parameter of Davis and Scott, assigned in BLOCK DATA and used in COMP.
COMPR
minus actual compressor ETAVLB slope of the fit of theoretical volumetric efficiency as a linear function of the correlating parameter of Davis and Scott, assigned in BLOCK DATA and used in COMP.
COMPR
efficiency of the compressor based on the shell ETAVOL volumetric inlet conditions, calculated and used in COMP, printed by OUTPUT.
PRNTl
ETP
heat exchanger effectiveness for the two-phase region of the calculated and used in EXCH, used and printed in condenser, OUTPUT.
PRNT4
ETPE
heat exchanger effectiveness for the two-phase region of the evaporator, calculated and used in EVAP, printed by OUTPUT.
PRNT5
EUATP
overall thermal conductance times heat transfer area (UA) for the two-phase region of the evaporator, calculated in EVAPR, printed in OUTPUT (Btu/hr-F).
UAS
thermal conductance times the heat transfer area (UA) EUATPA air-side for the two-phase region of the evaporator, calculated in EVAPR and printed in OUTPUT (Btu/hr-F).
UAS
thermal conductance times the heat transfer EUATPR refrigerant-side area (UA) for the two-phase region of the evaporator, calculated in EVAPR and printed in OUTPUT (Btu/hr-F).
UAS
EUAV
EUAVA
i
overall thermal conductance times the heat transfer area (UA) for the superheated region of the evaporator, calculated in EVAPR and printed in OUTPUT (Btu/hr-F).
UAS
air-side thermal conductance times the heat transfer area (UA) for the superheated region of the evaporator, calculated in EVAPR and printed in OUTPUT (Btu/hr-F).
UAS
L 1 -
.
113 UAS
tolerance parameter for the iteration on evaporator in the main program and in EVAPR, assigned a value DATA and used in the main program and EVAPR (F).
MPASS
superheat in BLOCK
heat exchanger effectiveness for the subcooled region condenser based on effectiveness-NTU correlations for flow heat exchangers, calculated and used in EXCH.
of the cross-
EFFARG
EXFRC
heat exchanger effectiveness condenser, assigned a value OUTPUT.
of the in
PRNT4
EXFRE
heat exchanger effectiveness for the superheated evaporator, calculated and used in EVAP, printed
region of the by OUTPUT.
PRNT5
EXFS
heat exchanger effectiveness for the superheated condenser based on heat capacity rates and delta ulated in EFFCT and used by EXCH.
region of the T's, calc-
EFFARG
EXFSC
heat exchanger the condenser, in OUTPUT.
' EXFR
3.
refrigerant-side thermal conductance times the heat transfer area (UA) for the superheated region of the evaporator, CalCulated in EVAPR and printed in OUTPUT (Btu/hr-F).
for the subcooled region in EXCH, used and printed
effectiveness for the superheated region of assigned a value in EXCH, used and printed
PRNT4
for the condenser, the FANEFC the combined fan-fan motor efficiency value is either read by DATAIN or calculated by FANFIT;used by COND.
CONDEN
for the evaporator, FANEFE the combined fan-fan motor efficiency the value is either read by DATAIN or calculated by FANFIT, used by COND.
EVAPOR
efficiency FANOUT combined operating calculated in FANFIT, printed
PRNT8
of the outdoor by OUTPUT.
fan-fan
motor,
FARC
ratio of the fin heat transfer area for the condenser to the total heat transfer area on the air side, calculated in CALC and used in COND.
CONDEN
FARE
ratio of the fin heat transfer area for the evaporator to the total heat transfer area, calculated in CALC, used by EVAP and EVAPR.
EVAPOR
FINTYC
switch to specify the type of fins in the condenser (smooth, wavy, or louvered), read by DATAIN, used by HAIR and PDAIR.
CONDEN
114 FINTYE
switch to specify the type of fins in the evaporator read by DATAIN, used by HAIR (smooth, wavy, or louvered), and PDAIR.
EVAPOR
FIXCAP
specified heating load for the heat pump, read by DATAIN and used by the main program (Btu/hr).
BLANK
FLMOT
full load motor power for the compressor, either read by DATAIN or calculated in COMP, used by COMPand printed by OUTPUT (kW).
COMPR
FLOCON relative tolerance parameter used in the refrigerant mass flow rate iteration in the main program, assigned a value in BLOCK DATA and used in the main program.
MPASS
FMOIST fraction of the evaporator where moisture removal from the air and used in EVAP, used in EVAPR and printed occurs, calculated in OUTPUT.
EVAPRS
CONDEN
FPC
the fin pitch for the condenser, by CALC and COND (fins/ft).
read by DATAIN and used
FPE
the fin pitch for the evaporator, by CALC and EVAPR (fins/ft).
FSCC
fraction of each circuit in the condenser which is subcooled calculated in EXCH, used by COND and EXCH, printed in OUTPUT.
CONDS
FSPC
fraction of each circuit in the condenser that is superheated, calculated and used in EXCH, used in COND, and printed in OUTPUT.
CONDS
FSPE
fraction heated, printed
EVAPS
FTPC
fraction of each circuit in the condenser which is in twophase flow, calculated and used in EXCH, used in COND, and printed in OUTPUT.
CONDS
FTPE
fraction of each circuit in the evaporator which is in twophase flow, calculated and used in EVAP, used in EVAPR, and printed by OUTPUT.
EVAPS
HAC
air-side heat transfer coeffiqient for the condenser, calculated by HAIR, used by COND, EFFCT, EXCH, and SEFF, printed by OUTPUT (Btu/hr-sq ft-F).
CONDS
air-side heat transfer coefficient for the evaporator, calculated and used by EBAPR, used by EVAP, and printed OUTPUT (Btu/hr-sq ft-F).
EVAPS
read by DATAIN and used
EVAPOR ' ,
HAE
of each circuit in the evaporator that is supercalculated and used in EVAP, used in EVAPR, and in OUTPUT.
by
L *
115 HAIR1
enthalpy of the moist air entering the evaporator, calculated and used in EVAP, used by TWISOL, and printed by OUTPUT (Btu/lbm).
TWI
HAIR0
enthalpy of the air leaving the evaporator, used in EVAP, printed by OUTPUT (Btu/lbm).
PRNT5
HAIR1
enthalpy of the air in the evaporator at the point where and used in EVAP, printed moisture removal begins, calculated by OUTPUT (Btu/lbm).
PRNT5
HAW
air-side heat transfer coefficient for the. wetted surface of the evaporator, calculated in EVAPR, used by EVAP and EVAPR, printed by OUTPUT (Btu/h-sq ft-F).
EVAPRS
.
calculated
and
A
. .
HCONTC contact conductance between the fins and tubes in the ft-F). condenser, read by DATAIN and used by COND (Btu/hr-sq
CONDEN
HCONTE contact conductance between the fins and the tubes evaporator, read by DATAIN, used by EVAij.~~iu/hr-sq'
EVAPOR
in the ft-F).
HIC
enthalpy of the refrigerant entering the condenser, calculated in COMP, modified by DESUPR, used by COMPand FLOBAL, printed by OUTPUT (Btu/lbm-F).
CONDSR
HIE
enthalpy of the refrigerant entering the evaporator, calculated in the main program, used by EVAPR, then calculated and used by FLOBAL, used by TXV, and printed by OUTPUT (Btu/lbm-F).
EVAPTR
HILOFC factor used to express the internal compressor heat transfer from the high side discharge line to the low side suction line as a fraction of compressor power, read by DATAIN and used by COMP.
COMPR
HINCMP enthalpy of the refrigerant at the compressor calculated and used in COMPor CMPMAP;'printed (Btu/lbm).
CMPRSR
shell inlet in OUTPUT
HOC
enthalpy of the refrigerant leaving the condenser, calculated CONDSR in COND and FLOBAL, used by the main program, COND, and FLOBAL, printed by OUTPUT (Btu/lbm-F).
HOE
enthalpy of the refrigerant leaving the evaporator, calculated and used by COMP, used by FLOBAL, and printed by OUTPUT (Btu/lbm-F).
at the compressor shell outlet HOUCMP enthalpy of the refrigerant calculated and used in COMPor 'CMP~P,"pn;ntea~in"OUTPUT (Btu/lbm).
EVAPTR
CMPRSR
116 refrigerant-side heat transfer coefficient for the subcooled region of the condenser, calculated in SPHTC, used by CHTC, COND, and EXCH, printed by OUTPUT (Btu/hr-sq ft-F).
CONDS
HSPC
refrigerant-side heat transfer coefficient for the superheated region of the condenser, calculated by SPHTC, used by COND, EFFCT, and EXCH, printed by OUTPUT (Btu/hr-sq ft-F).
CONDS .
HSPE
refrigerant-side heat transfer coefficient for the superheated region of the evaporator, calculated and used by EVAPR, also used by EVAP, printed by OUTPUT (Btu/hr-sq ft-F).
EVAPS
HTPC
refrigerant-side heat transfer coefficient for the two-phase CONDS region of the condenser, calculated in CHTC, adjusted by COND, used by COND and EXCH, printed by OUTPUT (Btu/hr-sq ft-F).
HTPE
refrigerant-side heat transfer coefficient for the two-phase calculated and used by EVAPR, also region in the evaporator, used by EVAP, printed by OUTPUT (Btu/hr-sq ft-F).
HSCC
HWALLI enthalpy of the air at the tube walls on the leading edge of the evaporator where moisture removal begins, calculated in EVAP, calculated and used in TWISOL, printed in OUTPUT (Btu/lbm-F). HWALLO enthalpy of the air at the tube walls on the leaving the evaporator, calculated and used in EVAP, printed (Btu/lbm-F).
.
EVAPS
TWI
PRNT5 edge of in OUTPUT
ICOMP
switch to specify which compressor DATAIN and used by FLOBAL.
read by
FLOWBA
IREFC
switch to specify type of flow control device, read by DATAIN, used by the main program, CAPTUB, FLOBAL, OUTPUT, and TXV.
FLOWBA
ITITLE
descriptive and printed
BLANK
title for labelling by OUTPUT.
model is used,
output,
read by DATAIN
constant used in correlating refrigerant heat capacity at a constant volume with temperature and specific volume, assigned a value in TABLES and used in SATPRP.
.
OTHER
* .
117 constant used in refrigeration equations of state correlating pressure with temperature and specific volume, assigned a value in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
LElO
multiplication constant used to convert to natural logarithms, assigned a value SATPRP and TSAT.
SUPER
LPRINT
switch to control levels of printed output, read by DATAIN, used by the main program, COMP, COND, EVAP, EVAPR, and FLOBAL.
MPASS
LlOE
multiplicative to logarithms SATPRP.
used to convert natural logarithms assigned a value in TABLES, used in
SUPER
for adding heat from the compressor read by DATAIN, used by COND,
MPASS
MFANFT switch to specify whether or not to use efficiency curve for the outdoor fan, read by DATAIN, used by the main program, COND, EVAPR, and OUTPUT.
MPASS
MFANIN switch to specify option for adding heat loss from the indoor fan to the air-stream, read by DATAIN, used by COND, EVAPR, and OUTPUT.
MPASS
MFANOU switch to specify option for adding heat loss from the outdoor fan to the air-stream, read by DATAIN, used by COND, EVAPR, and OUTPUT.
MPASS
MTRCLC switch ressor
COMPR
. .
*.
. c
constant base ten,
MCMPOP switch to specify option shell to the air-stream, EVAPR, and OUTPUT.
logarithms base ten in TABLES, used by
used to specify whether or not to calculate the compmotor speed, read by DATAIN, used by COMP.
MUNITC switch to specify whether the condenser is the indoor or the outdoor heat exchanger, assigned a value by DATAIN -nd used by COND.
CONDEN
MUNITE switch to specify whether the evaporator outdoor heat exchanger, assigned a value by EVAPR.
EVAPOR
is the indoor or the by DATAIN and used
118
NCAP
number of capillary tubes used, default value assigned in BLOCK DATA, optionally read by DATAIN if capillary tube submodel is used, used by CAPTUB and printed by OUTPUT.
FLOWBA
NCORH
switch to specify heating or cooling by DATAIN, used by the main program,
MPASS
NR
refrigerant number, assigned a value in BLOCK DATA, used by the main program, MUKCP, SATPRP, TSAT, TXV, and VAPOR.
mode of operation, read COND, EVAPR, and OUTPUT.
REFRIG
parallel refrigerant circuits in the NSECTC the number of equivalent condenser, read by DATAIN, used by OUTPUT, CALC, and COND.
CONDEN
parallel refrigerant circuits in the NSECTE the number of equivalent evaporator, read by DATAIN, used by EVAPR and OUTPUT.
EVAPOR
NTC
the number of tubes in the direction of' the air flow condenser, read by DATAIN, used by CALC and COND.
in the
CONDEN
NTE
the number of tubes in the direction of the air flow in the evaporator, read by DATAIN, used by CALC and EVAPR.
EVAPOR
tube pressure NZTBOP switch to bypass the nozzle and distributor drop calculations for the thermostatic expansion valve, assigned a value in BLOCK DATA or read by DATAIN, used by TXV.
ORIFD
diameter of the orifice or short tube, read by DATAIN or calculated by ORIFIC, used by ORIFIC and printed by OUTPUT (ft)
-
TXVDAT
FLOWBA
l
PA
the atmospheric pressure, assigned a value in BLOCK DATA, used by COND, EVAP, EVAPR, and XMOIST (psia).
AIR
PC
inside perimeter of each tube in the condenser, and used in CALC, used in COND (ft).
CONDEN
PCTFL
fraction of the full-lOad,compressor motor power, and used in COMP, printed in OUTPUT.
calculated calculated
pressure drop for the condenser, calculated PDAIRC air-side and used by COND and FANFIT, printed by OUTPUT (psi).
by PDAIR
PRNTl CONDS
h d
119
PDAIRE air-side pressure drop for the evaporator, initialized to 0 in the main program, calculated by PDAIR, used by EVAPR, and printed by OUTPUT (psi). PDC
CONDS refrigerant-side pressure drop for the condenser, calculated by PDROP in calls from FLOBAL. and COND, used by FLOBAL and COND, printed by OUTPUT (psi).
PDE
refrigerant-side pressure drop for each circuit in the evaporator, calculated and.used by EVAPR, printed by OUTPUT (psi).
EVAPS
PE
inside perimeter of each tube in the evaporator, and used by CALC, used by EVAPR (ft).
EVAPOR
-5
calculated
PFLODV pressure of the refrigerant entering the flow control device, calculated and used in FLOBAL, printed in OUTPUT (psia).
pressure of the refrigerant entering the condenser, calculated CONDSR in COMP, modified by DESUPR and FLOBAL, used by COMPand COND, printed by OUTPUT (psia).
PIE
pressure of the refrigerant entering the evaporator, calculated and used by EVAPR, also calculated and used by TXV, printed by OUTPUT (psia). shell inlet in OUTPUT
EVAPTR
CMPRSR
POC
pressure of the refrigerant leaving the condenser, calculated in COND and FLOBAL, used by COND and FLOBAL, printed by OUTPUT (psia).
POE
pressure of the refrigerant leaving the evaporator, calculated EVAPTR and used by COMP, used by EVAPR and TXV, printed by OUTPUT (psia).
POUCMP pressure of the refrigerant at the compressor calculated and used in COMPor CMPMAP, printed (psia).
l
PRNT2
PIG
PINCMP pressure of the refrigerant at the compressor calculated and used in COMPor CMPMAP, printed (psia).
.
EVAPS
shell outlet in OUTPUT
POW
compressor power consumption, in the main program, printed
POW2
rate of power consumption used by the main program,
PRA
Prandtl
PRINT
logical variable to control the output of intermediate calculations, set to TRUE or FALSE in the main program depending on the value of LPRINT.
number of air,
computed and used in COMP, used by OUTPUT (kW).
by the compressor, calculated printed by OUTPUT (Btu/hr).
calculated
and
and used by COND and EVAPR.
CONDSR
CMPRSR
COMPR PRNT8 AIR Al
120 QAC
air flow rate through by COND (cfm).
the condenser,
QAE
air flow rate through by EVAPR (cfm).
the evaporator,
QAIR
rate of heat transfer to or from,the air, includes heat transfered from the coil and the indoor fan, calculated and used by the main program, printed by OUTPUT (Btu/hr).
QC
rate of heat transfer from the condenser to the air, ulated in COND and used by the main program, printed (Btu/hr).
QCAN
rate of heat rejection from the compressor shell, read by DATAIN, optionally computed and used in COMP, used in COND and EVAPR, printed in OUTPUT (Btu/hr).
QDISLN heat loss rate in the compressor DATAIN, used by COMP (Btu/hr).
read by DATAIN, used read by DATAIN, used
discharge
-line,
read by
total rate of heat transfer for the evaporator, in EVAPR, used by the main program, and printed (Btu/hr).
QFANC
heat loss from the condenser printed by OUTPUT (Btu/hr).
QFANE
rate of heat loss from the evaporator and printed in OUTPUT (Btu/hr).
QHILO
internal heat tra,nsfer rate from the discharge line of the compressor to the suction line, read by DATAIN, optionally calculated and used by COMP, and printed by OUTPUT (Btu/hr).
QLIQLN rate of heat loss in the liquid line, the main program and FLOBAL (Btu/hr).
calculated by OUTPUT
calculated fan,
calculated
EVAPOR A PRNT8
_
calcCONDS by OUTPUT
QE
fan motor,
CONDEN
by COND,
COMPR
LINES EVAPS
PRNT3
-I
in EVAPR PRNT7
I
read by DATAIN, used by
COMPR
LINES
QMTPL
latent heat transfer rate for the two-phase region of the evaporator, calculated and used by EVAP, printed by OUTPUT (Btu/hr).
PRNTfj
QSCC
rate of heat transfer for the subcooled region of each circuit in the condenser, calculated in EXCH, used by COND and EXCH, and printed and used by OUTPUT (Btu/hr).
CONDS
QSPC
rate of heat transfer for the superheated region of each circuit in the condenser, calculated in EXCH, used by COND and EXCH, andp printed by OUTPUT (Btu/hr).
CONDS
QSPE
rate of heat transfer for the superheated region of each circuit of the evaporator, calculated and used by EVAP, also used by EVAPR, printed and used by OUTPUT (Btu/hr).
EVAPS
.
'
121 QSTPT
sensible heat transfer evaporator, calculated (Btu/hr).
rate for the two-phase region of the and used in EVAP, printed by OUTPUT ._... ...____-..".._. _ . I
PRNT5
.
4
QSUCLN rate of heat gain in the suctin line‘. ". between .,.., .,. the evaporator and the compressor, read by DATAIN and used by either'CMPMAPor COMP (Btu/hr).
LINES
QTPC
rate of heat transfer for the -. ‘.,. two-phase region (.*xIW-I ..__"~1 .^ "_- of each circuit in the condenser;“"' calculated and used in-‘%XCH*, used by COND, printed by OUTPUT (Btu/hr).
CONDS
QTPE
rate of heat transfer for the two-phase region of each circuit of the evaporator, calculated and used by EVAP, also assigned a value and used by EVAPR, printed by OUTPUT (Btu/hr).
EVAPS
R
constant used in refrigerant equations of state correlating pressure with temperature and specific volume, assigned a value in TABLES, used by SATPRP, SPVOL, and VAPOR.
STATEQ
RACKS
number of racks of resistance heaters, BLOCK DATA and used in PDAIR.
AIR
RAU
universal gas constant for air, assigned a value BLOCK DATA, used by COND and EVAPR (ft-lbf/lbm-R).
RESIST
amount of resistance heat-required to meet the house load, calculated and used in the main program, printed by.OUTPUT (Btu/hr).
PRNT8
RHI
the relative humidity of the air entering the evaporator calculated and used in EVAP, printed in OUTPUT.
PRNT5
RHIC
the relative humidity of the air entering the condenser read by DATAIN and used in COND and EXCH.
RHIE
the relative humidity of the air entering read by DATAIN, used by EVAP and EVAPR.
RHO
relative humidity of the air at the leaving edge of the evaporator, calculated and used by EVAP, printed by OUTPUT.
RHTERM constant in dehumidificati-on and used in TWOSOL. RPM
solution,
actual compressor motor speed, assigned COMP, printed by OUTPUT (rpm's).
assigned
a value
in
in
the evaporator
AIR
coil, unit, unit,
CONDEN EVAPOR PRNT5
calculated
in EVAP
TWO
a value
and used in
,PRNTl
122
RPMFAN motor speed for the outdoor fan, DATA, used by FANFIT (rpm's).
assigned
a value
in BLOCK
FANMOT .
RPMSLP slope of the linear fit to the fraction of no-load motor speed as a function of fraction of full-load assigned a value in BLOCK DATA and used in COMP. RTBCND the number of return bends in the condenser, scaled and used by CALC. RTBEVP the number of return bends in the evaporator, scaled and used by CALC.
compressor power,
CMPMOT
read by DATAIN,
CONDEN
read by DATAIN,
EVAPOR
for the superheated region in EFFCT and used by
EFFARG
RTOT
overall resistance to heat transfer of the condenser (l/UA), calculated EXCH (hr-F/Btu).
SEFFD
overall surface efficiency for hexagonal fins (staggered tubes) EVAPRS in the evaporator, calculated and used in EVAP, used by EVAPR, and printed by OUTPUT.
SEFFWA average surface efficiency for the portion of the evaporator where moisture removal occurs, calculated and used in EVAP, used in EVAPR, and printed in OUTPUT.
EVAPRS
SEFFXC overall surface efficiency for the condenser, SEFF, and used by COND, EFFCT, and EXCH.
CONDS
calculated
in
SIGAC
ratio total
of the free flow frontal frontal area, calculated
area of the condenser to the in CALC and used in COND.
CONDEN
SIGAE
ratio total
of the free flow frontal frontal area, calculated
area of the evaporator to the in CALC and used in EVAPR.
EVAPOR
ss
specific printed
STATIC
static superheat setting valve, assigned a value used by TXV and printed
STC
vertical spacing between tube passes in the condenser, by DATAIN, used by CALC and COND (ft).
STE
vertical spacing between tube passes in the evaporator, by DATAIN, used by CALC, EVAP, and EVAPR (ft).
speed of the outdoor by OUTPUT.
STRTOT ratio of the sensible evaporator, calculated
fan , calculated
in FANFIT and
PRNT8
for, the thermostatic expansion TXVDAT in BLOCK DATA or read by DATAIN, by OUTPUT (F).
to total heat transfer rates for in EVAP, printed by OUTPUT.
read read the
CONDEN EVAPOR PRNT5
123
STRTP
of the
sensible
to .total
heat
two-phase region of the evaporator, printed in OUTPUT.
.
e
ratio
!?zS?z.,rate?.
calculated
.,.,for.. the. ._
factor for suction gas heating in the compressor, SUCFAC correction assigned a value in BLOCK DATA and used in CMPMAP. SUPER
specified superheat of the refrigerant leaving the evaporator read by DATAIN, used by COMP, EVPTR, OUTPUT, and TXV (F).
SUPERB base superheat value for and used by CMPMAP(F).
the compressor
PRNT5
in EVAP and
map, read by DATAIN
MAPFIT FLOWBA MAPFIT
superheat for the thermostatic expansion SUPRAT rated operating valve, assigned a value in BLOCK DATA or read by DATAIN, used by TXV, and printed by OUTPUT (F).
TXVDAT
SYNC
synchronous or actual compressor motor speed depending whether or not the full load motor power is calculated, read in DATAIN and used in COMP (rpm's).
COMPR
TACI
temperature of the air at the point where moisture removal begins, calculated and used in EVAP, printed in OUTPUT (F).
TAIC
temperature heat losses and printed
TAIE
air temperature crossing the evaporator coil, inlet air temperature with any additions due to fan motor and compressor heat losses, calcul.ated in EVAPR and printed in OUTPUT (F).
PRNT7
TAIIC
air temperature entering the condenser, read by DATAIN, used by COND and EFFCT, printed by OUTPUT (F).
CONDSR
TAIIE
air temperature entering the evaporator which corresponds to EVAPTR the current values of TSATEI and TSATEO, the desired value for the ambient air temperature is read by DATAIN, current value is calculated and used by the main program, also used by EVAPR, EVPTR, and OUTPUT, printed by OUTPUT (F).
TAOC
air temperature after it passes over the condenser coil, ulated in EXCH, used by COND, and printed in OUTPUT (F).
calc-
CONDSR
TAOE
air temperature leaving the evaporator with any additions for compressor and fan motor heat losses, calculated and used by EVAP, used by EVAPR, and printed by OUTPUT (F).
EVAPTR
on
.
. L
of the air crossing the condenser, includes from the fan motor or compressor, calculated in OUTPUT (F).
PRNT5
any PRNT3 in COND
124
TAOOC
air temperature leaving the condenser, includes from the compressor and fan motors, calculated printed in OUTPUT (F).
any heat losses in COND and
TAOOE
air temperature leaving the evaporator, includes any additions due to heat losses from the fan and compressor motors, calculated in EVAPR and printed in OUTPUT (F).
PRNT7
TAOSC
outlet air condenser,
temperature calculated
from the subcooled region of the in EXCH, printed in OUTPUT (F).
PRNT4
TAOSPC outlet air condenser,
temperature calculated
from the superheated region of the in EXCH, printed in OUTPUT (F).
PRNT4
TAOSPE outlet air evaporator,
temperature calculated
from the superheated region of the and used in EVAP, printed in OUTPUT (F).
PRNT3
PRNT5
TAOTPC outlet air temperature from the two-phase region of the condendser, calculated in EXCH, printed in OUTPUT (F).
PRNT4
TAOTPE outlet orator,
PRNT5
air temperature from the two-phase region of the evapcalculated and used in EVAP, printed in OUTPUT (F).
TC
critical temperature for the refrigerant, assigned TABLES, used in SATPRP, SPVOL, and VAPOR (R).
TERAT
evaporating temperature used for rating the thermostatic expansion valve, assigned a value in BLOCK DATA, used in TXV (F).
TFLODV refrigerant temperature entering the flow control ulated and used in FLOBAL, printed in OUTPUT (F).
a value
device,
in
SUPER TXVDAT
talc-
PRNT2
TFR
constant used to convert from Farenheit to absolute temperature when using correlations for refrigerant properties (i.e. the value assumed by researcher deriving the correlation), assigned a value in TABLES, used by SATPRP, SPVOL, TSAT, and VAPOR (F).
SUPER
TIC
temperature of the refrigerant entering the condenser, initialized by COMP, calculated in DESUPR and COMP, used by COND, DESUPR, EFFCT, and EXCH, printed by OUTPUT (F).
CONDSR
TIE
temperature calculated
EVAPTR
of the refrigerant by EVAPR and printed
entering the evaporator, by OUTPUT (F).
TLQCOR correction factor to adjust tube and nozzle pressure drops in the thermostatic expansion valve for liquid temperatures entering the valve other than the rating temperature, calculated and used in TXV, printed in OUTPUT.
TXVDAT
125
of the.liquid refrigerant TLQRAT rating temperature the thermostatic expansion valve, assigned DATA and used in TXV (F). I
TOLH
TOLS
tolerance properties enthalpy,
t
MPASS
of the refrigerant at the compressor shell inlet and used in COMPor CMPMAP, printed in OUTPUT (F).
of the refrigerant at the compressor shell in COMPor CMPMAP, printed in OUTPUT (F).
TXVDAT
and
temperature of the refrigerant at the condenser exit, estimated in FLOBAL, calculated in EXCH, used by CNDNSR, COND, FLOBAL, and FLODEV, printed in OUTPUT (F).
TROCMP temperature calculated
.
parameter used to calculate superheated vapor when the two known properties are the pressure assigned a value in BLOCK DATA _....~_--..-..._ .,.,f used by COMP.
tolerance parameter used to calculate. .)_, superheated vapor properties when the two known properties are pressure and entropy, assigned a value in BLOCK DATA and used by COMP.
TRICMP temperature calculated TROC
at the inlet to a value in BLOCK
outlet
MPASS
CMPRSR CONDSR
CMPRSR
TROE
temperature of the refrigerant leaving the evaporator, calculated and used in EVAP, used by EVAPR and EVPTR, printed by OUTPUT (F).
EVAPTR
TRVDS
refrigerant temperature at the end of the desuperheating region in the condenser, initialized in COND, calculated by EXCH, and used by COND, EFFCT, and EXCH (F).
CONDS
TSATCI
saturation temperature of the refrigerant entering the condenguess read by DATAIN, calculated in the main ser, initial program, used by COMP, FLOBAL, COND, EXCH, and DESUPR, printed by OUTPUT (F).
CONDSR
TSATCO saturation temperature of the refrigerant leaving the condenser calculated in COND and FLOBAL, used by CNDNSR, COND, EXCH, FLOBAL, and FLODEV, printed in OUTPUT (F).
CONDSR
TSATEI
the evapand used
EVAPTR
TSATEO saturation temperature of the refrigerant leaving the evaporator, initial guess read by DATAIN, calculated and used by the main program, used by COMP, EVAPR, EVPTR, FLOBAL, and TXV, printed by OUTPUT (F).
EVAPTR
TSATFL saturation temperature of the refrigerant entering the flow control device, calculated and used in FLOBAL, printed in OUTPUT (F).
PRNT2
saturation temperature of the refrigerant orator, initialized in the main program, by EVAPR and TXV, printed by OUTPUT (F).
entering calculated
126
TSATRl
average saturation temperature in the evaporator, value in EVAP, used in TWISOL (F).
assigned
a
TWI
TSATR2 average saturation temperature in the evaporator, value in EVAP, used by TWOSOL(F).
assigned
a
Two
* 7
TSICMP saturation temperature shell inlet calculated in OUTPUT (F).
of the refrigerant at the compressor and used in COMPor CMPMAP, printed
TSOCMP saturatin temperature of the refrigerant at the compressor shell outlet calculated in COMPor CMPMAP, printed in OUTPUT (F). TWALLI
CMPRSR
tube wall temperature at the leading edge of the evaporator where moisture removal begins, calculated and used in EVAP, used by TWOSOL, printed in OUTPUT (F).
TWALLO air temperature at the tube walls on the leaving evaporator, calculated and used in EVAP, printed
edge of the in OUTPUT (F).
TWBI
wet bulb temperature of the air entering ulated in EVAP, printed in OUTPUT (F).
TWBO
wet bulb temperature of the air at the leaving edge of the evaporator, calculated in EVAP, printed in OUTPUT (F).
the evaporator,
expansion TXVRAT rated capacity of the thermostatic a value in BLOCK DATA and used in TXV (tons).
valve,
CMPRSR
calc-
assigned
TWO
PRNT5 PRNT5 PRNT5 TXVDAT
UATPAW air-side thermal conductance times the heat transfer area (UA) for the portion of the two-phase region of the evaporator which is wetted due to dehumidification, calculated in EVAPR, and printed in OUTPUT (Btu/hr-F).
UAS
UATPW
refrigerant-side thermal conductance times the heat transfer area (UA) for that portion of the two-phase region of the evaporator whose outside surface is wetted due to dehumidificalculated in EVAPR and printed in OUTPUT (Btu/hr-F). cation,
UAS
ULA
upper limit refrigerant,
V
constant in dehumidification solution, and used in TAOSOL and TWOSOL.
on the Reynold's number for laminar flow of the VAPHTC assigned a value in BLOCK DATA and used in SPHTC.
calculated
in EVAP
TWO
127
I
factor for the compressor volumetric efficiency, VOLFAC correction assigned value in BLOCK DATA and used in CMPMAP.
-
VR
compressor clearance used by COMP.
W
constant in dehumidification solution, and used in TAOSOL and TWOSOL.
WAIRI
humidity ratio ator calculated
WAIRO
humidity ratio of the air and used in EVAP, printed
calculated
PRNTS
WOUT
humidity ratio of the 'air at the leaving edge of the evaporator, calculated and used by EVAP, printed by OUTPUT.
PRNT5
WTC
spacing of the tubes in the condenser in the direction flow, read by DATAIN, used by CALC and COND (ft).
CONDEN
WTE
spacing of the tubes in the evaporator in the direction air flow, read by DATAIN, used by EVAP and EVAPR (ft).
. L ..
volume ratio,
read by DATAIN and
calculated
MAPFIT COMPR
in EVAP
of the air at the leading edge of the evaporand used in EVAP, printed in OUTPUT. leaving the evaporator, in OUTPUT.
of air of the
TWO PRNT5
EVAPOR
WWALLI humidity ratio at the tube walls on the leading edge of the evaporator where moisture removal begins, calculated and used in EVAP, calculated in TWISOL, printed in OUTPUT.
TWI
WWALLO humidity ratio of the air at the tube walls on the leaving edge of the evaporator, calculated and used in EVAP, printed in OUTPUT.
PRNT!j
X
constant used in correlating enthalpy of refrigerant vapor with OTHER temperature and specific volume, assigned a value in TABLES and used in SATPRP.
of the refrigerant entering XFLODV vapor quality in FLOBAL and printed device, calculated
the flow control in OUTPUT.
PRNT2
XIC
refrigerant vapor quality entering the condenser, set to 1.0 in FLOBAL, modified in DESUPR, used by COND and EXCH, printed in OUTPUT.
CONDSR
XIE
refrigerant vapor quality entering the evaporator, and used in EVAPR and printed in OUTPUT.
EVAPTR
calculated
128
XINCMP refrigerant calculated
CMPRSR -
shell inlet vapor quality at the compressor and used in COMPor CMPMAP, printed in OUTPUT.
XKFC
thermal conductivity of the fins in the condenser, DATAIN, used by CALC and COND (Btu/hrFt-F).
XKFE
thermal conductivity of the fins in the evaporator, DATAIN, used by CALC and EVAP (Btu/hr-ft-F).
read by
CONDEN r
read by
EVAPOR
factor for the pressure drop in distributor tubes XLCORR correction of the thermostatic expansion valve for tubes of lengths other than the rated length, calculated and used in TXV, printed in OUTPUT.
TXVDAT
XLLR
lower limit refrigerant,
XLEQDL equivalent reversing COMP (ft).
on the Reynold's number for laminar flow of the assigned a value in BLOCK DATAand used in SPHTC. length of the compressor discharge line from the valve to the condenser, read by DATAIN -nd used by
length of the high pressure XLEQHP equivalent discharge port to the reversing valve, by COMP (ft).
VAPHTC LINES
line from the compressor read by DATAIN and used
LINES
length of the liquid line connecting the indoor and XLEQLL equivalent outdoor heat exchangers, read by DATAIN -nd used by FLOBAL (ft).
LINES
length of the low pressure line from the compressor XLEQLP equivalent suction port to the reversing valve, read by DATAIN and used by COMP (ft).
LINES
length of the suction line from the evaporator XLEQSL equivalent the reversing valve, read by DATAIN -nd used by COMP (ft).
LINES
to
XMAC
mass flow rate of air condenser, calculated in OUTPUT (lbm/hr).
across each parallel circuit in the in COND, used by EFFCT and EXCH, printed
XMAE
mass flow rate of airy across each parallel circuit in the calculated in EVAPR and used in EVAP, printed evaporator, in OUTPUT (lbm/hr).
EVAPS
XMASP
mass flow rate of the air across the superheated region of calculated and used in EVAP, each circuit in the evaporator, printed in OUTPUT (lbm/hr).
PRNT5
XMATP
mass flow rate of the air across the two-phase region of each circuit in the evaporator, calculated and used in EVAP, printed in OUTPUT (lbm/hr).
PRNT5
XMM
rate of moisture removal from the air, printed in OUTPUT (lbm/hr).
calculated
in EVAP and
CONDS
PRNT5
129
XMR
guess or estimate read by refrigerant mass flow rate, initial DATAIN, calculated in COMP, used by the main program, COND, EFFCT, EVAPR, FLOBAL, and the flow control device subroutines, printed and used by OUTPUT (lbm/hr).
FLOBAL
XMUA
viscosity (lbm/hr-ft).
AIR
of air,
calculated
and used by COND and EVAPR
XNTUSC number of transfer units for the subcooled region of the condenser, calculated in EXCH, printed in OUTPUT.
PRNT4
XNTUSP number of heat transfer units for the superheated region of PRNT4 the condenser, calculated and used in EXCH, printed in OUTPUT. XNTUSU number of heat transfer units for the superheated region of PRNT5 the evaporator, calculated and used in EVAP, printed in OUTPUT. XNTUTP number of heat transfer units for the two-phase region of PRNT4 the condenser, calculated and used in EXCH, printed in OUTPUT. XNTUTP number of heat transfer units for the two-phase region of PRNT5 the evaporator, calculated and used in EVAP, printed in OUTPUT. . c c
xoc
refrigerant vapor quality leaving the condenser, initialized CONDSR in COND, calculated in EXCH, used by CNDNSR, COND, EXCH, FLOBAL, and FLODEV, printed in OUTPUT.
XOE
vapor quality of the refrigerant leaving the evaporator, calculated and used in EVAPR, used in EVPTR and COMP, and printed in OUTPUT.
XOUCMP refrigerant calculated
Y
vapor quality at the compressor shell outlet in COMPor CMPMAPand printed in OUTPUT;
constant used in correlating with temperature and specific TABLES and used in SATPRP.
entropy of refrigerant vapor volume, assigned a value in
EVAPTR
CMPRSR
OTHER
/
,
131
Appendix
E - Cross-Reference
The common blocks data marks.
to or from
the
listed
subroutines
of Common Blocks in
the
rows of Table
identified
by the
1 are columns
used to transmit with
the
"X1'
132
SUBROUTINES
AIR Al BLANK CMPMOT CMPRSR COMPR CONDEN CONDS CONDSR DENSIT EFFARG EVAPOR EVAPRS EVAPS EVAPTR FANMOT FLOWBA LINES MAPFIT MPASS OTHER PRNTl PRNTZ PRNT3 PRNT4 PRNT5 PRNT7 PRNT8 REFRIG
I 1x1 I lxlxl
SAT
I
lxlxl
STATE4
I
I
SUPER TWI TWO TXVDAT UAS VAPHTC
I iii/iii/
I I
X r.
I
XII
IIII
I
I
I
I
I IXIXI I I lxlxl
lxlxl I I I I4 II II 1x1 I I Ixlxlxl Ixlxlxl I I 1x1 1 I
I
I
4
I
I
X
x
1,x,,,
I
1
,,,,,
X ,,,.
I
I
133
Appendix
Subroutine
F - Subroutine Calls
Descriptions
and Cross-Reference
of Subroutine
Descriptions
serves as driving program for high- and low-side computaMain Program tions and contains iterative loop converging on evaporator inlet air temperature assigns BLOCK DATA parameters
geometric
to constants
CAPTUB
submodel
CHTC
calculates two-phase
CMPMAP
computes refrigerant mass flow from compressor map data
CNDNSR
serves as ,driving routine for COMP, COND and FLOBAL and returns the difference between the calculated and specified condenser subcooling or between the compressor and expansion device refrigerant mass flow rates
COMP
computes sumption
COND
calculates total condenser heat transfer rate, air ant properties and refrigerant and air-side pressure fixed inlet refrigerant conditions
capillary
tube
refrigerant-side region of the
for flow
both
heat
changed
calculates for
constants
and infrequently
CALC
exchangers
rate
heat transfer condenser
coefficient
and compressor
for
the
power consumption
refrigerant mass flow rate and compressor from efficiency and loss parameters
drops
for
DPFSOL
calculates the difference between the known saturation pressure of water vapor at the tube wall and the saturation pressure corresponding to an estimated, or given tube wall temperature
DPLINE
determines lines
EFFCT
computes the difference in condenser effectiveness values between the general effectiveness equation and the specific cross-flow effectiveness equation as a function of the fraction of the coil, fV, containing superheated refrigerant vapor
EHTC
calculates two-phase
refrigerant-side region of the
pressure
drop
heat transfer evaporator
in
for
and refriger-
determines an approximate dew-point vapor pressure of moist air
the single-phase
temperature
power con-
DPF
h 5
values
the
coefficient
a given
refrigerant
for
the
134
EVAP
determines the heat transfer, moisture removal, and outlet air temperatures and humidities from one circuit of the evaporator for given heat transfer coefficients and saturation temperatures at the beginning and end of the two-phase region
EVAPR
calculates properties, fixed exit
EVPTR
serves as a driving between the specified value
EXCH
determines the heat transfer and outlet temperatures for one circuit of the condenser for given heat transfer coefficients and saturation temperatures at the beginning and end of the two-phase region
EXF
determines using the
FANFIT
calculates pressure
FLOBAL
determines refrigerant conditions at the inlet control device and drives the expansion device
to the models
flow
FRICT
computes the general flow in tubes
s ingle
phase
GUESS2
brackets a solution prior to using end points by factors of 10
a root
finder
by shifting
GUESS3
brackets a solution prior to using end points by constant step
a root
finder
by shifting
INTER
interpolate
HAIR
computes air-side tube geometry
MUKCP
calculates viscosity, of 13 refrigerants
thermal
conductivity,
and specific
heat
MUKCPA
calculates of air
thermal
conductivity,
and specific
heat
OR1FIC
computes orifice
OUTPUT
prints
evaporator heat transfer rate, and refrigerant and air-side refrigerant conditions routine for evaporator
air and refrigerant pressure drops for
EVAPR and returns the difference superheat and the calculated
the effectiveness of a cross-flow effectiveness - NTU method combined fan-fan drop and volumetric
Moody friction
in a single
viscosity, refrigerant a detailed
motor efficiency air-flow rate
factor
dimension
heat
transfer
mass flow
using
summary of output
through data
exchanger
given
air-side
for
Lagrangian
coefficient
rate
heat
for
polynomial smooth
a short-tube
fin
and
135
a
PDAIR
calculates air-side fin and tube heat
PDRQP
determines singleand two-phase heat exchanger tubes
PVSF
calculates
SATPRP
evaluates refrigerant
SEFF
determines surface
SLAG
computes single dimensions from
SPFHT
calculates refrigerant volume and specific
SPHTC
computes transition, entrance
single-phase heat transfer coefficient for laminar, or turbulent gas flow from an abrupt contraction
SPHTCZ
computes developed
single-phase heat transfer liquid or gas flow
SPVOL
evaluates
SUPCOR
calculates power and mass flow correction factors for based compressor model to correct for superheat level
TABLES
assigns routines
TAOSOL
computes the exit air temperature evaporator where moisture removal
TRIAL
determines thermodynamic properties vapor given two known properties
TSAT
pressure exchangers
partial the
pressure
saturation
surface
drop
for
volume
for
'calculates saturation tion pressure
temperature
shaped
coefficient
thermodynamic from the occurs
of superheated of refrigerant
TXV
computes expansion
at the rate
and
fully
refrigerant given
exit
sub-
of the
moisture
through
map-
property
region
used to compute evaporator
mass flow
two
refrigerant
TWOSOL
the refrigerant valve
fin
pressure
for
used to compute the wall temperature at which begins on the leading edge of the evaporator temperature
in air
in
at constant
TWISOL
the wall
flow
of a specified
interpolation
of superheated
refrigerant
for
properties
heats
or louvered
in saturated
a hexagonal
Lagrangian data
specific heat ratio
wavy,
drops
vapor
thermodynamic
precision tabulated
constants
smooth,
pressure
of water
efficiency
specific
for
from
saturaremoval the
a thermostatic
1 3-6
VAPOR
determines thermodynamic vapor given temperature
WBF
determines
WTSFIT
computes coefficients to enthalpy of moist
for air
XMOIST
calculates wet-bulb
temperature, or relative
ZERO's
each of these routines solves for the function from two points which bracket section and Newton's method
wet-bulb
the
subprograms
"X'S".
temperature
dew point temperature
The subroutines listed
listed
properties and pressure
of moist
a quadratic
in the columns
in the
of superheated
rows of Table
refrigerant
air
fit
humidity humidity root, the
of Table F.1 which
of wall
temperature
ratio enthalpy, of moist air
and
or zero, of a solution using bi-
F.1 are called are
indicated
from by
137
. . +,
139
Appendix
G - Algebraic
A
frontal
a
intercept context
A A
a contact
area
of heat
fin contact
area
fin
frontal
area
A ft
air-side
fin
A heater
cross-sectional
A min
ratio of total free-flow area
A orifice
cross-sectional
A
total
A
filter
r
A side
heat
A
total
A
0
b b
fat
heat
transfer
of filter
transfer exposed
bypass
heat
and tubes
short
area
of each side
transfer
fit
to data
area
section
to the minimum
orifice area
for
each parellel
circuit
of each fin
of tubes
heat
for
heater
tube
transfer
transfer
factor
area
transfer
heat
area
unit
resistance
of the
surface
by
area
transfer
heat
heat
or bleed
fins
of the
area
of a linear
transfer
is determined
area
air-side
area
which
on indoor
and tube
refrigerant-side slope
heat
refrigerant-side
refrigerant-side
to data
between
area
A rh
tubes
fit
and tube
air-side
Af
exchanger
of a linear
air-side total
Notation
area which
in heat
exchanger
per unit
air
flow
is determined
a thermostatic
rate
by context
expansion
valve
C
compressor flow area
C
heat
C
ratio
C
max
maximum of the refrigerant
heat
capacity
rates
for
the
air
and the
C min
minimum of the refrigerant
heat
capacity
rates
for
the
air
and the
a fat
piston clearance volume ratio coefficient or Dittus-Boetter
capacity
rate
of humidity
of air ratio
or expansion exponent
device
(macPa) gradient
to the
temperature
gradient
140 C r
heat
C
constant
in Eq. 6.15
C P C pa C pm
specific
heat
of the
specific
heat
of air
specific
heat
of moist
C r
heat
C
correction factor to the air-side for wavy and louvered fins
D
total compressor piston displacement diameter depending on context
d
depth
D a
air-side
or outer
D ducts,
diameter
of air
D r
refrigerant-side
f
fraction of rated value, determined by context
F a
fin heat area
F moist
fraction of that part of the evaporator which flow that is wetted due to dehumidification
F P F sh
fin
f
combined
0
tot
capacity
rate
capacity
of refrigerant
rate
of heat
(mrcp)
refrigerant
air
of the
refrigerant heat
transfer
coefficient
or refrigerant
tube
exchanger diameter
of heat
exchanger
tubes
ducts or inner
transfer
area
diameter
of heat
fraction per unit
of coil, total
exchanger
tubes
or friction
air-side
heat is
in
factor transfer two-phase
pitch
compressor
suction
gas heating
wet and dry
F
W
fraction of the dehumidification
f
wet
friction wetted
friction
evaporator
factor that
factor for the region due to dehumidification
F
compressor
Fktt)
correlating transfer
volumetric function coefficient
factor
efficiency for
the
for
is wetted of the
the
evaporator
due to
evaporator
correction refrigerant-side
which
is
factor two-phase
heat
141
correlating
r F* .
. *
parameter
mass flux
G
gravitational
SC
acceleration
constant
enthalpy
of refrigerant
enthalpy
of saturated
air
ha, a, h fg h' fg h
enthalpy
of moist
at freestream
h
average
h h
a,s
ave
heat
driving
tP
enthalpy
difference
air or refrigerant heat subscript or context)
sensible
h* W
effective
condensation
h2
transfer
heat
transfer
heat
transfer
heat
coefficient
transfer heat
for
coefficient transfer
of tubes
in heat
exchangers
length
of tubes
in
exchanges
l
I
surface
using
Chaddock-
local evaporation heat transfer coefficient for the dry-out region using nonlinear interpolation between the Chaddockvapor Noerager value at xdo and the value for saturated
length
m
by
a wet surface
coefficient
e
ril
a wetted
for
equivalent
tubes
(determined
coefficient
L
c
condenser
coefficient
thermal
exposed
temperature
conditions
in the
k
L
surface
coefficient
local evaporation Noerager correlation
%
exchanger
average condensing or evaporating heat transfer coefficient (condensation value includes an adjustment due to Rohsenow)
h
W
air
at wall
of heat
of vaporization
mass transfer
‘d h
latent
or height
conductivity length
equivalent valve
length
mass flow
rate
mass flow
rate
of refrigerant
heat
of distributor
at "standardized"
pipes
exposed tubes
test
in
to the
air
thermostatic
or rating
expansion
conditions
,
.
.-.
142 N cap N circuits
number of parallel
tubes
number of equivalent heat exchanger
parallel
total
number of fins
in the
N Le
Lewis
number
N Pr N Re
Prandtl
Nf
refrigerant heat
circuits
in the
exchanger
number
Reynolds
number
NT
number of tube the air flow
f
N tot
total
/
NtU
rows in the
number of tubes
number of heat
heat
in the
transfer
exchanger
heat
in the
direction
of
exchanger
units
N
number of tubes
P
pressure
6
heat
transfer
6 can
rate
of. heat
6 cooling
rate of heat transfer to the suction gas
from
the
compressor
and compressor
motor
'hilo
rate of heat transfer the suction gas
from
the
compressor
discharge
to
V
volumetric
I
capillary
stacked
vertically
in the
air
rejection
flow
from
the
compressor
r
ratio
of UA values
used in Eq. 6.11
ratio
of UA values
used in Eq. 6.23
S
motor
speed
S
heat
exchanger
S oper
actual
S
specific
S sync
synchoronous
resistance
coil
or operating
shell
line
rate
thermal
S
exchanger
rate
R
rl
heat
per unit
of refrigerant
side-area
characteristic motor
speed of the
fan
speed for
the
speed of compressor
outdoor motor
heat
exchanger
143
3
vertical
T
air
T
aid
T T
aie v,ds
T
sat
temperature
bulk refrigerant heating region refrigerant
air
W
humidity
ratio
i
work
wa, m
humidity
i
fan motor
dehumidification
the evaporator at the
start
of the desuper-
temperature
performed ratio
occupied
of air
full
of heat
input
at freestream
load
exchanger
vapor
X tt
correlating
parameter
c1
total air-side volume
cx can
fraction from the
a
correction factor for inlet to a thermostatic fraction discharge
exchanger
conditions
power consumption
compressor-motor spacing
by heat
or energy
refrigerant
ahilo
by context)
conductance
X
corr
entering
saturation
of space
a
exchanger
(determined
vapor temperature of the condenser
volume
wT
in heat
temperature
of the
V hx
s,fl
tubes
evaporator inlet air temperature for which would occur at the finned tube surface
thermal
fan
between
or refrigerant
U
I;
.
spacing
shaft tubes
power in
the
direction
for
two-phase
transfer
area
heat per
transfer
unit
heat
of the compressor power consumption shell to the ambient air liquid refrigerant expansion valve
of compressor power consumption line to the suction line
r
total refrigerant-side exchanger volume
a
racks
correction factor number of resistant
flow
quality
heat
01
of air
heat
transfer
to adjust air-side heater racks
area
coefficient exchanger
which
is
rejected
temperatures at the other than 100°F transferred per unit
pressure
drop
from heat for
the
the
144
Bcorr
correction 30 inches
Y
ratio of specific heat at constant volume
6
fin
Ah
enthalpy
change
Af’
pressure
drop
pressure
drop
APJ3 0 2
factor
for
TXV distributor at constant
measured
in inches
level temperature difference, superheat depending on context
E
heat
E cf
effectiveness
n
compressor effectiveness
1-I
viscosity
V
specific
P
density
0
capillary
tube
free-flow exchanger
frontal -
a
pressure
lengths
other
to specific
than heat
thickness
AT
u
tube
exchanger
of water
of subcooling,
or level
of
effectiveness of a counter-flow
efficiency
or heat
heat
exchanger
exchanger
overall
surface
volume
flow
factor area
per unit
total
frontal
where
there
area
of heat
Subscripts a
air
actual
at modeled
ave
average
cm
compressor
d
dry
do
dry-out liquid
dp
dew point
operating
conditions
point in the evaporator film on the tube wall
is
no longer
a
145
eff
effective
evap
evaporator
f
fin
fl
compressor
frict
friction
ft
fin
i
heat
in
entering
inlet
compressor
isen
isentropic
R
liquid
lot
local
m
mean
map
from
max
maximum
mech
mechanical
mom
momentum
noz
TXV nozzle
0
heat
oper
at modeled
out
exiting
outlet
compressor
r
refrigerant
rated
at rating
rb
return
S
compressor
motor
full
load
and tube exchanger
inlet
shell
compressor
exchanger
inlet
map
outlet
operating
shell
conditions bends shaft
outlet
conditions
,
.
-
-,
,.~
._.
-^,“l
._..-
“__l^_._-“...
146 sat
saturation
sh
suction
static
static static
t
tube
tP
two-phase
tubes
heat
TXV
thermostatic
V
vapor
vol
volumetric
w
region
co
freestream
gas heating
L
fan efficiency or static expansion valve
exchanger
or TXV distributor expansion
wetted
valve
due to dehumidification
condition
superheat
tubes
value
of a thermo-
147 H - Derivation
Appendix H.l
of Dehumidification
Introduction A method
is
described
of an air-to-refrigerant The method
represents
Goodman [H.l], Threlkeld
[H.4],
H.2
Heat Transfer
be described,
evaporator
under
and Wiley from
as noted
The local
section
and McQuiston
derived
[H.7]
in this
an assimilation
McElgin
coefficients Myers
Solution
[H.2],
heat
in Section Equations
driving using
force
computing
the
dehumidification
and/or
[H.5,
dry
for
conditions.
expansion
Hiller H.61,
performance
of the work
and Glicksman
[H.3],
and uses wet heat
transfer
values
via
the
of
transfer
equation
of
6.
-Air
to Wall
for
simultaneous
simplifying
heat
assumptions,
and mass transfer
by Eq. H.l
can
[H.6]
db = h, (h a, to -ha , ,JdAs where a differential
quantity
of
local
dQ
-
‘d
-
ha, m
-
the
-
the enthalpy perature
of saturated
air
-
a differential
element
of wetted
h
a,s
dA Equation
S
H.l
the mass transfer
contains
enthalpy
heat
and mass transfer
coefficient
of humid
air
the assumption
at the
that
the
freestream
at the wall surface Colburn
condition surface
tem-
area. analogy,
given
by
Eq. H.2 CH.61,
h
W
’ pm ‘d
=
213
N
Le
(H-2)
148
where
hW
-
sensible
c
-
moist
heat air
transfer
specific
coefficient
heat,
c
pm
woo a, N Le holds is
humidity
-
Lewis
-
and that
neglected
(except wetted
and Elmahdy
assumptions; the
effect
transfer
rates
parallel
flat
equation
presented
the
+ .444 wa,co
pa
energy
content
film
thickness
the water sensible
heat
are analyses
in the
[~.8]
which
require
these heat
is
for
surface
and
There
simplified H.l
of
on the
however,
Equation
Also,
c: 1.
effect
surface).
[H.4] with
(MLe)2'3
the
of air,
= c
a wetted
number,
and the
for
ratio
pm
for
do not analyses
exchanger
strictly
plates.
such as humid air
by McQuiston
[H.5,
H.61
for
a
for
use
by Threlkeld of these
heat
pump model.
describing
heat
ignored
too complicated
in the for
condensate
coefficient
some or all
routines
For a fin-and-tube
is
literature
were judged
applicable
situations
transfer
of the
heat
flowing
and mass
between
exchanger,
two
an appropriate
is
di = hd nw,ft(ha,m- h,)dA,, where
the
subscripts
surface,
"t"
respectively.
is
subscript
a wet fin
similar meter
"m", m
fin
wet surface
= 1 --AAf ft
"f".refers
efficiency
to usual
to tube
and total
fin-and-tube
effectiveness,
nw ft, ,
as
nw,ft
the
refer
The overall
in Eq. H.3 is defined
where
and "ft"
(H-3)
which efficiency
where ,
to fin is
(1 - 'Iw,f)
surface.
defined equations
The term
by McQuiston except
that
nw f in Eq. H.4
[H.5]'in the
usual
a manner para-
149
is
replaced
by M, i.e.,
M = +)(l
+
C h fat fg)-jl/2 c pm
(H.5)
where h
h
W
h*
W
kf
heat
transfer
coefficient
for
a dry
sensible
heat
transfer
coefficient
for
a wet surface
effective
heat
thermal
conductivity
fin
6 h
sensible
coefficient of the
for
a wet surface
fin,
thickness
enthalpy
fg
transfer
surface
of vaporization
of the water
vapor,
and ratio of humidity ratio gradient to temperature (specifically defined in Appendix I).
C fat
C,__ in Eq. H.5 is
Since heat
gradients,
the
dehumidification H.3
Equation portion
tube
is
wall
-A ir
Equations
ht --
ratio
of
efficiency
latent
will
external
heat
when both
temperature).
to Mean Surface the
Since
temperature.
the mean temperature
surface
of wet fin
value
H.3 contains
of the
in determining
on the
to sensible
be dependent
on-the
rate.
Heat Transfer
at the
dependent
gradient,
transfer
the fin
enthalpy
of saturated
the
surface
tube
surface,
and tube
a better
start
of the fin-and-tube From Threlkeld
Conditions
[H.4],
is
air only
it
is
a small
variable
to condense
combination
evaluated
(the
shown that
for
water
use
vapor
effective for
a dry
surface
'd,f
T , = T:,,
- Tf,m - Tt
b-4.6)
150
where
the
subscript
wet surface, T
a,
03, it
"f,m"
since
follows
refers
McQuiston's
to a mean fin solution
for
Similarly,
value.
nm f is ,
in terms
for
a
of Tt and
that
T
rlw,f
a,03 - Tf,m
= T
(H-7)
a,03 - Tt
and T
a,a =T /Tt a,
nw,ft
Next,
using
a relationship
implied
-T
ft,m '
in the derivation
of nw f, ,
i.e.,
(H.9)
ha, co- hft ,m = (’ + Cfac)(Ta,az- Tft,m) the
overall
wet surface
rlw,ft
il
Substitution
of Eq.H.10
into
db
H.4
Heat Transfer Next,
and tube
=
hd(ha
Equations
a relationship temperature,
nw ft ,
effectiveness
is Tft
,
=
can also
-h ft,m a,m h co -ht a,
be written
as
h
(H.10)
'
Eq. H.3 yields
,
co
-
hft
,
(H.11)
l
m)dAft
-Mean
Surface
to Refrigerant
needed
between
conditions
m, and the
refrigerant
Conditions at the mean fin
temperature,
T,.
151
Using
McQuiston's
written
definition
in two ways,
of nw , ft,
,
do = h*w(Ta,m
of
Eqs.
total
heat
flow
can be
i.e.,
do = h*w 11, ft(Ta,m
Combination
the
- Tt)dAft
(H.12)
(H.13)
- Tft , m)dAft
H.12 and H.13 to eliminate
- VdAf
dQ = h*w "y,ft(Tft,m
T,,;
yields
(H.14)
t
- nw,ft The total
heat
flow
from
the
tube
dQ =
to the
refrigerant
is
given
by (H.15)
hr(Tt - Tr)dAr
t
Using
Eqs.
H.14 and H.15. to eliminate
Tt yields
(H.16)
M.5
Derivation Equations
ship
between
known as the
of the
Coil
Characteristic
H.16 and H.11 are next [ha,m coil
- hft
,]
characteristic,
and [Tft
and Total combined
m - T,].
s,'and
is
Heat Flow Equation
to arrive at a relationThis relationship is
given
here
by
(H.17) .
where
hd in Eq. H.11 has been replaced
by
hw/cpm.
152 Assuming
that
dAft dA=K-r
Aft r .
and using
the
relationship
C
h*
W
= h,(l + c fat
hfg)
,
pm
Eq. H.17 can be written
s
= Tft,m h a, m
At this
point,
we have a local
one cross-sectional gives
the
change
element
relation
in energy
content
i,
is
the
in Tft
air
local
of the
air
= -h,(h,,,
the mass flow
Equations equation
a,03
relationship
in the
between
r;ladh
where
as
flow
stream,
,
can then
from T,,,
Eq. H.17 and H.18 to solve for Tft 3m (x0) follows, Ta,03 ('0) is given.
after
flow
and the
is,
(H.19)
to give
be integrated
to Tft for
heat
H.19
,)dAft
th& exit of the heat exchanger ft,m at '0' evaluated, Eq. H.17 can be used to evaluate equation
Equation
of air.
T
transfer
that
H.18 and H.19 can be combined m which
ha,m and Tr at
direction.
cross-sectional
- hft
rate
between
T,,: which
coil. h
m can alsoaLe (x0). the
a differential
to solve
for
With co (x ).
Tft m (x0) A sinsible
wrytten
The derivation derivation
of the
the
value
of
heat
and used with of the
solution
solution
for
153
H.6
Derivation T ft,m ('0) First,
of the
Eq. H.18
dT
For a direct
Solution
is written
saturation
ventional
line
=
r
=
of T ft
since
is
and Wiley
calculations second
conditions if
instead,
are broken section hft
,
first
m is
Surface
form,
Temperature,
i.e.,
(H-20)
- dhft,m)
(H.21)
to - dhft,,)
m and T,,
m (between
chart)'is
into
related
ft,m
a* + b*Tft
is only
,
h and T on
required.
The con-
m
valid
has typically two sections of This
to Tf,
,
[H.3]) air
and "b"
sectionalization
,
small
where flow
derived
m by a quadratic
= a t bTft
a very
been recommended in the
"a"
(H.22)
within
and Glicksman
section.
m + CT& m ,
range
(see the
coil
direction from
and
the
exit
can be avoided expression,
i.e.,
(H.23)
of Eq. H.23 yields
dh ?
=
and Hiller
h
Differentiation
ft,m
uses new values
of the
Effective
that
solution
[H.2]
a,
hft
relationship
m, a sequential
McElgih the
this
s (dh
psychromet;ic
h
However,
S(dha,co
between
of the
assumption
Exit
dT, = 0, and Eq. H.20 becomes
ft,m
a relationship
the
in differential
coil,
dT
the
-dT
ft,m
expansion
Next,
for
ft,m
= dT
ft , mtb + 2CTft,m)
(H.24)
154 Using
Eq. H.24,
Eq. H.21 becomes . dT
Solving
=
ft,m
Eq. H.25 for
dh
a,
s [dh
H.19 is
replaced
by [Tft
,
now solved
for
m - TJ/s,
i.e.,
dh
Equating
Eqs.
where
9
m , )I
(H.25)
* 1 *
=-
,
+ b + 2cTft
,+-
dh a,oo and the
m) 9
term
[ha,m
(~26)
- hft , m1 is
-hd Oft m - Tr) dAft fi
(H.27)
S
a
H.26 and H.27 yields
(-+-I S Collecting
a,O"
o. - dTf m(b + 2cTft
o. yields,
dha,m = dTft
Equation
a,
like
b + 2cT ft,m)dTft,m
terms
on each side
= .$s a of
(Tft2;
the equation
-Tr)dAft
(~.28)
yields
u = 1 + bs and v = 2cs.
Integration are
the
of Eq. H.29 from
beginning
Tft
m (xi)
and end of the:wet
coil
to Tft
m (x ), where x. and x 1 0 0 regi& in the air flow direction,
yields
I
(u + vTr) Rn[Tf, , m(~) - T,l +
v(Tft
,
,(x)
- Tr)
I
T T
ft,mcxo)
= -hclAft
Fmoist ma
ft,mO(i) (H.30)
*
155
where
Fmoist
is
the
value
of Tft
m (xi)
fraction is
of the
calculated
coil
from
depth
that
is wetted.
Eq. H.18 and Fmoist
The
is
known from
v
the
of the
calculat!on
point
at which
dehumidification
first
occurs
L
I
(described T ft
in Appendix
m (x0).
of the
Rearrangement
equation
(u + vi,)
for
Ln[Tft,m(xo)
Eq. H.31 Tft the
is
the
transfer
term
knowns on one side
TJ
hdAft
-
(H.31)
Fmoist
6 a
it
must be solved
Eq. H.18 is
to calculate
flow
from
the
iteratively
h
air
used in conjunction
o3 (x ).
toaihe
With
h
rzfrigerant
a3 (x0) zin
be
the
the
ratio
with,
the dry
air
air
coil
heat
bulb
of
The solution bulb
transfer for
exit
either the
on,
do not rates
in
the
used in
the
and is preceding
air
total term
for
are not temperature
to obtain
Heat Pump Model derived
in
equations
the
region
ratio bulb
required
is
heat
a separation
the wet coil dry
(H.32)
region.
provide
and humidity is
the
second
in the wet coil
coil
temperature
and dependent
point
represents
and the
temperature
solution
dry
of Eq. H.31
section
to this
dry
at the
side
and latent)
equations
information. exit
right
(sensible
An additional
or humidity
on the
and latent exit
1 + inaCha,m(xi) - ha,,‘Xo’ 1
- T,,,(xI)
in
transfer
sensible
and thus
sistent
heat
(sensible)
the
of the
-
equation
has been evaluated,
~,[T,,,(n)
However,
extra
- T,l + v[Tft , m(xi)
routines
total
first
heat
known.
all
from
a=
of
in Eq. H.30 is
Trl + vDf t m(Xo) - Trl = ,
a transcendental
m (x0)
the
computed
total
unknown
of Eq. H.30 to collect
-
psy&rometric
evaluated,
where
only
m (x0)-
Onte Tft with
the
gives
(U + vTr) ln[Tft,m(xi)
Since
Therefore,
I).
is
a manner
this in terms con-
as follows.
'
156
H.7
Derivation T a,* ('0)
of the Solution
The sensible
heat
at a cross-section
for
transfer
in the
the
from
air
flow
Exit
the
Air
air
Dry Bulb
to the
direction
Temperature,
fin-and-tube
can be written
surface
as
dQ, = hwnw , ftcTa,m - Tt)dAft or,
alternatively,
using
Eq. H.7,
(H.33)
as
do, = $,$Ta,m - Tft ,,)dAft Equations
H.33 and H.34 contain
efficiency
for
transfer.
It
sensible
efficiency
for
for
total
transfer.
heat
slightly
to the
total
numerical
heat analysis
the
is
the
transfer
sensible
is heat
as well
local
sensible
heat
[H.9]
slightly
for
transfer to avoid
heat
the wet than
that
H.33 and H.34 rate
effect
total
that
lower Eqs.
this
and here
the wet fin
same as that
results,
However,
rate.
that
and Turner
[H.5]
in proportion
is the
ignored necessity
in of a
solution. change
the
in sensible
energy
macpdT
the
heat
assumption
Based on their
transfer
By equating to the
sensible
overestimate
McQuiston's
transfer
has been shown by O'Brien
fin will
heat
the
(H.34)
basic
governing
a,
equation
content
o3 = -h,(T,
for
the
,
gain
at the
of the
o. - T
air
a,
a,
=
-h,(T
a,
cQ- Tft
m ,)
stream,
sensible
surface
i.e.,
(H.35)
ft,mjdAft
transfer
When transforming Eq. H.35 into a integrable form, . . eliminate the term "dAftI' from Eq. H.35, i.e.,
macpdT
fin-and-tube
rate
is
obtained.
Eq. H.29 is used to
(~.36)
157
Cancellation
of terms
dTa,coz
Since
Tr is
yields
(Ta,m - Tft,m) - T (T ft,m r)
('
+ VTft,m)
Eq. H.37 can be rewritten
a constant,
dTft,m
(H.37)
in terms
of
(T - Tr)'s
to give
d(T, , a3- Tr) d(T
ft,m
-T
r
(Tap2 - Tr) Cu + v(Tft ,m - Tr) + vTrl )
T
-
ft,m
-T
r
-[u + v(Tft ,m - Tr) + vTrl Equation
H.38
is of the
=
(H.38)
form
+j$ + YPW = Q(z)
(H.39)
where '
-
Tft , meTr
p(z) -
‘(U + VTr + vz)/z
Q(z) -
‘(U + VTr + vz)
The solution
to Eq. H.39 can be written
yefdz
=pefdzdz
as:
t constant
(H.40)
158 Evaluating
the
innermost
1
Pdz
=
-
integral,
(u + vTr + vz)dz 1
1 Pdz =
- [(u
-
, and
2 + vTr)ln(z)+vt]
(H.41)
Then
z- (u + vT,),vz
Next,
the
outer
integral
(H.42)
becomes
(H.43) ./
1J Qe
PdZdz = -
J
(w + vz)z-'
eBVZ dz
(H.44)
159
Substitution the
of Eq. H.45
into
Eq. H.44 yields
an integrable
form,
with
result
f Qe
f
Pdz
(
"z(3'w) 3-w
wz(w
v
2-w
V2
2!
(
V3
3-r
Equation
i
H.46 can be written
+
wz(3'w) 3-w
wz(4-w) 4-w
+
1
+
vz(4-w) 4-w
+
vz(w 5 -w
>
-
+
(~.46)
--a
)
in summation
form
as
(H.47) i=l where the
'In"
series
is
the
number of terms which
expansion
solution.
gives
Substituting
acceptable
convergence
of
Eq. H.47 and H.42 into
Eq.
H.40 gives
y(z-’ eevz) =
t
,,C(i+U - WI (i
+ 1) -w)
+ C
(H.48)
i=l Equation
H.48 can be written
YF#)
in shorthand
= F&z) + C
form
as
(H.49)
160 The known conditions used to evaluate
at the the
beginning
constant
F1(zi)
in Eq. H.49 yields
z and y back to their
Ta,co(Xo)- Tr = fd
1
ft,m
0
the
exTiaairoW
Wa,,(xo)
=
The sensible the
coil,
values
is
a,m pa,&,,
heat
calculated
m , are obtained
(H.50)
for
exit
'.
conditions,
(H.51)
- F2(zi)l
yields
-Trl
+ *
- ‘r) - ‘2(‘Ft , m(x i ) from from
Eq. H.52,
the
-T,)>
(H.52)
humidity
ratio
' L?
of
the equation
+ 0.444*Ta,,.Jx
rate
4, and the rate
from
Eqs. H.54 and H.55,
and Wa,-
the
F2 CTft,m('o)
- cpTa ,(xo)]/[1061.2 , transfer
of TW a,m (x0)
Fl(Zi)
. i
representations
) - Tr'
(x ) has been evaluated
Once T
Yi
are
i.e.,
the solution
original
CT,, (Xi) - ‘,I ‘1(‘ft,m(~i’
region
- F2(zi)
(‘0) ’ y.=,&--r [FP Expanding
dehumidification
C in Eq. H.49,
c = yi
Use of Eq. H.50
of the
of water using
0
,]
(H.53)
condensation the
on
computed
(x0):
4
(H.54)
and
y
-5 ri
W
=1;1 a [wa,co'xi)
- wa,co(xo3
(H.55)
2
161
The calculations
.
heat
pump model
in a two-phase refrigerant of the
only
transfer References
for the
H.2
air
rates
and exit
air
for
Appendix
H
Engineering,
in
this
region
of
a fraction
in a two-phase
W. Goodman,
H.l
the
If
state.
vapor, coil
described
flow
the
coil
of the
rate
ia
for
state,
appendix
where
coil
the
use in in
the
the
of Air with p. 225-274.
and L. R. Glicksman,
No. 24525-96, of Technology
heat
Ref'igemrting Surface
Areas
Piping,
Improving Heat Pwnp Performance
- AnuZysis and Test, Laboratory, Massachusetts
Control
Heat Transfer (1976).
fraction
region.
"Calculation of Coil McElgin and D. C. Wiley, Air Cooling and Dehumidification," Heating, March, 1940, pp. 195-201. and Air Conditioning, C. C. Hiller
is
of the
Coils,"
J.
via Compressor Capacity
, the
calculati%
two-phase
the
superheated
by f
for
H.3
in
refrigerant
contains
is multiplied
conditions
"Dehumidification October, 1936,
are performed
Report
Institute
. l
L F
H.4
J.
H.5
F. C. McQuiston, "Fin Efficiency Transfer," ASHRAETransactions,
H.6
F. C. McQuiston
L. Threlkeld, Thermal Environmental Engineering, Hall, Inc., Englewood Cliffs, New Jersey, 1970.
and
and Air Conditioning
J.
Prentice-
with Combined Heat and Mass Vol. 81, Pt. 1, 1975.
Heating, Ventilating, and Design, Wiley, New York,
D. Parker,
Analysis
1977. H
f
l
7
of Dehumidification on the Air Side R. J. Myers, "The Effect Heat Transfer Coefficient for a Finned Tube Coil," Master's Thesis, University of Minnesota, 1967.
H.8
A. H. Elmahdy and R. C. Biggs, "Performance Simulation of Multi-Row Dry (and/or Wet) Heat Exchangers", Sixth IntemationaZ Heat Tpansfep Conference, Toronto, Aug 7-11, 1978, pp. 327332.
H.9
N. G. O'Brien and R. L. Turner, "Fin Simultaneous Heat and Mass Transfer," 1965, pp. 546-548.
Thermal
Efficiency
AICHE Joumal.,
During May
163
Appendix
I - Description
of the
The calculational the in
evaporator this
Appendix
(as shown in
The solutions
appendix. surface
respectively,
at H form
Algorithm
scheme used in the moisture
routine
fin-and-tube
Dehumidification
for
temperature
the the
trailing basis
Fig.
6.4
and air'dry edge of the
for
of the
m (x0)
Tft
the
removal text)
is
of
described
and Ta,,(xo)
bulb
(the
mean
temperature,
evaporator)
calculational
section
derived
procedure
in
discussed
here. A flow
1.1.
diagram
The algorithm
dehumidification exit
of the dehumidification is
would
air
temperature
calculated
assuming
occurs,
is
Taid
used only occur
T
compared
the
at the
a,m,senS only sensible
then
if
algorithm
(x0)
air
is
temperature,
fin-tube
surface
from
two-phase
Taid, is this
evaporator
inlet
air
is greater than Taid, moisture If Taie aie' the leading edge of the coil and the fraction
removal
than
region
If
T
at which
greater
transfer.
heat
to the
the
shown in Fig.
that
coil
is
condition temperature,
does not
of the
the
depth
occur
on
which
1
has only
*-
sensible
heat
transfer,
FSenS,
is
computed
-
TrJ] 'NtU,tp
from
the
equation
c w
F
sens
= en [ (Taie - T,)(T,,,
(1.1)
I
Equation Cmin'C max
I.1
is obtained
from
the
standard
effectiveness
equation
for
= 0, i.e.,
%
= 1 - e -N tu,tP
(1.2)
by using
T E=
T
ai,e aie
-T -T
aid r
0.3)
164 ORNL-DWG
81-4746~
Qa
t1
COMPUTE
Cfac (01,
7/w,ft
5 K%~/,
f+(O),
s(q)
CALCULATE
s(O)
ESTIMATE
Tft ,(x0),
COMPUTE
pyo AND
I
SET
L IEIRkJCH = 0
, 4
,
I
FMOIST
ha,a, (x0),
k;
= Sd
>
COMPUTE COEFFICIENTS FOR i, FIT -a +b Tft m+c T2 hf+,m , ft,m I , \ SOLVE FOR 7 c
Tft,m “J ITERATIVELY
I
NO
COMPUTE
hft ,Cx,‘,
I
ho ,(x,J
T~,oD(xb)~Wa,oD(:o),ct~~(~ol,
I
Fig.
I.1
Block
Diagram
of the
Dehumidification
Algorithm
4
165 The value
and rep1aci w NtU tp by Fsens N,,, tp. 1
at the
temperature,
Tf,
equal
dew'point
to the
i,
transition
point
temperature
of the
from
of the dry
entering
effective
surface
to wet coil,
air,
xi,
is
i.e.,
I
T ft,m('i)
-3.3
Values is
for
Cfa,(xi),
defined
r~,~~(x~),
by the
equation
(1.4)
= Tdp 9a , in
S(Xi)
and
are computed
next.
Cfac(X)
[I.11
wa,b) - Wft,b4 'faccx)
It
follows
Wa.,(xi)
from
Eq. I.4
= Wftam(xi)
appropriate
that
and,
to-use
is just
replace
h at xi,
the
because
of an increase
occur
in the where
predicted place
from
characteristic
from
Eq. H.5 that
transfer
coefficient.
was predicted Thus,
s(xi). W
Note that
It
is is
hw at xi Were
not since
into
temperature
the
Cfac(xi)
coil
due to the
A situation the
dry,
hw used to
further
effectiveness since
coil
when h was used,
h is used at
of fin
the
= 0.0.
shift
heat
calculation h*
would
surface
removal
of
point.
effective
when Fz~ was used.
coil
at this point
l
coefficient
in the
air-side
of kw in the
part
Cfac(xi)
transfer
beginning
(1.5)
' ft,JXJ
first
Eq. 1.5,
transition
moisture
- T
(x)
when the
the wet heat
dehumidification
increase
= T ' a,-
could
but
not
transition
point
nw ft(xi)
and the
='0,
it
in
follows
Therefore,
= h
W’
T.q,ftcXi)
= ~d,ft(Xi)
(1.6)
and
s (Xi)
at the
transition
point.
=
Sd
(1.7)
166 One feature
of this
procedure
needs to be explained
before
proI
ceeding
to discuss In the
edge e
in Appendix
the
case where
derivation
H, the
of the
integrations
dehumidification
solutions over
for
the
occurs Tft
coil
at the
m (x0)
depth,
leading d
and Ta,,(xo)
represented
by the .v-
variable (and is
"x",
were carried
from
Eq. H.16,
handled
by using
region
of
and exit f(xO)
h*w,
average
integration.
entrance
Returning
conditions
to Fig.
occurs
values
the 1.1,
on the
leading
values
entering
air.
iteratively
using
evaluated
at leading
(at
the
by f(O),
the entrance
f(xi),
and to calculate
is greater
the
and local
Tale,
In this
dehumidifi-
case,
the
m(O) is below the dew point The vilue of Tft m(O) must be found equation'from
edge conditions,
Appendix
l-l * *
i.e.,
1
a, m(O) - hft ,m(O) = nw, ,+$
than
T ft
characteristic
- Tr
simplified
parimeters
edge of the coil.
of the
which
assumption
and s over
as denoted
values
when Taid
temperature
h
That
s
more exactly.
surface
ft,nl(O)
of these
average
fin-tube
T
characteristic
was constant.
region)
at the
coil
the coil
&f hw, Jz*~, nw ft,
temperature
the
that
and nw ft)
removal
used to obtain and exit
assuming
Local
of the moisture
are
cation
h,,
out
hw(0)
w f&O)
i~~-~';c,~~tO), ,
Aft
+h c A r pm r
(1.8)
,'
is
s(O) = sw(o>
(I.9
For the
first evaluation of Tft m(OL s(O) , a-. and h (0) in Eq. I.8 are not known where f I3 and hfg(0) is given by
hf,(0)
= 1061.2 + 0.444 Tf,
is not known because C fat(O) C is given by Eq. I.5 facto)
, I?)O)
+ h&N
(1.10)
and T ft,m-32 hQ(0)
=
I0
, if Tft,m 1
32
,
if
T
ft,m
< 32
(1.11)
*
167
The initial
evaluation
of Tft
,
m(O) is made assuming
s(o)
Equation
I.9
then
used with
(0)] 7-l. ft m(O) to Tft m
to solve
ini&ial are is
is
value
calculated until
iteration
The value is heat
If
transfer
caused
the
entering the
dry
Cfac(o) an error
smaller is
humidity
air.
the
enthalpy
of the
found
fraction air
dehum idification
ft,mtx) conditions
coil
transition
occurs
at the
for
the
start
of the
The coefficients
are determined
CTft
and an equation
,
m(xi)
h ft,m(~)
- 51,
= 0.24Tft,,(x)
point
point, leading
from
using
T
ft,m('i" relating hft
+ Wf,,,(x)
region
LlC61.2
,
If
of calculations, stopped. coil
or
to wet conand the
= 1 and xi
= 0.
as a function
once the
of
initial
have been obtained.
an estimate m to Tft
by
are evaluated.
,(x)
H) can be estimat;d dehumidification
replaced
FmoiSt,
,(x,),
of hft
has
of the
dry
edge~'fmoist fit
coil
edge of the
is wetted, h
a dry
repeated. is
leading
that
to surface
s(O)]
series
execution
that
quadratic
(as assumed in Appendix at the
depth
one
is
m(O) is
seco;d
at the
from
from
dew point
program
process
freestream
switch
of Tft
This
positive,
above the
transition
at the
The coefficients T
and the
of the
is
[in
end of the
at the
it
the wet coefficient
m(O) has been evaluated
) has'been
the
,(O)
coefficient
calculation
message is printed
and s(O)
tilerance.
from
to zero,
hf,(0),
of Tft
relate
Once an
computed.
difference
to rise
at the
j:
to see if
to a wet coil
positive
is
[to
m(O)a
checked
is equal
and the
Tft
a prescribed
such has occurred,
coefficient is not
ft,m"' in the value
then
temperature
If
for
T
than
ratio
Cfa,(0)
surface
ditions, If
is
coefficient
Once Tft
Tft,m (xi
next
the
positive.
for
routines
found, Cfac(
m(O) has been
of Cfac(0)
to see if
psychrometric
iteratively
the difference
to the
(IJ2)
Sd
the
and a Aew value
repeated
is,
T ft
f&
=
of T
ft rn"o' m and Wft ,,'i.e., , J
+ O*444Tft,m(x)]
(1.13)
168 The saturation Wft
curve
on the
m as a function
by i quadratic derived
of Tft
equation;
from T ft
,
standard m.
This
hiwever,
m(xi)
and Tft
,
Wft,m(x)
-' ASHRAE psychrometric
curve
is most closely
we will
assume that
m(xO) will
be adequate,
=
chart
represents
approximated
a linear
equation
. @. a' + b'Tft,,(x)
(1.14)
-
(1.15)
b" Tft,m(Xj,)
and
b'=
Substitution
[Wft,mo(i)
- Wft,m(Xo)l'CTft
of Eq. I.14
hft,mb)
into
,
Eq. I.13
m('i)
- Tft,m(Xo)l
(1.16)
yields
= a + b Tft,mb)
+ c T:t,mb)
(1.17)
where
Note that
the
a
= 1061.2
b
= 0.24
C
= 0.444
assumption h
ft,mtx)
d
i.e.,
where
a' =W ft,m(xi)
.
a'
+ 1061.2
(1.18) b'
+ 0.444
b'
a'
(1.19) (1.20)
of = a* + b*T
ft,mcx)
(1.21)
169
made in simpler
analyses
implies
1
that
wft,m(x)
(1.22)
= a*
t 4 Once estimates culated, H is
the
of the
coefficients
transcendental
evaluated
solution a.., .,.."~," ^ .
iteratively
until
on T
Tf,
is
to obtain
h ft
a value
The exit the
values
humidity
air
dry
of the
‘a ,m(‘c)) = Values using
averages calculate until
from
new value
and c
b, until
more
convergence
Tr' After
of s(xo)
routines
[s(O)
or sd]
of Eq. H.17,
exit
Ta,,(xo)
s(xO),
b
Ta mu
air,
are
is
used
i.e.,
(X.23)
and
c
is computed and the
and h,,,(xd)
ia,_
next
series
using
solution
are evaluated,
can be calculated
- cpa -&,ho)l/Wj~.2
the
surface
computed
are then
from
the
a rearrange-
exit
computed
for
conditions is
obtained
of Cfac(xo) ,(x0).
on T is
If
hft(x,),
and air
T-I, ft(xo),
condition:.
hW and ?z; from which from
nW ft
ave and aie
This
iteration
a positive
value
used to continues
co(xO).
checakid
Cfac(xO)
coil.
and
Arithmetic
of s ave and s(xdi the
(L24)
+ 0.444 Ta,Jxo)l
Cfac(~o)9
The new estimates
new exit
a,
psychrometric
temperature,
for
on T
a,
are repeated
a form
can be calculated
convergence
iteration
co(xo)
bulb
ba,,(xo)
The value '1 -_
This
i.e.,
Save are obtained. t
in Appendix
ig obtained.
coefficients
and an estimate
a,
of T
ratio
derived
= Uft,m(xO) - T,l/s(xo) + hft,m(xo)
ment of Eq. 1.13,
s(xO)
m(xO)
ft lnW~ in Appendix H.
derived
m(~o)
has been evaluated,/the +or h
ha/&)
-*i
and c have been cal-
achieved.
m(~o)
used to obtai;
Tft
m(xo)
T
ft ,I& xo) Once T ft
b,
for *_",
is used to evaluate thk ft ,JXo) accurately, upon which the calculations for
a,
is
for
equal
to zero,
the
following value
of
the
I
170 J?.Jxd)
is
solution
set of T ft
bega"ns near from
equal
m(~O)m
the
and latent
two-phase
Appendix
edge of the
could
Upon completion region
iterations
Such a situation
the'trailing
h to k,(x,)
sensible
to h and the
-
cause of the
heat
might coil
where
the wet region calculation
transfer
of the
repeated
rates
evaporator
starting
occur the
if
of the exit and the water are computed
dehumidification
discrete
to revert
at the jump
to a dry air
surface.
conditions, removal
as described
the
rate
in in
Ma
References “for Appendix I 6.1
F. C. McQuiston, "'Fin Efficiency Transfer," ASHRAETransactions,
with Combined Heat and Mass Vol. 81, Pt. 1, 1975.
171 J - Ratios of Geometric Calculations
Appendix
*
The correlations
. -4
for
make use of ratios The need for repeated are
of air-
calculated;
This
manipulations
section
and they
A
-
NT
-
transfer
coefficients
geometric
parameters.
discussed
is-devoted
elsewhere
to deriving
the
equations
were obtained
and is
how these degrees
The input
data
not
parameters
used are not
by various
of terms,
area,
number of tubes
In the
*= number of equivalent, distance
between
horizontal
D
outer,
or air-side,
D r
inner,
or refrigerant-side,
FP
fin
pitch,
6
fin
thickness,
suff icient
distance
altogether
of algebraic required
for
between
diameter
the
h
height
of the
heat
exchanger
e
length
of the
heat
exchanger
d
depth
geometric
calculated:
of the
heat
number of tubes basic ratios
of tubes,
and condenser
nor explicitly
geometry-related
some basic
evaporator
as input
These remaining
(center-to-center),
of tubes,
there
V
c'ircujts,
(center-to-center),
tubes
diameter
(horizontally),
refrigerant
tubes
However9
N
of ajr-flow
and
to describe are
direction parallel
wT a
c
is
cross-sectional
vertical
are
and heat
the Heat Exchanger
exchanger:
N eLircl.lit
I* c
terms
and cancellation
each heat
drops
Used in
and refrigerant-side
in many instances
straightforward,
*
pressure
and use of these
here.
Parameters
parameters
parameters (see
Fig.
completely,
that
are not
required
J.1)
exchanger in
a row (vertically) are not
can be computed
computed by careful
because
all
use of the
of the given
172
h
--_ I w I J-1 I
I I I
Fig.
J.l
Sample Tube-and-Fin
Heat Exchanger
173
input
The following
parameters.
geometric
ratios
N,
number of fins
=
used in the
equations Design
N,=R
N tot
total
=
are
used to compute
Model.
(J+U
Fp
number of refrigerant N tot
the
tubes =
NvNT
d
=
h
= NV s,
(J-2)
NT!!,
(J-3)
CJ.4)
(J-5)
V
hx
volume
=
V
hx
A
side =
A,
of heat
exchanger
=
=
lhd
~(NVST)(NTWT)
the area of each side
Aside
=
=
fin
of each fin
Ntot~Da2 = N
hd-
NtotRSTW,
(5.6)
is
TID '
totSTWT - NtotT-
4 heat
=
transfer
area, (J-7)
Af
=
2NfAside
=
21FFpNtot
.
174 A
=
contact
contact
area
A
A a
=
total
heat
A a
area
A
,where l
exposed be written using
is
=
contact
air-side
Atubes =
between
transfer
of exposed
=N
the
=
XFpNtOtnDa6
(5.8)
area
tubes
tube
surfaces
totnDa
length
and tubes
NN f totnDaS
=A,+A
tu'bes
fins
' exposed
of each tube
not
' covered
by fins,
which
can
Eq. J.l:
&
= 1 - Nfcs = L(1
exposed
- Fp6)
(J 09)
rrD ' A a
A r
=
= 2Ntot~Fp~S,W,
refrigerant-side
- -+I
heat
A r
A free-flow
=
free-flow
A
free-flow
frontal
+
NtotD~aR(l
transfer
=
NtotRIID
P
(J.10)
area
r
area
= A - A*tUbes
- F 6)
-A*,
(J.11)
175 where
A”
and using
total
=
tubes
obstructing
air
flow
NvDa 'exposed
free-flow
The ratios
that
=
NvD, le(l
Eqs.
=
area
Nfh6
l(NVST)
=
RNV[S, - D,(l
=
lNV(S,
fin
J.l,
-
=
=
are needed
c,
Substituting
=
cross-sectional
A*f A
=
tubes
of tubes
Eq. J-9
A"
PIQ,
' area
= "projected
AXtubes
of the
fins
obstructing
NVD,R(l
- Fps) - R(NVS,)Fp6
- Fp6) - STFpS]
-- Da)(l
(3.72)
- FpS)
by the heat
pump model
=
A
contact
and 5.8
2RF Ntot(STWT -4 =
are: Af
area
lIDa Cf
air-flow
RFp"(NVST)
heat transfer contact area
5.7,
- FpS)
rrD 6 RF N p tot a
)
=
m,w,
l-IDa2 - -q-1 __riDas
(5.13)
176
a
Substituting
a
= total air-side volume of heat
Eq. J.6
a
a
and J.10: llDaL (S,W, ---q-)
2NtotLF
=
A a = V hx
area exchanger
+
NtotKDal(l
- F 6)
=
Ntot%WT
2F (S,W, -+
)
+
r~D,(l
-F
6) (5.14)
%Jm
a
= refrigerant-side volume
r
of
heat transfer heat exchanger
area
c b
Substituting
Eq. J.6
and J.11: f
a
N toPD N,,pTw;
=
nDr s,w,
= fin heat transfer area total air-side heat transfer
F a Substituting
r
=
5.7
(5.15)
area
and J.10:
lIDa 2&F Ntot(S,W, F a
= 2LN totFp(S,W,
-4
)
Trn 2 - "T ) + N&D
(J.16) a
J
(1 - F a) P
6
177
u
Substituting
I free-flow frontal frontal area
a
u
=
a
+,(S,
7r
Eqs.
( NVNT) ((q) 0
=
rrDr
GaAaa
(5.18)
area
A +zGA a 5.2, flDr
5.5,
(3.19)
N#nDr =
refrigerant-side heat transfer mass flow rate of air
Substituting
transfer ci,rcuits
llDrL N tot Ncircuits
N circuits
0
- F 6)
ST
=
A (NAN,) G
=
A
(ST - Da)(l
=
and then Eqs. 5.2 and 5.5:
A r =N circuits
A
Q
- F 6)
(5.17)
= refrigerant-side heat number of parallel
Eq. J.11
Arh
=
a
R(NvsT)
Substituting
0
Afree-flow A
=
- Da)(l
A rh
A
and 5.12:
Eq. 5.4
u
area
(5.20)
'TNcircuits
area
A r a free-flow J.11,
and 5.17
(5.21)
into
5.21:
178
A min
=
total total
air-side heat transfer area free-flow frontal area A A
min
=
A
a
(5.22)
free-flow
Note that
Aa = ($)A, r
and then
A
min
=
=
'2"ArhNcircuits' r
substituting
'aArhNcircuits LXrAo
a
Eq. J.. 17 into
5.22:
179
AVAILABILITYAND DISTRIBUTIONOF
APPENDIX K:
THE HEAT PUMP MODEL
4
\ The ORNL Heat Pump Model heat
pump research
form
of six
(1)
is
available
It
development.
files
on a nine
FORTRAN source
track
to organizations
is customarily
magnetic
of the model
set
tape
performing
distributed which
up for
in
consist
batch
the
of:
processing
on a
computer,
(2)
sample
input
(3)
FORTRAN source evaporating
data
for
the
of a program
temperatures
I
4 I
I
compressor
(4)
sample
input
(5)
FORTRAN source batch
the
a fuction
in Section
the map fitting
version
program
that
of the
(some of these
addition
in file
1 and others
FORTRAN source set
up a file
for of
1 that
have the
an interactive input
data
for
program
to use the 4.3),
in file
Heat Pump Model
program
in file
and
3,
can be used to convert
to an interactive to those
(this
required
described
subroutines
of condensing
map data
coefficients
model
for
for
processing
subroutines
(6)
data
to fit
to compressor
can be used to generate map-based
Heat Pump Model,
subroutines
1 are
in
are used to replace
same name),
program the
in file
the
batch
that version
and can be used to of the
Heat
Pump Model. These files the
FORTRAN source The tapes
e
are
listed
in Append ices
L-Q (with
files).
are written
on an IBM computer
and:
cross-references
of
180
o
can be at
800, 1600,
e
can be in
ASCI,I or EBCDIC,and
e
can have either
After
a tape
Model
(file
on the
tape
The printed from
these they
logical
It
compressor
are executed
using
of the
who wish contain
is
not
recipient
of the
programs
The interactive into
seven
compressor to fit
on the
shown in Table
input
user
is
screen
for
value.
Data
giving
the
e save the
data
for
save the
data
so that
to override
computer
can vary related
for
indoor
the
returned.
features.
in files
The
one computer will is
be different. left
5 and 6 is
related
changed
to each
exchanger)
by first the
version
can be read
by the
default and get
data
designed
an item which Once the
he can: of the model, interactive
values,
output
:
pages are
specifying
results
batch
(e.g.,
and is
of these
I( B
organized
information
number and new value.
the
the
output
need to be aware
and it
heat
with
Heat Pump Model
terminal.
is
subroutines
Examples
item
use with
or replace
tape
are
2 and 4).
on a page is assigned
satisfied
it
3) that
himself.
programs
is
Heat Pump
and the
from
differences
entry
and is
o execute
I/O
the
tape
programs
Each page contains
Each data
on and then
the
(ffle (files
dependent
of a CRT terminal.
K-1.
data
when the
interactive
the
has made changes
l
generated
these
for
parameters
number and a default page it
that
to do this
pages of data. data,
appropriate
and output
to identify
of the
program
the
and installation input
fiies. copies
map fitting
and some of the
difficult
and the
are enclosed
to use the
used for
to another
job
runs
system
or non-labeled
is rewound
1) and the
units
system
it
two computer
People that
labeled
is written
listing
6250 BPI,
or
listed
model
and/or on the
i
I !
181 The user 0 C
needs to be aware that
(e.g.,
map-based
control
device)
when he selects
or loss-and-efficiency the
The programs
last
that
provided
as a convenience
formally
documented
to be understandable.
compressor
one chosen
is what will
are distributed to the
although
there
from
with
users
model,
type
of flow
be used.
the Heat Pump Model
of the model.
should
two or more alternatives
be sufficient
are
They are not comments
in them
4
P .
t
183 Table K.l
I
Sample Pages of Input Data for Interactive of Heat Pwnp Model
Version
c I
. I
PAGEO-1: COMPRESSOR DATA (EFFICIENCY & GEOMETRY: 1. DISPLACEMENT (CU IN) 4.520 2. CLEARANCE-VOLUME RATIO .060 14OTOR: 3. SPEED(RPM) -. 3450. 4. OUTPUTAT FULL LOAD (KW) 2.149 EFFICIENCIES: 5. MAXIMUMMOTOR .820 6. ISENTROPIC .700 7. MECHANICAL .800
LOSSMODEL) LOSSES: 8. HEATREJECTIONRATE FROM COMPRESSOR SHELL (BTU/HR) 0.0 OR POWER 9. FRACTIONOF COMPRESSOR REJECTEDFROMTHE SHELL -350 HEATTRANSFERRATE FROMHIGH SIDE TO LOWSIDE (BTU/HR) 300.0 OR 11. FRACTIONOF COMPRESSOR POWER TRANSFERED FROMHIGH SIDE TO LOWSIDE .ooo 10.
NOTES: 1. ITEM 8 OR 9 IS SPECIFIED, NOTBOTH 2. ITEM 10 OR 11 IS SPECIFIED, NOTBOTH 3. IF ITEM 4 IS SET TO 0, THE NECESSARY OUTPUTAT FULL LOADWILL BE CALCULATED AND ITEM 3 IS THE ACTUALMOTORSPEED 4. IF ITEM 4 IS NOT 0, ITEM 3 IS THE SYNCHRONOUS MOTORSPEEDAND THE ACTUALSPEEDIS CALCULATED TYPE "RETURN"TO CONTINUE
PAGE#2: COMPRESSOR DATA (MAP-BASEDMODEL) GEOMETRY: LOSSES: 1. DISPLACEMENT (CU IN) 4.520 6. HEAT REJECTIONRATE FROM 0.0 2. BASEDISPLACEMENT FOR COMPRESSOR SHELL (BTU/HR) 4.520 THEMAP (CU IN) 7. FRACTION':F COMPRESSOR POWER MOTOR: 3. RATEDSPEED(RPM). REJECTED FROM SHELL .350 3450. SUPERHEAT AT SHELL INLET: 4. BASESUPERHEAT FOR MAP (F) 20.00 CURVEFITS FROMCOMPRESSOR MAP (AS FUNCTIONSOF CONDENSING & EVAPORATING TEMPERATURES): 8. .REFRIGERANT&SSFLOWRATE 5. POWERCONSUMPTION -1.509E-04 -2.6753-02 CONDENSING **2 CONDENSING l *2 4.089E-02 4.633E+OO CONDENSING CONDENSING -1.3383-04 4.7033-02 EVAPORATING l *2 EVAPORATING l "2 5.8603-04 9.64OE+OO EVAPORATING EVAPORATING CONDENSING*EVAPORATING 3.6383-04 CONDENSING*EVAPORATING -1.8683-02 1.207E-04 CONSTANT 9.7593-05 CONSTANT 'NOTE: 1. ITEM 6 OR 7 IS SPECIFIED, NOT BOTH TYPE,"RETURN"TO CONTINUE
x.
I.,
(2
. .
.\..,a.
/.
._
.L^
..“”
,
.”
..-.,
. .
”
.._.
PAGE#3: INDOORHEATEXCHANGER GEOMETRY: 3.167 1.. FRONTALAREA (SQ FT) 2. NUMBEROF PARALLELCIRCUITS 3.0 3. NUMBER OF TUBE ROWS 3.0 4. NUMBER OF RETURNBENDS 72.0 FAN - FAN MOTOR: 5. COMBINEDEFFICIENCY .200 INLET AIR: 6. FLOWRATE (CFM) 1200.0 7. TEMPERATURE (F) 70.0 8. RELATIVE HUMIDITY .500
FINS:
TYPE (1 SMOOTH,2 WAVY, OR 3 LOUVERED) PITCH (FINS/INCH) THICKNESS(IN) THERMALCONDUCTIVITY (BTU/HR-FT-F) CONTACTCONDUCTANCE (BTU/HR-FT2-F) TUBES: 14. O.D. (IN) 15. I.D. (IN) 16. VERTICAL SPACING(IN) 17. HORIZONTALSPACING(IN)
2. 14.0 .00636 128.0 30000. .400 -336
1.ooo 0.875
TYPE "RETURN"TO CONTINUE
OUTDOOR HEAT EXCHANGER PAGEf4: GEOMETRY; 1. FRONTALAREA (SQ FT) 5.040 2. NUMBEROF PARALLEL-CIRCUITS 4.0 3. NUMBEROF TUBEROWS 4. NUMBER OF RETURNBENDS 62:: FAN - FAN MOTOR: 5. COMBINED EFFICIENCY .160 INLET AIR: / 6. FLOWRATE (CFM) 2300.0 7. TEMPERATURE (F) 47.0 8. RELATIVE HUMIDITY .700
FINS: 9. TYPE (1 SMOOTH,2 WAVY, OR 3 LOUVERED) lo. PITCH‘ (FINS/INCH) 11. THICKNESS(IN) 12. THERMALCONDUCTIVITY (BTu/HR-FT-F) 13. CONTACTCONDUCTANCE (BTU/HR-FT2-F) TUBES: 14. O.D. (IN) 15. I.D. (IN) 16. VERTICAL SPACING(IN) 17. HORIZONTALSPACING(IN)
2. 14.0 .00636 128.0 30000.
NOTES: 1. IF ITEM 85 IS SET TO 0, THE FAN VS STATIC EFFICIENCY CURVEIS USED WITH A 55% EFFICIENT MOTOR TYPE "RETURN"TO CONTINUE
.400 -336 1.000 0.875
Table K.1
(continued)
DEVICE PAGEP5: FLOWCONTROL FIXED CONDENSOR SUBCOOLING IS BEING USED THERMOSTATIC EXPANSIONVALVE: SPECIFIED CONDENSER SUBcooLING:. '5. VALVE CAPACITY RATING (TONS) .ooo 45.00 1. DESIREDSUBCOOLING(F) 0.0 STATIC SUPERHEAT (F) 6. OR 7. RATEDOPERATING CAPILLARYTUBES: o*o SUPERHEAT(F) 0 2. NUMBEROF CAPILLARYTUBES .ooo BLEEDFACTOR 3. CAPILLARYTUBE FLOW . NOZZLE& TUBE PRESSURE ;: 0.000 FACTOR 0 DROP(1 INCLUDE, 0 EXCLUDE) OR SHORT-TUBEORIFICE: .ooooo 4. DIAMETEROF ORIFICE (IN) OR NOTES: 1. CATEGORIESARE MUTUALLYEXCLUSIVE, I.E. SPECIFY 1, OR 2 & 3, OR 4, OR5-9 2. IF CONDENSER SUBCOOLING IS SPECIFIED, THE VALUESOF ITEMS 2-9 THAT WOULDPRODUCE THAT LEVEL OF SUBCOOLING ARE COMPUTED TYPE "RETURN"TO CONTINUE
I
..I
..L.....h*.-.-wU
184
PAGEP6:
CONNECTING REFRIGERANT LINES
GEOMETRY: LIQUID LINE SUCTIONLINE DISCHARGELINE VAPORLINE VAPORLINE
FROM: INDOORHX REVERSINGVALVE COMPRESSOR OUTLET REVERSINGVALVE REVERSINGVALVE
TO: OUTDOOR HX COMPRESSOR INLET REVERSINGVALVE OUTDOOR HX INDOORHX
ID (IN) 1. 0.190 3. 5. 7. 9.
0.000 0.000 0.680 0.556
EQUIVALENT LENGTH(FT) 2. 30.4 4. .6. 8.
10.
RATESOF LOWSIDE HEATGAINS: RATESOF HIGH SIDE HEAT LOSSES: 13. IN SUCTIONLINE (BTU/HR) 11. FROMDISCHARGE LINE (BTU/HR)2000. 200. 12. FROMLIQUID LINE (BTU/HR)
0.0 0.0 6.0 30.4
> +
300.
NOTE: 1. GAIN & LOSSESIN ITEMS 11-13 C!N BE FORVALIDATION RUNSOR TO SIMULATEHIGH TO LOWSIDE HEATEXCHANGERS ANDREVERSINGVALVES TYPE "RETURN"TO CONTINUE
PAGE#7: MISC. DATA 1. TITLE: EFFECTSOF HEAT LOSSESON INLET AIR: 2. HOUSELOADAT OUTDOOR 7. INDOORFAN (SEE NOTEBELOW) 15000. 2 AMBIENT(BTU/HR) 8. OUTDOOR FAN (SEE NOTEBELOW) 10.0 i SUPERHEAT(F) 3. EVAPORATOR 9. COMPRESSOR (SEE NOTEBELOW) 8.0 1 4. DUCTSIZE (IN) ESTIMATESTO BEGIN ITERATIONS: 5. MODEOF OPERATION TEMPERATURE (F) 2 10. EVAPORATING (1 COOLING,2 HEATING) 30.0 11. CONDENSING TEMPERATURE (F) 115.0 6. LEVEL OF OUTPUT(0 SUMMARY, 12. REFRIGERANT MASSFLOW 1 COMPRESSOR, 2 COILS, RATE (LBM/HR) 1 400.0 3 STANDARD) NOTE: 1. HEATLOSSESFROMTHE FANS AND COMPRESSOR CANBE ADDEDTO THE AIRSTREAM EITHER BEFOREOR AFTER‘IT CROSSES THE COIL (SIMULATEBLOW-THROUGH & DRAW-THROUGH FANS) SPECIFY 0 TO NEGLECTTHE HEATLOSS THE COIL 1 TO ADD THE HEAT To THE AIR BEFOREIT CROSSES HEAT To THE AIR AFTER IT CROSSES THE COIL 2 TO ADD'THE TYPE "RETURN"TO CONTINUE
'-
f*
‘(I
tr
185 Appendices
L, M, N, 0, P, and Q are
J
8 ;
1 -li
i :
zr
to the
back cover.
.
on microfiche
in an envelope
attached