Thermal Analysis Research Program Reference Manual
Thermal Analysis Research Program Reference Manual
March 1983
Sponsored by
National Bureau of Standards
1
and
George N. Walton
March 1983
Sponsored by
National Bureau of Standards
and
The Thermal Analysis Research Program (TARP) has been developed as a research
tool for the thermal analysis of buildings. It especially aims to study the
interactions of many complex heat transfer phenomena. TARP uses the detailed
heat balance method for the simultaneous calculation of the energy requirements
of multiple rooms. Interroom conductive and convective processes are simulated
in detail. This program reference manual describes the algorithms, input, out-
put, and program structure of TARP. The program is written to be portable and
modifiable. It is written in FORTRAN 77 and has run on CDC and UNIVAC compu-
ters. Future expansions of the program are anticipated, particularly for the
simultaneous simulation of the equipment performance and building thermal
response.
iii
PREFACE
This report is one of a series documenting NBS research and analysis efforts in
developing energy and cost data to support the Department of Energy/National
Bureau of Standards Measurements Program. It was prepared by the Thermal
Analysis Group, Building Physics Division, Center for Building Technology,
National Engineering Laboratory, National Bureau of Standards (NBS). This work
was jointly sponsored by NBS and the Department of Energy (DoE). The develop-
ment of multiroom air flow modeling was supported by DoE/NBS Task Order A008
under Interagency agreement No. EA-77-A-0 1-6010. This report describes a
computer program which was written as part of an effort to develop a comprehen-
sive modeling technique for predicting the simultaneous transfer of air, mois-
ture, and heat in and through multiroom buildings.
The author wishes to acknowledge the support and direction given by Dr. T.
Kusuda, Thermal Analysis Group Leader.
The computer program described in this report is intended for research purposes
only. It is furnished by the government and is accepted and used by any recipi-
ent with the expressed understanding that the United Stated Government makes no
warranty, expressed or implied, concerning the accuracy, completeness, reli-
ability, usability, or suitability for any particular purpose of the information
and data contained in this program or furnished in connection therewith, and
the United States shall be under no liability whatsoever to any person by reason
of any use made thereof.
iv
TABLE OF CONTENTS
Page
I. INTRODUCTION 1
A. BACKGROUND 1
B. OBJECTIVE 1
C. SCOPE 2
III. ALGORITHMS 13
A. INTRODUCTION 13
1. The Heat Balance Method 13
2. Table of Symbols 13
B. ENVIRONMENT 16
1. Time 16
2. Psychrometrics 18
3. Temperatures 21
4. Wind 22
5. Solar radiation 22
C. TRANSIENT HEAT CONDUCTION 26
1. Use of Conduction Transfer Functions 26
2. Calculation of Conduction Transfer Functions 27
D. WINDOW OPTICS 29
1. Optical Coefficients 29
2. Shading Coefficients 30
3. Uncoated Pane Optics 31
4. Reflective Coatings 33
5. Multi-Pane Windows 36
6. Diffuse Radiation 37
7 . Interior Shading Devices 38
E. SOLAR GAINS AND SHADOWING 46
1. Surface Geometry 46
2. Shadow Projection 47
3. Homogeneous Coordinates 48
4. Overlapping Shadows 49
5. Solar Gains 31
F. INTERNAL HEAT GAINS 58
1. Sources and Types of Gains 58
2. Distribution of Radiant Gains 59
G. NATURAL AIR MOVEMENTS 63
1 . Wind Forces * 63
2. Thermal Forces 63
3. Flow Equations 64
4. Solution of Flow Equations 65
H. ROOM RADIANT INTERCHANGE 68
1. Detailed Radiant Interchange 68
2. MRT Network Algorithm 68
v
Page
IV. INPUT 94
A. OVERVIEW 94
B. SYNTAX 94
1. BNF Description 95
2. English Description 95
3. Error Processing 97
C. INPUT DETAILS 98
1. Project Section 98
2. Library Section 102
3. Building Section 108
4. Zone Section Ill
D. INPUT SUMMARY 124
E. ZONE TEMPERATURE CONTROL 127
F. GEOMETRIC DATA 136
1. Surface Coordinate System 136
2. Window and Door Subsurfaces 137
3. Shadowing Subsurfaces 137
4. Zone Coordinate System 138
5. Building Coordinate System 138
6. Detached Shadowing Surfaces 139
G. TABLES 150
1. Shadowing Dates 150
V. OUTPUT 151
A. INPUT REPORTS 151
1. Default Reports 151
2. Run Control Modifications 156
3. Numbered Reports 160
4. Library Reports 170
vi
Page
vii
LIST OF FIGURES
Page
viii
LIST OF FIGURES (cont.)
Page
ix
LIST OF FIGURES (cont.)
Page
x
LIST OF FIGURES (cont.)
Page
xi
,
I. INTRODUCTION
I. A BACKGROUND
I.B OBJECTIVES
1
computing room loads and temperatures by simultaneously considering conduction
and air movements between all rooms in a building. This program provides the
framework for the integration of many more advanced simulation features.
I.C SCOPE
Chapter III reviews the algorithms that the program uses to compute the
building energy and air flows.
References are indicated by brackets and are listed at the end of each
section.
REFERENCES
[1] Kusuda, T., "NBSLD, the Computer Porgram for Heating and Cooling Loads in
Buildings," NBS Building Science Series 69, (July 1976).
2
II PROGRAM OVERVIEW
II .B PROGRAM STRUCTURE
3
The TARP is divided into two subprograms: the input processor (IP) and
the loads calculator (LC). This structure is illustrated in figure II. 1. The
IP checks the user input, prepares input files for the LC, and reports the
array sizes that the LC will require. The LC performs the building loads
simulation and reports the results. This separation of subprograms allows a
reduction in computer core requirements. The IP uses minimal array sizes
because it processes the input one element (surface, schedule, material, etc.)
at a time. The LC uses large arrays which consider these elements simultan-
eously. The IP and LC may be run as separate jobs when It is necessary to
recompile the LC for correct array sizes, or they may be run as a single job
if the LC is known to have sufficiently large arrays.
The loads calculator is divided into three sections. The first reads the
input file created by the IP and fills the LC arrays. It checks for sufficient
array space and performs some checks which require that all data for the pro-
blem description be present in the arrays. The second section uses the geo-
metric and optical descriptions of the building surfces to compute solar gain
factors for various periods of the year. The third section does the actual
heat balances necessary to compute the multi-room loads. It does this on an
hour-by-hour basis.
4
The TARP input for the example is shown in figure 3. The input begins by
describing project data in the first nine lines between "PROJECT [" and A
project description is given; the location is described; ground temperatures
are given. Then there is a description of the weather between "DDY(" and ")"
on the following line. This sets up a mild, late winter condition when the
solar and internal gains nearly offset the infiltration and conduction losses.
The low temperature (LOW) for the day is 30F and the high (HIGH) is 50F. A
clear sky (CLR) for 21 February (DATE) is assumed. A 15 mph (WS=1320 ft/min)
wind blows from the northwest (DIR^IS).
There are three "ZONE" sections to describe the three zones of the example
building. Each zone is identified by a user supplied name, e.g., "NORTH".
Activities related to the entire zone are the temperature control (CTRL) and
people (PEO), lights (LIT), and equipment (EQP) loads. There are three people
maximum (CAP=3) in the zone. The hourly number of people is determined by that
maximum and the previously defined week-schedule (WS=OCCN). The other activi-
ties have similar capacities and schedules. The rest of the zone description
defines the surfaces, both envelope and internal, of the zone. The mass (MASS)
of furnishings is handled by a single surface of appropriate density and area.
"SRF" then describes most of the north wall of the zone as an exterior (EX)
opaque (OP), base surface (BS). It is made of the previously defined WL-I
construction. The wall faces north (AZM=0) and is 41 ft long by 8 ft high
(SIZE). It is located relative to the other surfaces of the building by the
origin (ORG) command. The next line of input specifies a window on the wall.
The window is an exterior (EX), transparent (TR) subsurface (SS) made of the
,
WDW2 construction. Its size is made to represent all windows on the north
wall. Openings in the wall for identification are combined into a single area
under the "AFS" command. The next surface command describes the remainder of
5
the north wall: the portion at the studs which Is construction WL-S. Four
lines later there is a description of the massive partition described as an
interior (IN), opaque (OP), base surface (BS). It is thermally connected to
the center zone by the command "ZONE-CENTER".
Figure 4 shows the computed loads for the sample run. The report is headed
by the project, location, and environment descriptions. Then the following
values are presented for each zone: Total heating and cooling requirements,
peak heat and cooling loads, maximum zone air temperatures, hours at peak
heating capacity (HX), other heating hours (HT), ventilating hours (V), cooling
(below peak capacity) hours (CL), and hours at peak cooling capacity (CX).
Building-wide values are summarized by the line labeled "TOTAL:". These values
include the total heating and cooling requirements for all zones, the peak
coincident heating and cooling loads , and maximum and minimum of all zone
temperatures. The hour and day are reported for the occurance of the peak
loads and temperatures.
This example was not intended to provide complete instruction. The user
should read the remainder of the program documentation before attempting to run
TARP.
REFERENCES
[1] Hittle, D. C., "The Building Loads Analysis and System Thermodynamics
(BLAST) Program, Version 2.0: User Manual", CERL Technical Report E-153,
(Construction Engineering Research Laboratory, 1979).
[2] Herron, D. Walton, G., Lawrie, L., "Building Loads Analysis and System
,
6
Figure II. 1 Structure of the Thermal Analysis Research Program
Figure II. 2 House configuration for example problem
8
1 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/07, 05.38.40. .713
1 PROJECT C
2 RC(DEM/UIN=ENGLISH/UOUT=ENGLISH/
3 DESC= THERMAL MASS TEST - SIMULTANEOUS ZONES* DIRECT SUN')
'
4 RPT ( RNOS = 9)
5 LOC<DESC=' COLUMBIA* MO. '/LATD=38.?/L0NG«92.3/TZ*6>
6 GRND(GRT=40*42*48*55*62*67*70*69*63*55*48»42)
7 DDY DESC= MILD WINTER ' /HI GH=50/L0W=30/UB=40/DATE-21FEB/
< '
8 CLR=1./M0N/US=1320/DIR=315)
9 3
10 L I BR AR Y C
11 DS(NAME=NONE/DESC='ZERO ALL HOURS ' /FFC=24*0
12 DS<NAME=FULL/DESC='ONE ALL H0URS'/FFC=24*1
13 DS<NAME=ON/DESC=' NORTH OCCUPANCY'/
14 FFC=6*1, .5* .16*8*0*4*. 16 *2*. 33* .67*1)
15 DS(NAME=OC/DESC*' CENTER OCCUPANCY'/
16 FFC=16*0*3*. 17* .67* . 5 * . 67 * 17 0
. »
48 MATL(NAME=R19/DESC='R19 INSULATI0N'/R=19.
49 MATL < NAME = R30/DESC= R30 INSULATION /R=30 .
' '
9
) )
> ) > ) )> ) )
>
PAGE 2 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/07. 05. 38. 41. 1.050
83 1
84 BUILDING C
85 SIM(SLDS=3/HTB=0/FL0W/CNVG=.l. .01. .34.12)
86 BSS(AZM=180./TILT=153.4/0RG=0. 12. 16/SIZE=44.2)
87 DSS< AZM=180./TILT=206.6/0RG=0.0.8/SIZE®44.2)
88 DSS< AZM=0./TILT=116.6/0RG=0.36.8/SIZE=44.2)
89 1
90 ZONE C
91 GEOM(NAME=NORTH>
92 CTRL(CWS=HEATN/HEAT=12)
93 PE0<WS=0CCN/CAP=3>
94 LIT <WS = LITN/CAP = 1 ,4)
95 EQP WS=EBPN/CAP=1 . 0
<
96 MASS(C0NS=PTN-A/AREA=480)
97 SRF(EX/0P/BS/C0NS=WL-I/AZM=0/0RG=44.36.0/SIZE=41.8)
98 SRF(EX/TR/SS/C0NS=WDW2/0RG=12.5.-,25/SIZE=20.2)
99 ‘
AFS(AREA=0. 2/EXP®. 65)
100 SRF(EX/0P/BS/C0NS=WL-S/AZM®0/0RG=3.36.0/SIZE=3.8)
101 SRF(EX/CP/BS/CONS=RF-T/AZM=0/0RG=44.36.8/SIZE=41.8,95/TILT=26.6)
102 AFS(AREA=0. 1/EXP®. 65)
103 SRF(EX/0P/BS/C0NS=RF-I/AZM=0/0R6=3.36»8/SIZE=3.8.95/TILT=26.6)
104 SRF( IN/0P/BS/C0NS=PTN-M/AZM=180/0RG=0.24.0/SIZE=44. 12/Z0NE=CENTER)
10
PAGE 3 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/07. 05.38.42. 1.467
105 AFS(AREA=0.5/EXP=.50)
106 SRF(IN/0P/BS/C0NS=PTN-M/AZM=270/0RG=0,36,0/VRTS=12,0, 12,12, 0.8)
107 SRF < IN/0P/BS/C0NS=PTN-M/AZM=90/0RG=44 . 24 » 0/VRT8=l2 , 0 . 12.8. 0.12)
108 SRF(EX/0P/BS/CQNS=SLAB-C/AZM=180/0RG=0,36»0/SIZE*44»12/TILT=180)
109 OSC TGP=1
(
110 3
111 ZONE C
112 GEOM(NAME=CENTER)
113 CTRL(CWS=HEATC/HEAT=12)
114 PE0(WS=0CCC/CAP=3)
115 LIT(WS=LITC/CAP=1.4)
116 EQPCWS=E0PC/CAP=1.0)
117 MASS(C0NS=PTN-A/AREA=480)
118 SRF(EX/0P/BS/C0NS=FRAME/AZM=180/0RG=0»12»12/SIZE=44,3)
119 SRF<EX/TR/SS/C0NS=WDU2/ORG=2,0»-.25/SIZE=4O,3)
120 AFS(AREA=0.25/EXP=.65)
121 SRF<EX/0P/BS/C0NS=RF-T/AZM=0/0RG=44»24,12/SIZE=41»8. 95/TILT=26 . 6
122 AFS(AREA=0. 1/EXP=.65)
123 SRF(EX/0P/BS/C0NS=RF-I/AZM=0/0RG=3»24»12/SIZE=3»8.95/TILT=26.6)
124 SRF(IN/0P/BS/C0NS=PTN-A/AZM=180/0RG=3,12,0/SIZE=41,12/Z0NE=S0UTH>
125 AFS(AREA=0.5/EXP=,50)
126 SRF IN/OP /BS/C0NS=PTN-S/AZM=180/0RG=0 .12 »0/SIZE=3.1 2/Z0NE=S0UTH
(
127 SRF<IN/0P/BS/C0NS=PTN-M/AZM=0/0RG=44»24»0/SIZE=44,12/Z0NE=N0RTH>
128 AFS(AREA=0.5/EXP=.50)
129 SRF(IN/0P/BS/C0NS=PTN-M/AZM=270/0RG=0,24.0/VRTS=12.0. 12.16, 0,12)
130 SRF(IN/0P/BS/C0NS=PTN-M/AZM-90/0RG=44»12»0/VRTS=12»0» 12,12, 0,16)
131 SRF(EX/OP/BS/CONS=SLAB-C/AZM=18O/0RG=O»24»O/SIZE=44»12/TILT=18O)
132 OSC TGP= 1
<
133 3
134 ZONE C
135 GEQM(NAME=SOUTH)
136 CTRL(CUS=HEATS/HEAT=12)
137 PE0(WS=0CCS/CAP=3)
138 LIT (WS=LITS/CAP=1 .4)
139 EQP(WS=EQPS/CAP=9.4/L0ST=,4)
140 MASS(C0NS=PTN-A/AREA=480)
141 SRF EX/OP/BS/CONS=WL-I/AZM= 180/0RG=0 ,0»0/SIZE=41»8)
(
143 AFS(AREA=0,2/EXP=.65)
144 SRF(EX/0P/BS/C0NS=WL-S/AZM=180/0RG=41,0»0/SIZE=3,8)
145 SRF (EX/0P/BS/C0NS=RF-T/AZM=180/0RG=0,0»8/SIZE=4 1,8. 95/TILT =26, 6)
146 AFS(AREA=0.1/EXP=.65)
147 SRF(EX/0P/BS/C0NS=RF-I/AZM=180/0RG=41 0 , 8/SIZE=3 » 8 95/TILT=26 6
, . .
11
PAGE 6 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/07. 05.39.04. 5.977
COLUMBIA* MO.
MILD WINTER
12
Ill ALGORITHMS
I I I. A INTRODUCTION
where
CE = energy convected from internal sources
Si = area of surface i
HC^ = convection coefficient of surface i
T^ = temperature of surface i
TZ = room air temperature
Fj = mass flow of air through opening j
Cp = specific heat of air
TSj = temperature of the air through opening j
L = energy from the air handling system
The summations apply to all surfaces/openings which are part of the room.
The values of the surface temperatures, T^, are obtained by evaluating a
heat balance at every surface of the room using an assumed room temperature.
A better estimate of the room temperature is then obtained by solving
equation A. 1.1 for TZ. Repeating this process leads to successively better
approximations of the surface and air temperatures. The process is rapidly
convergent. Multiple rooms can be handled within the iterative scheme given
an appropriate data structure and proper heat balances on the surfaces between
rooms. Multiple rooms are thermally connected by conductive, convective, and
radiative heat transfer. The following sections describe the algorithms for
the various heat transfer mechanisms and thermal driving functions which
culminate in multi-room heat balances.
I I I. A. 2 TABLE OF SYMBOLS
13
Subscripts
Subscripts may be omitted if context is clear.
Symbol Meaning
a ambient air
g ground surface
i surface i
j surface j (usually meaning opposite side of surface i)
k order of conduction transfer functions (CTF)
m number of CTF
m zone m (usually meaning adjacent to zone n)
n zone n
t time
s sky
sg surface to ground
ss surface to sky
Summation limits •
•
Parameters:
Symbol Units Meaning
<*i absorptance of surface i for short wave radiation
ei emittance of surface i for long wave radiation
g m/s 2 gravitational constant (9.80)
Pi reflectance of surface i
0 W/m 2 K4 Stefan-Boltzman constant (5.6688E-8)
Ti transmittance of surface i
0i angle of incidence of sun’s rays on surface i
<t>i
surface tilt angle
o
surface azimuth angle
Ai area of opening i
Wh/kgK specific heat (often refers to air - 1004.)
cl„ W convective gain from all internal sources in zone n
«1 W/m^ convective flux from zone air to (inside of) surface i
CO,
1
W/m* convective flux from environment to surface i
D kg/m3 density
F kg/s mass flow rate
F ij o radiation interchange factor from surface i to surface
HA, W/m;K surface i outside convection film coefficient
HC i W/m 2 K surface i inside convection film coefficient
HG^ W/m2 K surface i radiative conductance
r\
(outside surface to ground radiant temperature)
H0 i W/m z K overall outside surface conductance
HR^ W/m2 K surface i radiative conductance
(inside surface to zone mean radiant temperature)
HSi W/m 2 K surface i radiative conductance
(outside surface to sky radiant temperature)
14
Ib W/m 2 intensity of solar beam (direct) radiation
x
d
W/m2 intensity of solar diffuse radiation
*S
W/m2 intensity of sky diffuse radiation (sky brightness)
W/nr intensity of ground reflected radiation
g
J i,k flux referring CTF
KI i,t W/m 2 inward conductive flux at inside of surface i at time t
KIP i,t W/m2 portion of inward conductive flux based on previous
temperature and flux history terms
KO i,t W/m 2 outward conductive flux at outside of surface i at time 1
15
III.B ENVIRONMENT
III.B . 1 TIME
Currently, TARP uses a one hour time step. The time from midnight to
1:00 AM is hour number 1; from 1:00 to 2:00 AM is hour 2; and from
11:00 PM to midnight is hour 24. The hour number (t 8 ) is used as an index in
the program to select array elements. Most driving functions, such as solar
radiation, wind, and scheduled activities, use their average values for the
hour. Temperatures are for the mid-point of the hour. The one hour time
step was chosen because it is the shortest for which weather data is commonly
available, and it is possibly the longest for significant detail in the thermal
analysis. The time step can be changed by slight modification of the program,
primarily in the handling of day (24 hour) schedules.
LM(M0=1) = 31
LM(2) = 28 + LY
LM(3) = 31
LM(4) = 30
etc.
MO-1
DY = DM + E LM(m) (B.1.1)
m=0
16
The Julian date (JD) is the most powerful date representation. It is
a continuous count of days from 1 January 4713 B.C. Any number of days
added to or subtracted from a JD is a valid JD. The following algorithms
[1] convert between conventional and Julian dates:
L » (MO - 14) / 12
JD - DM - 32075 + 1461 * (YEAR + 4800 + L) / 4
+ 367 * (MO - 2 - 12*L) / 12
- 3 * ((YEAR + 4900 + L) / 100) / 4 (B.1.2)
J JD + 68569
83
N - 4 * J / 146097
K = J - (146097*N + 3) / 4
YEAR' = 4000 * (K + 1) / 1461001
L - K - 1461 * YEAR’ / 4 + 31
MO’ - 80 * L / 2447
DM = L - 2447 * MO’ / 80
M - MO’ / 11
MO - MO' + 2 - 12 * M
YEAR = 100 * (N - 49) + YEAR 1
+ M (B.1.3)
The cyclic date (CD) is used in TARP to represent time periods which
begin and end certain functions. They are given by
CD - 32 * MO + DM • (B.1.4)
MO - CD / 32 (B.1.5a)
DM = MOD (CD, 32) (B.1.5b)
where
M0D(A,B) * the remainder of A/B
= A - (A/B)*B (MOD is a FORTRAN intrinsic function)
If(CD2 > CD1), then (CD > CD1 and CD < CD2) for CD to be in the period.
If(CDl > CD2), then (CD > CD1 or CD < CD2) for CD to be in the period.
(B.1.6)
17
Building activities typically differ on holidays from normal schedules.
The standard U. S. federal holidays are:
If any of the four holidays with dates fall on a weekend, the nearest Friday
or Monday becomes a holiday. Five tests can describe the federal holidays
(given that MO is correct).
A test is also needed for the New Year's holiday on Friday, 31 December.
III.B.2 PSYCHROMETRICS
18
1
Range of T A0 A1 a2 A3
Range of P A0 A1 a2 a3
a4 a5 a6
0.2617 to -80.481222 60.488509 -104.02603 116.27016
1.08 -70.240148 17 .403927
1.08 to -71.907528 17.059641 -7.6219223 2.2337690
3.94 -.35643002 .023469826
3.94 to -62.643477 4.8931575 -.58140203 .047441082
22.35 -.0023302718 6. 2036235E-5 -6. 8559241E-7
22.35 to -50.166219 1.0796773 -.025433803 .00041771901
103.26 4. 1788472E-6 2. 2878800E-8 -5. 2415828E-1
103.26 to -36.261166 .23970771 -.0010835426 3.3697165E -6
617.56 -6.2986874E-9 6. 3667218E-12 -2 .6646050E-15
617.56 to -21.4716785 .052658152 -3.8998869E-5 1.9227789E -8
2364.9 -5.1661347E-12 5 .7225294E-16
2364.9 to -6.2548222 .017442339 -3.7406571E-6 5.4267058E -10
7471.5 -4.3379212E-14 1. 4427502E-18
7471.5 to 9.3624054 .0065309333 -4.6769042E-7 2.4048849E -11
31605.7 -7.6034040E-16 1 .3232127E-20 -9 .6773935E-26
31605.7 to 33.328934 .0017898571 -2 .7520024E-8 2 .8970085E-13
102862. -1.6828604E-18 4.0679479E-24
19
The polynomials are faster than the ASHRAE detailed algorithm which
involves a logarithm and an exponential. The other relationships are based on
the ideal gas law and are generally accurate to better than 0.7 percent. The
humidity ratio may be obtained from one of the following algorithms [3].
W (T d , PB)
Pw = PS(T d )
W - 0.62198 * P
w / (PB - Pw ) (B.2.2)
W (T a , T w , PB)
P* = PS(T )
W
W (T a ,
RH, PB)
P w = RH * PS(T a )
W = 0.62198 * P
w / (PB - Pw ) (B.2.4)
W (T a ,
H)
The following values may be computed when the humidity ratio is known.
T d (W, PB)
Pw = PB * W / (0.62198 + W)
T d = TS(P W ) (B.2.6)
RH (T a , T d )
V (T a , W, PB)
H (T a , W)
T w (T a , W, PB)
20
III.B.3 TEMPERATURES
+
T at - T a - f * (T
+ ' T <B - 3a >
t a a*>
where f t is given by
t: 1 2 3 4 5 6 7 8 9 10 11 12
f .85 .90 .94 .98 1.0 .99 .95 .89 .78 .64 .48 .31
t: 13 14 15 16 17 18 19 20 21 22 23 24
f: .17 .06 .00 .01 .07 .16 .23 .41 .53 .63 .72 .79
The net thermal radiative flux between a horizontal surface and the sky
(Q ss ) is given by
4 - 4
°ss = H * (a * T0 i e
s
* a * Ta ) (B.3.2)
4 - 4
Q ss = H * a * (T 0l T
s )
(B.3.4)
so that
(B.3.5)
21
III. B. 3. GROUND TEMPERATURES
The user specifies monthly ground temperatures. Useful advice for the
selection of ground temperatures is a subject of ongoing research.
III.B.4 WIND
V = VQ * * Z / Z a (B.4.1)
3 (
Q )
where
Class Description a 8
1 ocean or other body of water with at 0.10 1.30
least 5km of unrestricted expanse
2 flat terrain with some isolated obs- 0.15 1.00
tacles (buildings or trees well
separated from each other)
3 rural areas with low buildings, 0.20 0.85
trees, etc.
4 urban, industrial, or forest area 0.25 0.67
5 center of large city (Manhattan) 0.35 0.47
Wind speeds are usually measured at class 2 sites (airports). Therefore, the
reported velocity is V Q . Eq (B.4.1) can be rearranged to compute V Q for V
measured at non-standard sites.
22
]
AU 2 for the years 1800 through 2100. This is more accurate than is needed
for building thermal analysis. Equations (B.5.3a-c) can be truncated with
some loss in accuracy. That loss can be estimated by the sum of the
magnitudes of the truncated coefficients. TARP uses the first 4, 4, and 5
terms of these equations, respectively. With the shortened equations, only
G and L are evaluated.
Then compute the solar declination angle, equation of time, and the
earth-sun distance.
e -arcsin( <f>
/ /p-0 2 ) / 15 [hours (B.5.4b)
23
2 ,
R - /p (B.5.4c)
where
The solar direction cosines can be determined from these values plus the
latitude, X (north positive, south negative), of the site.
CS 3 is the cosine of the solar zenith angle which also equals the sine of
the solar altitude angle. The solar azimuth angle, (relative to north), t|j
(B.5.7a)
(B.5.7b)
These angle calculations apply for the northern and southern hemispheres.
I
b =
Ae“ B/CS 3 (B .5 8 a)
Is = C * I b (B.5.8b)
kT = I / I0 (B.5.9)
where
( - Ib * CS 3 + I8 )
24
Iq the solar radiation on a horizontal surface above the earth's
atmosphere [ - ( sc /R 2 )*CSo
I - the solar constant (1353 w/m 2 )
9C
Weather tape data is used by means of the BLAST program weather file
processor. For transportability, a version of this processor should be
incorporated into TARP. This version should include the most current solar
radiation correlations.
REFERENCES
[1] Fliegel, H. and Van Flandern, T., "A Machine Algorithm for Processing
Calender Dates," Communications of the ACM Vol. 11, No. 10
,
Computers for the year 1978 (United States Naval Observatory, 1978)
,
pp All-17.
[9] Orgill, J. F. and Hollands, K. G. T., "Correlation Equation for
,
25
III.C TRANSIENT HEAT CONDUCTION
The conventional equation for calculating steady state heat flux from the
inside to the outside of a wall is
This equation is valid only when the wall is non-massive or when the inside or
outside temperatures are constant in time. In the more common and complex
situation where transient effects cannot be safely ignored, heat conduction
calculations which properly account for the dynamic effects must be utilized.
The methodology implemented in TARP is based on the standard response factor
technique [1] developed by Stephenson and Mitalas [2] and Kusuda [3] modified
for efficiency according to a suggestion by Peavy [4]. The modified response
factors are called "conduction transfer functions" (CTF). The heat conduction
at (and toward) the inside surface at time t is given by:
NT NT NF
KI t - E Y m*TO t _ m - E Zm *TI t -m + Z Jk*KI t -k (C.1.2)
m=0 m°0 k*l
NT NT NF
K0 t - E Y m*TI t _ m - E Xm*TO t _ m + E J k *KO t _ k (C.1.3)
m=0 m=0 k=l
TARP uses the convention that the X CTF refer to outside (or exterior) surface
temperatures, and that the Z CTF refer to inside (or Zone side) temperatures as
in [1]. The opposite convention is often used in other programs.
The current inside and outside conductive fluxes are functions of the
current inside and outside temperatures plus a factor containing the previous
temperatures and fluxes. TARP initializes the temperature and flux histories
by starting from a steady state and repeatedly calculating the building thermal
performance for the first day of the simulation period until essentially equal
results are achieved. Equation (C.1.2) can be rewritten to emphasize the
unknown values at time t.
where
NT NT NF
KIP t = E Y m *TO t _ m - E Zm*TI t _ m + E Ji*KI t _ U (C.1.5)
m=l m=l k=l
NT NT NF
KOP t Z Y^TI^- Z Xjn*TO t _ m + Z Jk *KO t _ k (C.1.7)
m*l m«l k"l
For lightweight surfaces, such as sheet metal or thin glass, there are no CTF
terms for m > 1 or k > 0. In that case Xq Yq “ Zq U, and equations (C.1.4)
and (C.1.6) both reduce to equation (C.1.1). Air cavities are handled by
massless 'slabs' which have only thermal resistance. This resistance cannot
be described as temperature dependent. The conductive properties of a multi-
pane window are described in this manner and appropriate conduction transfer
functions are computed.
One useful property of the CTF is that the steady state conductance of
the wall is related to the CTF by
NT NF
U - Z Xm / (1 - Z Jk )
m=0 k=l
NT NF
- Z Ym / (1 - Z Jk )
m=0 k=l
NT NF
- E Zm / (1 - Z Jk ) (C.2.1
m=0 j=l
This provides the method for determining the number of CTF terms to produce
the desired accuracy.
REFERENCES
[4] B. A. Peavy, "A Note on Response Factors and Conduction Transfer Functions"
ASHRAE Transactions, Vol 84 (ASHRAE, 1978), pp 688-690.
27
[5] Hittle, D. C., "The Building Load Analysis and System Thermodynamics
(BLAST) Program, Version 2.0: Users Manual", CERL Technical Report E-153.
(Construction Engineering Research Laboratory, 1979).
[6] Hittle, D. C., "Calculating Building Heating and Cooling Loads Using the
Frequency Response of Multilayered Slabs", CERL Technical Manuscript E-169
(Construction Engineering Research Laboratory, 1981), pp 20-70.
28
Ill .D WINDOW OPTICS
Q = I b *cos0 + Id (D.1.1)
The beam component (I b ) comes directly fom the sun and is incident at angle 0.
The diffuse component (I d ) comes uniformly from all directions. Some of the
incident radiation is transmitted according to
QT = T b * Ib * cos 0 + T d * I d (D.l. 2)
where
T b is the beam transmittance at 0 and
T d is the transmittance for diffuse radiation.
QA = A b * Ib * cos0 + A d * Id (D.l. 3)
where Ab and A d are the beam and diffuse absorptances . If an internal shade
is present, it also absorbs radiation with absorptances AS b and AS d .
QS = AS b * I b * cos 0 + AS d * I d (D.l. 4)
The energy absorbed in the shade is convected directly to the room air. The
energy absorbed in the glass is convected both into the room and to the outside
according to proportions (neglecting the thermal resistance of the glass which
is small compared to HI and HO) defined by
29
QA in - QA * HI / (HI + HO) (D.1.5a)
and
where
and
where
The angle dependent properties of T^, AO^, AI^, and AS(j are represented
by polynomial equations in terms of cos0. Two ranges of cos9 are used for
greater accuracy, so each property (X) has two equations
The coefficients c-^ and d^ are determined by curve fits to optical values
computed by methods described in the following sections. TARP also computes
transmittance and absorptances for diffuse radiation.
30
For single glazing this is
where
QT » radiation transmitted
QTd S a = Q T of DSA glass
QA = radiation absorbed
QAdsa " QA of DSA glass
HO “ outside surface heat transfer coefficient ( 22.7 W/m^k)
HI » inside surface heat transfer coefficient ( 8.3 W/m^K)
The reference glazing material is double-strength sheet glass (DSA) with 0.86
transmittance, 0.08 reflectance, and 0.06 absorptance at normal incidence.
where n^ and n 2 are the indexes of refraction of the two media. The refractive
index of air is effectively equal to unity. The refractive index of glass is
dependent on both the composition of the glass and the wavelength of the light.
The variation is small for typical wavelengths between 0.3 and 3.0 microns,
but it is enough for a glass prism to separate a ray of white light into the
spectrum of colors.
Reflectance is the ratio of the energy reflected from the interface to the
energy incident upon it. The electromagnetic description of light shows that
reflectance depends on the polarization of the incident light. The two compo-
nents of polarization are perpendicular to the plane of incidence (also called
transverse electric or TE) and parallel to the plane of incidence (called trans-
verse magnetic or TM) . A simple formula [3] for reflectance at the interface
shown in figure D.2 is
31
Pi = cos0i * ni for TE light (D.3.2a)
and
Pi **
cos0i / ni for TM light (D.3.2b)
Unpolarized light may be considered to consist of one half TE and one half TM
polarized light. (D.3.2) is a computationally simpler form of the more familiar
Fresnel equations.
2
t = (p 2 /Pi )* [ 2*p i / ( p i + p 2 )] (D.3.3)
"kL/cos0 2 (D.3.5)
i = e
g
A ray of light incident on a sheet of glass will be both reflected and
refracted at the front surface. Part of the refracted component will be
absorbed in the glass before reaching the back surface where it is again
reflected and refracted. This continues until all the light is either absorbed
or escapes from the front or back surfaces of the glass. Figure D.3 shows the
first several reflections of the incident light ray. Summation of the trans-
mitted, reflected, and absorbed components of the incident ray give the total
transmittance (t), reflectance (r), and absorptance (a) of the sheet of glass.
These values must be determined separately for both polarization components,
p, x (= 1 - p), ig, and cig (= 1 - Tg) are functions of the angle of incidence
and wavelength of the light. They are identical at both sides of the glass.
t = (l-p) 2 T 2 2 3 2 5
+ (l-p) p Tg + (l-p) p ^Tg + ...
g
= (
- 2 Tg
O- + 2
T 2 + ^
T g^ + **•
p p p ^
g
2
= t Tg / (1 - p 2 Tg 2 ) (D.3. 6)
32
Single pane reflectance:
r p + (l-p) 2 pT 2
+ (l-p) 2 p 3 T ^ + ...
g g
» p + (l-p)2pr
g
2
(1 + p
2
T
g
2 + ^ Tg + ••• )
2 2 2 2
“ p + t p T / (1 - p T ) (D.3.7)
g g
a » [
1— P )— 1— P Tg ] + [ (l-p)pTg-(l-p)pTg 2 ] + [(l-p)p 2 Tg 2 -(l-p)p 2 Tg 3 ] + ...
2 2 +
(l-p)(l-Tg)(l + P Tg + p T ... )
g
B t ag / (1 - p Tg) (D.3.8)
2
p = [(nj - n )/(n + n 2 )] (D.3.9a)
2 1
which combine to give a quadratic equation for Tg. The solution of that
equation is
2
t » ( -B + / B + 4C ) / 2 (D.3.10)
o
k = - ln( / L (D.3.11)
g)
This provides a means of describing the glass pane in terms of simply measurable
optical properties and makes it unnecessary to use the less exact shading
coefficients when those properties are known.
33
1
n
2
= n
2
" lk
2
[
The amplitude (not energy) reflectance at the interface between media i and j
for light incident from medium i is given by
=
Pil (Pi “ P-j)/(Pi + P-i) (D.4. 3a)
where
jk = cos * n^ [TE polarization] (D.4. 3b)
pi = n^ / cos
^ [TM polarization] (D.4. 3c)
where
D = thickness of the film and
\ = wavelength of the light
so that
~ - e -id
_ /"=
t [i = ]
(D.4 ,4b)
£
2 ? 2
P
n (P
12 23
T
f
) / (1
12 23
) (D.4. 6a)
and
34
t
13
t
12 / (1 + p 12 * p 23 * ) (D.4.6b)
Host of the above expression are complex. Rather than expanding them to
get solutions in terms of real numbers, it is easier to let the computer
evaluate the expressions by using complex arithmetic.
Coated glass has the same multiple reflection property as uncoated glass.
Derivation of the total reflectance, transmittance, and absorptance is compli-
cated by the fact that reflectances at the coatings depend on the direction of
incidence. Again consider a light ray incident from the left (as was done in
figure D.3). Properties (p, t) for both the left and right surfaces of the
glass are dependent on whether the light is incident from the left or right.
Consider the left side to be the "front" and the right side to be the "back"
of the glass. The following subscripts specify these conditions:
f 2 2
r =
pj + p^ t t
2
T / (1 - p
2 p^
r ) (D.4.8a)
1 g g
and
f = 2
t t t T / (1 - p T ) (D.4.8b)
1 4 g 2 p^ g
A simple exchange of subscripts (1:3 and 2:4) will give the expressions for
light incident from the back (superscript b).
b = 2
- p 2
(D.4.9a)
r + t t
4 Tg 4 p 2 Tg )
p3 p2 / ( 1
3
and
b = - 2
(D.4.9b)
t t t Tg / (1 p4 p2 Tg )
3 2
Since light is absorbed in the glass and in the coatings on either side, it is
easiest to compute absorptance from the energy balance: If a=l-r-t.
both sides of the glass are uncoated, Pi = P2 = P3 = P4> an<* T 1 = T 2 = T 3 = T 4
and equations (D.4.8a and b and D.4.9a and b) reduce to (D.3. 7) and (D.3. 6).
35
III.D.5 MULTI-PANE WINDOWS
Figure D.4 shows two panes of glass with fluxes of light going to both the
right (I) and to the left (J). The radiation fluxes are made up of all trans-
mitted and inter-reflected components. Let
then
I, * r0 r + Jo * t- (D.5 . la)
J
0
= x
0
* r + J1 * Z (D. 5. lb)
l l
- b
h II * t
2 + J2 * r2 (D.5. lc)
I - I„
L * t, + J,
'1
* r (D.5. Id)
x o 1 1
The optical properties for light incident only from the front are given by
setting J 2 “ 0, and substituting (D.5.1a-d) into the following definitions.
Tf - I I
2 / 0
= B * t2 (D.5. 2a)
R £ - j0 / 1
0
f *
= r/ + B * r
2
t
1
(D.5. 2b)
f .= * IL n D *
l
l (3^ o
+ ai J,) / I,
= a, f +
L
~ OB * r^
1 (D.5. 2c)
1 2
f = f *
a0 I / I
,
1 0
= B * a f
2
(D.5. 2d)
where
36
b * f
B » tx / (1 - r x r2 ) (D.5.2e)
Tb = J J = C * tl (D.5.3a)
0 / 2
b *
Rb = I
2 / J
2
= r2
b + C * rj t
2
(D.5.3b)
b b *
A2 b = (a b * J + a2^ * 1^) / J = a2 + C * r^ a
^
(D.5.3c)
2 2 2 2
Al b a b * J
l 1
/ J 2 * C * ai b (D.5.3d)
where
b
C = t / (1 - r ^ * r^ ) (D.5.3e)
2 2
Note that = Tb . All other front and back coefficients may be different.
Similar heat balances can be set up for triple glazing. They lead to the
following equations for front radiation.
f -
T e * D * t3 (D.5 .4a)
Rf - r f + E * t * (r 2
f
+ D * r^ f * t ) (D.5. 4b)
x l 2
Al f = a l f + a: b * E * (r
2
f
+ D * r f * t
3 2)
(D.5. 4c)
A 2 ^ = E * (a ^ + a b * D * r^) (D.5.4d)
2 2
f » a, f * E * D
Ao (D.5 .4e)
where
D - t? (1 - r ?
b * r (D.5.4f)
/ 3 )
and
E - t / (1 - r^ * r f - D * r^ * r * C 2> (D.5.4g)
x 2 3
There is a similar set of equations for radiation incident from the back. TARP
also uses expressions for quadruple glazing, which are not repeated here.
where
37
dQ / dAj ° energy per unit area
= 2 tt (r 8in0) r d0 / r2
tt/2 tt/2
Xd = / X(0)2 tt I sin0 cos 0 d0 / / 2 ir I sin0 cos0 d0 (D.6.2)
0 0
As - 1 * a ' *0 -
lf * c * a - r * R b> D-7 -2 >
s s s s I s d <
b
n
= A J
n
+ J*s * nd
,
- A L + A * Tf * r b
/ (1 - r e * Rj ) (D.7.3)
n nd‘
REFERENCES
38
[3] Born, M,, and Wolf, E., Principles of Optics, 4th ed, (Pergamon Press,
1970), pp 36-42.
[4] Heavens, 0. S., Optical Properties of Thin Solid Films , (Academic Press,
1955), pp 51-58.
[5] Born, M., and Wolf, E., Principles of Optics, 4th ed, (Pergamon Press,
1970), pp 627-632.
[6] Hottel, H. C., and Sarofim, A. F., Radiative Transfer, (McGraw-Hill,
1967), p 197, fig 5-12.
[7] Kreith, F., Principles of Heat Transfer , 3rd ed. (Intext Educational
Publishers, 1973), pp 230-231.
39
1.0
0.8
0.6
TRANSMITTANCE
0.4
0.2
A = DSA glass
B = 6mm (0.25 in) clear glass
C = 6mm (0.25 in) heat absorbing glass
40
Figure III.D.2 Reflection and refraction at an interface
41
2
(1-p) r l
f
2 2 3
|1-p) P r |l
2 4
(1-p) p r*l
42
Figure III.D.4 Two-pane window system
43
Intensity
in
44
shade
45
.
x zo = x br * cos
^b
“ Y br * sin ^b (E.l. la)
z zo = z br (E.l.lc)
The user is also permitted to rotate a zone z degrees relative to the i\>
x zo + X
A zr * cosi|> z - v
L * sinij;. (E.l. 2a)
so zv
so Y zo + x zr * sini|>z T
+ yI zr * cos4>
;
(E.l. 2b)
so z zo + Z zr (E.l. 2c)
46
.
= + ^z + (E.1.3)
^ 'I's ’J'b
Y = Y so + X sr * - Y * (E.1.4b)
sinij/
sr * cosi|> cos<J>
Z = Z so + Y sr * sin<{> (E.1.4c)
X* = X - X (E. 1 .5a)
so
Y’ = Y - Y (E. 1 .5b)
so
Z' = Z - Z (E . 1 .5c)
so
47
9 .
enclose a solid, and only those surfaces which are not sunlit at a given
hour are considered shadowing surfaces.
CW 3 = cos<{> (E 2 . lc)
The cosine of the angle of incidence of the sun's rays on the surface are
given by the dot product of surface and sun (B.3.3) direction cosines.
48
The sequence in the cross product is a convention to determine sign. The
condition that a point (x, y, w) lie on a line (a, b, c) is that
the point is to the left of the line, and if it is less than zero, the point
is to the right of the line. The intercept (x, y, w) of line (a^ , b^ , c^)
and line (a 2 , b 2 , 02 ) is given by
After transforming the shadows onto the plane of the receiving surface,
the basic job of the shadow algorithm is to determine the area of the
overlap between the polygons representing the shadows and the polygon
representing the receiving surface.
The vertices which define the overlap between two convex polygons, A
and B, consist of:
49
Coordinate transformation retains the order of the vertices of a
polygon, while a projection reverses the order. The sequence of vertices of
the receiving polygons should be reversed so it and all shadow polygons will
have the same sequence.
Once the vertices are determined, they must be sorted into clockwise
order for the area to be computed. Given a plane polygon of n sequential
vertices (x ls y^, (x 2 y 2 ), •••, (*n yn ), its area is given by [4]
, ,
If two shadows overlap the receiving surface, they may also overlap
each other as in figure E.6. The vertices of this overlap can be computed.
The areas of all overlaps can be computed. The total sunlit area can be
expressed as the sum of all polygon areas given a proper sign on each of the
areas. The following convention was adopted:
and so on through multiple overlaps where the sign of the overlap area is
the product of the signs of the overlapping areas.
It is easy to determine the sunlit area of a window once all the shadow
and overlap vertices on the wall have been computed. Consider wall 2 of
figure E.l. First the wall is considered a simple rectangle and the window
on it is ignored. The shadow overlapping is performed and the sunlit portion
of the gross wall area is computed. Then the window rectangle is overlapped
with the shadow to determine its sunlit area. The sunlit area of the window
is subtracted from the gross wall sunlit area to determine the net wall sunlit
50
area. During this calculation it is not necessary to recompute the shadows,
because they were precisely determined on the wall.
where
a = solar absorptance of the surface
9 = angle of incidence of the sun's rays
S = area of the surface
Ss = sunlit area
Ib = intensity of beam (direct) radiation
Is = intensity of sky diffuse radiation
Ig = intensity of ground reflected diffuse radiation
F S g = angle factor between the surface and the sky
Fgg = angle factor between the surface and the ground
F ss = (1 + cos<j>) / 2 (E.5.2a)
and
REFERENCES
51
building
shaped
L~
an
on
Shadows
,15.1
II
I
Figure
52
Sun’s rays Sun’s rays
53
I
54
polygons
convex
Overlapping
III.E.4
Figure
55
segments
line
of
Intercepts
III.E.5
Figure
56
Figure III.E.6 Multiple shadow overlaps
57
Ill .F INTERNAL HEAT GAINS
Internal heat gains from lights, people, and equipment of various types
are usually significant elements in the zone thermal balance. TARP allows
the user to specify heat gains from equipment of various types, lights,
occupants, and baseboard heat controlled by outdoor temperature. The total
heat gain is made up of convective, radiant, and latent gains in various
proportions from each of these sources. Convective gains are instantaneous
additions of heat to the zone air. Radiant gains are distributed on the
surfaces of the zone, where they are first absorbed and then released back
into the room (and some fraction conducted through the surface) according to
the surface heat balances (section I). Latent gains must be handled by venti-
lation or air conditioning equipment. Recommended heat gains for many items
of equipment are given by ASHRAE [1]. These recommendations include the sen-
sible (convective plus radiative) and latent proportions. Sensible gains from
equipment are predominately radiant. The user can specify the heat gains and
proportions appropriate for any type of equipment. The gains from lights,
people, and baseboard heat are slightly more complicated.
III.F.l. 2 PEOPLE
58
S “ 6.461927 + .946892*M + .0000255737*112
This function is compared to the original data points in figure F.l. The
radiant fraction of the sensible gain is determined by the user.
Some fraction of the beam solar radiation transmitted into the zone is
directly absorbed by the interior surfaces according to the solar distribution
algorithm (3.1.5) selected by the user. The beam radiation not directly
absorbed, plus the diffuse sky and ground-reflected radiation, plus the short
wavelength radiation from lights are combined and distributed over the surfaces
of the zone according to:
59
NS
OSIi - QS n * ot / E S ± * (1 - Pi) (F.2.1)
i=l
REFERENCES
60
people
from
gains
heat
Sensible
at
data
III.F.l
Published
Figure
61
OUTDOOR TEMPERATURE
62
Ill .G NATURAL AIR MOVEMENTS
2
PW 1 ’
= 0.5 * D
a * V (G. 1.1)
D a = 0.0034838*PB/T a (G.1.2)
The wind speed is evaluated at the highest point of the surface using the
wind speed to height correlation (B.3.1).
The pressure difference (AP) across the surface is Pj - Pi. It is more useful
to know this AP in terms of zone pressures Pm and Pn:
63
:
Pj = p
m + Dm * g * (Zm - Zj ) (G.2.1a)
The sane formulae apply to a vertical surface as in figure G.2. They also
apply to an envelope surface where Pj is the pressure of the outside air
(figure G.3). Ambient air variables are substituted for the zone ra
variables
Pj = P + Da * g * (Z a - Zj) (G.2.2)
a
AP = Pj - Pi
- Pm - Pn + ’ D *g*(Zn-Zi)
n
= P - P + PSi (G.2.3)
m n
where
PSi - D m *g*(Zm-Zj) - D n *g*(z n -z1 )
F=C*A*D*/2*AP/D (G.3. 1)
where A is the observed opening area and where the flow coefficient, C, is
a function of the Reynolds number and the ratio of the opening size to the
entire surface. When the opening is small, C equals 0.6 for a wide range of
Reynolds numbers [4]. TARP assumes this value of C. The user may modify C
by using an effective flow area (FA) which differs from the actual proportional
opening area. Experimental results indicate that flow through real surfaces
is a combination of orifice flow (Fa /AP ) and capillary flow ( F a AP )
with an exponent of about 0.65 [5]. In order to retain the concept of opening
area and provide a simple grouping of all other terras in (G.3.1), C is multi-
plied by /2 and D divided by /D to give:
64
Xi
F
±
« 0.848 * FA £ * /d * (AP) (G.3.2)
where FA-^ is the effective flow area of the openings in the surface, D is
the density of incoming air, and is the flow exponent. The pressure
difference for envelope surfaces is
X1
F
a
- 0.848 * FA a * /Dm * (AP) (AP > 0) (G.3.4a)
or
F
i
= -.848 * FA
t
* /D n * (~AP) Xi (AP < 0) (G.3.4b)
Fs + I Fi = 0 (G.3.5)
where F g is the net mass flow into the zone from the air handling system and
the F-i are the flows through openings in the zone envelope surfaces. These
mass balances form a system of non-linear algebraic equations which can be
solved for the zone pressures.
—
3F,
3P n
i = -X,
1
* F,
1 / AP (G.4.1a)
—
3Pm
3F
i = Xj
1
* Fj
1 / AP (interzone surface) (G.4 . lb)
or
3F
1 =
0 (envelope surface) (G.4. lc)
3Fm
65
These equations have been successfully solved by a modified Newton’s
method in all cases tested. A standard Newton's method was found to sometimes
give very slow convergence because of overcorrections during the iteration.
The modification reduces the Newton correction when this problem appears to
occur.
REFERENCES
[4] ASHRAE Handbook - 1981 Fundamentals (ASHRAE, 1981), p 4.13, Fig 18.
[5] ASHRAE Handbook - 1981 Fundamentals (ASHRAE, 1981), p 22.8.
66
Xi
F
1
= 0.848 * FA t * /D * (AP) (G.3.2)
where FA-^ is the effective flow area of the openings in the surface, D is
the density of incoming air, and is the flow exponent. The pressure
difference for envelope surfaces is
AP = PWi + ?S ± - P n (G.3.3a)
FAt Xi
F
t
= 0.848 * * /Dm * (AP) (AP > 0) (G.3 .4a)
or
= -.848 * FA ± * /D n * (-AP) Xi (AP < 0) (G.3. 4b)
Fs + Z Fi = 0 (G.3.5)
where F s is the net mass flow into the zone from the air handling system and
the F-^ are the flows through openings in the zone envelope surfaces. These
mass balances form a system of non-linear algebraic equations which can be
solved for the zone pressures.
3F,
—
8P
1 = -Xj
1
* F, / AP (G.4. la)
n
—
3F •
i = Xj
1
* Fj
1
/ AP (interzone surface) (G.4 . lb)
or
3F
—
3p
1 = 0 (envelope surface) (G.4. lc)
m
65
These equations have been successfully solved by a modified Newton's
method in all cases tested. A standard Newton's method was found to sometimes
give very slow convergence because of overcorrections during the iteration.
The modification reduces the Newton correction when this problem appears to
occur.
REFERENCES
[4] ASHRAE Handbook - 1981 Fundamentals (ASHRAE, 1981), p 4.13, Fig 18.
[5] ASHRAE Handbook - 1981 Fundamentals (ASHRAE, 1981), p 22.8.
66
Zm
ROOMm
Zj
zzzzzzzZzzzzzzzz
Z|
Pi
ROOM n
z„
*Pn
Zm
Pm
ROOM n I I
Z|,Zj
ROOMm
•Pn I Z„
Pi*^*pj Zi,Z,
ROOM n AMBIENT
P.
1 Z,
67
Ill .H ROOM RADIANT INTERCHANGE
4 -
QI i = E F
ii
*(T.j T^) + QLI ± + QS^ (H.1.1)
j
This expression is not usually substituted into the surface heat balance
equations since the T 4 terms would make the equations nonlinear and difficult
to solve. The radiation interchange factor, F-jj (which includes geometry and
radiocity factors) is difficult to calculate for rooms with arbitrary geometry.
Since each surface interacts with all other surfaces of the room (total of N
surfaces) there are N^ interactions to be considered. This requires a time
consuming solution technique. Quicker solutions are possible by linearizing
and assuming constant radiative and convective coefficients for all surfaces
[1]. Alternate methods are available which use mean radiant temperature
instead of the detailed radiant interchange. These methods allow efficient
treatment of the problem without assuming constant coefficients. They so
significantly affect the solution time for problems with many surfaces that
they make simultaneous heat balances for many rooms feasible [2].
Since the MRT includes the effect of all surfaces, coupling the surface
to the MRT couples it to itself. The MRT network method includes a factor
F^, the "MRT view factor", which compensates for the self weight of surface
i in the MRT.
68
HR t = A*a*T 3 /[l/F 1 + ( l-e ± )/e ± ] (H.2.3)
3. Then evaluate
This method scales each HR^ by two factors nearly proportional to T** 3 and
,
estimates the radiant fluxes to within about 1 percent for surfaces between 0C
and 50C.
Shortwave radiation, such as from the sun and lights, should not be handled
by this modification of the MRT because the room surfaces can have entirely
different emissivites for short and long wavelength radiation. Therefore,
H.2.1 should be written:
The MRTN method of Carroll is slightly more accurate than the "balanced
mean radiant temperature" (MRTB) method of Walton. In addition it can
easily handle multiple reflections and air emittance effects. However, when
MRTN is used in the TARP iterative heat balance, it is slower than MRTB A.
single MRTN iteration is slightly faster than MRTB, but more are required.
This is because MRTN updates TR n once per heat balance iteration of all
69
.
surfaces while MRTB updates TR n for each surface during the heat balance
iteration. The effect is the same as Jacobi versus Gauss-Seidel iteration [4]
The MRTN method was chosen over MRTB because of its accuracy and potential for
future expansion. It thus became important to find a method to reduce the
number of iterations to regain speed similar to MRTB. This method is described
in section L.
REFERENCES
70
I
The convective plus radiative heat flux from a surface to the environment
may be represented by:
The surface coefficient, HO, accounts for both convection and longwave radiant
interchange. QSO is the absorbed shortwave radiant flux. A simple correla-
tion for HO uses a polynomial expression [1] based on an ASHRAE relationship
[2] which relates the surface coefficient (W/m^K) to wind speed (m/s) and
surface roughness.
H 0 - *0 + A l* V a + V v
a
2
0.1.2)
71
.
V v
aZ = a * b
* < z/z o> a 0.2.1)
H If - 0.86*D
a
*C*P L " 2 / 3 *(.25*V cLZ *C*v/A)!/ 2 (1.2.3)
p
where
For air Pr~ 2/3 = 1.245, C - 1004 J/Kg K , and D a iTv 5 0.00472 Kg/s m 5/2
p
which combine to give
The local value of Hf was found to be greatest near the edges of the surface
where the thermal boundary layer is thinnest. Therefore, to avoid computing
too large a value for small subsurfaces (windows and doors), TARP uses the
same C/A for the subsurfaces as for the wall they are on.
A simple surace roughness multiplier, Rf, was developed from the ASHRAE
graph of surface conductances [6]. At 0 mph all lines on this graph intercept
at about 1.9 Btu/hrft 2 F. The correlations for the six roughnesses are nearly
linear. Therefore, a simple roughness multiplier was constructed by subtrac-
ting 1.9 from all values of H (at 25 mph) and dividing by the quantity H [
72
Roughness H (at 25 mph) Rf
1 13.6 2.1
2 10.8 1.67
3 10.1 1.52
4 8.0 1.13
5 7.9 1.11
6 7.3 1.00
1 1 1. 1. 2. 2 NATURAL CONVECTION
Heat flow up: H 9.482* 3 >/| TO-Ta| /(7 .238- cos (I. 2. 5a)
n - |
4> |
There is radiant interchange between the surface and the sky, and between
the surface and the ground surface or adjacent structures. The effective
temperature of the ground surface is assumed to be equal to T a (pending further
algorithm development). The effective temperature of any point in the sky
depends on the optical path length through the atmosphere. Section B.3 des-
cribes the algorithm which computes the effective sky temperature which would
be seen by a horizontal surface. For tilted surfaces the radiant interchange
coefficient (H r ) must be divided into three parts: one (HS) to sky temperature,
one (HG) to ground surface temperature, and one (Ha) to air temperature.
4
HS = a * e * F * cos(<j>/2)* (TO 4 - T o )/|TO - T |
(1. 2. 6b)
^o
Pending further research, TARP assumes the ground surface temperature, Tg,
equals the air temperature, T a This three component split of H r must still
.
73
4
Since the expressions for Hn and Hr are used to evaluate TO, the most
recent value of TO is used in their calculation. The three factors combine
to give the following expression for combined convective and radiative flux.
KO + CO + QO = 0 (1.3.1)
In contrast to the internal surface heat balance which treats all surfaces
simultaneously, the external thermal balance for each surface is performed
independent of all other surfaces. This implies that there is no direct
interaction between the individual surfaces.
The TARP includes four possible representations for the basic outside
surface heat balance. The first two depend on which of the optional surface
conductance algorithms the user selects. The simple outside surface conduc-
tance (1.1) which includes both the convective and thermal interchange between
the surface and the environment in a single coefficient, is represented by
the thermal network in figure 1.1. (1.3.1) can be expressed in terras of
(C.1.6) and (1.1.1) as
74
The third and fourth representations occur when the outside surface has
been covered with a movable insulation. The insulation has a conductance of
UM. This case is represented by the thermal network in figure 1.3. The
insulation must be massless because it is not generally possible to perform a
correct thermal balance at the juncture of two surfaces each modeled by CTF.
The equation for the thermal balance betwen the surface and the insulation is
In this case the values of HA, HS, and HG must be found by using an estimated
value of TM in place of TO in (1.2).
75
] ] ] ]
Fj_ = Y
0 / (Zq + HI + HR) (I. 4. 2a)
F2 = UM / (UM + HO) (I. 4. 2b)
F3 = UM / (UM + HA + HS + HG) (I .4 .2c)
/ (X 0 + HA + HS + HG) (1.4.5)
76
CASE 7: Yq small, detailed conductance, with movable insulation. From (1.3.7)
and (1.3.9):
REFERENCES
[1] T. Kusuda, "NBSLD, the Computer Program for Heating and Cooling Loads in
Buildings," NBS Building Science Series 69, (July 1967), p 66a.
[2] ASHRAE Handbook - 1981 Fundamentals (ASHRAE, 1981), p 23.2, Figure 1.
[3] R. J. Cole and N. S. Sturrock, "The Convective Heat Exhcange at the
External Surface of Buildings", Building and Environment, Vol 12 (1977)
pp 207-214.
[4] E. M. Sparrow, J. W. Ramsey, and E. A. Mass, "Effect of Finite Width
on Heat Transfer and Fluid Flow about an Inclined Rectangular Plate,"
Journal of Heat Transfer , Vol 101 (May 1979), p 204.
[5] R. J. Cole and N. S. Sturrock, "The Convective Heat Exchange at the
External Surface of Buildings", Building and Environment, Vol 12
(1977), p 212.
[6] ASHRAE Handbook - 1981 Fundamentals, (ASHRAE, 1981), p 23.2, Figure 1.
77
Figure III.I.l Thermal network for simple outside surface coefficient
Tl
78
III.J SURFACE INSIDE HEAT BALANCES
The convective heat flux from the zone air to a surface is given by:
The coefficient for natural convection heat transfer is based on data from
ASHRAE [1], The ASHRAE values are for combined convection and radiation, so
a radiative component of (1.02 times emissivity of 0.9) is subtracted from
the tabulated values which are then converted to metric. When warm air is
below a cool horizontal surface (or cool air above a warm surface), there is
more convective mixing than for air next to a veritcal surface. Conversely,
when cool air is below a warm surface (or warm air above a cool surface)
there is reduced convection.
When the air handling system is moving air through the zone, a single value
of 6.14 W/m 2 K (= 2.0 - 1.02 * 0.90 Btu/h*ft 2 *F) should be used for HI.
ASHRAE has the following expressions for natural convection from large
surfaces [2].
1.31* 3
/| TZ-TI |
Horizontal (J.2. lb)
It may be noted that these expressions are inconsistant with the values
given for the simple algorithm.
To account for all surface tilts, the following expressions are used.
HC = 9.482 * 3
/| TZ-TI |
/ (7.238- |
cos<|) |
) (J.2. 2a)
79
)
The value of TO^ must have been previously calculated by the methods in
section 1.4.
The following two partial derivatives are useful in improving the rate
of convergence of the heat balance iterations.
3TI t (J.3.4a)
= HCi / (HCi + HR i + Zo.l)
3TZ
n
3TIi (J.3.4b)
= HRi / (HCi + HRi + Z0,i)
3TR
n
3TIi ( J.4 . 2a
HCi / (HCi + HRi + Z 0j i Y 0 ,i)
3T Z
n
80
) )
and
“
3TX
= HR
i 1 (HC 1 + HR 1 + ZQ.i " Y 0,i)
(J.4.2b)
HC i *(TZ n -TM 1 ) + HR i *(TRn -TM1 ) + UM1 *(TI 1 ~TMi ) + QS!^ = 0.0 (J.5.1)
A similar set of equations exist for the far side of the surface (described
as surface j adjacent to zone m) . The far side heat balance (without
movable insulation) is
KIPj + Y 0) j*TOj - Zq
^ j
*TIj + UMj *(TMj-TIj ) + QSIj =0.0 (J.5.5)
- X
Zo.j 0>1 (= X 0 ); Y0>j - Y 0>1 (= Y 0 ).
81
Tit = [HCi^Zn + HRi*TRn + QSIi + KIPi +
and
3TI,
)Yf - HCj
,
/
,
<HCi + HR i + ZO - W (J. 5 .8a)
3TL
^ = H RjL / (HCi + HR i + Z0 Fi*Y 0 )
(J.5.8b)
3TX
—i = HCi / (HCi + HRi + UMi)
(J.5 . 10a)
n
and
—
3TI
9T R
i
n
- HRi / (HCi + HRi + UMi)
(J.5. 10b)
The partial derivatives are given by (J.5. 8a) and (J.5. 8b).
CASE 4: Movable insulation on both sides. QSIi and QSIj are both zero
because the insulation is assumed opaque. From equations (J.5.1), (J.5.2),
(J.5.4) and (J.5. 5):
82
}
The partial derivatives are given by (J.5. 10a) and (J.5. 10b).
REFERENCES
83
TZ
84
Ill .K ZONE AIR HEAT BALANCE
The summations apply to all surfaces which are part of the zone.
The system performance model (III.K.l) may be used with the detailed
natural air movement model (III.G) if the supply and exhaust air flows of
every zone are equal (that is, Fs in equation (G.3.5) must be zero). The
simple air movement models (III.K.2 through III.K.5) may not be used
siraulatneously with the detailed air movement model.
The energy from the system is based on the user defined control profile
which approximates the performance of the actual system. The control
profile consists of a set of control points which define the heating or
cooling capacity of the system as a function of the control temperature as
in figure K.l. Control range p is between temperatures T and Tp+j . The
p
control profile is linear between the control points. Simple
proportionality gives
L = LT * T + LC (K.l. 2a)
where
The control temperature consists of both the air and the mean radiant
temperatures according to the formula
T = 0 * TR + (1 - 0) * TZ (K.l. 3)
85
III.K.2 SIMPLE INFILTRATION MODEL
An equation for the mass flow of infiltration air can be developed from
measured values in similar structures [1]. The TARP assumes a correlation
of the form
FI - (A0 + A l *|TZ - T a |
+ A 2 *V + A 3 *V 2 ) * D
a
(K.2.1)
LV = FV * C (T a - TZ) (K.3.1)
p *
86
are intimately coupled, it is best to seek an exact simultaneous solution.
The description of that solution is simplified by rearranging equation
(K.2.2) to
B n * TZn - ^ (K.4.1)
where
Bn = E HCi + FI n *C pn - (1 - B)*LT n
A similar equation exists for zone m. The effect of an interzone air flow
is given by the two equations
where MM = FM * C .
p
There are three possible conditions for the interzone air flow. First
is that TZjh - TZ n < TD and therefore MM = 0. Second is that the maximum
flow rate still leaves T - TZ n > TD. Simultaneous solution of equations
(K.4. 2a) and (K.4. 2b) gives
where Fj = MM / (B m + MM)
The third case occurs when no flow gives TZq - TZn > TD and full flow would
cause TZ m - TZ n < TD. It is then necessary to solve for the partial flow
during the hour which will cause TZtn - TZn = TD. This gives
A model similar to the zone mixing model (III. K.4) can be developed
which describes the effects of a fan moving air only from zone n to zone m
with movement of air to and from ambient completing the mass balance. Such
a case is the exhaust of house air into its attic. A simple control scheme
allows a mass flow of up to FF whenever TZ n > TD > T a . The effect of the
air flow on the zone energy balances is given by
87
(B n + MF)*TZn - An + MF*TA (K.5.1a)
and
where
MF = FF * C
p
There are three possible conditions for the interzone air flow. First
is that TD < T a or TZ n < TD which gives MF=0. Second is that the maximum
flow leaves TZjB > TD. Direct solution of equations la and lb gives
The third case occurs when no flow gives TZn ^ TD anc^ full flow would cause
TZ n < TD. It is then necessary to solve for the partial flow during the
hour which causes TZn = TD.
REFERENCES
88
Figure III.K.l Zone temperature control profile
89
Figure III.K.2 Ventilation control profile: T a higher than TV
90
Figure III.K.5 Ventilation and cooling control profiles
91
Ill .L HEAT BALANCE CONVERGENCE
In any zone the air temperature, mean radiant temperature, and the
surface temperatures are intimately related. A change in any one will cause
changes in all of the others. Those secondary changes should be predictable.
It was observed that, with the zone air temperature held nearly
constant by a large heating capacity, there was a regular pattern in
successive radiant temperature values during the heat balance iterations for
a given hour.
where
with i for all surfaces of the zone. The partial derivatives are described
in section J.
The following logic led to a method for improved convergence for fixed
or floating zone air temperature.
AT Zq = ATR 0 * E HC i *S i /
9TI
ATR X = ATZq * E HC *S *
i i
i / E HC^S 1
8TZ n
92
aTI
= ATRq * E IIC
i
*S i * i /
8TZ n
[Z HCi^Si + E Fj*Cp - (l-e)LT ]
= Y * A TR 0 (L.1.4)
(6) Combine the changes in radiant temperature in (3) and (5) to give the
total change in TR for one iteration.
= ATRq / (1 - a - y) (L.1.6)
2
ATZ = [ATRq * E HC i *S i )
* [ (crHy) + (cc+y) + ... ] /
[E HC i *S 1 + Fj *C (1-S)*LT (L.1.7)
p -
E ]
(8) The current estimates of TR and TZ are therefore increased by ATR and
ATZ, respectively. Then the surface temperatures and resulting new zone air
temperature are calculated in the normal manner. For the convergence test
the new TZ is compared to the predicted value of TZ instead of the previous
iteration value. If the two temperatures (or associated heating/cooling
load) are sufficiently close, the solution has converged.
93
IV. INPUT
IV. A OVERVIEW
This input method was chosen from the following four basic
alternatives.
94
.
an English description of the same features for those who are unfamiliar
with BNF.
9. <REAL> <INTEGER> |
<INTEGER> . . <UNS INTEGER> |
|
15. <LETTER> :
A |
B |
... |
Y |
Z
18. <C0MMENT> : := $ |
) |
] {<CHAR> |
<SYMB0L> |
<NULL>} END-OF-LINE
1. The TARP input deck consists of several sections. There are presently
four sections: PROJECT, LIBRARY, BUILDING, and ZONE.
2. Each section consists of a section name and one or more keys enclosed in
brackets
95
3. A key is a key name and one or more subkeys enclosed in parentheses.
5. There are five different data types. Only descriptions may be longer
than 10 non-null characters.
9. There are five distinct representations for real numbers. Note that a
decimal point is not necessary and that exponential notation is possible.
13. There are several generally recognized symbols in FORTRAN 77. These
symbols are used as punctuation in the input language. Most systems also
provide other symbols which may be used in those language elements which use
"symbols" such as ", %, @, etc.
14. The blank and plus sign are nulls. They are ignored by the input
reader. Therefore, elements of the language may be separated by blanks to
improve readability. The plus sign does not affect the value or character
count of a number.
18. Comments begin with a dollar sign or at the end of a key or end of a
section and continue to the end of the line. They may be included on lines
with input or on separate lines.
96
Other features:
.
Successive identical subkey values " . . . <DATA> <DATA> <DATA>
, , , , . . " may
.
be replaced by " . . . ,<UNSIGNED INTEGER>*<DATA> , . . " .
The final subkey slash before ending a section may be dropped. That
is, )" may be replaced by ")" .
"/ The slash may be considered a separator
rather than a terminator.
The TARP handles incorrect input by reporting the error and then
skipping to the next readable input. It makes no attempt to correct the
input error. It is more important for the user to correctly type
punctuation than names or values. The approximate order of importance of
punctuation is: '$[]()/=,. Similarly, section names are more
important than key names which are more important than subkey names which
are more important than data. For example, an incorrect data value would
cause the program to not process only that value, but an incorrect subkey
name causes the program to skip all values following it (because the program
cannot know to which subkey they apply) and go to the next subkey, which
will be after the next "/"• Similarly, an incorrect key name causes the
program to skip the entire key including all its subkeys. The key name will
have to be corrected and another run made to find any errors in the subkeys.
The punctuation is critical because it determines where the input skipping
stops. The program must then find a recognizable name or another skip will
occur. Thus the brackets and parentheses are more important than slashes
because an error in them causes more skipped input. The apostrophe and
dollar sign are most critical because descriptions and comments can include
punctuation symbols which could totally confuse the input reader. Although
this process sounds complex, the consistancy of the syntax makes it very
easy for the user to check punctuation before submitting the input to the
computer. The input deck can be debugged in a small number of runs because
the entire deck will be processed during every run instead of quitting at
the first error.
97
IV. INPUT DETAILS
1. The section, key, and subkey numbers (e.g. 2.4.1). These serve in
cross-referencing to other input. Several related subkeys may be grouped
under one subkey number.
2. The subkey name, "", data types. The data types mean: n * name, r =
real, i = integer, d = date, desc ° description.
4. The units. means dimensionless. English units are given below the
metric units.
5. The status, "reqd" means a value must be given, "opt" means a value is
optional; otherwise the default will be used, "cond" means conditional.
One of several values must be chosen or a value must be chosen if some other
input has been set.
98
1.1.4 2SPC no - opt
This will cause double spacing. A blank line will be inserted between
each user input line which is echoed in the input processor report.
This enhances readability.
OLDW no - opt
This subkey may be used if a environment file (ENVTFL) has been prepared
1 and saved in a previous run. It causes keys 1.3, 1.4, 1.5, and 1.6 to
1.1.7
be skipped, saving the time that would be spent preparing a new
identical environment file.
99
1.3 LOC ( Location.
Ground temperatures for the 12 months starting with January. They are
used only with other side coefficients (4.3).
100
.
SNOW no - opt -
SNOW causes the ground reflectance to be increased from its default of
0.2 to 0.7. This gives a much higher diffuse solar gain to the zones.
101
IV. C. LIBRARY SECTION
NDS = i 16 - opt -
NWS = i 16 - opt -
NMATL = i 16 - opt -
NCONS « i 16 - opt -
NCPFL = i 16 - opt -
NCDS = i 16 - opt -
NCWS i 16 - opt -
These subkeys control the number of library elements which can be
defined. Library names are stored in an array which is 256 elements
long. The total number of named library elements cannot exceed 256.
The sum of these subkeys (including those left to default) cannot
exceed 256.
i
ALL no - opt -
DS no — opt -
WS no - opt -
MATL no - opt -
CONS no - opt -
CPFL no - opt -
CDS no - opt -
CWS no - opt -
These eight subkeys control the printing of library data. The eight
individual names will each produce reports for that class of library
data. ALL will cause all eight classes to be reported.
2.2 DS ( Day-schedules
A day-schedule is a set of 24 numbers used to provide an hourly
variation in some building parameter.
102
2.2.2 DESC desc’ * blank - opt
A description will more fully describe the purpose of the schedule. It
will appear in the day-schedule report.
2.2.3 FFCr,r,...<24>
=* 24*1.0 - opt 0.<“ROl.
Fraction of full capacity. 24 values are required: one for each hour
of the day. They may range from zero to one. A zero value can be
replaced by OFF and ON can be used in place of one. This may be more
intuitive for day-schedules which are controlling the on/off operation
of something.
2.3 WS ( Week-schedules.
A week-schedule is a set of 12 day-schedules: one for each day of the
week, one for holidays, and one for each of 4 special day types (SP1,
SP2,
SP3, and SP4). It is referenced by name in subsequent building
input
2.3.2 -
2.3.1 NAME = n blank - reqd
Each week-schedule must have a name different from all other
week-schedule names.
2.3.3
DESC = 'desc' blank - opt -
A description is useful in more fully describing the schedule, but its
use is optional.
WD = n blank - cond
Set weekday (Monday through Friday) pointers to the named day-schedule.
That day-schedule must have been previously defined (2.2). Other
options are:
WE = n Set weekend (SAT & SUN) pointers.
ALL = n Set all 12 pointers.
SUN = n Set Sunday pointer only
MON = n Set Monday pointer only,
and similarly for TUE = n, WED = n, THU = n, FRI = n, SAT = n, SPl = n,
SP2 = n, SP3 = n, and SP4 = n. At least one day must be defined (other
pointers will be set by default rules).
103
2.4.4 CP - r .8368 kJ/kg*K cond 0.<r
0.2 Btu/lbm*F
Specific heat.
104
.
GLASS - - cond
This subkey indicates the glass panes in a transparent construct. A
GLASS material must have conductive properties and optical properties
chosen from: (TRNS, IR) or (TRNS, IR, FTRN) or (SC). The program
assumes that multiple panes are separated by air layers for the optical
calculations
This subkey describes the material layers which form the construct.
The values are pairs of (thickness, name of material) from the outside
to the inside layer. Layers which are described by R-value (2.4.6)
must be given a thickness to maintain the proper sequence of values.
The value of that thickness is not used in the conductivity
calculations, so even a zero is adequate.
105
2.5.6 ORGH - i - - opt l<-i<-6
This describes the roughness of the outside surface. It has the values
of subkey 2.4.14. If no value is given, the default is the roughness
of the outer layer.
106
2.8 CWS ( Control week-schedules.
A control week-schedule is a set of 12 control day-schedules: one for
each day of the week, one for holidays, and one for each of 4 special
day types (SP1, SP2, SP3, and SP4). It is referenced by name in
subsequent building input.
WD = n blank - cond -
Set weekday (Monday through Friday) pointers to the named control
day-schedule. That control day-schedule must have been previously
defined (2.7). Other options are:
WE = n Set weekend (SAT & SUN) pointers.
ALL = n Set all 12 pointers.
SUN = n Set Sunday pointer only
MON = n Set Monday pointer only.
and similarly for TUE = n, WED = n, THU = n, FRI = n, SAT = n, SP1 = n,
SP2 = n, SP3 = n, and SP4 = n. At least one day must be defined (other
pointers will be set by default rules).
107
3 .
108
.
3.1.6 HO = i 0 - opt i = 0 or 2
Compute exterior convection coefficient using simple (0) or detailed
(1) algoritrah (1.1 and 1.2).
3.1.7 HI = i 0 - opt i = 0 or 1
Compute interior convection coefficient using simple (0) or detailed
(2) algorithm (J.l and J.2).
X = r 0. m opt -
0. ft
Y = r 0. m opt -
0. ft
Z = r 0. m opt
''
-
0. ft
ORG = r,r,r
Coordinates of lower left corner of surface relative to building
coordinate system. ORG provides a simpler way of describing all 3
coordinates at once (ORG = X, Y, Z) . See IV. F on geometric data
3.2.4 W = r 0. m reqd -
0. ft
H = r 0. ra reqd -
0. ft
SIZE = r,r
or
VRTS = r,r,r,r or r,r,r,r,r,r
The first two forms specify the width and height of a rectangular
surface (SIZE = W, H). The last form is for triangles or quadri-
laterals. The first pair of numbers are the surface relative x-y
coordinates of the first vertex counter-clockwise from the lower
left vertex. The subsequent pairs of numbers are the coordinates
109
. F
The surface may apply for only part of the year, Any dates may be
used, but for greatest efficiency dates should be chosen from the
shadowing surface dates table.
110
IV. C. ZONE SECTION
X = r 0. m opt
0. ft
Y = r 0. m opt
0. ft
Z = r 0. m opt -
0. ft
ORG = r,r,r
4.1.4
Coordinates of the lower left corner of the zone relative to the
building coordinate system. ORG provides a simpler way of describing
all 3 coordinates at once (ORG = X, Y, Z) . See IV. F on geometric data.
Zone origin is necessary only for the more detailed shadowing
4.1.5
calculations (SLDS > 0; sec 3.1.2).
MIRX no - opt
This key causes all surface x-vertices to be mirrored about the y-axis
(north-south axis).
MIRY no - opt -
This key causes all surface y-vertices to be mirrored about the x-axis
(east-west axis).
Ill
.
112
4.3.7 FROM » d 1JAN - opt
THRU = d 31DEC - opt
The lighting week-schedule may be applied to only part of the year.
The schedules used by other LIT keys may overlap. This allows a simple
way to combine lighting from different sources as might be necessary to
describe area and task lighting.
4.4.1 WS = n - - reqd -
The equipment is selected from the library of week-schedules (2.3).
The schedule refers to clock time.
EL cond
electricity
GS gas cond
ST steam cond
HW hot water cond
OTHR other cond
These simple subkeys identify the energy source of the equipment. They
allow summing of the building energy requirements by energy type. One
of them must be given. OTHR allows a load to occur in the zone without
adding to the energy requirements.
113
speed dependencies. It may not be used with the detailed air flow
option (3.1.4).
4.5.1 WS = n - - reqd
The infiltration schedule is selected from the library of
week-schedules (2.3). The schedule refers to clock time.
C = r .606 - opt -
T = r .03636 1/C opt -
.0202 1/F
V = r .1177 s/m opt -
.000598 min/ ft
>> ii
j-i 0. (s/m) 2 opt -
0. (min/f t) 2
114
4.6.4 MAX = r 90. C opt
This causes ventilation to occur whenever the zone air temperature is
greater than this value. If this subkey is used, do not use MIN
(4.6.3).
4.7.1 WS = n - - reqd
The mixing schedule is chosen from the library of week-schedules (2.3).
The schedule refers to clock time.
4.8.1 WS = n - - reqd
The fan schedule is chosen from the library of week-schedules (2.3).
The schedule refers to clock time.
115
.
4.10.1 WS - n - - reqd -
The fraction of beam sunlight which is usable for lighting is modified
by the week-schedule (2.3) element for the hour (standard time).
116
4.10.2 BEAM “ r 0 opt O.Or<“l
This is the fraction of beam sunlight entering the zone which can be
used to offset the visible energy from the lights.
4.11.1 WS = n - - reqd
The peak heating capacity avalable at any hour (clock time) is modified
by the baseboard heating schedule selected from the week-schedule
library (2.3).
4.11.6 ST cond
HW cond
These refer to the source of heating energy: steam or hot water. One
or the other must be given.
117
4.12 1 CONS = n - - reqd
Select the composition of the mass from the constructs library (2.5)
(2.5).
4.13 EX - -
cond -
IN - -
cond -
OP - -
cond -
TR - -
cond -
BS - -
cond -
SS - -
cond -
The user must specify the following three factors: Is this an EXterior
or INterior surface? Is the surface OPaque or TRansparent? Are its
coordinates defined as a Base Surface or a Subsurface? (The
coordinates of a subsurface are defined relative to the most recently
described base surface.) There are eight combinations of these
factors. Three of those are not permitted because of various
algorithmic restrictions in the program. The following table explains
the three factors in terras of the common English names for such
surfaces (at three basic tilt angles).
118
4.13.5 W » r 0. m reqd 0.<r
0. ft
H « r 0. ra reqd 0.<r
0. ft
SIZE ® r,r
VRTS *= r,r,r,r or VRTS r,r,r,r,r,r
The first two forms specify the width and height of a rectangular
surface (SIZE = W, H). The last form is for triangles or
quadrilaterals. The first pair of numbers are the surface relative x-y
coordinates of the first vertex counter-clockwise from the lower left
vertex. The subsequent pairs of numbers are the coordinates of the
remaining vertices in counter-clockwise sequence. See IV. F on
geometric data.
X - r 0. m opt -
0. ft
Y - r 0. m opt -
0. ft
Z - r 0. m opt -
0. ft
ORG = r,r,r
For base surfaces (BS - 4.13.2) these are the coordinates of the lower
left corner of the surface relative to the zone coordinate system. The
base surface origin is not used in the simple shadow calculation mode
(3.1.2). For subsurfaces (SS - 4.13.2) these are the coordinates
relative to the lower left corner of the base surface. The Z
coordinate is the amount of the setback. ORG provides a simpler way of
describing all three coordinates at once (ORG = X, Y, Z) . See IV. F on
geometric data.
119
4.14 AFS ( Air flow surface - Algorithm G
AFS describes the crack area and other openings in the base surface
which are used in the detailed air movement simulation. FLOW (3.1.4)
must have been specified. The base surface must be either an exterior
surface or an interzone surface (4.13.2 and 4.13.9). The AFS assumes
the azimuth of the base surface for computing wind direction effects.
When the AFS applies to an interzone surface, there must be an
identical AFS specified for the appropriate base surface in the
adjacent zone.
4.14.3 ZS = r - m opt
ft
This is the height of the centroid of the surface or opening relative
to the zone origin (G.2). The default value is the average of the
heights of the vertices of the base surface.
4.14.4 ZT = r - m opt
ft
This is the height at which the wind speed is evaluated for the opening
(G.l). The default value is the highest vertex of the base surface.
120
4.15 3 X = r 0. m opt ~
0. ft
>* it
0. m opt
0. ft
N II
0. m opt 0 ,<=r
0. ft
ORG = r,r,r
The positioning of the shadowing subsurface is given relative
base surface coordinate system. See IV. F on geometric data.
121
.
4.16.4 CT = r 0. C opt
0. F
Constant temperature.
4.17.2 DS = n - - cond
SNUP - - cond -
SNDN - - cond
One of three methods must be chosen for control of the EMI. The first
is reference to a day-schedule (2.2). A value of zero means no EMI for
that hour (standard time) and a one means the EMI is in effect.
Fractional values in the day-schedule modify the R-value. Any schedule
value greater than zero causes sunlight to be absorbed on the EMI
rather than on the surface. Second is SNUP which means the EMI is in
effect whenever the sun is above the horizon. SNDN means the EMI is in
effect when the sun is below the horizon.
122
4,18 IMI ( Interior movable insulation.
Interior movable insulation applies to interzone surfaces only. That
is, the most recently defined surface should be interior (IN; 4.13.2)
and have an other side zone name (4.13.9).
4.18.2 DS = n - - cond
GAIN - - cond
LOSS - - cond
One of three methods must be chosen for control of the IMI. First is
reference to a day-schedule (2.2). A value of zero means no IMI for
that hour (clock time) and a one means the IMI is in effect. Fraction
values in the day-schedule modify the R-value. The second method is
GAIN which causes the IMI to be in effect whenever the surface
temperature is greater than the zone control temperature (4.9.4). LOSS
is the inverse of the GAIN operation.
123
IV. D INPUT SUMMARY
Data types:
n = name, r » real, i » integer, d date, desc description
1 . PROJECT [
2 . LIBRARY [
124
2.6 CPFL ( NAME - n / DESC - 'desc' / TMP = r,...<8> / CAP - r,...<8> )
3 . BUILDING [
4 . ZONE [
FROM = d / THRU = d )
125
4.12 MASS ( CONS - n / AREA - r )
126
E
Point 2 might occur on a mild day in the spring or fall. At this point
(73 F [22.8 C]), no room heating or cooling is accomplished by the building
fan system (although heating or cooling energy may be consumed to condition
ventilation air). For the heat balance to occur at point 2, convective heat
gains caused by such factors as lights or occupancy must be balanced by
convective heat losses to room surfaces or by infiltration. The load at
point 2 is zero.
Point 3 shows a cooling condition. Here heat gains to the room air are
balanced by cooling capacity available at the balance point temperature.
The tempertures and capacities of the control points must be given in the
sequence shown. Very steep profiles are possible. The user may set two
successive control temperatures equal to each other to hold the room
temperature constant. TARP will make the second 0.01 C greater than the
first to permit the air heat balance calculations.
127
,
System type can subtly affect room temperature control profiles. For
example, assume that a package DX (direct expansion) condensing unit with
electric heat has been selected to serve the zone whose control strategy
is shown in figure E.l. Futher assume that the unit delivers 3850 cfm
(1.82 m^/s) to the space, that the cooling coil was selected to allow 55 F
(12.8 C) air to be delivered at the design sensible load of 100 kBtu/hr
(29.3 kW), and that the heating coil was selected to deliver 86.2 F (30.1 C)
air at the design heating load of 80 kBtu/hr (23.4 kW). In this case,
figure E.l closely approximates the achieved room temperature control. The
room thermostat will modulate the condensing unit from "off" to "full" capa-
city through as many capacity control steps as are available on the unit.
This modulation will occur in a roughly linear fashion as the balance point
room temperature varies from 77 to 79 F (25 to 26.1 C). The heating will be
similarly modulated between 69 and 67 F (20.6 and 19.4 C) and both heating
,
and cooling will be "off" between 69 and 77 F (20.6 and 25 C). Since the
condensing unit's maximum heating and cooling capacity are approximately
constant, the unit adds or removes fixed quantities of heat whenever the
room temperature is below 67 F (19.4 C) or above 79 F (26.1 C).
Figure E.3 shows the equivalent control profile that should be used if
a three-deck multizone unit is selected to serve several zones with the same
thermostat throttling ranges and deadband. The rationale for selecting this
profile is:
2. In range B-C, air from the hot deck and the bypass deck of the fan
system are mixed to meet zone heating requirements. This is the heating
throttling range for the thermostat. At its lower end (67 F [19.4 C]), all
air to the zone comes through the hot deck and is delivered to the zone at
the hot deck terapeature (86.2 F [30.1 C] for example).
128
3. If heat losses are severe, the room temperature may fall into range
A-B. In this range, air entering the room comes from the hot deck only. As
long as the fan system heating coil is not overloaded, the air is supplied
to the zone at constant temperature. However, the amount of heat added to
the room air continues to increase as the room air temperature falls. The
amount of heat added per hour can be calculated from the following formula:
L = F * Cp * (Ts - TZ)
where
F » mass flow rate of air supplied by the fan system to the room
Cp “ specific heat of air
Ts supply air temperature to the room (the hot deck temperature)
TZ =* air temperature in the room
For this system, F and Cp are constant and Ts is constant for range A-B.
Therefore, as room temperature (TZ) falls, the heat added to the room
increases linearly as shown in figure E.3.
4. Ranges D-E and E-F are analogous to ranges B-C and A-B, but are for
cooling.
CPFL ( NAME-C2 / TMP=60 ,67 ,69 ,77 ,79 ,85 / CAP“1 .73
, ,0 ,-0 ,-.8 ,-l )
When defining a control profile like the one shown in figure E.3, users
should select arbitrary low and high room temperatures to establish the full
heating and cooling (+1 and -1) points on the profile. In the example
above, 60 and 85 F (15.6 and 29.4 C) were selected. When specifying
corresponding heating and cooling capacities for each zone, the "design"
capacity generally should not be used. Capacity should correspond with the
selected low and high temperatures.
The same reasoning used for the three-deck multizone system can be used
to construct control profiles for spaces served by conventional multi-zone,
dual duct, or reheat systems (see figure E.4). Ranges A-B and C-D
correspond to ranges A-B and E-F of figure E.3, respectively. Range B-C
covers the temperature range where the mixing box, zone dampers, or reheat
coils are modulated by the room thermostat. The TARP description of such a
profile is:
129
CPFL ( NAME=C3 / TMP = 60,73,77,85 / CAP = 1 , .58 ,-.74 ,-l )
Figure E.5 shows a typical control profile for a VAV (variable air
volume) system with reheat. It was constructed based on the assumption that
3850 cfm (1.82 m-Vs) of 55 F (12.8 C) delivery air is to be supplied to the
zone with the VAV dampers fully open at 79 F (26.1 C). The dampers close to
20 percent of full flow at 77 F (25 C). Reheat is to operate between 69 and
67 F (20.6 to 19.4 C), and the reheat coil will deliver air at 140 F (60 C)
when fully energized. The following temperatures and capacities define the
control profile:
CPFL ( NAME=C4 / TMP=60,67 ,69 ,77 ,79 ,85 / CAP=1 , .91 .-.09 ,-.15 ,-.8 ,-l )
REFERENCES
[1] Hittle, D. C., BLAST, Version 2.0: Users Manual, Vol 1, U. S. Army
Construction Engineering Research Laboratory, Technical Report E-153
(Champaign IL, 1979), pp 41-46.
130
CAPACITY
SYSTEM
131
+ 1.0
CAPACITY
SYSTEM
NORMALIZED
132
CAPACITY
SYSTEM
NORMALIZED
133
Figure IV. E. 4 Control profile for multizone system
134
Figure IV. E. 5 Control profile for VAV system
135
F
To accurately account for shading, solar gains, and the effects of wind
on exterior walls, windows, and roofs, a fairly detailed description of the
building's geometry is required. Every building surface is represented by
a plane polygon of three or four sides. Within the program, each polygon is
represented by the locations if its vertices in a simple Cartesian coordinate
system. In that system, the X axis points due east, the Y-axis points north,
and the Z-axis points upward. The orientation of each surface is maintained
by two conventions. First, the vertices are stored in counter-clockwise
sequence as the surface is viewed from outside the zone it encloses. (Looking
at the surface from the opposite side would give a clockwise sequence.)
Second, the orientation is described in terms of the "azimuth" and "tilt"
angles.
The simplest level (0) of shadowing detail (3.1.2) requires that only
the azimuth and tilt angles be specified. Detailed air flow analysis (3.1.4)
and the detailed exterior surface coefficient algorithm (3.1.6) require
correct heights for every surface. The first and second levels of shadowing
detail require that all surfaces of the building envelope be correctly posi-
tioned in the building coordinate system. The third level of shadowing
detail requires that internal surfaces also be correctly positioned and that
the zone is completely enclosed.
136
Non-rectangular surfaces can also be described by giving the
coordinates of the vertices. For example, VRTS“10,0,5,30 describes the
triangle in figure F.2. Negative coordinates like VRTS=12 ,-3 ,7 ,7 ,-5 10 are
,
The location and size of these features are described in the coordinate
system of the base surface (wall, roof, etc.) on which they lie. The size
of each window or door is given by its width and height (for example, W=4 /
H=6 or equivalently SIZE“4,6 ). The subsurface is located on the base
surface by giving the coordinates of its lower left corner relative to the
origin of the base surface (for example X es 4 / Y“1.5 or 0RG“4, 1.5,0 ).
Coordinates are given relativve to the base surface coordinate system.
Figure F.4 shows the window thus described placed on the wall from figure
F.l. The subsurface must lie completely within the base surface. Non-
rectangular subsurfaces must be described as equivalent rectangles.
It is common for windows and doors to be set back slightly behind the
plane of the base surface. This causes shadowing of the subsurface,
especially as the sun's rays become nearly parallel to the surface. The
setback is described by the Z coordinate ( Z=-.5 ) of the location of the
subsurface. (A "set forward" is not permitted.)
Figures F.5 through F.8 show how the size and position of these
shadowing features are described in terras of W, H (or SIZE) and X, Y, Z (or
ORG).
A shadowing subsurface will shade the base surface and all its window
or door subsurfaces. It will not shade any other base surface.
137
4. )
The surfaces of a zone are located relative to each other in the zone
coordinate system. It is a Cartesian coordinate system with a vertical
Z-axis and the X and Y axes oriented parallel to major zone features to
simplify placement of surfaces. A surface is located in the zone coordinate
system by giving the coordinates of the surface origin relative to the zone
origin ( X=r / Y-r / Z-r or ORG»r,r,r [4.13.6]) and the azimuth ( AZM - r
[4.13.3]) and tilt ( TILT-r [4.13.4] angles of the surface.
The tilt angle is the angle from the Z-axis to the normal to the
surface. A vertical surface has a tilt of 90 degrees. A horizontal roof
has a tilt of 0. A floor has a tilt of 180 (the outward normal points
down)
Azimuth is the clockwise angle between the Y-azis and the horizontal
projection of the outward normal of the surface. For vertical surfaces this
is simply the direction the wall faces. Figure F.9 shows the locations,
azimuths and tilts of three vertical and three tilted surfaces. For
horizontal surfaces the true azimuth angle is undefined, but one must be
given to properly locate the surface vertices in the zone coordinate system.
This is most easily done by imagining the horizontal surface rotated about
an edge until it is vertical to allow identification of origin, width, and
height as in figure F.10 for a flat roof and figure F.ll for a floor.
Because of the general orientation of the building, the Y-axis may not
line up with true north. The angle from true north to the Y-axis of the
building coordinate system ( N«r [3.1.1]) must be given to allow TARP to
compute the true coordinates of the surface vertices.
138
Zone 4 in figure F.13 is very similar to zone 1, but its geometric
description cannot be obtained by any combination of rotations and transla-
tions. It is a mirror image of zone 1. Its description must use one of
the mirror image commands ( MIRX [4.1.5] or MIRY [4.1.6]). Figure F.15
shows how zone 1 is first mirrored in the Y direction (Y coordinates in the
zone system are replaced by their negatives) and then translated to a new
origin to produce the description of zone 4. Figure F.16 shows that the
same results could have been obtained by mirroring in the X direction,
rotating 180 degrees, and translating the origin. Mirror imaging works for
all surfaces and subsurfaces except those cases where a subsurface is on a
non-rectangular base surface which was described with the first vertex not
lying on the X-axis (such as the base surface in figure F.3).
139
Y
Y
140
Figure IV. F. 5 Overhang Figure IV. F. 6 Left fin
• Subsurface origin
• Base surface origin
141
1: ORG = 0, 18.2, 0 / SIZE = 10.7, 8 / AZM = 270 / TILT = 90
2: ORG = 0, 7.5, 0 / SIZE = 10.7, 8 / AZM = 225 / TILT = 90
3: ORG = 7.5, 0, 0 / SIZE = 10.7, 8 / AZM = 180 / TILT = 90
4: ORG = 0, 18.2, 8 / SIZE = 10.7,8 / AZM = 270/ TILT = 20
142
ORG = 0, 0, 8 / SIZE = 12, 12 / AZM = 180 / TILT = 0
143
ORG = 0, 12, 0 / SIZE = 12, 12/ AZM = 180/TILT = 180
144
3: ORG = 0,0, 10
4: ORG = 15, 0, 10
145
Figure IV.F.13 Portion of building with similar zones
146
Figure IV.F.14 Zone 3 described by rotation and translation of zone 2
147
GEOM (MIRY/ORG = 0, 55, 0/
148
Figure IV.F.16 Zone 4 described by mirror imaging, rotation and
translation of zone 1
149
IV. G TABLES
Shadow patterns are normally computed for the 14 time periods listed
below. These periods were chosen so that the error in the solar position
would be less than 2 degrees on any day in the period. This causes the
periods to be shorter around March 21 and Sept 21 when the solar declination
angle changes rapidly. Shading features which do not use these dates will
increase the number of time periods and therefore the amount of time
necessary to calculate shadowing. Most shading features will occur over
several time periods, e.g. FROM - 16MAY / THRU = 30SEP.
150
V OUTPUT
The input processor produces an "echo" of the user input. This echo is
organized by numbered lines on numbered pages as shown in the report for the
sample problem from II. C shown in figure 1 on the following pages. Each page
heading includes a page number, title, date, time of day, and central proces-
sor execution time. This last number is based on a non-ANSI function and may
not be available on all computers. The input lines are numbered primarily to
aid in editing the user input file. If user lines are not over about 72 char-
acters long, the output will appear on most terminals without "wrap-around".
The last page of the input processor report tells the unused library
elements and the required array dimensions (parameters) for the load calcula-
tor. The report of unused library names can help to identify typing errors
which cause library elements to not be found when they are referenced later
in the input. The report also identifies elements which can be removed from
the input deck (possibly by commenting - $) to use less array space. The
array dimensions report tells the minimum parameter sizes (see VI. C for a
complete list of parameters) necessary to run the simulation.
151
1 ) ) )) )) )) ) ) )
1 PROJECT t
2 RC<DEM/UIN=ENGLISH/UOUT=ENGLISH/
3 DESC* THERMAL MASS TEST - SIMULTANEOUS ZONES DIRECT SUN')
' »
8 CLR= 1 ./M0N/WS=1320/DIR=315)
9 3
10 LIBRARY C
1 DS(NAME=NONE/DESC='ZERO ALL HOURS /FFC=24*0 '
52 MATL(NAME=DIRT/K=.1/D=65./A=,70)
152
PAGE 2 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/17. 06.07.56. 1.048
79 CDS NAME = HS/DESC= 'SOUTH HE AT I NG /CPFL = 6*P 1 2*P2 8*P 1 5*P2 3*P 1
( ' , , , ,
83 3
84 BUILDING C
85 SIM(SLDS=3/HTB=0/FLOW/CNVG=,1» .01, ,34.12)
86 DSS< AZM=180./TILT=153.4/0RG=0, 12, 1 6/SI ZE=44 , 2
87 DSS(AZM=180,/TILT=206,6/0RG=0»0»8/SIZE=44»2)
88 DSS( AZM = 0./TILT = U6.6/ORG = 0,36»8/SIZE = 44,2)
8? 3
90 ZONE C
91 GEOM( NAME-NORTH)
92 CTRL(CWS=HEATN/HEAT=12)
93 PE0(WS=0CCN/CAP=3)
94 LIT <WS = LITN/CAP=1 ,4)
95 EQP (WS=EQPN/CAP=1 ,0)
96 MASS(C0NS=PTN-A/AREA=480)
97 SRF(EX/0P/BS/C0NS=WL-I/AZM=0/0R6=44»36»0/SIZE=41»8)
98 SRF(EX/TR/SS/C0NS=WDW2/0RG=12,5,-.25/SIZE=20,2)
99 AFS(AREA=0.2/EXP=,65)
100 SRF EX/0P/BS/C0NS=WL-S/AZM=0/0RG=3,36,0/SIZE=3,8)
(
101 SRF(EX/0P/BS/C0NS=RF-T/AZM=0/0RG=44,36»8/SIZE=41,8.95/TILT=26.6)
102 AFS AREA-0 1/EXP= .65)
( .
153
3 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/17. 06.07.58. 1.491
108 SRF(EX/0P/BS/C0NS=SLAB-C/AZM=180/0RG=0.36»0/SIZE=44»12/TILT*180)
109 OSC < TGP= 1
no 3
1 1 ZONE C
112 GEOM<NAME=CENTER)
113 CTRL(CWS=HEATC/HEAT=12)
114 PE0(WS=0CCC/CAP=3)
115 LIT (WS=LITC/CAP=1 .4)
116 EGP(US=EQPC/CAP=1 .0)
117 MASS(C0NS=PTN-A/AREA=480)
118 SRF (EX/O P/BS/C0NS=FRAME/AZM*180/0RG=0» 12. 12/SIZE=44»3)
119 SRF(EX/TR/SS/C0NS=UnU2/0RG=2.0.-.25/SIZE=40.3)
120 AFS( ARE A=0. 25/EXP*. 65)
121 SRF <EX/0P/BS/C0NS = RF-T/AZM = 0/0RG = 44»24» 12/SIZE = 41 »8.95/TILT = 26.6)
122 AFS(AREA=0. 1/EXP*. 65)
123 SRF(EX/0P/BS/C0NS=RF-I/AZM=0/0RG=3.24.12/SIZE*3.8,95/TILT*26.6)
124 SRF(IN/0P/BS/C0NS=PTN-A/AZM=180/0RG=3»12»0/SIZE=41.12/Z0NE=S0UTH)
125 AFS(AREA=0. 5/EXP*. 50)
126 SRF IN/0P/BS/C0NS=PTN-S/AZM=180/0RG=0. 12»0/SIZE=3» 12/Z0NE=S0UTH)
(
131 SRF EX/OP/BS/CONS=SLAB-C/ AZM* 1 80/0RG=0 24. 0/S I ZE=44» 12/TILT* 180)
( .
132 0SC(TGP=1)
133 1
134 ZONE C
135 GEOM(NAME=SOUTH)
136 CTRL(CWS=HEATS/HEAT=12)
137 PE0(WS=0CCS/CAP=3)
138 LIT <WS=LITS/CAP=1 *4)
139 EQP(US=EQPS/CAP=9. 4/LOST*. 4)
140 MASS(C0NS=PTN-A/AREA=480)
141 SRF(EX/0P/BS/C0NS=WL-I/AZM=180/0RG=0.0»0/SIZE=41»8)
142 SRF(EX/TR/SS/C0NS=WDU2/0RG=12.3.-.25/SIZE=20»4)
143 AFS< ARE A=0. 2/EXP*. 65)
144 SRF(EX/DP/BS/C0NS=WL-S/AZM=180/0RG=41»0»0/SIZE=3»8)
145 SRF(EX/0P/BS/C0NS = RF-T/AZM=180/0RG = 0.0.8/SIZE = 41 8 95/T I L.T = 26 6
» . .
154
PAGE 4 MBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/17 06.07.59 1.884
zon: 3
srf: 32
hts : 32
afs : 10
cos: 12
gtp: 1
tfh : 283
ctf: 140
gsc: 3
dsh: 11
wsh: 10
con: 11
pfl: 2
cds: 3
cws: 3
Figure V. 1 (concluded)
155
V.A.2 RUN CONTROL MODIFICATIONS
The 2SPC command (1.1.4) causes a blank line to be printed before every
input line echo as shown in figure 3. Note that double spacing does not
begin until the command has been recognized, so there is no skip between
lines 1 and 2.
1 PROJECT C
2 RC(DEM/2SPC/UIN=ENGLISH/U0UT=ENGLISH/
4 RPTCRNOS = 9)
5 LOC(DESC='COLUMBIA» MO .
' /LATD=38 9/L0NG=92 3/TZ=6
. ,
6 GRND<GRT=40»42»48>55»62»67»70»69»63»55»48»42)
7 DDY DESC= 'MILD WINTER /HIGH=50/L0W=30/WB=40/DATE=21FEB/
( '
8 CLR=1 ./M0N/WS=1320/DIR=315)
9 I
10 LIBRARY C
The XREF command (1.1.5) allows the user to match the data in the
reports and, especially, the dumps to his input. In general, the input is
read and organized into various arrays. The XREF reports the array
locations. In figure 4, line 87 describes a detached shadowing surface
(3.2). Information about this surface will be stored as element 2 in arrays
having a MAXSRF parameter dimension. The next DSS (line 88) is stored as
element 3 in those arrays. Zone controls and scheduled loads (lines 92-95)
are stored in arrays with MAXCGS dimensions. The MASS is the first heat
transfer surface, so its data goes into the MAXHTS dimensioned arrays at
location 1. The surface description (line 97) has two cross references
because both geometric data (in MAXSRF dimensioned arrays) and heat transfer
data (in MAXHTS arrays) must be stored. The zone number appears after the
zone description is completed (line 110).
157
87 DSS<AZM=180./TILT X 206.6/0RG=0»0»8/SIZE*44»2)
SRF = 2
88 HSS<AZM = 0./TILT = U4.6/0RG = 0>36»8/SIZE = 44*2)
SRF = 3
89 3
90 ZONE C
91 6E0M<NAME=N0RTH)
92 CTRL<CWS=HEATN/HEAT=12)
CGS = 1
93 PE0<US=0CCN/CAP=3>
CGS = 2
94 LIT<US=LITN/CAP=1 .4)
CGS = 3
95 EGP<WS=EGPN/CAP=1 *0)
CGS = 4
96 MASS<C0NS=PTN-A/AREA=480)
HTS = 1
97 SRF EX/0F‘/BS/C0NS = UL-I/AZM = 0/0RG = 44>36>0/SIZE = 41 » 8)
<
SRF = 4
HjS — 2
98 SRF<EX/TR/SS/C0NS=UDW2/0RG=12>5>-.25/SIZE=20»2)
SRF = 5
HTS = 3
99 AFS<AREA=0,2/EXP=,65>
AFS = 1
SRF = 12
HTS = 10
109 03C TGP= 1 <
SSC = 1
110 ]
ZON = 1
111 ZONE C
112 GEOM(NAME=CENTER)
113 CTRL(CWS=HEATC/HEAT=12)
CGS = 5
114 F'E0<WS = 0CCC/CAP = 3)
CGS = 6
158
BUILDING LOADS FOR
COLUMBIA* MO,
MILD WINTER
TOTAL TOTAL PEAK PEAK MAX MIN
ZONE HEATING COOLING HEATING COOLING TEMP TEMP HX HT V CL CX
(KWH) (KWH) (KU) < KW <C> <C>
DESCRIPTION OF CONSTRUCTION
159
V.A.3 NUMBERED REPORTS
This provides summaries at the end of the PROJECT and LIBRARY sections
The DUMP OF PROJECT CONTROL DATA shown below lists the values associated
with different variables which can be set in the PROJECT section. Refer to
VI. C for a description of the individual variables.
COMMON /TITLES/?
F'TITLE = THERMAL MASS TEST - SIMULTANEOUS ZONES * DIRECT SUN
LTITLE = COLUMBIA* MO.
ETITLE = MILD WINTER
COMMON /RUNCTL/J
EMU I = 2
SYNTAX = F
CALCS = T
DEM = T
OLDW = F
COMMON /EFLHDR/:
l.ATD 38.900
LONG 92.300
it NR 6.000
ALTD 0.000
TRRAIN = 3
SINL .628
CQSL .778
MERID = 2.300
TGRND = 4.4 5.6 8.9 12.8 1 A 7 19.4 21.1 20.6 17.2 12.8 8.9 5.6
.
TMKUP = 12.8 12.8 12.8 12.8 12.8 12.8 12.8 12.8 12.8 12.8 12.8 12.8
NRENVS = 1
NRDAYS 1
WSTATN = 0
COMMON /REPORT/*.
EMUO
T T T T F F F
RFLAGS = T T T
F F F F F F F F F
F
F F F F F F F F F F
160
The DUMP OF LIBRARY CONTROL DATA shown below gives the number and
names of all library entries. First, summary data is given for each of
the seven libraries. This data is taken from common block LIBCTL. It
begins with the day-schedule (DS) by giving the 0 location (DSO) the ,
LIBRARY NAME =
DS — 0 11 16
WS = 16 10 16
MATL = 32 13 16
CONS 48 11 16
CPFL = 64 2 16
CDS 80 3 16
CWS = 96 3 16
ELEMENT names:
NONE FULL ON OC OS LN LC LS
EN EC ES
OCCN OCCC OCCS LITN LITC LITS EQPN EGF'C
EQF'S CONST
GYP WD-S SHTH C-HW B R5 R 1 R30
SF'DS ASH DIRT AL CPT
WL-I WL-S PTN-A PTN-S WDW2 FRAME RF-T RF-I
F'TN-M SLAB-C SLAB-T
PI P2
HN HC HS
.3.2 RNOS = 2 ;
dump of DDY and WTAP data
This is a dump of all weather data for every day described in the DDY
(1.6) and WTAP (1.5) keys. It is a dump of data in the EVDAT common block.
The data names are summarized as follows:
161
.. 1
8 CLR= 1 ,/M0N/WS=1320/DIR=315)
TD TU BP H US UD TS QB QD QR RN SN
.33 -.87 100000. .0031 6.7 315. -16.16 0.0 0.0 0.0 0 0
- .22 -1.19 100000. .0031 6.7 315. -16.69 0.0 0.0 0.0 0 0
- ,67 -1.45 100000. .0031 6.7 315. -17.10 0.0 0.0 0.0 0 0
-1.00 -1.64 100000. .0031 6.7 315. -17.42 0.0 0.0 0.0 0 0
-1.11 -1.71 100000. .0031 6.7 315. -17.52 0.0 0.0 0.0 0 0
-.89 -1 .58 100000. .0031 6.7 315. -17.31 0.0 0.0 0.0 0 0
-.33 -1.25 100000. .0031 6.7 315. -16.79 0,0 0.0 0.0 0 0
.67 -.67 100000. .0031 6.7 315. -15.85 307.4 18.6 10.2 0 0
2.11 . 09 100000. .0031 6.7 315. -14,49 734.5 44.5 51 . 0 0
3.78 1.05 100000. .0031 6.7 315. -12.93 875.8 53.1 88.0 0 0
5.67 2.12 100000, .0031 6.7 315. -11.15 937.7 56.8 116.1 0 0
7 . 44 3.09 100000. .0031 6.7 315. -9.48 965.4 58.5 133.2 0 0
8,78 3.80 100000 .0031 6.7 315. -8.23 972.2 58.9 138.0 0 0
9 . 67 4.27 100000. .0031 6.7 315. -7.39 961.0 58.2 130.2 0 0
10.00 4.44 100000. .0031 6.7 315. -7.08 927.0 56.2 110.4 0 0
9.67 4.27 100000. .0031 6.7 315. -7.39 852.7 51 .7 80.0 0 0
8.89 3.86 100000. .0031 6.7 315. -8.12 675.5 40.9 41.5 0 0
7.67 3.21 100000. ,0031 6.7 315. -9.27 107.0 6.5 2.6 0 0
6.22 2.43 100000. .0031 6.7 315. -10.63 0.0 0.0 0.0 0 0
4.78 1 .62 100000. .0031 6.7 315. -11.99 0,0 0.0 0.0 0 0
3.56 .93 100000. .0031 6.7 315. -13.14 0.0 0.0 0.0 0 0
2.44 .28 100000 .0031 6.7 315. -14.18 0.0* 0.0 0.0 0 0
1 .56 -.17 100000 .0031 6.7 315. -15.02 0.0 0.0 0.0 0 0
.89 - . 55 100000. .0031 6.7 315. -15.64 0.0 0.0 0.0 0 0
162
V.A.3.3 RNOS “ 3; dump of conduction transfer function calculations
N ROOT RATIO
1 .013546079 .516548286
2 .027011651 ,072319215
3 ,040288196 .002899261
4 ,053198718 .000037605
5 065446912 ,000000201
6 .076662966 .000000001
7 .036835154 .000000000
B .096795305 .000000000
9 ,107590765 ,000000000
N 0 0 0 1 0 2 0 3 4 0 5
163
V.A.3.4 RNOS = 4; dump of glass calculations
The polynomial coefficients are the least squares fit values for
(D.1.7b) and (D.1.7a). T refers to transmittance, A to absorptance, 0 to
outside, I to inside, and S to shade.
PANE= 1 T 1 = .858 Rll* .079 R 12= .079 T2= .861 R2 1 = ,075 R22* .075
PANE* 2 T 1 = .858 R 1 1 = .079 R 1 2= .079 T2= .861 R2 1 = ,075 R22* ,075
•
PANE* 1 T 1 = .828 R 1 1 = .106 R12= .106 T2= .882 R21* .051 R22= .051
PANE* 2 T 1 = .828 Rll* ,106 R 1 2= .106 T2= .882 R21* .051 R22* .051
OPTICAL VALUES AT COS THET A = .938
< )
F'ANE= 1 T 1 = .845 Rll* .090 R 1 2= .090 T2= .871 R21= .064 R22= .064
PANE* 2 T 1 = .845 Rll* .090 R12= .090 T2= .871 R21* .064 R22* .064
OPTICAL VALUES AT COS THET A = 1.000
< )
FANE* 1 T 1 = .860 Rll* ,077 R 1 2= .077 T2= .860 R21* .077 R22* .077
PANE* 2 Tl= ,860 Rll* .077 R 1 2= .077 T2= .860 R21= .077 R22= .077
164
POLYNOMIAL COEFFICIENTS
I T2 A02 AI2 AS2
1 3.8544537E-01 427646 1 E + 00
2 » 5 , 75769 1 2E-0 -I
2 1 .0308017E+01 -2 2363210E+01
, -3.6428286E+00 -I
3 -3.3949733E+01 1 .0021852E+02 1 .4986929E+01 -I
4 5.0154231E+01 -2.3554957E+02 -3.4400820E+01 -I
5 --3.9590559E+01 2 7725274E+02 3 » 984 1 064E + 0 -I
6 1 .4092558E+01 -1 2867078E+02
. -1 .8204343E+01 -I
I T A0 All AS 1
1 -3.8380648E-01 1 .2150528E-01 .2238673E-02
1 -I
2 4.402 1 943E+00 -1 .3248196E-01 2 4052260E-0
, -I
3 -6.7817725E+00 1 .6526246E-01 -4 1181869E-01
. -I
4 4 835 1 028E + 00
» -1 .2219971E-01 2 9334734E-0
. -I
5 -1.3277812E+00 3 52326 1 9E-02
. -7.9688593E-02 -I
This dump begins with the surface number (same as SRF number in
V.A.2.3). Then the coordinates of the vertices of the surface are presented
as X, Y, and Z values (in meters) arranged in the same order as the surface
is viewed. For example, (-.000 3.113 5.150) are the X, Y, and Z coordinates
respectively of the upper left vertex. (13.411 3.658 4.877) refers to
the lower corner. Several individual values are then presented:
CA = cosine of the axirauth (global coordinate system) angle
SA = sine of the aximuth angle
CT = cosine of the tilt angle
ST = sine of the tilt angle
TRNS = surface transmittance
ON = start date in cyclic date (B.1.4) form
OFF = end date
SI ,S2 = variables SHDSRF(I,1) and SHDSRF(I,2). See VI. C.
85 BUILDING C
86 SIM(SLDS=3/HTB=0/FL0U/CNVG=.1» .01 .34.12) »
surface: 1
VERTICES :
- ,000 3.113 5.150 13.411 3.113 5.150
0,000 ro >0 in CO 4.877 13.411 3.658 4.877
CA SA CT ST TRNS ON OFF SI S2
-1 ,000 -.000 -.894 448 0.000 33 415 0 0
SURFACE 2
VERTICES
- . 000 -.545 2.165 13.411 -.545 2.165
0.000 0.000 2.438 13.411 -.000 2.438
CA SA CT ST TRNS ON OFF SI S2
-1 .000 -.000 -.894 -. 448 0.000 33 415 0 0
165
V.A.3.6 RNOS “ 6; dump of ZONE data
The zone data dump produces output as shown in figure 12. There are
different data formats corresponding to the many different elements of zone
data. Controls and scheduled loads produce similar output as shown for
lines 92 through 95. First the schedule number (corresponding to CGS in
V.A.2.3) is given. Then the contents of CGSZON (zone number), type of
load), GGSON (on date), CGSOFF (stop date), CGSDX (week-schedule pointer),
and elements 1 to 5 of RCGS (see VI. c) are given order.
Lines 96 and 101 produce dumps of heat transfer surface data. They
begin with the heat transfer surface number (corresponding to HTS in
V.A.2.3). On the next two lines the following values are given:
The surface described in line 101 is part of the envelope of the zone.
Therefore, it also has shadowing surface data identical to that described in
V.A.3.5.
The air flow opening at line 102 produces another dump beginning with
the XREF number (corresponding to AFS). Then following variables are dumped
in order:
Other side coefficients (line 109) also produce a data dump beginning with
the XREFF number (corresponding to OSC). Then the following variables are
dumped in order:
166
90 ZONE C
91 GEOM<NAME=NORTH)
92 CTRL<CWS=HEATN/HEAT=12)
schedule: i
93 PE0(WS=0CCN/CAP=3)
schedule: 2
ZON TYP ON OFF DX R1 R2 R3 R4 R5
1 1 33 415 1 3.000 131.800 .500 0.000 0.000
schedule: 3
ZON TYP ON OFF DX R1 R2 R3 R4 R5
1 2 33 415 4 410.025 0.000 .200 .200 0.000
95 EGP<WS=EGPN/CAP=1 .0)
schedule: 4
ZON TYP ON OFF DX R1 R2 R3 R4 R5
1 3 33 415 7 292.875 0.000 .300 0.000 .700
96 MASS(C0NS=PTN-A/AREA=480)
HEAT TRANSFER SURFACE: 1
AREA INV ONV FSS FS6 K
44.59 1 -1 0.000 0.000 3
AZM CTILT HOC THE THI QHE GHI
0. 0.000 0.000 0 0 1 1
101 SRF<EX/0F/BS/C0NS=RF-T/AZM=0/0RG=44»36»8/SIZE=41>8.95/TILT=26.6)
HEAT TRANSFER SURFACE: 5
AREA INV ONV FSS FSG K
34.09 1 -3 .947 .053 7
AZM CTILT HOC THE THI GHE GHI
0. .894 3.339 20 26 32 37
surface: 7
VERTICES .
13 .411 8.534 3.660 .914 8.534 3.660
13 .411 10.973 2,438 .914 10.973 2.438
CA SA CT ST TRNS DN OFF SI S2
1 ,000 0 .000 894 .448 0.000 33 415 5 5
109 0SC(TGF'=1)
110
Figure V.12 Dump of zone data
167
.
The air flow opening at line 105 also created an entry in the AFO
table. The count and zone names are similar to those used for the IZP. The
number are the opening area, flow exponent, and opening height, respectively.
The input processor continues to not find matches (lines 124, 125, 126,
148, and 149) or to match IZP and AFO (lines 127, 128, and 150). When a
match is found, the reduced number of items IZP or AFO tables is reported as
well as the HTS or AFS of the match. In this particular example, there is
a typing error in the name of the adjacent zone at line 148. This error
prevents a match of both IZP and AFO. At the conclusion of input processing
168
this failure to match is reported as SEVERE errors. The data accompanying the
error messages (produced by default) should be sufficient to determine the
incorrect or missing surfaces and openings.
124 SRF(IN/0P/BS/C0NS=F'TN-A/AZM=180/0RG=3»12>0/SIZE=41*12/Z0NE=S0UTH)
IZP - NO MATCH*. 2 SOUTH CENTER 45.708 4.8998 1 1
125 AFS(AREA=0.5/EXP=.50)
IZF - NO MATCH: 2 SOUTH CENTER .04645 .500 1.829
126 SRF< IN/0P/BS/C0NS=PTN-S/AZM=180/0RG=0> 1 2 0/S I ZE = 3 , 1 2/Z0NE = S0UTH
>
169
V.A.4 LIBRARY REPORTS
ON l NORTH OCCUPANCY
3 5
1 2 3 4 5 6 7 8 9 1 0 1 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 19 20 21 22 23 24
1.0 +XXXXXXXXXXXXXXXXXX XXX+
+XXXXXXXXXXXXXXXXXX XXX+
.8 +XXXXXXXXXXXXXXXXXX XXX+
+XXXXXXXXXXXXXXXXXX XXXXXX+
6 +XXXXXXXXXXXXXXXXXX XXXXXX+
+XXXXXXXXXXXXXXXXXXXXX XXXXXX+
.4 +XXXXXXXXXXXXXXXXXXXXX XXXXXX+
+XXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXX+
+XXXXXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXXXXXXXXXX+
+XXXXXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXXXXXXXXXX+
0 +XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX+
100100100100100100 50 16 0 0 0 0 0 0 0 0 16 16 16 16 33 33 67100
0C : CENTER OCCUPANCY
4 »
1 23456789 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
1.0 + +
+ +
.8 + +
+ XXX XXX +
.6 + XXX XXX +
+ xxxxxxxxx +
.4 + xxxxxxxxx +
+ xxxxxxxxx +
.2 + xxxxxxxxxxxxxxxxxxxxx +
+ xxxxxxxxxxxxxxxxxxxxx +
0 +XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX+
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 17 17 17 67 50 67 17 0
The first two lines of the week-schedule report give the user defined
name, description, and the schedule number. These are followed by a list of
12 days: SUN through FRI on one line and SAT through SP4 lines below. Below
each day name is the name of the day-schedule in effect on that day. Below
that is the day-schedule number. Figure 15 shows three week-schedules.
170
0CCN *, NORTH OCCUPANCY
i
SUN MON TUE WED THU FRI
ON ON ON ON ON ON
3 3 3 3 3 3
The first two lines of the materials report give the user defined name,
description, and the schedule number. The thermal properties of the material
are given below the names (using the input subkey names) of those properties.
Optical properties are also reported when they apply.
MATERIALS S
K CP D R E RGH
.420 .300 45.000 0.000 .7! .90
WD-S : SOFTWOOD
OTHER
CP D R A E RGH
K
.330 32.000 0,000 .75 .90
,068
171
V.A.4.4 Constructs (CONS)
The first two lines of the construct report give the user defined name,
description, and the schedule number. The name and thermal properties of
each layer are given. The number and values of the conduction transfer
functions are reported. "Internal" refers to the Z functions, "cross" to
the Y functions, "external" to the X functions, and "flux" to the J func-
tions. The thermal conductance and surface properties are reported as shown
below in figure 17.
DESCRIPTION OF CONSTRUCTION
172
WDU2 : 2 PANE WINDOU
5 *
DESCRIPTION OF CONSTRUCTION
OPTICAL VALUES
COS< THETA) .0625 .1250 .1875 2500 .3125 .3750 .4375 .5000
COS ( THET A ) .5625 .6250 .6875 7500 .8125 .8750 .9375 H oo DIFFUSE
transmittance: .6743 ,6962 .7117 7228 .7308 .7366 .7409 .7440 .6554
reflectance: .1832 .1639 .1511 1429 .1380 .1352 .1341 .1341 .2089
ABSORPTANCE: ,0811 .0788 .0767 0747 .0728 .0709 .0691 .0673 .0782
ABSORPTANCE: ,0615 .0611 .0605 0595 .0584 .0572 .0559 .0546 .0575
173
V.A.4.5 Control profiles (CPFL)
The first two lines of the control profile report give the user defined
name, description, and the schedule number. The next line is a list of the
user defined control temperatures followed by a line of matching fractions
of full capacity. The number of control points is determined by the user
subject to array dimension constraints. Figure 19 shows this report for two
control profiles.
CONTROL profiles:
PI : HEATING SETBACK
: STANDARD HEATING
After the two line of standard heading information the CDS report gives a
list of 24 hours (on four lines) with the profile name and index reported below
each hour. Figure 20 shows a CDS report for two control day-schedules.
The first lines of the control week-schedule report give the user defined
name, description, and the schedule number. These are followed by a list of
12 days: SUN through FRI on one line and SAT through SP4 four lines below.
Below each day name is the name of the control day-schedule in effect on that
day. Below that is the control day schedule number. Figure 21 shows two
control week-schedules.
174
CONTROL day-schedules:
HN : NORTH HEATING
i
i: 2 3 4 5 6
PI PI Pi Pi Pi PI
1 1 1 l l 1
7 8 9 10 li 12
P2 P2 PI PI Pi PI
2 2 1 1 1 1
13 14 15 16 17 18
PI PI PI PI P2 P2
1 1 1 1 2 2
19 20 21 22 23 24
P2 P2 P2 P2 PI PI
2 2 2 2 1 1
HC J CENTER HEATING
2
1 2 3 4 5 6
PI PI PI PI PI PI
1 1 1 1 1 1
7 8 9 10 11 12
PI PI PI PI PI PI
1 1 1 1 1 1
13 14 15 16 17 18
PI PI PI PI P2 P2
1 1 1 1 2 2
19 20 21 22 23 24
P2 P2 P2 P2 PI PI
2 2 2 2 1 1
CONTROL UEEK-SCHEDULES:
175
V.B LOADS REPORTS
The next default report is shown as page 2 in figure 22. It describes the
calculational effort used to determine the building loads. First, the tempera-
ture and flux history initialization period is reported. The execution time
for initialization can be determined by comparing the completion time with the
start time (which was given with the page 2 header). Initialization proceeds
until load and temperature extreme values for two successive days are within
the convergence limits. In this example that required seven days. There were
635 surface and air heat balances, or about 3.78 iterations per hour. This
number can be reduced by widening the convergence limits. That will also
result in fewer days for initialization. Wider convergence limits tend to
change peak values more than they change the total loads. There were 2181
Newton's method iterations to compute air flows. That is about 3.43 iterations
per heat balance iteration. There were no heat balance iterations which failed
to converge. Failure to converge may be due to an inadequate iteration limit
(3.1.3) or to cycling between points on two different segments of a control
profile. The loads for that hour may or may not be reasonable. Check the
peak loads (below) for impossible values. Similar information on time, days,
and iterations is given for the simulation period.
176
!
loads, and the maximum and minimum of all zone temperatures. The hour and day
are reported for the occurrence of the peak loads and temperatures. In the
example, there Is no time reported for peak cooling since there was no cooling
load.
zqn: 4
SRF 36
HTS J 36
AFS 5 24
cgs: 16
GTP 5 4
TFH *
512
CTF! 192
ssc: 16
dsh: 16
wsh: 16
con: 16
pfl: 8
cds: 8
ows: 8
zon: 3
srf: 32
hts : 32
afs: 10
cgs: 12
gtp: 1
tfh: 283
ctf: 140
ssc: 3
dsh: 11
wsh: 10
con: 11
pfl : 2
cds: 3
cws: 3
REQUIRED CMB = 82
REQUIRED HCS = 10
177
) ) )
COLUMBIA* MO.
MILD WINTER
178
.
V.B.2. 1 RNOS =»
8; Dump of SINPFL as read by loads program
This is a dump of every value transferred from the input processor to the
loads calculator by the simulation input file (SINPFL). It includes all project
data except the weather and can, therefore, be rather long. Figure 23 for a
SINPFL dump of the sample problem data is four pages long.
The SINPFL dump begins with a dump of variables which can be described as
run control data. Refer to VI. C for a description of the individual variables.
Next comes a dump of library data beginning with the day-schedules. The data
for each day-schedule consists of the array index, the schedule name, and 24
hours of capacity fractions expressed as percents. The week-schedules are
dumped in a similar manner. Each line describes one schedule and consists of
an array index, the schedule name, and 12 pointers to day-schedule array
indices
There is no data for materials. All materials have been used in constructs.
The headings for the constructs data are:
The conduction transfer functions are listed directly from the CTF array
at ten numbers per line. The construct pointer values (PE, PC, and PI) point
to the location immediately before CTF location. For example, the X(l) value
of WL-S is -6.3074297 (PE = 7; therefore CTF index = 8). The X(0), Y(0), and
Z(0) values are stored separately. The window pointer, WTYPE points to the ,
array index of the window coefficients. In this example, there is only one
window type, but each type has four line of data: one each for transmittance
179
and outside, inside, and shade absorptance. Each line of data consists of
the window type index, data abbreviation, value for diffuse radiation, and
polynomial coefficients for angular properties.
Control profile data consists of the array index, profile name, number of
control points, and control temperatures and capacity fractions. In the
example, the capacities for PI are 1.0 at 1=15.56 and 0.0 at T-16.67. The
control day-schedule data consists of the array index, schedule name, and 24
pointers to control profile array indices. The control week-schedule data
consists of the array index, schedule name, and 12 pointers to control day-
schedule array indices.
Inputs (4.2) through (4.11) create scheduled loads and activities. Each
line of the scheduled loads dump consists of the array index, zone number
(CGSZON), type (CGSTYP), start and stop dates (CGSON and CGSOFF), schedule
pointer (CGSDF), all 5 values of RCGS, and zone names (CGSZON). See VI. C for
descriptions of these variables.
Next the surfaces involved in shadowing and solar gains (SRF) are described.
This consists first of the array index, number of vertices, and X, Y, and Z
coordinates of the vertices for each surface. Then comes the subsurface data
previously described in V.A.3.5.
The heading for the columns of data for the heat transfer surfaces (HTS)
are:
The special coefficient data consists of the array index and nine values
from ISSC or RSSC (see VI. C) depending on the type of data being described
(see also V.A.3.6). Special surface coefficients are stored at the bottom of
the array and movable insulation data are stored at the top.
180
30©0©-«00 © © o o
KNNNNNN f^NNfj
0'©©^aoo‘ 0 ‘ m rs rs O'
rs o o » ruj >o co © in ^
©CNOOCNO'COIf) in -r in
~
ooooooooooo
oo ro
U*) CN QO
©orsps©©©©©©© oo©o*oo©ooo©
u.if)ininmroinminininin
©< ** <r ro o ro ©rsrsrsrx^rsrs rs rs rs rs
oowisooooooo
o n n o
ro >o ro
©©«CNWrt©COO©
O h h M M N N C©
00'CNWK0©©©0 ©©o©©©© ©©©©
o©©©©©©© ©©©©
Ch H«n N o© in UO'OO'O'O'OO' O' O' O' O'
©©c©oo©o©
©
©o ro
3.085
©0©©©0©0©©0 ©©©©©©©©©©©
oinmifiininmin © © in in
© K ©rsrsfsrsrsrsrs rs rs rs rs
©©©©©©©COCO
e ro
co<r<r<rco^ij )0'^nro
©O©©©©©©©©©
,
<raococoinincDroro(Nro
©©©©©©©©©©©
© ro
ro <r in >o n co o© ro ro
c in ro o o© © ©
©©©©©CO©©©©
© ro
ro^ifi'Orscoo'© ^ cn
©
r.
CN
o ch
© in ©
cc
O'
rv c
(f)©
in lo
h«oo
in «* ©
©
^
62/05/17.
©OOCOO©©©©©
© ro
L*)'Orvooos©^CN
L*)
©©©©©©©©©©© M^-ll'OfVOOOsO^M
'OrO^’TCDHMM^CO
©OOCLCsiCDO-rsCCOO
©CDlf)'T«COo-'C CD CC ro ro
(N<raicoif)if)'Orororsrs
©©>c©©©©©©©©
©^ © lo rs rs
M^LI'ONOOO'OhCN ro
L.'COTNO'inTiniiT-'O
PROGRAM
©©©©©©©©©©©
© ©
in in iv rs
ro^jiTsorsCOOsO-iCNj Cl. ~«(\roro«r'Ofscc-<ro
©©©©©©©©©©©
©© ro
ro <r in>0 NCD 0'©-.(\ U. CN — CN © ^ ID
RESEARCH
©©©©©©©©oo©
©© ro
ro^if)sOrsCDOs©s-<CN
M^bTronoocc ro cn ro O'
u. csi ro ro «r n is co rj ©
©©©©©©©©©©©
©© ro
ro«rmorscoo'©»-<(N
(j (\l ^ CN O' 0s ^ (N fH O' ^ f 1
x
^(Mr'jro^WNNO'N
©o©©©©©©©©©
©© ro
ro o* in >o rs co oo ** cn
a.
ANALYSIS
ujosHioo>0"CO"ONin
©©©©©©©©©©©
oo ro
ro <r in o rs coo*© cn U_ CNCNfOrO^ <0 N CC -•
>0 © © © ©
©©©©©©©©©©©
©© ro
ro<r ut 'O rs co O' o cm rj«r»-'«TOCN'CcNrocDrs
THERMAL in © © O ©
UD ^ Os ^ © ©
© © © © © ©
• CN © <
NBS
ii ii h u h n ii n ii u ii n ii
h © CO UJ 2 I
3 2 WOO
cn ^ CN cn ro III
I 2
UJ J 20 CO O Z 2
UUL)l-KWiiIu.2
CO O COCO jo
22 3 ©
i
2© ©
i rj r r- h i x a
I
3 O CD CD
CD H- _J CD CD -J
Jjr*wwiiiU,wjj
I
©ou
I
oo
l |
181
1 o
•0 CD ^ h
inN M I
(\o ^
O'N O
O o CM
mcNi'OtNOrvCNo cm *i rs O' 00 GO
nNOo^in’OCD o^fS'Oin -4 CM^
ian oo« ro Nw © © * ro ro ^
-.1701522
m h f inin o on
rsjM'O^oaoin'O o ogo 10
hoomn *CO <
I
I
**
ooorso'O'O&oO'O o © m oo cm inn o
NrtNOIO'OODOD 000^00 O (N H
i
O' U*5 O^
(N K MH I I Ps I in CM GO oo
I ^ ««* I GO oo
I 0^ rv os
KN I")
o- o- v (h o n *o <
cd <n in ifi I CM
'ooo -o rg co © O' © © rs
oo in
«o < in O' in n w >o o do
rs 9-«a, PO fM CM Hi
'Oo ro in rs o- o © rg
o* o- O' *-« CM in ffl I
2.1913321
oo
oo Nin
m^ino'Oco o © © © © in
>o
oo'Oocn
t: ^ O'
in o-
o
©
00 O (N CD O H f O © *h © PO hino*
1 I<o I O* 1 0 in o- oo
1 CM CM inn w oc
1
CM I
Mh
I
PO PO rg rv rg ro ro o- O' oo 0 oo
o O' H O' o ~ O' ^H 00 oo
-6.3074297 oo CM O' © © © rv O' rg O'
cc o CM po m © © 00
PO o © ro © © © rs MOT
no h
in 1 rs rg 1 'C 1 oo
CM i oo CC 0C CC CC
i CM
OCJ
HMH
CD I
XX
p- >-
XXLj Id Id Id
»-->- p— »— »— p—
X XX X
I— >— »—
o- <r in O' ©00 OC PO CM Ctcrz Z Z Z X =5 X D
O' O'
O
CD
PO o- o- ro
CO 94 rs
00 oo in
C *0
ro ro •
*
oo
co
Z2oou
ZZ
uiiducoao
(J (J (J (J 05 05 05 05
^©
00 O' O' 94 O' rg oo
.0005898
oo m © o CM -4 ooooo ooooooo
ONO ooooooo
O'
9^
00
©
o
© o o^
rs rH
O © © Tl CO
rg I
oo
OO
ooo ©© © ©© rs PO
*o
rl
*0
PO
| in in
T'fl
nw
N I
oo
oo
ooo ooo ooo
1 1 m pm in
oo CM PM
rg O' ro rv O' •o ©*H ro <r in
oo ro rg ©
rg po rg
oo O'
CD
ro
00
©
^
'O o' in
©
©
fs
O'
N CM
oo ooo
oo ooooooo
ooo ooooooo
oo (moo o rg o o o cn t
.0023383
r* in 4T © © rg ©
©
PO in O © ro «->
©©
O'
©T © Ps O'
^
oo ooo ooo
co m* O'
~
ro i CO n pg *o i oo
^H i 1
*4 pv in O' oo
I I
— ©
1
wooN
CM P'5
oo OOO ooooooo
om ro © mooooooo
O' rv so ro rg rv
© in ro O O'© ^
© ^
fS *H I oo ooo
<, O' o- rg 9i O' oo cm ro © w
cm cm po
O' ro (M ro tr© in o' oo
-3.9713113 PO O' 00 o© © > U")
©
ro c © rs \n oo m *-•
©o 00 © 9H © ro © hn n i
in rg ro
94 1 1 in i PO CM O' oo
00 *h CM oo
<r
4+0
4+0 004+++ 0004+000
004 0004+000
o rg 94 Pv O' IM 00 "C
i
O oo oo cm rs © o rg n o o eg rg
o com oooomooo
i 1 i I
ino
©
CM *h rg ©
«r rH <r ro o cm t n o rg t n o n
©
co «r 00 co © ro
oo O' 00 CM hi o ^ 0 in ' 9-4
o
*T O' *T CD © o-
«r h t m i m o- cm in t rg in <r rs
ro© rs © 9H in © rg cm rg o PO P'5 PO CM
© rg o o © o ro
.0585821
O' OM T
©© o o© ©©^ t h rg
o © ** © ©©© i '-t rs cm rg in co ro po *o O'
o* i
i 'O i rg oo
ro i oo
1
i
oo H CM c w cn _j _J J « <t «
3 p- <r <t <t O 94 CM OZCJ05 Q Ld Id Id
TRANSFER
o
a
CC Cu ul
P—
iZXXX
—
ccxxx
P"
1 z z r-4 rg zh cm m z n cm ro i HCMMTingJNODCpOHrg
irgfOTin'ONoooHrgro hhh hh o o o
u u o
182
, 0 0
HP in m
IS rO is ro po CM CM CM N S3 -h GO 00 '0 M> CM M3 CM 'O © ro o «r o- M3 -0 © «o ©
** in cm pH o* in ro ro O' PO m >0 O' 00 00 >0 >0 PO >0 PO
3i
'O o^ © o- ^>0 <0 © * o
PO 1 pH
pH rH
ono pH
GO 00 rs rs 10,
7,
PO o- <r po ro
7. 7.
PO | CM I CM CM PO PO ro i
MD © rs rs cm rs rs rs CM CM rs rs M3 © CM CM M3 M3 CM *C CM CM © ^4 © ©© 'C <1 © M3 M3
M3 © O' O' in O' O' O' ro ro O' O' M3 © PO ro '0 M) ro ms ro ro © O' © ©© 'O '0 © M3 'O
PO 10, 10,
1 10, 10, 10, rs rs 10. 10, PO
0, 7.
rs ro ro
7, 3, 7.
rs i i i
3. 3. 0, 3. 3.
i
ph —4 rH rH in o rH © rH © rH rH rH pH © © rH rH ©© r4 rHO VO ,H ©
rH rH ©© rH rH
<r O' O' O' rs © O' © *T © <r «r o- 00 O' © O' O' ©o p- C in rs O' in o’ O' ©© O’ O*
CD © © © © © © © © © © ©© © © © © ©© © © © © © © © ©© © ©© © ©
Z © © © © o © © © © © ©© © © © © ©© © © © © © © © ©© © ©© © ©
© © © © o © © © © © ©© © © © © ©© © © © © © © © ©© © ©© © ©
© o* <t © co © o- o* © © © © M3 CD O M3 © © © © © © © 00 © <r t ©©© ©o cc
00 O* «r ©©© o- O' © © © © O © 'O >0 ©© ©©© © ©© © «r nr © ©© ©© r— © © © © © © © © © © ©© © © © © ©© © © o © © © © ©© © ©© © ©
3, 0. 0, 0, 0, 0, 0, 0, 0,
O* CM CM © © i CM CM © © © © P0 ro 1 ©© © 1 CM CM ©©
u. in in in ©©© © ©© ©©© ©©© ©©© © ©© © ©© © ©o © ©© ©©
L. rH rH
o o- O' O'
M3 © rs rs cm rs rs N CM IS cm rs M3 © CM CM S3 M3 rs CM 'C CM © rH © ©© M3 •O M3 © M3
M3 © O' O' m O' O' O' PO O' ro O' *0 © ro ro '0 M3 ro ro S3 P0 © O' © ©© M3 S3 M3 © M3
ro P0 ro © ©© © ©© ©©© o©© ©©© © ©© © o© © ©© © ©© © ©
PO © © © rH © © ©
rH —4 H rH rH
rs ©
rH
is ©
rH
PO © is is ro ro 773 rs © © ©© ro ro ro © ro o ro ro ro
© 0 0. P0 ro
*“ c
ro
rH
© © ro
rH
© © ro
H
© MOM
*4 rH
o© ro CM
rH
© CM
rH
ro © ro ©
pH
1
in in rs ro is ro ro rs rs CM rs M3 •h CD 03 S3 M3 cm rs *o M3 © ro © o* O' 'O S3 M3 © © ©© © © © © ©© © © ©© © ©
rH in cn O' rH O' in ro O' ro ro m 'O O' CO 00 S) M3 ro ro >o >0 © ~4 © O' 4T MD S3 M3 ©© o* «r oo ©©© O' O' © © © © © © O' O' O' O'
© ©© © ©©© ©oo ©© © ©© © ©© © © ©© ©©
©U
O' O' O' O' O' O' O' O'
ro 1 -pH © CM © co 87 ©77 ro o- 433 773 ro i rs i 22 ro 3 3 1 00 CD o' ©©© CD 00 © © © © © © 00 00 ©© © ©© o ©© © CD CD © o© ©©
rH rH
1 1 1 1 1 1 1 1 1 pH 1 1 1 1 1 1 1 1 1 • 1 i i •
1 i
^
i
© © ©
i i
rH
i i
SP
pH pH ** 1 o ©
1 © © © © 1 1 1 1 © 1 1 1 1 1 i 1
1 i
1 1 i
SURFACES
1 1 r 1 1 1 1 1 1 1 1 1 1 i
r-4 cm ro "T in sd 00 O' o rs ro o' in M3 rs 00 rH rs ro O' in M3 IS GO O' CM © rs ro O' in S3 is 03 O' o pH rs ro 4T in S3 is CD O' © rs ro «r in s> rs 00 O' © pH rs
rH rH rH rH 4
pH rs rs 4 rs rs rs rs rs rs rs rs ro ro ro *4 pH pH rH pH rH rH pH tH rs rs CM rs rs rs rs rs rs
183
« n
onnconoooonoiflintPJOooooMOMncDMOoooort
HOC
owrojnM'OOoo^oininw'Oooooo^owwinw'Ooooo^
OOOOO'O'OOOOOOO^O'OOOOOOOOOOO-^'OOOOO
OOOOCOGDOOOOOOO GOGOOOOOOOO 000000300000
CTLT
O I I
I I I I ^ O I I I I I I I O I I I I I I 1 ~
I I I
CL CL
z z
o rg o «r Oho co ro in rs *h O' *r -h O' ^rg»H -ooo'Org'O^-O'O'O-©
UJ UJ
^ -o co o oo UJ H L- L- —
FHE ^m^min'Orscoo'O' ©rg~rgmro«-'OrsrsiscD©~<<Nrgrgroin o oo 2 2 3 CL 2 0.0000 0.0000 0.0000
h h M rg hh^cnnhhmmcvi o• oo <c UJ O O UJ
VCONV
• • 2 CJ 03 2 o
o oo
o
o o rg o o fO -oo rg
co rg 0 <r co ro co -h rs *h in h to m in *h >o ^ <r o o-cocN'0 <rcoO'0 ^ro CL
THI THCN^r-cininrscDO'O' o ^ rg rg ro ro m rs rsrsaoo'HHTHrgrgro*- I-
~ r-« ^ *h *h *h hh h rg rj rg rg rg rg ooo CL
.0056 .0047 .0055
oo©
o o o ooo
VMULT
o o• o• © ©o o
• rg o o o o o © O ©o ©
o rgcoco o mh m ©m rs O' in inrsoo
^
o ©
«
h >oo iNcoooorgroin'O^'O tH HH W CI3o oo oo O o ©o o
rg o- ro in in 'O oo co O' O' *h rg ro <r <r rs
o» rg mis rs rs *h rg ro ro o
O ooo o o
THE
•“* »h rg rg HHHHHHHHtNHnrgrgrg 1 © i© i
li.
o oo 1.500 1.500 1.500
o oo ZZ
K ro th in rg rs co O' O' O' o ro in is oo ro <r cs cs o © ro **mrgrscoroo ,
0' 0'»H o• o* o
• o CO o O CD © © © 00 o
o oo tH o «r o © *T © © © <r ©
03 o «r o © «* © © © o* o 3
O 10
— o tH O o tH © 7
SSC
ooooooooo'oooooooooooinoooooooooo^- o
o
oo
oo
CJ
li. 1 I
1
LAFS
o oo 1A 8
o oo O' oo ro o* «r © © © ©
'O rs rs n
rs rs rs 'O rs rs FAFS
FSG omooo
o
in in
id
o o
inr- o oo ©
©© m
ininoomin©©
in fr- l-OOO© OOOlfllilhl-OOO
minooinininoo mm© OOO ©
N 10 21 32
OOO
ooo LHTS
FSS
ooooinmoooooooominoo
o m O'O'^inmo o
in in in in O' O' rg rg
£>» rs o o o o ooomin'orsooo
ininooininino'&'HHiny3o
ooo rg i/in hs ro ro rg in ro
1 11 22
rgoco^rgoDOorgoco
o o o oo o FHTS
ooo
ooo cn HMH«r^HHHrOH
V HMnMWMrgrgMg HMMMnnMH(NN ,
rorororororgrgrgrg ooo isl
22 32
ON i i i i i i i i i i i i i i i i i i i i i ill •
ooo
• •
LSRF
12
O o o o o o ooo o
03 cl. in in © m m © o in in o A 13 23
Z -« hhh ^ ,h th th ^ rg rg rg rg (N rg rg eg rg rg rg ro rororororororororo ro k ooo x 'O >o in -o o in in 'O >o in
*••••*
SURFACES?
2 000 •• uj FSRF
u ooo 03 U.
VALUES?
UJ
u ooo u
<t
O' o in m in in O' o in O' o-
rs rg
ro ck m H NMrochO'HViniDin c h u.
rgo^ NroHHo
Li.
ii3<rNNOvOHHOiDNHO’rNroo'0'OOin^g , ,
U. cl g >o ro in rg ro i in <o ro in
O
AREA
uj ro ro ro 3 uj C0O''<3ro0''O go O' 'O
h h o c rg h ^ osM^rgyironHo. H o t rg O <T
TRANSFER
o- ro
o ro rg «r
ro
n O'
^ HH«r <thh «
«r rg in
<r
ro ow in
«r h H
O'
<r <r ro ro <r *h
o
u
03 CL
3
<L O oO OoO oooo
<r «r
o RELATED
NORTH
CENTER
SOUTH
O NAME
h u 12 3
rg ro v m >0 IS ^Hrgro<rb*3'0fsco0'0’-‘rg rOg*b3>ON©0'OH u o in <r rg ro ^ Li >o rs ao O'
HEAT
hhh h h h n hh rg rg rg
rg
rg rg rg rg rg rg rg ro ro ro 0.
03
HHH CL
HH
<x
ZONE
184
The headings for the air flow opening data are:
The zone related values consist of the array index, the zone name, the
indices of the first and last surfaces which are part of the zone, the first
and last heat transfer surface numbers, the first and last air flow surface
numbers, the height of the zone for pressure calculations, and multipliers for
visible and convective gains.
This report gives a description of each zone (one zone per page). It
begins with the zone name, for example, NORTH in figure V.24. Next is a
listing of the scheduled loads and controls. The data is presented using
the input variable names. The surfaces of the zone are described next. The
first zone includes a description of the detached shadowing surfaces. The
reported U-values do not include surface coefficients. Finally, the air
flow openings are described. ZONE is the name of the zone on the other
side of the opening.
PAGE 4 NBS THERMAL ANALYSIS RESEARCH PROGRAM 82/05/17. 05.38.07. 3.234
DESCRIPTION OF ZONE NORTH
REP=0 00
.
EQP/EL EQPN 1 JAN 31DEC CAP= 1.00E+03 L AT=0 00 RAD* .30 L0ST=0.00
.
SURFACES
This report gives the fraction sunlit area for all heat transfer surfaces.
It begins with the date and solar declination and equation of time values used
in computing the solar position. The values of the cosine of the angle of
incidence on each surface and the sunlit fraction are presented for 24 hours on
the half hour. Nine values from ISSC or RSSC (see VI. C) depending on the They
are in the HTS number sequence. The sunlit fractions apply only to the outside
of envelope surfaces, so partitions and such cannot have sunlit fractions but
are reported to maintain the sequence of surfaces.
186
V V N N V V NV
o o ©© o © © O © © 'O ©
o © ©© © © © © © O' ©
© © ©© © O >0 © © © m©
0 0» © ©
» at *•
©
a.
©» © © © ©
t-i >— in 0^ O' >0 m sd >0 S3
1 1
1 1 ss oo CD in CN «r cn
<r o S3 in cn rs IN o
•
1 i i 1
\ V VV v\ V NV
© © © © © © \n © © © (N ©
© © © © © © «r © ©© ©
© © © © © © in © ©© m©
0» 0* 0 0» © ^ © © -* ©
•> *• *
NH •— t-t O' CD O' in 0 o >o in
'
| 1 1 1 CD CD a >0 <r cn <r cn
S3 -s >0 «r rs cn IN CN
* * • • • • • •
n s X\ \ V \\ \N s\
© © h< © © © © ©© ^©
© © © 1 © © © i ©© 1 ©
© © © ©• © © •
©• ©• ©
o © © © a-4 «-*© «-4 ©
o a. a.» •> a. a. •t at at
LJ V s N V \\ \\ NN
z
*-* © © © © © O' © © © CD © ©
© © © © © CN © © © rc O ©
© © © © O O' © © © O' © ©
C © © © © © © a-t ©©
a. •»
\\ NN \ Nt \\ \N \N
ro
CN
CO
o \\ \ \ \N \\ NN
\ ©© © © © © © © ** ©©
II ©© © © o O' © O O O' ©O
z ©© © © ©© © © ©© ©©
o
-1 *+ 0* 0 00 © © © © ©©
(J b- •*
u 1-
H \\ \\ \\ \ \ s\ NN
-1
z © ©© •“«
© © © *-« ©« ©©
3 © ©© 1 © © © 1 © ©© i
cn © ©© © © ©• © ©©
in
rv * ©. * © © © © i-i
•»
© ©
CD UJ at
«-« CJ HH M Ok »-• in lO S3 in «r
• z 1 1 1
1 CD 1 S3 © <T 1 CJ O'
cn i UJ *o «T CD IN CN CD
z n Ci • • • . • •
o 1 1 1 1
HH CJ
b~ -j z s \ \ N \ N N % N N N N
<x o M
-J LJ © © © © © © © © © © © o
3 a U- © © © © © © © © © © © ©
U w o © © © © © © O © © © © ©
z •
<X HH UJ
z
© © ©. © © ©* © *« 00 ©
u cn
H
* at «t at •t •t at at at at
\
1
o a
Ot UJ II \\ \N II \\ N N 11 NN \N
<z u.
I H •• Ui ©•-* o © UJ -« © © S3 UJ HH o © IN
cn CN <x z © 1 © © z © i © O z © O CN HC
M
1
b- >— © o © ©• © FN* 1— © © IN
U- II <r b- • b~
o 3 o © © o © * © ©
LJ a at at at at a
cn b- U. ct 1— 1 >-< DC m in O' O' cn hh m S3 C
- <T o 1 1
1
1 <X 1'O CD CD <z 1 c-t <r
REQUIRED
Q Ci a <r* S3 o* a CN (N IN
3 z z z . z • •
cn lC <x <i <i
Lu o b- i— -
iZ u. cn \\ \\ cn N\ \ \ cn N \ NN
187
V.B.2.4 RNOS = 14; Report of shadowing combinations
REQUIRED CMB = 82
188
V.B.2.5 RNOS **
15, 16; Dump of shadowing calculations
The shadowing calculation dumps are intended for checking the shadow
calculation process. They produce a large amount of output. Figure V.27
shows part of the output for only one hour. Report 15 provides summary
results such as a record of which subroutines where called and the net
results of the calculations. Report 16 gives the details of the homogeneous
coordinate calculations. The output from 16 is printed as the calculations
occur, and so it is intermingled with report 15. The output consists of
subroutine names (see VI. B) and variable names (see VI. C) with values. To
use this report, follow the calculations through the program listing and
simultaneously through the report.
189
n 1 1
\ \N\
o o o©
O o -h o
0 0-00
o o o
in O' O' 'O
S3 CO CO in
>0 O ©
I I
\ \ \s \
GO
o o in o O'
©
'O
o o <r o PO
o o in © I
O ^ O
O' 00 O' 10
QC CO CD -O
-o h so^r o o
• a • * H © O' ro rN © © ©a M O «h ©© m
II
H
CM CN © • a* • o
©
\ Ns Ns Ns II A II II II II II
X CJ X X X X h*
o o Hh © CJ o o u CJ
o o 1 o © © © © ©
o• o• ©• CJ O u CJ CJ L) © oo
o ©a o o o o o
T-i aH © -I —J -J -1 -J CN o
at at at a» © a
© © *— in II Ns CM \ H\ \H \ Ns aH
© © 1 S3 1-4 © © © s-4 © © 1 00«* ** oo <r 00 «T
A
1
rH «r m 00 I in ^ 00 I
• • • II <X II II II II © in O' in
1 © © © © © 1 1
A A A A A
Ns \X © in © © © © II it
z rs z z z z «
o O' o© • UJ UJ
o cj o© CC cc
o• O' © ©• N II © © © © <x
• • o
tH ©© II A II II II II ©
•> a at at © © © © © 'O <T »-H
© O' S3 O' 2 2 2 2 2 II !l **
© © in © A © A A A A i © ©m
aH S3 i S3 Z © Z Z z Z © ©
• • • * cj © i frH o ooo o hh h CM CM
1 1 S3
H
a O' X z
CN s-H © vH © © M© II II II II ii II II
S X s Ns S3 © oo u u U
a
©
II
11
X
A
II
©
II
©
II
©
II
©
©
II
XIII
CJ
©
cn co co ;i
I
cn
<
;i
1
cn cn
CJ
© II 1- © © © © Ns \m 222 2 m v2 Ns K> N, P0
'O CD © © © © © ©© aaa a © ncc
:> <r © =>
Ns Ns \ Ns «r Z M z z z Z II ro z z zz z Z hZ © z H ZZ © Z
IN Z © © © O' m O'
© T-i ©© a © © H <• i
in rs X fN CN CN • fN © CN II II II II II
o * © ©a CN a -J © • © * ^ a CN '* m z U U CJ U - ©
CJ
a. •> CN 1 © © 1 m a © H-t 0 I XXX X X © HI X X © H-4
1 AII
© © II © II © II © II ©© - CN II en
II
H* © II s-H rH II T-i II
~ CD ~
i
H CD CD ~ CD CD n ro cn «r tH *H at
\ N \ Ns «s* CD © II © © © © ©© © S3 H Z M- <r © S3 S3 Z ©
v-f z © >0 z
© *w IICJ 1- ~ h- h- ~ h- © CN © w O' O' z cm »* rj fN H^ CM CN CJ H
o hh © © o <t XX «
— X « X <r x <l \n CJ O Z© 00 CM CN fN fN fN z CM © z CN S3 rs z
o ©© 1 CJ h~ M ©w w 1- ^ h- ^ H- CN a a ©•
©• ©• ©a z UJ Z A <Z Id <1 Ul «I Ul <1 UJ a 1
* I
©©
I
© <x <t
NWS i cn
O
• o
II
cn
cn uj
H- cc
I <x
<r uj
I-
UJ
2
^ II UJ UJ
<1 •! <1 <1 UJ
K H- |- h- ©
UJ UJ UJ UJ
II <x liJ
Id
“ P-
Ul UJ II
o o oo cn 1 CJ x x ~ x x x x x XX CN
o o oo <r h- II X oX fN X ©X ro K © 1- K © h- H- © CN
o o oo i «“V a-l <r u <r Z CJ U i-l- o u u cj u <r CJ u
© «r S3 fN © II II X it ©
o o o ** © II II II ii:> ^ z> •w © ©
K © © © © z «T Z T-t rs <r Z
O O' -o O' in X
w CD II Z II Z II z || z ro H W CN CN UJ
PO CD © CD O © CD ©2 ©2 © 2 II
•• *o in -o co © II «Z ©Z ©Z ©Z ©^ CN ^II <C CM II ii ii ii II II aH || II II
©© Ul © ©© © © ©© © 2 ©2 2 © © ©© © © © © ©
z ©© © © © IIZ Z2 z Z Z z2 Z
II II z z z z
© <1 © © © © © ©© © © © © ©© © © © © ©
© •* 2 * * * 2 z 2 2
«* CD CO © A A A A A A A
U 1 w Oo 'O © ©
x
h
i © o s: © A •• © © © © © CD © © .. •• a. aa <1 o © •• -
o...
o rs A ••X ©2 A 2 A 2 A 2 A© © aa -J •• A © •* © •• -I
X © © A© X © X © X © X 2 O © © CJ 2©©!©©©©©2CJ © U © 2 CJ © CJ © © 2 0©
© AX©C
h-
o o © A ©A © A © A © A © © X © CD X © © ©©2 2 222AX => x z AX© X ©
© A X X X X © CJ OO© A A A A CJ O A CJ O A A O
cd
<r
>-< in O' O'
o co oo
Z A Z A©
X Ul Z
z
HH
© z
HH
© z
HH
© zx
HH
CJ
AX
X
AAAAAAAAX
xxzxxxxx
ill
A
X
A X
X
© ©UJ©©©©©
i
a o o• © © Ul A CD A © A © A © w © ©
Z * • Id Z Ul 2 UJ 2 Ul z Ul
<X A UJ A Ul A Ul A UJ A
h-
Ifl \ \ S\
190
NWS
w v u
O O OO
o o no
o o ©
in 'O >o >o
cn *<N
cn in rs >o
s S
© O
w(N| ©
© © «r ©
© © in ©
© o M t cn eg
© ro © N rg in eg
• © * ©
© © © ©© ©
© © © ©© ©
© © © ©© ©
© © ^O
© © © ©© © © © © I
© © ©
o*
© ©© ©© ©
O' 'S © © © © «* © *a -o h in
'O
© ©© ©© © cn eg I CN
CN CN CN
II
S S \ N
<c © ©© © cn ©
UJ © © © ©A © © A © ©
CL ©©© ©© © © eg © ©
<Z ii
o © O' © ©
eg A
«T © <r © ©
II A ii ^ 'C 'C 'O 'C
cn eg cn *h M © ©© ©© © eg «r eg «r
cn *h in eg © © rv ©© © eg rs >o is
S3 rs ©© o ©© ©
I i
\ s \\
ii
\
ro O' © *H ©© ©
© C ©© ©
=> <r
2H
O'
2
<r
*h
S3
TH ©
IN
>0 ©© © NSW
O H © »H © ©
OUh II
© O' © ©
X .
© © ©©
— i tH
ro ©
in cc <1 M3 bn <
N3 2 O' O' rv © © © «r <r eg <r
eg ro 'O N3 © © © © in in eg in
© *
© © ©
co © CN
th eg
2H
II
ii cn
CJ i-
© © ©
NWS
~ x x © >H © ©
n cn v o O' ro © © rv © © ©
2A©
I
cn 2 © N3 ro © © eg © ©©
cn a ^ cn uj © -TH © ©© ©
H- 2 cn 2 iX
\ X A ©
<T cn H bn «r
^ © eg cn ’-i i eg O'
O UJ 2o eg a
u: ii ck ro eg o eg ro © *h © fN oO
©u ro * ^ © on © © o- © O' ro ro
cn i © © © ro ©©
II CC
cn <x
cn
• o un <r o*
NSW
© © o©
II II x ii © © ©©
x 2 © © ^ ©©©©in© © © © ©
ii eg a <r
© © *4 © o x in
<r
ii cn
ii
ro
o X
u ^ i
in
ii
u ^ cn ro
cn h
©ho
©bn© © o o ©ro©
CNunbO
«r ro ro
O O © tH»
O' 2 A 2 00 -I n <r ^ u ii
i
tH o tm ecu u ii _j o -J A © © ©© © ro «r eg <r o-
x in 2x X X •
••IN 'O in a
cj cn u
cn cn cn cn
02 •* o
i
© © © o © © © o© © a
*2
•
cn
22 OHO?
1-2
©©©©©©©©©©-H^H
©oo©©©©©©©egeg
cn
» \ \ SN
A
Cl. J h ** •*
OOCOQOOOOOhh UJH
x
O
©
ONeg
0-1 <L OA a? ©
— -i _j •• X H- _J <C i
O -I •* X O -i — x x cn cc H © © fN
u © © cj cn z i- cc o i- cc u cn A CD
a cn I cccc X — O X A — © x a
t-i cn A
J
cn u A CJ A A O O 2 0 O 2 A X
ax XX 2X X UJ X UJ 2 A cn bn VD *0
cn u cn UJ 2 eg «r o-
UJ A eg ns
2
Ul
191
V.B.2.6 RNOS = 21; Hourly dump of primary zone loads data
DB dry-bulb temperature
WB wet-bulb temperature
BP barometric pressure
SPD wind speed
DIR wind direction
TS sky temperature
HB direct normal solar flux
HD sky brightness
HR ground brightness
TG ground temperature
Z LU Ui z © UJ Z© ©
<r © 0- <1 O' © <L CN ©
O ID © o \n TH O CN CN
' rx
•
>0 ID
• •
O’ o >0
•
©
ID ID ID «
* © o • <c 08 • © O'
ID ID ID
II II II
O © O
K iXi t- A t- A
A A A
«»
ID O id in ID >0
© © • • • © .
o o © © © 00 CN 00 ©
II © ii <r II
Ui LJ UJ
H II i- ii »- II
CN
rx
• • » •
O'
• • • •
•* ©• •
CN cn • •
O © © © © 00 'C 00 o
II ID ii «r Ii
0. Cl a.
> II > II >- ii
K a H*A 1- A -
>- Z h- >-ZK >- z 1
<C o © O <L O
O «•>
A A
*•>
IS
, ,
ID © .
CN CN O 00 CN ID O ©© cn rx 0 Oi
©
II ID ii rx ii ©
3 © 3 ND 3 ~4
u. II U. II U II
o A O A O A
> Z3 >- Z 3 >Z 3
© Z
a Z Z
<x <X
a A
O' CN rx
H * o © © t-4 * © 00 tH • 00 ©
CN IS CN© CN O'
II 1 II 1 II 1
r Z z
Li. II U II u II
O © H» O © - 1 o © r-
>- K © >- h- © >- h- ©
<c © <1
A «» A A
(»
• « • * • •
ID O © © *•* ID © 00 «* ID 00 ©
cn CN *-i CN ^
ro ro ro
ii ii ii ii ii ii
x CL Z CL Z CL
H- - UJ r- hh uj - >H ©
Z a©
CJ
ZAO
c
1
ZA O
O
Z z z
rx • • • rx • • • IS • * •
«•> • o © © «K . © 00 «•* • 00 ©
CN o CN O CN <t
ID ID ID
ii CN II CN cn ro II CN CN ro
11 a ©- II A © ©© A ©© ©
II
U IS U O' u rx
o © ID CN m «T is rx O' ro CN
>0
+ CL* in ©• • CL <T rx • cl rx rx ©
«•*© i- « • . ** is H- - •
©
•
4 T
'O CO CN ro rx O' CN <T O'
T-i tH CN CN •rH T-i CN CN tH CN CN
II II II 11 11 ll
x CO Z A Z A
t— CO N CN rx O' w u N rx © o *“ A
1—
1
N rx <r 'O
»- ID CN ro
»- 1— 'O ro rx 1- t- o* IS CN
A • • • A — • * * A • - •
1— •* rx CN ro 1— o CN If) 1- .. OCN ©
CO cc ^ CN CN © -* cn CN CN
CfC © iZ CN CN CN
—LU UI © •*
Z© X© I ©
H- A ^ CN ro cn ro •• h- A t-i CN ro
© <T <Z LU © © © <r <t
z UI o X LU O z© O
HH 3 _J *—<
1—
3— *-•
1—
3 ©
1-
193
V.B.2.7 RNOS “ 21, 22, 23, 24; Dump of zone loads calculations
Reports 22, 23, and 24 give details of the loads calculation process. All
are in metric units. These dumps can be quite lengthly. Figure 29 shows part
of the output generated for one hour. Report 21 should be specified to report
times for the other reports.
Report 22 describes the control profiles and solar gains for the hour.
The report of control values describes the control profile (IV. E) for each
zone. It consists first of the zone number, the number of control temperature
ranges (NRCTR) and the temperature range which will be used for the first air
,
heat balance (EFSTR) .The second line is a list of control point temperatures.
It is followed by the QFSTDC and QFSTIF variables for the control ranges between
the indicated temperatures. The solar gains data consists first of the total
short wavelength (including from lights) flux (W/m^) absorbed by each heat
transfer surface on the outside (QSO), the inside (QSI), and on any movable
insulation (QSM). Then the radiant gains converted directly to convection
(QC), the total diffuse solar gain (QD), the total diffuse shortwave radiant
gain (QS), the total beam sunlight into the zone (QBV), and the total diffuse
sunlight into the zone (QDV).
The following dumps (23) describe the heat balance calculations. They
begin with the airflow calculations. The following data is presented for each
airflow opening (I).
Then the partial outside (KOP) and partial inside (KIP) conductive fluxes
(C.1.7 and C.1.5) are presented for each heat transfer surface. These values
are constant for the hour. The convection coefficients (HC) are also presented
for each surface followed by the outside surface temperatures (TO) and then the
inside surface temperatures (TI). The results of the zone air heat balances
are summarized under ZON:
194
:
The FLOW, TI and ZON values are dumped for each heat balance iteration.
HC and TO may also be dumped if new values are computed.
TIME! STDT I M= 1 7 CLKTIM=17ff D A YOFS= If BAYOFY= 52f HONTH= 2f PAY0FM=21I DAY0FW = 2 f DAYT YP= 2f ci
WEATHER: DB = 8.89J WB = 3,86. BP=100000.J SP0 = A. 7 > DIR=315.f TS= -8.12f HB = 675 . 5 1 HD= 40,9 t i
QC QD QS QBV QDO
i 0 . 1 274E+02
, 7.130E-01 0. 1 .003E+02
2 0. 8 084E+02
. 3.801E+00 2.053E+03 3 ,009E+02
3 0 . 4.971E+02 2 737E+00
. 1.05 2E+03 2 .006E+02
flow: I N M PS PW PZ DP FLOW
1 l 0 18.06 2.57 19.75 4.25820 161 .67
2 1 0 18.72 -1.14 19.75 -.11059 -7.42
3 1 2 .06 0.00 19.75 -.98394 -153.97
4 2 0 19.60 -3.30 20.79 - 2.10995 -125.08
5 2 0 19.69 -1.27 20.79 -.16838 -9.67
6 2 3 .02 0.00 20.79 -.01496 -18.83
7 2 1 -.06 0.00 20.79 .98394 153.97
8 3 0 17.99 -2.57 20.83 .27151 27.01
9 3 0 19. 10 -3.61 20.83 - 1 .88684 -46.41
10 3 2 -.02 0.00 20.83 .01496 18.83
kop: 91,796 28.722 .000 54.964 81 .925 46.633 771.327 771.327
771 .327 46.261 91.796 59.301 ,000 82 829 . 46.871 91 .796
199.303 771.327 771 .327 771.327 46.923 91 .796 41 .915 .000
85.136 133.525 76.771 91 .796 199.303 771 .327 771 .327 49.593
kip: 71 .786 71.192 -.000 138.661 131.319 57 395 . 698.483 603.115
603.115 39.761 91 .968 70.383 -.000 162.679 73.416 93.171
197.591 760.825 747.766 765.739 44.098 96.059 95.667 -.000
192.235 174.818 77 400 94.865
. 202.337 765.428 761 .470 794.945
195
OAHFL 4.559E-02 8.418E-03
©
© © © ©© © ©© © ©© © © © © © © © o© © ©© © © 0.
4- 4- b b b b b 4* 4* 4- b 4* 4- 4- b b b 4- b +
i
4* 4-
Id Id Id id Id Id Ld Id Id Id Ld Ld Id Id Ld Ld Ld Id Ld Ld Id Id Ld Id Ld
co rv rg rv 44 44 <3 00 CM <r rv CN m O' >0 © <T in 44 rv CN X >0 O (N
rv © © rg *o in m O' 44 ro cm in •4 rg m 44 rg O' >0 CN ro iv in >0
co © rv © co in «o in <r ro rg CM rg in -• CN cn rg rg V CM m rg
BB
O' CM - CD CN - in ^ - m - CN CM in rg - CN CM CN in CN
i 7 1
©o 44 O © ^ © © ©^ © 4- o 0. 0.
©© © © © © ©o o©© © © © © © © © ©© © ©© © ©© O.
b 4- +- b b b 4- 4- 4- 4- b 4- 4- 4- 4- 4- 4- 4- 4- 4- 4-
i 4- 4- 4-
id Id Id UJ Id Ld Id Ld Id Id Ld Id Ld Ld Ld Ld id Id Id Id id Id Ld Id U. Id
<r rv *0 *r in © >o rv ro <r CM 00 O' <3 <3 <3 O' <r in iv 44 4T 43 X O'
© O' 44 © O' <r in o O' O' <3 44 rv ro m rv CM © CN X X O' m -X
vO © rv cd m rg m rv CN ro © CN CN m 43 O' CN in «r r. LI ro in
RA
00 in - 44 ^ + in 44 CN - 44 - CN CN in CN m X CN ro rg in rg ro
1 i i i 1 i
i i i
© o 44 44 © ^ © 44 o 44 © 1-4
0. 0. 0.
© © © © © © ©© © ©© © ©o © © © o ©© © ©© © ©©
4- 4- b -b 4- 4- b 4- 4- 4- 4- 4- 4- b 4» 4- 4- 4- 4- 4- 4- b 4- 4-
id Id id Id u Ld Ld Ld Id Id Id Ld Id Ld Ld Ld Id Ld Id Ld Id Id Id Id Id Id
rv o CM 44 «r >0 ro © O CO CM 4“ rg < <r 44 rv X O' CN © <3 *0 LI
in ao in cn ro 'O in O' <r <r o rg X <r m X 'O 43 iv m rv <r ro LI >o rg
m >o co O' © rv in 'O CM CN O 00 CN CN m © © O' ro *0 -• ro LI ro 44 0T
m O' - in ^ * 44 CN CN CN m CN ** 44 CN ro CM CN in rg ro
i i 7 1 1 i
•4 44 44 .4 © -• o H 44 © 44 44 44 44 o 0. 0. 0.
©© © ©o o ©© ©©© © © © © © © © ©© © ©© © ©©
4- -» 4- 4- b 4- 4- 4- b b 4- b 4- 4- 4- 4* 4- 4- 4- b 4- b
4 4- -b
-
Id Id U Id Ld Ld Ld Id Ld Ld Id Ld Id Ld Ld Id Ld Ld id Id Id Id Ld LU Id Ld
co ro ro 4 4 rv •o rg 44 rv o 44 >0 ro 44 rg <r rv 44 X © — O' X 44 43
rv ro © in O' O' in O' rg go m X CM X <3 O' X © rg rg 44 o rg © ro
rv >o CN ro rv rv m rv 'O o © rv CN O' O' ro ro rg ro <r >0 ro
HW
«4 H -
i
d m <r - CN CN CN o- 2. 6. rg (N rg CM CN * cm rg
i 7 1 1 1 1 7 i
CO GO co ©m o co rv <r rg rv ro 44 44 © 44 © -h © © © -1 44 © 44 O — -4 0. 0. 0.
rv rv o* >0 CO O' LI 43 CD N MN
<T CN O ©© © ©© © ©© ©© © © © © © © © © ©© © o© © ©o
6.140 © o o CN *0 *0 LI CM o o
u. rv rv oo 4- b b b b b b 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4- 4-
Id UJ id Id Ld U Id Ld Ld Id Id Id Id Id Ld Ld Ld Id Id Id Ld Ld Ld Id Id Id
roro<rrvroouiivrg O' ©© O' © IV O' <, O' >0 O' ro ro ro <r rv X «r <r O' rg X X m •4 rg <T V 44
-4 cn *h rg CN IN CM <3 <3 rv «r O' in rg CN in ro IV © co X «r rv in li ro rv rv iv rg — ro
i i i O' 'C ro ro O' rv in rv O' ro m in X ** ro <r cn in O' «r rv rv <r © ST
44 *4 - m 44 <4 in ^ <r CN CN o- 44 ro CN rg CN CN rv rg rv rg rg rg
rococo o0 o- <r rv ^ T-< »-4 i 1 1 1 i i i i i
xrvfvcNrsuix-^rvrox O'NOO
6.140
rg © © O' -« © rg ro iv in _j 'OO'O' 44 -« © 44 © ^ © 44 44 ©© 44 44 44 44 -©
rgrororg©roro&.rgroro
a*.*
© co co
©© © ©© © ©© © © © © © © © © © © ©© © ©© c o©
4- 4- b b + b 4- 4- 4- 4- 4- 4- 4- 4- 4" 4- 4- 4- 4- 4- 4- 4- b 4* b b
0, 0. 0.
CM •* ih rg CN C4 cm li -• Id Id Id Id Ld Id Id Id Id Ld Id Id Ld Id Ld Id Ld Ld Ld Id Id Id (d LU LU Id
*H O' li o- ro rv o 44 rg «r CM rg <3 m X <3 X ro rg < in rg x — CN
O' CO ro iv m in o- <r m <3 O' ro © ro in >0 © O' <3
111 43 O' 43 ro rv
© *o 'O rg 00 rv ro © ** cm ro <r CN in ro >0 rg in
u*) © rg CN <r >
K GO CD rv oo O' CM «r rg O' CM GE
rv rv in -o O' 'O cm ro O' 43 rg cm oo
^
-- - «r 44 ** in 44 O' ro rg rg CM CN in >0 - rg cn - rg rg rg rg x
6.140
co © © H HO> © rv <r O' ©U co rv i i
• • £
ro ro ro © © ro ro in ro ro z) o- co ro © *4 44 44 44 44 © 44
O' © *4 © © 44 44 44 44 ©© 0. 0. 0.
»-«** (N CM *h C4 cm cm co ©^ ro ©©
li © c© O o© ©© © © © © © ©© © ©© © ©© © o©
ro © 4-
CM 4- b b b b 4- 4- 4- 4- 4- 4- b 4- 4- 4- b 4- 4- 4- 4- 4*
1
ii Id u Ld U Id Ld Ld Id Id Ld Id Id Id Id Ld Id Ld Id Id Id Ld Ld id Ld LU
Id
X O' ro <r >0 IV © CM O' o ro 44 O' ro rv >C Hi 4? ro t in rv cc ro >0 «r
O'
GO ©
GO O' oo o ro < O' oo rg O' CD 44 44 <r in CM in -0 © X >0 ro in 43 © < © in rg ro X Li 43
Id 43
rv cm rv «r li in 43 «r 43 ro < cd ro ro © O' rg rv rv ©
H- rv CM in ro ro in ro ro ro © li ro rg LI O'
6.140
© O' © rg rv li >0 <3 GO CM 'O ©Hfi EE
O' 44 - -- 4* <r 44 in rg rg rg CN m >0 - rg cn ro 44 rg rg LI X
1.353E102 4.979E+02 3.040E+03
© © -«
ro ro © © in cm in co ^
ro _i 1 i
X «0 44 - ^ rv
2.701E101
eg CN 44 44 CN CM CD * • • in rg rg
«r © ©
S.915ET01
ID 1 i i i i i
1-
m i4 *4 ©
44 44 © © © © © 44 *4 44 44 ©©
co CD ro ro «r ro <r *T «T ro rv ro ©© © © © © ©© ©©o © o © © ©© © ©© © ©© © ©©
rv rv CD o <3 «T rv 43 43 <r *4 •• b 4- 4- b 4- 4- 4- 4- b 4- 4- 4- b 4- 4* b 4- 4- 4- b 4- 4- 4-
6.140
o© CN rv •o ro CN <T 43 ro O' CD Id Id Id Ld Ld Ld Id Ld Id Id Id Ld Id Ld Id Ld Ld Id Id Ld Id Ld Ld LU Id Id
Id © o* in ro 44 © IV >0 ro © <T © <r <T rg <3 T 44 CO © «r 44 <r <3 •O Li
ro ro CM O' CM CM m ro cm CN ro r- ^4 uo in rg rv ro rv O' in 44 ro © © O' CN >0 m © X LI LI X O' CN ro li ro
CM CN *4 CM CM CM © O -4 © <3 rg O' rv m rv rg rg ro CM <3 >0 m © «T O' Li o- <r © LI o
SN 2.040EL03
h- v4 ro ro »- 44 44 - 44 oo 44 in m - - 44 CN 44 in in CN rg w- rg 44 rg LI O'
u CD i 1 1 i
i I
O 43 >0 a »4 h c mn c ui n -< -* rg ^ iv
140
©good <r x o o cz li rg rv
co co
<rvvo>-o-oiiiO"0'Oirw • • •
44 © © go cm 4> cn co cm 43 ©
«r «r
li rg ro
z
<£
o pg ri in h h
II iii
III (III!
t-i in rg*-4
CM CM Id
ro ro ui rv ro «r in N N ro
TZ
-O -<3 X
H ^ ^ ^ ^ eg (>j X 20.47 22.78
25.20
o ro ^
N V* CCD•
h* •
••
co rg
^ rg cm
1 2 3
LOAHS!
cj o o
x h- M ^4 cm ro
196
VI COMPUTER PROGRAM
The input processor (IP) is the first of the two TARP subprograms. It
checks and interprets the user input, creates the input files for the loads
calculator (LC), and writes various reports. This processing involves all
of chapter IV and primarily sections B, C, and D of chapter III. Data about
the location and weather goes to the LC in the ENVTFL (environment file).
The description of the building, including library data, goes in the SINPFL
(simulation input file).
AFS subroutine
Evaluate the subkeys of AFS (4.14). The AFSDAT common block contains
data to match interzone openings. Data is output in the AFSIL common
block.
197
. .
BBH subroutine
BEAM subroutine
BUILDG subroutine
Evaluate the keys of BUILDING (section 3). This is the root subroutine
of the building segment.
CDS subroutine
CHKCP subroutine
COMPMS subroutine
Compute the optical properties for the pane-shade system for both
polarizations of light incident from the outside. Compute overall
transmittance and reflectance and individual absorptances
198
CONS subroutine
CPFL subroutine
Evaluate the subkeys of CPFL (2.6). Call CHKCP to check the control
profile. Data is output in the CNTPFL common block.
Transfer constructs reports (2.1.3) from unit SCRFL to the output file.
CTRL subroutine
Evaluate the subkeys of CTRL (4.9). Data is output in the CGSIL common
block.
CTF(I,CONVG) subroutine
CWS subroutine
DCHXYZ subroutine
DDY subroutine
DER subroutine
199
.
DIFFUS subroutine
DIFSHD subroutine
DLT subroutine
Evaluate the subkeys of DLT (4.10). Data is output in the CGSIL common
block.
DMPAFS subroutine
Dump air flow opening data in the AFSIL common block. Dump initiated by
RN0S=6 (1.2.1).
DMPCGS. subroutine
Dump zone control and general schedule data in the CGSIL common block.
Dump initiated by RN0S=6 (1.2.1).
DMPHTS subroutine
Dump heat transfer surface data in the HTSIL, HTSIS, and HTSISL common
blocks. Dump initiated by RN0S=6 (1.2.1).
DMPSRF subroutine
Dump surface geometry data in the SRFIS common block. Dump initiated by
RN0S=6 (1.2.1).
DMPSSC subroutine
Dump special surface data in the SSCISL common block. Dump initiated by
RNOS=6 (1.2.1).
DS subroutine
200
DSS subroutine
Evaluate the subkeys of DSS (3.2). Coordinates are converted from the
building to the global coordinate system by (E.1.1) and (E.1.4). Data
is output in the SRFIS common block.
DUMPGL(MODE) subroutine
DUMPLB subroutine
DUMPNV subroutine
DUMPPJ subroutine
Dump project data: common blocks TITLES, RUNCTL, EFLHDR, and REPORT.
Initiated by RN0S=1 (1.2.1).
DUMPRF(MODE) subroutine
DUPLIB(LO,LEN) subroutine
EMI subroutine
201
: )
EQP subroutine
Evaluate the subkeys of EQP (4.4). Data is output in the CGSIL common
block.
FAN subroutine
Evaluate the subkeys of FAN (4.8). Data is output in the CGSIL common
block.
FINISH subroutine
GEOM subroutine
Evaluate the subkeys of GEOM (4.1). Values are set for the subsequent
geometric processing of zone surfaces.
202
.
Skip NSKIP days of weather data from the BLAST 3.0 weather file
(WTHRTP, unit 3). Transfer the following NDAYS of data to the ENVTFL.
GETLIN subroutine
Get the next line on user input from the INPUT file. Echo the line in
the output file. Print error message if data is being skipped.
GLASSX subroutine
GRND subroutine
Evaluate the subkeys of GRND (1.4). Data goes into ENVTFL header
record.
ICTF subroutine
Use modified false position method to compute root (RT) in the interval
bounded by A and B. Convergent solution is found the upper bound -
lower bound < TOL. If solution is not found, set REPEAT to true.
IMI subroutine
INF subroutine
Evaluate the subkeys of INF (4.5). Data is output in the CGSIL common
block.
203
.
INITCF subroutine
Compute the units conversion factors. Most factors are derived from
CFL(2)“.3048 CFM(2)=. 45359237 and CFC(2)«4 . 184 to give a consistant
, ,
INITGL subroutine
INITRF subroutine
I SHADE subroutine
If JFLAG = 1, compute the Julian date (JD) from the year (YR), month
(MO), and day-of-month (DM) (B.1.2). Otherwise, compute YR, MO, and DM
from YR (B. 1 .3)
LIB subroutine
204
LIBRRY subroutine
LIT subroutine
Evaluate the subkeys of LIT (4.3). Data is output in the CGSIL common
block.
LOC subroutine
Evaluate the subkeys of LOC (1.3). Data is stored in the weather file
header (common block EFLHDR)
Transfer library reports (2.1.3) from unit SCRFL to the output file.
Each report element is NLINES long and is given a TITLE. Close the
SCRFL when finished.
205
LYF(YR) integer function
Compute the leap year indicator from the four digit year number (YR)
If YR is a leap year, set LYF to 1. Otherwise set LYF to 0.
MASS subroutine
Evaluate the subkeys of MASS (4.12). Thermal properties are read from
UNIT 5. Data are output in common blocks HTSIL and HTSISL.
MATL subroutine
Evaluate the subkeys of MATL (2.4). The input data are checked against
appropriate ranges. The material properties are written to UNIT 4 for
reference by the CONS subroutine. The materials report (MATL - 2.1.3)
is written to UNIT 23.
MATRIX(W) subroutine
MIX subroutine
Evaluate the subkeys of MIX (4.7). Data is output in the CGSIL common
block.
OSC subroutine
206
PAGE (I ACT, NEED) subroutine
Write output report page headings. This routine can be used to count
the number of lines being printed on issue page eject commands at
appropriate intervals according to the use of IACT and NEED.
IACT“0: Initialize page routine. Set page length to NEED.
IACT=1: NEED lines are to be printed. If there is insufficient space
remaining on the page, eject page and print heading.
IACT=2: Return page number in NEED.
IACT=3: Return number of lines left on page in NEED.
IACT=4: Eject page and print heading if NEED lines not left on page.
Do not increment line count.
The page heading includes a report of cumulative program execution time
determined from the non-ANSI function SECOND. This function can be
removed without otherwise affecting the program.
PANE(N) subroutine
PEO subroutine
Evaluate the subkeys of PEO (4.2). Data is output in the CGSIL common
block.
PROJCT subroutine
PSYCH subroutine
207
PSYDPW(W,PB) real function
Compute the dew point temperature (TD) from the humidity ratio (W) and
the barometric pressure (PB) (B.2.6).
Compute the enthalpy (H) from the dry bulb temperature (TA) and
humidity ratio (W) (B.2.9).
Compute the relative humidity (RH) from the dry bulb (TA) and dew point
(TD) temperatures (B.2.7). Check for 0 < RH < 1.
Compute the wet bulb temperature (TW) from the dry bulb temperature
(TA) ,
humidity ratio (W) and the barometric pressure (PB). TW is
,
Compute the specific volume (V) of moist air from the dry bulb
temperature (TA) , humidity ratio (W) and barometric pressure (PB)
,
Compute humidity ratio (W) from dew point temperature (TD) and
barometric pressure (PB) (B.2.2). Check for W > 1.61E-6 (W for TD =
-70C).
Compute humidity ratio (W) from dry bulb temperature (TA) and enthalpy
(H) (B.2.5) . Check for W > 1.61E-6 (W for TD = -70C).
Compute humidity ratio (W) from dry bulb (TA) and wet bulb (TW)
temperatures and barometric pressure (PB) (B.2.3). Check for W >
1.61E-6 (W for TD = -70C).
208
.
Compute humidity ratio (W) from dry bulb temperature (TA) relative
,
humidity (RH), and barometric pressure (PB) (B.2.4). Check for W >
1.61E-6 (W for TD = -70C)«
QTF subroutine
RC subroutine
Evaluate the subkeys of RC (1.1). Data are stored in the RUNCTL common
block and used in subsequent processing of user input.
READER subroutine
Decode the user input (in UINREC) according to input syntax rules
(IV.B.1). Output is current status (UISTAT) and individual words in
the AWORD array. The longest element of input processed in one call to
this routine is a subkey with data.
REPTGL(LSC) subroutine
REPTRF subroutine
RESPNS subroutine
I
RFCOMP subroutine
209
RPT subroutine
Evaluate the subkeys of RPT (1.2). Data is stored in the RFLAGS array.
SEARCH subroutine
Determine upper and lower bounds within which a root must exist. Call
ILLINI to converge to root. Trigonometric identities are used (with
APPROX) to evaluate AS and BS as PT is incremented. The search for
roots of BS is helped by the fact that roots of BS must lie between
roots of AS.
SETUPV subroutine
SHD subroutine
Evaluate the subkeys of SHD (4.15). Data is output in the SRFIS common
block.
210
SIM subroutine
Evaluate the subkeys of SIM (3.1). Data is stored in the RUNCTL common
block and output as one of the last steps of input processor operation.
SLR program
SLR stands for "simplified language reader". This is the root program
of the input processor. It initializes program variables and opens
files, directs the processing of input sections, prints final reports
and closes files. The variable SEQ is used to insure that sections are
not used out of sequence.
Compute sky temperature from dry bulb (TA) and dew point (TD)
temperatures and tenths cloud cover (N - integer) (B.3.1 through B.3.4)
SPCTRM subroutine
SRF subroutine
211
SUBKEY (KEYS LKEYS ,NKEYS ,LOC)
, subroutine
This subroutine evaluates the subkey name (syntax element [IV>B>1] 4).
The subkey name is always the first word in ALIST. KEYS is a single
string of all subkey names, each no more than LKEYS characters long.
The number of names in KEYS is NKEYS. LOC is the word number location
of the subkey name in KEYS. The INDEX function is used in determining
LOC. If the name is not found, an error message is printed. All valid
subkey names are optionally (DEM- 1.1.3) printed.
SUDDEN subroutine
This routine fills a single day of weather data in common block ENVDAT
from the user's dedign day description (DDY - 1.6). Calendar values
are computed. A dew point temperature (TD) for the day is computed.
It is used in computing the hourly wet bulb temperatures, humidity
ratios, and sky temperatures. Solar intensities are computed.
SUN2(DAYYR,A,B,C,AVSC) subroutine
SUN5(H0UR) subroutine
VENT subroutine
Evaluate the subkeys of VENT (4.6). Data is output in the CGSIL common
block.
212
WDWPOL subroutine
WS subroutine
WTAP subroutine
Evaluate the subkeys of WTAP (1.5). Get weather from the BLAST weather
file unless SYNTAX (1.1.2) or OLDW (1.1.6) are in effect.
WXREF(A,I) subroutine
Write cross reference data consisting of the array name (A) and the
location (I) in that array. WXREF is initiated by XREF (1.1.5). It
relates the input keys to the arrays reported by DUMPIN in the loads
calculator. DUMPIN is caused by RN0S=8 (1.2.1).
ZERORE subroutine
ZONE subroutine
213
This diagram (continued over several pages) shows the general relationship
of the input processor subroutines to each other. A subroutine calls other
subroutines which are to the right and below it on the diagram. Multiple calls
to the same subroutine are not shown, only the first. Circled letters are con-
tinuation points for often used groups of subroutines and for large groups of
subroutines. Not all subroutines and functions are included. The intrinsic
FORTRAN 77 routines have been left out. Several simple utility routines which
are called from many places are not always shown. These include: ERROR,
INTCON, I ZERO, LOCLIB, MODY, PAGE, RELCON, WXREF , and ZERO.
BDSLR
BDPROJ
BDLIB
SLR —
INITCF
PAGE
READER
-BLANK
-GETLIN
PAGE
ERROR
-ERROR
-INTCON
-VALID
— PROJCT
-READER
-VALID
-RC
(B). l
-READER •(A) I
-VALID I
PAGE
j
-ERROR I
-RPT J
-SUBKEY (B)
-LOC
•SUBKEY (B)
—GRND
-SUBKEY -—(B)
— WTAP
—(B)
~k SUBKEY
ETB3W
DUMPNV
l—FERROR
—DDY
1
l-DUMPPJ
Figure VI. A. 1
SUBKEY
SUDDEN
DUMPNV
—(B)
(C)
214
-LIBRRY—
-READER- •(A)
-VALID
-LIB.
-SUBKEY (B)
-DS.
SUBKEY. (B)
DUPLIB
-WS
-SUBKEY- (B)
-DUPLIB
— MATL-
SUBKEY- (B)
-DUPLIB
-CONS.
-SUBKEY- (B)
DUPLIB
•INITRF- (D)
•INITGL- (E)
— CPFL-
•SUBKEY -
(B)
DUPLIB
CHKCP
-CDS.
•SUBKEY- (B)
•DUPLIB
-CWS.
SUBKEY- (B)
•DUPLIB
1-DUMPLB
— BUILDG
-READER •(A)
-VALID
-SIM
*— DSS
SUBKEY- — (B)
SUBKEY-
DMPSRF
— (B)
-ZONE
I— READER (A)
VALID
AFS
SUBKEY- (B)
DMPAFS
—BBH
SUBKEY-
DMPCGS
— (B)
•CTRL
SUBKEY -
(B)
DMPCGS
215
216
DYS
JDS
IDSTF
SUN2
PSYWTP
t SATUPT
PSYWTR
SATUPT
PSYDPW
SATUTP
PSYTWD
PSYDPW
t PSYWTP
SATUPT
1 SATUPT
PSYWTR-
SATUPT
PSYWPD
“LSATUPT
•SKYTF
SUN 5
•RESPNS —
-ZERORE
-SEARCH
1 MATRIX
ILLINI
“L,MATRIX
DUMPRF
-DER
-RFCOMP
ZERO
ICTF
ERROR
1 CTF
QTF
•DCHXYZ — DUMPRF
ICTF
•REPTRF
1 CTF
QTF
—DUMPRF
217
•i * »»
•EXTCF
.FINDFD—
Upane-
-FILM
Lerror
•SPCTRM—
V—DUMPGL
DIFFUS
COMPMS—
PANE-
FILM
r
Uglassx
LisHADE
-DIFSHD—
t -TSBF
RSBF
U*TSBF
Lrsbf
Lbeam-
.COMPMS —
Upane-
•FILM
L•GLASSX
Lishade
-WDWPOL
t RSBF
TSBF
LSQFIT
•DUMPGL
-REPTGL
1 P0LY2F
POLY IF
The loads calculator (LC) subprogram of TARP computes the zone loads
based on the algorithms in chapter III, primarily sections E through L. It
is divided into three sections. The first reads the SINPFL created by the
IP and fills the LC arrays. It checks for sufficient array space and
performs some checks and calculations which require that the entire building
description be present in the arrays. The second section section uses the
geometric and optical descriptions of the building surfaces to compute solar
gain factors for various periods of the year. The third section computes
the zone loads by simultaneous, hour-by-hour heat balances.
AIRMVO subroutine
AIRMV1 subroutine
Compute infiltration and interzone air flows using the detailed air
movement algorithms in III.G.
CALOAD subroutine
Calculate values for the building loads report including system status
hours, zone loads, zone peaks, building loads, and building peak
values
219
.
CHACS subroutine
CHKBKS ( NRS NB S )
, subroutine
Check that a convex zone has been defined. No vertex of surface number
NBS may lie in front of surface number NRS.
CHKSBS subroutine
CISAF subroutine
Compute the fraction absorbed by each surface of the total beam solar
energy transmitted into the zone. This is for the simple solar
distribution which assumes all beam solar gains are incident on the
floors of the zone. The fraction absorbed is proportional to the solar
absorptance of each floor.
CIVAF subroutine
CLIP subroutine
’Clip* the shadowing surface polygon so that none of it lies below the
plane of the shadow receiving surface (E.2). This prevents the casting
of false shadows.
CROSS(A,B,C) subroutine
220
CTRANS(NS,NR) subroutine
DAILY subroutine
Perform all operations which can be done once per day. Read the weather
data and set date values. Set the solar gains file group number. Check
external movable insulation and set pointers if the date is in effect.
Set the schedule and control pointers according to the date. Set the
monthly report flags.
DUMPAF subroutine
DUMPEH subroutine
DUMPHW subroutine
DUMPIN subroutine
Dump all values read from the simulation input file (RN0S=8 - 1.2.1).
DUMPQB subroutine
DUMPSC subroutine
DUMPSG subroutine
DUMPST subroutine
221
:
DUMPSl(I) subroutine
Print the solar incidence angles and sunlit fractions for all heat
transfer surfaces for hour number I (RN0S a 13 - 1.2.1).
DUMPS2(M0DE) subroutine
DUMPS3(M0DE) subroutine
DUMPTC(NZ) subroutine
DUMPTH subroutine
DUMPTM subroutine
DUMPVC(NZ) subroutine
DUMPZL subroutine
222
GAUSSY(A,B,X,N,MAX) subroutine
HBAIR(CNVG) subroutine
HBOUT subroutine
HBSRF subroutine
IBLRPT subroutine
IMRTNW subroutine
IMVINS(I,N) subroutine
223
INITBL subroutine
Direct the set-up of the simulation data. This is the root routine in
the initailization segment. The segmentation is helped by temporary
storage of data on units 4 and 5. The report of array limits is the
first action. It will occur even if there is no SINPFL or ENVTFL
attached. Next, those files are opened and the environment header is
read. Then the SINPFL is read [READER] and the required array space
reported. Some calculations are made that require data simultaneously
for all zones and surfaces [SETIZP, CISAF, and CIVAF] . Dumps [DUMPEH,
DUMPIN] ( RN0S=8 , 1.2.1) and reports REPTZD] (RN0S»9, 1.2.1) are
[
INITSG subroutine
INITTH subroutine
INITZP subroutine
INTCPT subroutine
Determine all intercepts between the sides of figure NS2 and the sides
of figure NS1 (E.3.4 and E.4).
If JFLAG = 1, compute the Julian date (JD) from the year (YR) , month
(MO), and day-of-month (DM) (B.1.2). Otherwise, compute YR, MO, and DM
from YR (B.1.3).
224
MRL( OUTPUT SINPFL ,ENVTFL)
, program
Multi-room loads program. This is the main (or root) program of the
loads calculator. It calls the three segments of the loads
calculation.
MRTNET subroutine
Compute the zone mean radiant temperatures for the second and
subsequent heat balance iterations. Use the simple (H.2.3) or the
detailed (H.2.4) algorithm. This routine includes equations to improve
heat balance convergence.
MRTNT1 subroutine
Compute the zone mean radiant temperatures for the first heat balance
iteration of the hour. Use the simple (H.2.3) or the detailed (H.2.4)
algorithm. This routines includes equations to improve convergence of
the MRT on the first iteration. Development of those equations is not
complete.
MULTOL(NNN,LOCO,NRFIGS) subroutine
Evaluate the heat exchange coefficients HA (=H0), HS, and HG for heat
transfer surface number I. Input also includes the U-value of movalbe
insulation, if present, the solar absorptance of the surface/movable
insulation, and the most recently computes surface outside temperature
TO. Use the algorithms in (1.2). The surface temperature of movable
insulation is estimated from (1.3.8) with HO as input computed from the
simple algorithm.
ORDER(NS3,NV3) subroutine
OVRLAP subroutine
225
PAGE (I ACT, NEED) subroutine
Write output report page heading. See PAGE description in the list of
IP subroutines.
QINSRC subroutine
QSDIST subroutine
QSUN subroutine
Determine solar gains. Use the solar gains multipliers and the solar
radiation intensities to determine the diffuse solar gains in the zones,
the solar gains on outside surfaces, and the beam soalr gains on inside
surfaces. Also evaluate placement of movable insulation for the hour
[IMVINS and OMVINS].
READIN subroutine
Read the simulation input file (SINPFL) and determine the array space
used. Data which was prepared in single element arrays in the IP is
read into multiple element arrays. The reading is directed by
two-element records that tell which arrays will be filled when the next
record is read. The unformatted read statements must match the IP write
statements. Note that each data read is protected to prevent array
overflow. Also note that a simulation control information (created in
FINISH) must be the last record read. Otherwise, an end-of-file is
encountered which terminates the run.
226
REPTBL subroutine
REPTSM(KEY) subroutine
REPTZD subroutine
SETIZP subroutine
SGACTN subroutine
SGBEAM subroutine
Direct the calculation of solar gains multipliers for a day. For each
hour, compute the angle of incidence of the sun's rays on every surface
(E.2.2); compute the solar multipliers [SHADOW]; store the results on
unit 3 for later heat balance calculations.
SGCOMB subroutine
Determine which surfaces can shade each other. This is constant for
the entire simulation and can, therefore, be used to avoid (recomputing
such data at every hour. Data is stored in the SHDCMB array and is
used to direct the calculation sequence of the SHADOW subroutine.
Shadowing combinations are computed for heat transfer, exterior, base
surfaces. All surfaces which can shade the base surface are determined
CHK.GSS .
[ ] All heat transfer subsurfaces are determined [CHKSBS]. For
227
.
Compute the sunlit areas and solar gains multipliers for all building
envelope surfaces at time IHOUR. From the shadow combinations array,
determine the number of shadowing surfaces (NGSS), inside (back)
surfaces (NBKS), and heat transfer subsurfaces (NSBS). If the surface
faces away from the sun, its sunlit area (and all subsurfaces) must be
zero. If it faces the sun and there are no shadowing surfaces, it must
be totally sunlit. Otherwise, compute the shadow patterns on the
surface [SHDGSS]. For window subsurfaces and the detailed solar
distribution call SHDBKS. Then process the subsurfaces [SHDSBS] to
compute their sunlit areas and solar gains multipliers and to determine
how they modify the gross sunlit area of the base surface. Finally,
compute the solar gain multiplier of the base surface.
Project the zone inside surfaces along the sun’s rays to the plane of
the current general shadow receiving surface. This is done for the
detailed internal solar distribution (3.1.2) if there are windows on
the base surface and if the sun might shine through the window onto the
inside surface. This is not possible if the inside surface faces the
sun. The projected figure of the inside surface is reduced to fit
within the base surface by the overlap algorithm (E.4) [OVRLAP]
SHDINT(L) subroutine
228
projections of the internal surfaces along the sun’s rays to the plane
of the base surface and computing the sunlit area of each internal
surface. Tf the internal surface is transparent, its diffuse
transmittance is used to pass sunlight into the adjacent zone. The
simple distribution puts the sunlight on the floor(s) of the zone
according to ISABSF [CISAF]. In either case, that portion of the
sunlight which is not directly absorbed by the internal surface becomes
a diffuse gain for the entire zone.
SHDRVL subroutine
SHDSBS(IHOUR) subroutine
SIMBLD subroutine
SIMDAY subroutine
229
.
SMVINS(K,IHOUR,OABS) subroutine
SOLVZP subroutine
Solve the detailed air flow equations for the zone pressures. On the
first iteration of the hour approximate the zone pressures (PZ) by
solving a linear set of equations (flow exponents set to one). This
gives a reasonable set of starting values for PZ. Then use a modified
Newton’s method to solve the simultaneous non-linear algebraic
equations (G.3.5, G.4.1a to c). This method always produces a
diagonally dominant matrix which can be solved by GAUSSY, a simplified
Gauss elimination routine.
Compute the sine of the solar declination angle (SINDEC) and the
equation of time value (EQTIM) for a given day-of-year (DAYYR)
SUN4(H0UR) subroutine
UTHRMH subroutine
Update the surface temperature and flux history values. All histories
are stored contiguously in the TFH array. New values of temperature
and flux are added to the array by overwriting old values (that will
not be needed next hour) of other surfaces. This causes the histories
230
to "crawl" slowly through the TFH array. Compute the inside (C.1.4)
and outside (C.1.6) conductive fluxes and store in the TFH array.
Compute KIP (C.1.5) and KOP (C.1.7) for the next hour. Shift the
history values back one hour by incrementing TSHIFT. If the histories
have reached the bottom of the TFH array, shift all values back to the
top of the array and reset TSHIFT.
VARTMP subroutine
231
•
This diagram (continued over several pages) shows the general relationship
of the loads calculator subroutines to ech other. A subroutine calls other
subroutines which are to the right and below it on the diagram. Multiple calls
to the same routine are not shown, only the first. Circled letters are contin-
uation points for often used groups of subroutines. Not all subroutines and
functions are included. The intrinsic FORTRAN 77 have been left out. The
simple utility routines ERROR and PAGE which are called from many places are
not always shown.
BDMRL
BDSIM
MRL —
PAGE
INITBL
-READIN
j-IZERO
I— ZERO
Udumpin
•-ERROR
RPTASU
-ERROR
— SETIZP
— CISAF
T ERROR
LOCRAN
— CIVAF
-DUMPEH
-DUMPIN
— REPTZD
— INITSG
•SGCOMB —
— CHKGSS
— CHKSBS-
T CROSS
DOT
-CTRANS
— HTRANS mm$
— OVRLAP (A). r*
i
-INCLOS 1
i
— INTCPT •
i
-ORDER
i -HTRANS ;
i -DUMP S3
l
L-error
i j
•—ERROR i
i 1
— CHKBKS
CROSS
DOT
ERROR
—DUMPSC
-ERROR
232
GACTN —
r
L-s
t GBEAM —UdDUMPS
rsSUN3
ZERO
z
r
L-sHADOW-
-DUMPS2
-CTRANS
-HTRANS
-DUMPS3
-SHDGSS
-DUMPS2
-CTRANS
-CLIP
-HTRANS
-DUMPS3
— OVRLAP (A)
-MULTOL (B)
j|— DUMPS3
t-OVRLAP—
[
—SHDBKS—
I-DUMPS2
[-CTRANS
[-HTRANS
I-DUMPS3
Lovrlap— (A)
— SHDSBS*
-DUMPS2
-HTRANS
-DUMPS3
-MULTOL •(B)
— SMVINS
-SHDRVL "I
I-DUMPS2
I— HTRANS
-POLY IF
I-DUMPS3
Lmultol — (B)
-P0LY2F
-SHDINT ~
1_DUMPS2
-ERROR
L-ERROR
234
c
VI. VARIABLES
short
parameter name description
ALEN maximum number of words decoded at one time
DESLEN - maximum number of characters in a description
MAXAFS AFS maximum number of air flow openings
MAXCDS CDS maximum number of control day-schedules
MAXCGS CGS maximum number of control/general schedule descriptions
MAXCON CON maximum number of constructs
MAXCMB - maximum number of shadowing combinations
MAXCTF CTF maximum number of conduction transfer functions (CTF)
MAXCWS CWS maximum number of conrol week-schedules
MAXDSH DSH maximum number of day-scheduled
MAXENV - maximum number of environments
MAXGTP GTP maximum number of optical constructs
MAXHCS - maximum number of homogeneous coordinate figures
MAXHCV - maximum number of vertices in an HC figure
MAXHTS HTS maximum number of heat transfer surfaces
MAXLNM - maximum number of library names (entries)
MAXNL - maximum number of layers in a construct
MAXNRF - maximum number of response factors
MAXNRT - maximum number of roots for CTF calculation
MAXNWL - maximum number of wavelengths
MAXPFL PFL maximum number of control profiles
MAXPNT - maximum number of user defined points in control profile
MAXRNG - maximum number of control ranges (=MAXPNT+1)
MAXSSC ssc maximum number of special surface coefficients
MAXSRF SRF maximum number of surfaces (ex MASS, inc shadowing)
MAXTFH TFH maximum number of temperature and flux histories
MAXTMP - maximum number of control points (=MAXPNT+2)
MAXWSH WSH maximum number of week-schedules
MAXZON ZON maximum number of zones
RFGLEN - number of report flags
UINLEN - number of characters in one user input record
WRDLEN - maximum number of characters in a word
Nearly all significant variables in the TARP are in common blocks. The
list below describes all such variables by giving the name, common block,
subprogram using the variable (IP, LC, or BOTH), type, dimensions (including
parameters), and units on one line followed by a verbal description.
235
A /GLAS3/ IP real (2,4,2)
Absorptances of single panesof glass
(polarization, pane number, side of incidence).
236
AD /RF4/ IP real
Upper left element of the layer derivative matrix.
AL /RF4/ IP real
Upper left element of the layer matrix.
237
.
AN /GLAS9/ IP real
Index of refraction of air (=1.0).
AP /IZPDAT/ IP real - m2
Area of interzone partition.
AS /RF4/ IP real
Upper left element of the total construct matrix.
AT /RF4/ IP real
Upper left element of the total derivative matrix.
238
BAREA /BASES/ IP real - m2
Net area of base surface (** gross area minus subsurface areas).
BD /RF4/ IP real
Upper right element of the layer derivative matrix.
BL /RF4/ IP real
Upper right element of the layer matrix.
239
BMXH /BLDRPT/ LC real - W
Peak heating load for the building.
BS /RF4/ IP real
Upper right element of the total construct matrix.
BT /RF4/ IP real
Upper right element of the total derivative matrix.
C / SHDCTL/ LC real
X displacement of a shadow in the surface coordinate system (E.2.3a).
240
CC /RF/ IP real - W/m2 K
Layer thickness * density * specific heat (= thermal capacitance).
CD /RF4/ IP real
Lower left element of the layer derivative matrix.
241
CFS /REPORT/ BOTH real
Speed conversion factor: ft/min * CFS(2) = m/s
Type of action:
1 = people
2 = lights
3 = electric equipment
4 = gas equipment
5 = steam equipment
6 = hot water equipment
7 = other equipment
8 = infiltration
9 = ventilation
10 = controls
11 = day lighting
12 = steam baseboard heat
13 = hot water baseboard heat
242
CL /RF4/ IP real
Lower left element of the layer matrix.
243
CPAIR /AIR/ LC real Ws/kgK
Specific heat of air.
CS /RF4/ IP real
Lower left element of the total construct matrix.
CT /RF4/ IP real
Lower left element of the total derivative matrix.
244
CTZ /RF1/ IP real (0:5,0:MAXNRF) W/m 2 K
Z (internal) conduction transfer functions (orders 0 through 5)
DD /RF4/ IP real
Lower right element of the layer derivative matrix.
245
DDMON /DDENVT/ IP integer
Design day month number (1.6.5).
DL /RF4/ IP real -
Lower right element of the layer matrix.
DS /RF4/ IP real
Lower right element of the total construct matrix.
246
DSTIND /ENVDAT/ BOTH integer
Daylight savings time indicator.
DT /RF4/ IP real -
Lower right element of total derivative matrix.
BOTH real - hr
EQTM / SLPSSH/
Equation of time.
247
.
248
FHTS / ZONISL/ BOTH integer (MAXZON)
Number of the first heat transfer surface in the zone description.
249
)
250
HCT /HCFS/ LC real (MAXHCS)
Transmittances of homogeneous coordinate figures.
251
.
JJ / GLAS4/ IP integer
Inner or outer surface indicator.
252
KOP /HTSL/ LC real (MAXHTS) W/m2
Constant part of outside conductive flux (C.1.7).
253
LOCSCM / SHDCTL/ LC integer
Current location in the shadowing combinations (SHDCMB) array.
254
MDHMXT / BLDRPT/ LC integer
Mo nth/ day /hour at which maximum temperture (BMXT) occurs
255
MXMATL /LIBCTL/ IP integer
Maximum number of materials (2.1.3).
NC /IZPDAT/ IP integer
Number of temperature CTF to determine interzone surface match.
256
NDSH / SIMCTL/ LC integer
Number of day-schedules.
NL /RF/ IP integer
Number of layers in construct.
257
.
NQ / IZPDAT/ IP integer
Number of flux CTF coefficients for matching interzone surfaces.
NS SC /RUNCTL/ IP integer -
Number of special surface coefficients (RSSC).
/SIMCTL/ LC
258
. .
259
.
260
OTABS /CONDAT/ BOTH real (MAXCON)
Surface outside thermal absorptance/emittance (2.4.8 and 2.5.8).
261
QFSTDC /TCDESG/ LC real (MAXRNG, MAXZON) W/K
LT in (K.1.2a).
262
s
2 = activity level W
3 = fraction radiant
Type = 2 (lights - 4.3)
1 = design lighting load W
2 = fraction of load to return air
3 = fraction radiant (long wavelengths)
4 = fraction visible (short wavelengths)
5 = fraction of load replaceable *by daylight
Types = 3, 4, 5, 6, and 7 (equipment - 4.4)
1 = design load W
2 = fraction latent
3 = fraction radiant
4 = fraction lost
5 = fraction convected
Type = 8 (infiltration - 4.5)
1 = design volumetric flow m-Vs
2 = constant part of Achenbach-Coblenz equation
3 = temperature part
4 = velocity part s/m
5 = velocity squared part s^/m^
Type = 9 (ventilation - 4.6)
1 = design volumetric flow m^/s
2 = minimum temperature C
3 = maximum temperature C
4 = temperature difference C
Type = 10 (two-zone mixing - 4.7)
1 = design volumetric flow m^/
4 = temperature difference C
Type = 11 (whole house fan - 4.8)
1 = design volumetric flow m^/
4 = desired temperature C
Type = 12 (controls - 4.9)
1 = maximum heating capacity W
2 = maximum cooling capacity W
3 = fraction of control based on mean radiant temperature
Type =13 (daylight - 4.10)
1 = fraction beam usable
263
REFRAC /MATDAT/ IP real (MAXNL)
Index of refraction (2.4.10).
264
RRD / GLAS5/ IP real (2,2)
Reflectances of a window system for diffuse light
(polarization, side of incidence).
S /SHDCTL/ lc real - -
Y displacement of a shadow in the surface coordinate system
265
SA /SRFIS/ BOTH real (MAXSRF)
Sine of surface azimuth angle (4.13.3).
266
SIND / SLPSSH/ BOTH real
Sine of the solar declination angle.
267
.
2 = in heaing mode.
3 = in ventilation mode.
4 = in cooling mode.
5 = cooling capacity exceeded.
6 = no heating, cooling, or ventilating.
268
SUMX /RF1/ IP real (5)
Sum of X CTF to test convergence (C.2.1).
269
TMKUP /EFLHDR/ BOTH real (12) C
Make-up water temperatures (not used).
270
TTD /GLAS5/ IP real (2)
Transmit tances of a window system for diffuse light (polarization).
271
UM /HTSL/ LC real (MAXHTS) W/m 2 K
U-value of movable insulation. 0 means no movable insulation
272
WFSDAY /WFLHDR/ IP integer
Day-of-month of first day of BLAST weather file.
273
.
274
YVT / SHDVRT/ LC real (5) m
Y-coordinates (surface relative) of the shadowing surface.
27 5
. . .
27 6
VI. FILES
277
U.S. DEPT. OF COMM. 1. PUBLICATION OR 2. Performing Organ. Report No 3. Publication Date
REPORT NO.
BIBLIOGRAPHIC DATA NBSTR 83-2655
SHEET (See instructions) February 1983
4. TITLE AND SUBTITLE
5. AUTHOR(S)
G. N. Walton
6. PERFORMING ORGANIZATION (If joint or other than NB S, see instructions) 7. Contract/Grant No.
9. SPONSORING ORGANIZATION NAME AND COMPLETE ADDRESS (Street, City, State, ZIP)
National Bureau of Standards and U.S. Department of Energy
Building Energy Sciences Branch
11.Office of Building Energy Research and Development
The Thermal Analysis Research Program (TARP) has been developed as a research tool
for the thermal analysis of buildings. It especially aims to study the interactions
of many complex heat transfer phenomena. TARP uses the detailed heat balance method
for the simultaneous calculation of the energy requirements of multiple rooms.
Interroom conductive and convective processes are simulated in detail. This program
reference manual describes the algorithms, input, output, and program structure of
12.TARP. The program is written to be portable and modifiable. It is written in FORTRAN
77 and has run on CDC and UNIVAC computers. Future expansions of the program are
anticipated, particularly for the simultaneous simulation of equipment performance
13.and building thermal response.
KEY WORDS (Six to twelve entries; alphabetical order; capitalize only proper names; and separate key words by semicolons)
r~] Order From Superintendent of Documents, U.S. Government Printing Office, Washington, D.C.
20402. 15. Price
$22.50
* [_XJ Order From National Technical Information Service (NTIS), Springfield, VA. 22161
USCOMM-DC 6043-P80