Smoke Management Systems - WTC
Smoke Management Systems - WTC
Smoke Management Systems - WTC
P 2
CA
K
= Q
0
(B1)
where:
Q =volumetric flow rate (ft
3
/min)
C =dimensionless flow coefficient (0.6)
A =flow area (leakage area) (ft
2
)
P =pressure difference, in m.g.
=density of air (0.075 lb/ft
3
)
K
0
=unit conversion factor (776)
A number of smaller openings were also included in the model. These consisted of elevators, stairs, doors,
walls, and floor leakages. Both elevators and stairs were modeled as ASCOS openings, governed by
equation 1. Flow though smaller leakage areas was modeled using the power law equation:
) p c( = Q
n
(B2)
In this equation, p is the pressure difference across the opening in in. of H
2
O, and Q is the mass flow
across the opening in ft
3
/min, and n is the power law exponent. The value of c, the flow coefficient, is
determined using the following equation:
) p ( )
2
( L
C K
= c
n - 0.5
r
0.5
D 0
(B3)
where:
c =flow coefficient
n =power law exponent
p
r
=reference pressure at which leakage factor was determined (in. H
2
O)
K
0
=unit of conversion factor (776)
A =density of air (0.075 lb/ft
3
)
Appendix B
134 NIST NCSTAR 1-4D, WTC Investigation
L =leakage area determined at reference pressure (ft
2
)
C
D
=discharge coefficient
In order for CONTAM to calculate flow through an opening using the leakage area approach, a flow
coefficient, reference pressure difference, and power law exponent are required to be input for each
opening type. For different types of openings, the values listed in Table B2 were used.
Table B2. Leakage area parameters.
Opening Type Size Coefficient
Reference P
(in. H
2
O) Exponent
Floor 0.008 in.
2
/ft
2
0.65 0.1 0.5
Exterior wall 0.03 in.
2
/ft 1.0 0.3 0.65
Interior wall 2 in.
2
/ft 0.65 0.3 0.5
Doorway see Table B3 0.65 0.3 0.5
Large opening 0.82 0.3 0.5
Entrances to openings can be characterized as intruding, flush, or well-rounded inlets. Exits can be
characterized in the same manner. Entrance conditions can be characterized as either plenum or pipe-
flow, self-generated eddies caused by stagnation or separation can occur based on the angle of attack of
the airflow. Virtually all of the opening areas in the WTC complex can be reasonably approximated by
use of an orifice discharge coefficient, C
D
, equal to 0.65. This coefficient is well documented and
applicable to virtually all of the flow conditions where the Reynolds number, Re, is fully turbulent, i.e.,
greater than 10,000.
Some of the very large areas with very low velocities are not well characterized by this same C
D
. In these
instances, it is more appropriate to use values observed by F.C. Lea (Lea 1942) and those of Newell.
These two works lead to the conclusion that for these conditions a number lying between 0.74 and 0.91 is
more appropriate. Given that the experimental data reported has a spread of nearly 15 percent, the value
of 0.82 was selected for C
D
to characterize these openings. The parameters for the various leakage paths
in the model are shown in Table B2. Leakage areas for the various doors in the model are shown in
Table B3.
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 135
Table B3. Door leakage areas.
Description/Use Type
Average
Measured
Leakage Area
(ft
2
)
Corrected
a
Leakage Area
(ft
2
)
Side hinged single leaf door
(stairs, tenant space, etc.)
A
0.1688 0.20
Center parting double leaf doors (tenant space,
escalator cut-off, etc.)
A
0.3376 0.081
Center parting single speed elevator doors
(freight elevators)
B
0.4017 0.4017
Center parting double speed elevator doors
(freight elevators)
C
0.7885 0.7885
Center parting single speed elevator doors
(local passenger elevators) B 0.5673 1.00
Center parting double speed elevator doors
(shuttle passenger elevators) C 1.2014 1.2014
Freight elevator door D 1.2500 1.25
Tenant doors A 0.3866 2.60
Elevator pit doors A 0.63 0.63
Mechanical room door A 1.00 1.00
Automatic closing dock door 1.00 1.00
Rolling steel dock doors 6 in.
2
/ft
2
6 in.
2
/ft
2
Rubber door at dock 1.00 1.00
a. 1996 WTC model.
B.1.3 1996 CONTAM HVAC System
Air handling systems in the 1996 CONTAM model were simulated by using CONTAMs Simple Air
Handling System feature. A simple air handling system consists of two implicit airflow nodes (supply and
return), three implicit flow paths (recirculation, outdoor, and exhaust air), and multiple supply and return
points, each with specified airflows, in the building zones. Outside air supply to a system is set either by
definition, as a fixed percentage, by profile, or by calculation. If the sum of supply flows exceeds the sum
of return flows, the balance is made up from outside air. Excess return flow is exhausted. The model does
not require that both supply and return points be defined for a system.
To facilitate the various scenarios and cases that were simulated, four simple air handling systems were
constructed for each HVAC zone in the tower: Interior space supply, interior space return, core supply,
and core return. (The perimeter and interior systems were combined into a single system in the model).
Additional systems were constructed for the elevator machine room (EMR) exhausts, the mechanical
equipment room (MER) exhausts, the sub-grade transfer fans, and the garage and truck dock systems.
Supply and return points for each system were added to the appropriate zones in the model. The flow for
each point was derived from a careful study of the HVAC plans for the building. The fan schedules and
riser diagrams were reviewed, and the supply and return flow for each zone was determined. To simulate
actual operating efficiency (based on building engineer estimates), the simulations were conducted using
Appendix B
136 NIST NCSTAR 1-4D, WTC Investigation
85 percent of the design flow to each floor. Tables B4a through B4e show the HVAC flows used in the
WTC 1 model.
Table B4a. HVAC system flows in tower model.
Fan Service Schedule Flow
Design Flow
(85 % of Schedule)
Supplies to Sub-Grade
VS-S5-3 MER Ventilation 30,030 25,526
VS-S5-7 MER Ventilation 38,200 32,470
Supplies to B-1B-6
ACS-S5-1 B-1 6,100 5,185
B-2 5,015 4,263
B-2 5,470 4,650
Fan Total 19,340 16,439
ACS-S5-2 Electrical sub station on B-1 13,800 11,730
ACS-S5-3 B-1 12,940 10,999
ACS-S5-5 B-1 7,145 6,073
B-1 2,400 2,040
B-2 4,430 3,766
B-2 3,340 2,839
B-2 6,000 5,100
Fan Total 37,625 31,981
VS-S5-4 Tenant Storage B-4 4,940 4,199
Tenant Storage B-3 4,670 3,970
Tenant Storage B-2 3,320 2,822
Tenant Storage B-4 6,980 5,933
Tenant Storage B-3 6,410 5,449
Tenant Storage B-4 7,265 6,175
Tenant Storage B-3 6,580 5,593
Tenant Storage B-4 5,290 4,497
Tenant Storage B-3 4,670 3,970
Tenant Storage B-2 7,740 6,579
Fan total 49,720 42,262
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 137
Table B4b. HVAC system flows in tower model.
Fan Service Schedule Flow
Design Flow
(85 % of Schedule)
ACR-S5-1 B-1 5,875 4,994
B-2 7,740 6,579
Fan Total 13,615 11,573
ACR-S5-2 Electrical sub station on B-1 13,800 11,730
ACR-S5-5 B-1 8,100 6,885
B-1 5,500 4,675
B-2 8,200 6,970
B-2 3,400 2,890
Fan Total 26,185 22,257
MER Exhaust to Garage
E-S5-7 To Garage on B-4 24,000 20,400
E-S5-8 To Garage on B-4 44,000 37,400
Transfer Fans from Tower A to Garage
T-S3-1 From B-4 thru north wall 24,470 20,800
T-S2-1 From B-3 thru north wall 7,000 5,950
T-S2-2
From B-3 thru north wall on
B-4 14,910 12,674
T-S1-1 From B-2 thru north wall 3,320 2,822
Garage Exhausts
E-S3-S3 90,500 76,925
E-S2-S4 92,000 78,200
E-S2-S5 92,000 78,200
E-S1-S1 90,500 76,925
E-S1-S2 90,500 76,925
E-S1-S3 90,500 76,925
E-S1-S4 101,000 85,850
E-S1-S5 101,000 85,850
E-S1-S6 101,000 85,850
E-S1-S7 101,000 85,850
Appendix B
138 NIST NCSTAR 1-4D, WTC Investigation
Table B4c. HVAC system flows in tower model (100 percent).
Floor Core Supply
Core
Return
Interior
Supply
Interior
Return LMRS LMRX SMRS SMRX
B-4 55,630 68,000 0 0 0 0 0 0
B-3 28,794 28,765 0 0 0 0 0 0
B-2 26,271 21,985 0 0 0 0 0 0
B-1 26,486 12,420 0 0 0 0 0 0
Service 204,451 23,245 0 0 0 0 0 0
1 8,200 2,500 219,000 200,100 0 0 0 0
2 985 985 0 0 0 0 0 0
3 3,820 3,880 0 0 0 0 0 0
4 3,820 3,880 0 0 0 0 0 0
5 3,880 3,880 0 0 0 0 0 0
6 1,680 1,620 0 0 0 0 0 0
9 2,140 1,785 26,060 23,320 0 0 0 0
10 2,140 1,785 26,060 23,320 0 0 0 0
11 2,140 1,785 26,060 23,320 0 0 0 0
12 2,140 1,785 26,060 23,320 0 0 0 0
13 2,140 1,785 26,060 23,320 0 0 0 0
14 2,140 1,785 26,060 23,320 0 0 0 0
15 2,140 1,785 26,060 23,320 0 0 0 0
16 2,245 1,595 26,060 23,320 0 0 0 0
17 2,000 1,675 26,060 23,320 4,530 4,530 0 0
18 2,005 1,675 26,060 23,320 2,270 2,270 0 0
19 2,380 1,785 26,060 23,755 0 0 0 0
20 2,380 1,785 26,060 23,755 0 0 0 0
21 2,130 1,785 26,560 23,425 0 0 0 0
22 2,130 1,785 26,560 23,755 0 0 0 0
23 2,130 1,785 26,560 23,755 0 0 0 0
24 2,335 1,565 26,360 23,755 0 0 0 0
25 2,105 1,565 26,365 23,585 5,930 5,930 0 0
26 2,175 1,565 26,060 23,320 2,970 2,970 0 0
27 1,855 1,565 27,500 24,570 0 0 0 0
28 1,855 1,565 27,500 24,570 0 0 0 0
29 1,855 1,565 27,500 24,570 0 0 0 0
30 1,855 1,565 27,500 24,570 0 0 0 0
31 1,855 1,565 27,500 24,570 0 0 0 0
32 2,105 1,565 27,190 24,570 0 0 0 0
33 1,855 1,565 27,190 24,300 6,670 6,670 0 0
34 1,855 1,565 27,190 24,300 3,330 3,330 0 0
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 139
Floor Core Supply
Core
Return
Interior
Supply
Interior
Return LMRS LMRX SMRS SMRX
35 1,845 1,565 28,030 25,025 0 0 0 0
36 1,845 1,565 28,030 25,025 0 0 0 0
37 1,845 1,565 28,030 25,025 0 0 0 0
38 1,845 1,565 28,030 25,025 0 0 0 0
39 1,845 1,565 28,030 25,025 0 0 0 0
40 12,950 12,350 24,945 21,661 0 0 0 0
43 1,600 2,315 17,390 11,650 0 0 0 0
44 6,960 765 21,475 11,150 0 0 0 0
Key: LMRS/X, Elevator Machine Room Supply/Exhaust; SMRS/X, Machine Room Supply/Exhaust.
Appendix B
140 NIST NCSTAR 1-4D, WTC Investigation
Table B4d. HVAC system flows in tower model (100 percent).
Floor Core Supply
Core
Return
Interior
Supply
Interior
Return LMRS LMRX SMRS SMRX
45 2,585 2,210 25,710 23,010 0 0 0 0
46 2,830 1,785 25,710 23,010 0 0 0 0
47 1,925 1,785 25,710 23,010 0 0 14,300 14,700
48 2,100 1,785 25,500 22,830 0 0 0 0
49 1,895 1,785 27,900 24,760 0 0 0 0
50 1,895 1,785 27,900 24,920 0 0 0 0
51 1,895 1,785 27,900 24,920 0 0 0 0
52 1,895 1,785 27,900 24,920 0 0 0 0
53 1,895 1,785 27,900 24,920 0 0 0 0
54 2,130 1,785 27,900 24,920 0 0 0 0
55 1,760 1,785 27,745 24,620 3,160 3,160 0 0
56 1,760 1,785 27,745 24,620 1,590 1,590 0 0
57 1,760 1,785 28,575 25,505 0 0 0 0
58 1,760 1,785 28,575 25,505 0 0 0 0
59 1,835 1,785 28,640 25,270 0 0 0 0
60 1,670 1,785 28,315 25,280 0 0 0 0
61 2,310 1,825 27,890 25,500 0 0 0 0
62 1,765 1,785 28,225 25,200 4,430 4,430 0 0
63 1,760 1,785 28,575 25,215 2,220 2,220 0 0
64 1,665 1,565 29,290 26,130 0 0 0 0
65 1,665 1,565 29,310 26,145 0 0 0 0
66 1,665 1,565 29,175 26,145 0 0 0 0
67 2,010 1,565 38,755 37,195 0 0 0 0
68 1,580 1,565 28,910 25,760 4,750 4,750 0 0
69 1,510 1,565 29,085 26,445 2,400 2,400 0 0
70 1,525 1,565 29,905 26,660 0 0 0 0
71 1,525 1,565 29,905 26,660 0 0 0 0
72 1,525 1,565 29,800 26,570 0 0 0 0
73 1,525 1,565 29,650 26,570 0 0 0 0
74 1,525 1,565 24,005 26,240 0 0 0 0
77 2,460 2,335 26,795 23,115 0 0 0 0
78 6,350 1,320 37,040 32,660 0 0 0 0
79 2,535 1,735 26,845 23,995 0 0 0 0
80 2,470 1,735 27,125 24,240 0 0 0 0
81 1,780 1,735 27,125 24,250 0 0 19,000 19,000
82 2,205 1,785 26,915 24,060 0 0 0 0
83 1,600 1,755 30,730 27,375 0 0 0 0
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 141
Floor Core Supply
Core
Return
Interior
Supply
Interior
Return LMRS LMRX SMRS SMRX
84 1,600 1,755 30,730 27,375 0 0 0 0
85 1,600 1,755 30,730 27,375 0 0 0 0
86 1,600 1,755 30,730 27,375 0 0 0 0
87 1,605 1,755 30,385 26,845 3,160 3,160 0 0
Key: LMRS/X, Elevator Machine Room Supply/Exhaust; SMRS/X, Machine Room Supply/Exhaust.
Appendix B
142 NIST NCSTAR 1-4D, WTC Investigation
Table B4e. HVAC system flows in tower model (100 percent).
Floor Core Supply
Core
Return
Interior
Supply
Interior
Return LMRS LMRX SMRS SMRX
88 1,605 1,755 30,365 26,825 1,590 1,590 0 0
89 1,605 1,755 31,300 27,870 0 0 0 0
90 1,605 1,755 31,300 27,870 0 0 0 0
91 1,605 1,755 31,300 27,872 0 0 0 0
92 1,605 1,655 31,300 27,870 0 0 0 0
93 1,575 1,585 31,300 27,870 0 0 0 0
94 1,445 1,585 30,525 27,270 4,430 4,430 0 0
95 1,445 1,585 30,525 26,995 2,220 2,220 0 0
96 1,445 1,585 31,275 27,850 0 0 0 0
97 1,445 1,585 31,275 27,850 0 0 0 0
98 1,445 1,585 31,275 27,850 0 0 0 0
99 1,445 1,685 31,275 27,850 0 0 0 0
100 1,665 1,685 31,275 27,850 0 0 0 0
101 1,555 1,685 31,275 27,850 4,750 4,750 0 0
102 1,555 1,685 31,275 27,850 2,400 2,400 0 0
103 1,465 1,765 31,500 28,145 0 0 0 0
104 1,465 1,765 29,585 25,335 0 0 0 0
105 1,465 1,685 28,215 28,065 0 0 0 0
106 1,640 1,825 34,011 30,513 0 0 0 0
107 0 0 5,250 0 0 0 0 0
110 5,000 5,000 0 0 0 0 14,000 14,000
Key: LMRS/X, Elevator Machine Room Supply/Exhaust; SMRS/X, Machine Room Supply/Exhaust.
B.1.4 1996 Site Survey and Physical Testing
As part of the 1996 study, a site survey and physical testing of the WTC complex was conducted. The site
survey was performed on J une 28J uly 2, 1996 with a follow up visit on J uly 23, 1996. The physical
testing was conducted on June 30, between 1:00 a.m. and 3:00 p.m. Temperatures early in the morning
were nominally 60 F 5 F, and the wind was fairly constant at 2 mph (HAI and DCE 1996). The
National Oceanic and Atmospheric Administration reported a New York City temperature range of
60 F75 F on J une 30, 1996.
The purpose of the testing that was conducted was to determine the performance of the buildings air
movement pathways as well as the various barriers to air movement. Pressure measurements were taken
at various locations through the building under various HVAC modes. The HVAC modes included:
Mode 1All fans on
Mode 2Fans set on summer normal
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 143
Mode 3Fans set for core pressurization for the full height of WTC 1
Mode 4Floors 59 through 91 set on core pressurization and all other floors set on
100 percent purge
Site survey pressure measurement results, for HVAC Modes 1 through 4, are listed in Tables B5
through B8.
Measurements of the leakage areas of several types of door assemblies were also taken during the site
survey. Survey measurements are presented in two tables:
1. Table B9 is the raw data recorded during the site survey.
2. Table B10 is the adjusted data with outlying measurements discarded.
Table B11 lists the adjusted average leakage areas proposed to be used for the various types of leakage
paths by the 1996 HAI/DCE study (HAI and DCE 1996).
A
p
p
e
n
d
i
x
B
1
4
4
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
Table B5. Measured pressure differentials in Mode 1.
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23 W X Y Z
1 -0.02 -0.02 -0.05
29 0.01 0.01 0 0.01 0 0.01 0.01 0
44 -0.006 0 0 0.014 0.016 0.006 0.01 0.006 0.03 0.018 0.016 0.038 0.03 -0.018 0 0 0
56 0 0 0 0 0 0 0 0
60 0.006 0 0.006 0.01 0 0 0 0
78 -0.015 -0.02 -0.01 -0.01 -0.01 0.02 0.01 0.01 0.015 0.02 0.015
83 0 0 0 0 0 0 0 0 0
92 -0.015 -0.01 0 -0.01 -0.01 0 0 0 0
104 -0.04 -0.005 0 0 0 0 0.01 0.01 0.01
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 -0.02 -0.01 -0.02
29 0.02 0.01
44 0 0 0 0 -0.012 -0.006 -0.006 0
56 0 0
60 0.006 0.004
78 -0.02 -0.02 -0.01 -0.01 -0.01 0.001 -0.02 0.001
83 0.001 0.001
92 0 0
104 -0.01 0.001
Local Elevator Banks
I-C II-A II-B Floor
Floor
III-D
Tenat Doors Shuttle & Frieght Elevators Stair
II-C II-D III-A III-C III-B
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate air flowing out of a space.
Source: HAI and DCE 1996.
C
O
N
T
A
M
M
o
d
e
l
C
o
n
s
t
r
u
c
t
i
o
n
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
4
5
Table B6. Measured pressure differentials in Mode 2.
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23 W X Y Z
1 0.02 0.01 0
29 -0.19 -0.25 -0.23 -0.24 0.03 0.03 0.01 0.06
44 -0.014 0 0.004 -0.05 -0.066 0 0 0 -0.06 -0.048 -0.066 -0.038 -0.038 -0.074 -0.058 -0.056 -0.06
56 -0.02 -0.07 -0.07 0 0 -0.11 0.03 0.03 0.02 0.04
60 -0.024 -0.08 -0.08 -0.1 -0.09 -0.11 0.034 0.03 0.032 0.016
78 -0.05 -0.07 -0.04 -0.075 0.01 -0.03 -0.02 -0.01 0.01 -0.005 -0.01 -0.015
83 -0.04 -0.03 -0.04 -0.06 -0.028 0 0.03 0.026 0
92 -0.17 -0.11 -0.075 -0.13 -0.075 0 0 -0.015 -0.01
104 0 -0.03 -0.01 -0.035 0.02 0.05 0.03 0.03 0
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 0.03 0.01 0.12
29 -0.2 -0.19
44 0.024 0.044 0.072 0.068 0.07 0.062 0.044 0.062
56 -0.01
60 -0.02 -0.02
78 0.001 0.001 0.01 0.01 -0.005 -0.005 -0.09 -0.05
83 -0.008 -0.008
92 -0.01 -0.01
104 0.01 0.001
III-B Floor
Local Elevator Banks
I-C II-A II-B II-C II-D III-A III-C III-D
Floor
Stair Shuttle & Frieght Elevators Tenat Doors
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate air flowing out of space.
Source: HAI and DCE 1996.
A
p
p
e
n
d
i
x
B
1
4
6
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
Table B7. Measured pressure differentials in Mode 3.
1 2 3 17 49 50 W X Y Z
1
29 -0.01 0.01 0.02 -0.02 0.02 0.01 0.03 0
44
56
60 0.03 0.006 0.1 0.07 0.05 0.05 0.008 0 0.03 0.01
78
83
92 -0.01 -0.025 -0.01 0.03 0.01 0.01 0.01 0.01
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.04 0.03
44
56
60 0.032 0.038
78
83
92 -0.01 -0.01
104
Floor
Stair Shuttle & Frieght Elevators Tenat Doors
Local elevator Banks
Floor I-C II-B III-C
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values
indicate air flowing out of a space.
Source: HAI and DCE 1996.
C
O
N
T
A
M
M
o
d
e
l
C
o
n
s
t
r
u
c
t
i
o
n
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
4
7
Table B8. Measured pressure differentials in Mode 4.
1 2 3 5 6 17 49 50 W X Y Z
1
29 -0.18 -0.23 -0.26 -0.25 -0.32 -0.04 -0.03 -0.01 -0.08
44
56
60 0.006 0 0.076 0.018 0 0.078 0.054 0.024 0.02 0
78
83
92 -0.16 -0.13 -0.055 -0.055 0.03 0 0 0 0
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 -0.19 -0.19
44
56
60 0.114 0.11
78
83
92 -0.11 -0.11
104
Shuttle & Frieght Elevators
Local elevator Banks
Tenat Doors
Floor I-C
Floor
Stair
II-B III-C
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate
air flowing out of a space.
Source: HAI and DCE 1996.
Appendix B
148 NIST NCSTAR 1-4D, WTC Investigation
Table B9. Door leakage path measurements (raw data).
Source: HAI and DCE 1996.
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 149
Table B10. Adjusted door leakage path measurements.
Source: HAI and DCE 1996.
Appendix B
150 NIST NCSTAR 1-4D, WTC Investigation
Table B11. Door leakage areas.
Description/Use Type
Average
a
Measured
Leakage Area
(ft
2
)
Corrected
b
Leakage
Area (ft
2
)
Side hinged single leaf door
(stairs, tenant space, etc.) A 0.1688 0.20
Center parting double leaf doors
(tenant space, escalator cut-off, etc.) A 0.3376 0.081
Center parting single speed elevator doors
(freight elevators) B 0.4017 0.4017
Center parting double speed elevator doors
(freight elevators) C 0.7885 0.7885
Center parting single speed elevator doors
(local passenger elevators) B 0.5673 1.00
Center parting double speed elevator doors
(shuttle passenger elevators) C 1.2014 1.2014
Freight elevator door D 1.2500 1.25
Tenant doors A 0.3866 2.60
Elevator pit doors A 0.63 0.63
Mechanical room door A 1.00 1.00
Automatic closing dock door 1.00 1.00
Rolling steel dock doors 6 in.
2
/ft
2
6 in.
2
/ft
2
Rubber door at dock 1.00 1.00
a. Door leakage areas used in 2004 WTC CONTAM model.
b. Door leakage areas used in 1996 WTC CONTAM model.
B.1.6 1996 WTC CONTAM Model Calibration
The pressure measurements obtained during the site survey were used to validate the 1996 WTC
CONTAM model. In order to fine tune the model, the leakage path areas for stair doors, passenger
elevators, and selected other doors were increased. These increases lumped together as part of the stair
and elevator door openings other construction leakages that may have existed, such as penetrations in
elevator/stair shaft walls. Therefore, the leakage areas used may appear higher than that typically
associated with stair doors or elevator doors. The corrected leakage path areas used in the 1996 WTC
CONTAM model are listed in Table B11.
B.2 2004 WTC CONTAM MODEL
The construction process of the 2004 WTC CONTAM model was guided by efforts to provide a better
calibration of the model to the actual pressure measurements taken within the building. Details of this
calibration process are discussed in Appendix C. The following section outlines the final version(s) of the
2004 WTC CONTAM model, including aircraft impact damage models for WTC 1 and WTC 2.
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 151
B.2.1 Leakage Areas
The leakage areas used in the 2004 WTC CONTAM model for the various building construction elements
and their associated tightness classification are listed in Table B12, for model configurations A/B and C.
As is discussed in Appendix C, there were several potential building configurations, representing
combinations of assumed leakages and ventilation system parameters that provide relatively good
agreement with the calibration data.
Configuration A represents a best approximation of the various leakages and ventilation system
parameters that existed in the towers on September 11, 2001. Configuration B uses the same leakage areas
assumed for configuration A. The only difference between model configurations A and B is that
configuration A assumes 85 percent fan efficiency, while configuration B assumes 60 percent fan
efficiency. Tightness classifications are based on typical leakage area measurements of walls and floors in
commercial buildings (see Table B12) (Klote and Milke 2002).
Several architectural features were added to the 2004 WTC CONTAM model. Wall transfers grills, which
were observed during the 1996 site survey, were added to the model. A perimeter leakage was added,
intended to model floor slab penetrations for the high-pressure perimeter air system, discussed in
Chapter 4. Given the buildings approximate 800 ft perimeter, a perimeter leakage of 8 ft
2
and 16 ft
2
correlates to a 1/8 in. and 1/4 in. perimeter gap respectively.
Table B12. Leakage area parameters used in 2004 WTC model.
Model Configuration A/B Model Configuration C
Construction Elements Leakage Area
Tightness
a
Classification Leakage Area
Tightness
a
Classification
Floor 0.00095 in.
2
/ft
2
Tight 0.00095 in.
2
/ft
2
Tight
Exterior walls 0.02 in.
2
/ft
2
AverageTight 0.03 in.
2
/ft
2
Average
Stair walls
b
0.6 in.
2
/ft
2
Loose 1.0 in.
2
/ft
2
Very Loose
Shaft walls
b
1.5 in.
2
/ft
2
Average 3.1 in.
2
/ft
2
Loose
Wall transfer grills 6 in.
2
/ft
2
NA 10 ft
2
NA
HVAC louvers 0.3 in.
2
/ft
2
NA 0.3 in.
2
/ft
2
NA
Floor slab perimeter 8 ft
2
NA 16 ft
2
NA
Doorway
c
Average Measured Areas (see Table B11)
a. Tightness classifications are based on typical leakage area measurements of walls and floors in commercial building
(Klote and Milke 2002).
b. Listed wall leakages are per linear ft of wall based on a 12 ft slab-to-slab height.
c. Modeled doorway leakage areas were based on average measured leakage areas.
Key: NA, not applicable.
In the 1996 CONTAM model only one wall leakage type was used, termed interior wall. Variations in
wall leakage were taken into account by lumping extra leakage area into stair and elevator leakage paths.
This methodology presents problems in trying to keep track of the effect of the various leakage
parameters used.
The 2004 CONTAM model added several new leakage path types to provide unique leakage paths for the
various interior wall types (stairs, elevators, tenant/core separations). Leakage areas for doors (elevator,
stair, tenant) were reset to correspond to the actual measured values, as shown in Table B11.
Appendix B
152 NIST NCSTAR 1-4D, WTC Investigation
B.2.2 2004 CONTAM HVAC System
The basic framework of the HVAC modeling approach used in the 1996 WTC CONTAM model
(discussed in Sec. B.1.3) remained the same for the 2004 WTC CONTAM model. However, several
additions were made in the 2004 model, which had a sizeable impact on simulated pressures. The
additions to the 1996 HVAC modeling approach included:
Added schedules to the vertical shaft openings, which allowed the user to easily open and
close the shafts when necessary. When both supply and exhaust fans are turned off, the
ventilation shafts are modeled open in CONTAM. This is done to model the passive pathway
created when no air is being mechanically driven within the ventilation ducts. Conversely,
when both supply and exhaust fans are turned on, as in the case of summer normal mode,
ventilation shafts are modeled closed. This is done to model the mechanically driven air,
which eliminates the passive pathway for air transport via the ventilation ducts.
Dampers were added to core and tenant ventilation shafts. Damper dimensions were taken
from WTC mechanical riser diagrams.
Outdoor air louvers, connected to each ventilation unit/shaft, where added on each MER
level. Modeling the outdoor air louvers provided additional pressure relief to the building.
Outdoor air louver dimensions were taken from WTC mechanical riser diagrams.
In the 2004 WTC CONTAM model, configurations A and C assume that air is supplied and returned at
85 percent of listed capacities, while configuration B assumes that air is supplied and returned at
60 percent of listed capacities (see Tables B4a through B4e for details). It is important to note that in all
three configurations (AC) the ratio between supply and return air remains constant (as listed in
Tables B4a through B.4e). The ratio between supply and return air can have a substantial effect on
building pressures.
Zoned Smoke Control With Stair Pressurization
As discussed in Sec. 8.1.5 of this report, the state of the art smoke management approach referred to as
Zoned Smoke Control With Stair Pressurization was analyzed. This smoke management approach
consisted of pressurizing Stairs 13, while exhausting air within the tenant (interior) space on the fire
floor. Theoretical fans sizes for stair pressurization and tenant space exhaust were determined based on
providing 0.05 in. H
2
O across the stair and tenant doors under the following HVAC configuration:
Supply and exhaust fans are turned off in the multi-floor ventilation zone containing the fire.
For example, if a fire occurs on the 80th floor (served by MER 75), interior and core (supply
and exhaust) fans on floors 59 through 91 are turned off. As previously discussed, when both
supply and exhaust fans are turned off (within a ventilation zone), the ventilation shafts are
modeled open in CONTAM. This is done to model the passive pathway created when no air
is being mechanically driven within the ventilation ducts.
Supply and exhaust fans in ventilation zones not containing the fire are set to summer normal
mode. In summer normal mode, both supply and exhaust fans are running and therefore a
passive pathway for air movement is eliminated. Under summer normal mode, ventilation
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 153
shafts are modeled in CONTAM as being closed to address the mechanically driven air
within the ventilation ducts.
Class-II smoke dampers activate in the multi-floor ventilation zone containing the fire, as
smoke dampers would typically be provided in a state-of-the-art smoke management system.
Class-II leakage rated smoke dampers (UL 555S) were modeled in CONTAM by multiplying
the flow area of each ventilation shaft damper by a leakage area ratio of 0.6 in.
2
/ft
2
(Klote and
Milke 2002).
Some smoke control designers limit the height of single-injection stairways to 8 stories; however, other
designers feel this limit can be extended to 12 stories (Klote and Milke 2002). Given the height of the
WTC stairways (40 or more stories within each stair section), a multiple-injection point system was used.
A multi injection point stair pressurization system was modeled in CONTAM using the Simple Air
Handling System feature. Supply (injection) points were placed approximately five stories apart within
each stairway. CONTAM schedules were added to each supply point, which allowed the user to adjust
the amount of air supplied to each stair section. Stair sections were determined by the location of the
horizontal transfer corridors.
Stair 1 and Stair 2 consists of three sections:
Section 1: Floors 2 through 42
Section 2: Floors 43 through 75
Section 3: Floors 75 to Roof
Stair 3 consists of two sections:
Section 1: Floors 1 through 75
Section 2: Floors 75 through 107
Through an iterative simulation process using CONTAM, the stair pressurization fan size required to
achieve 0.05 in. H
2
O across the stair doors for each stair section was determined. Stair pressurization fan
sizes were determined for both model configurations A/B and C (as discussed in Sec. 8.3.3). Model
configurations A and B assumed identical leakage areas for the various building construction elements
and therefore the stair pressurization fan sizes required were the same. Modeling results for the stair
pressurization fan sizes required to achieve 0.05 in. H
2
O across stair doors are shown in Table B13. The
listed fan capacities in Table B13 assume that all stair doors, stair transition doors, elevator doors, and
tenant doors were closed.
Appendix B
154 NIST NCSTAR 1-4D, WTC Investigation
Table B13. Stair pressurization system fan sizes.
Stair 1 Stair 2 Stair 3
(cfm) (cfm) (cfm)
1* 19,200 17,600 21,000
2 10,500 9,600 7,800
3 13,300 15,400 -
Total 43,000 42,600 28,800
1* 22,400 20,800 27,000
2 14,100 12,600 10,200
3 17,500 19,600 -
Total 54,000 53,000 37,200
Model
Configuration
A/B
C
Stair
Sections
The listed fan capacities for Stairs 1 and 2 in stair section 1 (floors 2 through 42), do not achieve 0.05 in.
H
2
O across the stair doors on floors 2 through 6. Floors 2 through 6 are unique in that a stair vestibule is
located between the core space and stairway. As shown in Table B13 the fan size needed to pressurize
the stairway in model configuration C are considerably higher due to the relative looseness of the building
construction elements compared to model configuration A/B.
In conjunction with stair pressurization, an iterative simulation process was used in CONTAM to
determine the interior zone exhaust fan size required to achieve 0.05 in. H
2
O and 0.1 in. H
2
O across the
tenant doors for both model configurations A/B and C. Modeling results for interior exhaust fan sizes
required to achieve 0.05 in. H
2
O and 0.1 in. H
2
O across tenant doors are shown in Table B14.
Table B14. Calculated interior exhaust fan sizes.
0.05" H
2
0 0.1" H
2
0
(cfm) (cfm)
A/B 29,100 39,500
C 51,000 73,000
Tenant Doors
Model
Configuration
The listed fan capacities in Table B14 assume that all stair doors, stair transition doors, elevator doors,
and tenant doors are closed.
B.3 DAMAGED WTC CONTAM MODELS
A damaged CONTAM model, for both WTC 1 and WTC 2, was developed from the 2004 WTC
baseline CONTAM model (discussed in Sec. B.2). The damaged models reflect the aircraft impact
damage to a towers exterior walls, concrete floor slabs, core walls, and shafts (elevator and stairs) due to
the events on September 11, 2001.
For the damaged models, a northerly wind was modeled in CONTAM at a velocity of 11.2 mph (5 m/s)
and a direction of 328 degrees, depicting wind conditions on September 11, 2001. In the CONTAM
model, the wall azimuth angle (reference point) for the top wall on the SketchPad is defined as 0 degrees.
A 0 degree wind modeled in CONTAM would be directly perpendicular to the north face of WTC 1. As
shown in Fig. B2, a northerly wind is approximately 32 degrees off the perpendicular, equaling a wind
direction of 328 degrees. The outside temperature was modeled in CONTAM as 70 F, depicting
temperature conditions on September 11.
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 155
Source: Adapted from McAllister 2002.
Figure B2. Illustration of a northerly wind modeled in CONTAM at 328 degrees.
B.3.1 Damaged WTC 1 CONTAM Model
The top of the CONTAM SketchPad represents the north face of the building; where as the bottom of the
SketchPad represents the south face of the building. The impact damage to each exterior face of the
building, discussed below, was modeled accordingly. This was to done to ensure that the effects of a
northerly wind across the large exterior (damage) openings were modeled properly.
Table B15 shows the estimated impacted damage to the exterior walls of WTC 1. The impact damage
was added in CONTAM to the external walls as two-way flow openings. The estimated impacted damage
to exterior walls is based on the World Trade Center Building Performance Study performed by the
Federal Emergency Management Agency (FEMA) (McAllister 2002).
Appendix B
156 NIST NCSTAR 1-4D, WTC Investigation
Table B15. Estimated impact damage to exterior walls, WTC 1.
North South East West
(ft
2
) (ft
2
) (ft
2
) (ft
2
)
92 743 0 1,572 0
93 958 0 1,356 0
94 592 54 1,163 0
95 1,055 54 0 420
96 797 151 0 1,518
97 926 151 0 1,798
98 1,335 0 0 0
Floor
Table B16 shows the estimated impact damage to the concrete floor slabs and core walls of WTC 1. The
impact damage was added in CONTAM to the floors in the tenant area and to the core walls as a two-way
flow opening. The estimated impact damage is based on the National Institute of Standards and
Technology (NIST) preliminary impact estimates shown in Appendix A. The floor slab damage is based
on the solid lines up to the core area in Figs. A1 through A6. The core wall damage is based on the
dotted lines across the core area walls in Figs. A1 through A6.
Table B16. Estimated impact damage to floor
slabs and core walls, WTC 1.
(ft
2
) (ft
2
)
94 2,800 216
95 1,650 432
96 2,830 588
97 870 0
98 330 0
Floor Slab
Damage
Core Wall
Damage Floor
Table B17 shows the estimated impact damage to the stair and elevator shafts in WTC 1. The impact
damage was added in CONTAM to the associated shaft wall as a two-way flow opening. The estimated
impact damage is based on the NIST impact estimates shown in Appendix A. The stair and elevator
damage is based on the dotted lines in Figs. 71 and 72.
Table B17. Estimated impact damage to stair/elevator shafts, WTC 1.
1 2 3 50 6-7 90-92 93-95
(ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
)
94 192 0 0 0 0 0 0
95 480 120 120 468 0 0 0
96 288 576 576 708 672 864 864
Local Elevators Express & Freight Elevators Stairs
Floor
B.3.2 Damaged WTC 2 CONTAM Model
Similar to the damaged WTC 1 model, the damaged WTC 2 model considers the north face of the
building as the top of the SketchPad. The impact damage to each exterior face of the building, discussed
below, was modeled accordingly. This was to done to ensure the effects of a northerly wind across the
large exterior (damage) openings were modeled properly for the different orientations of WTC 1 and
WTC 2. It should be noted, however, that the WTC 2 model does not consider the impact of WTC 1 on
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 157
wind driven flow around the building, which is important since WTC 1 was directly in the path of the
prevailing wind on September 11, 2001.
The impact damage to the core walls of WTC 2 were rotated by 90 degrees with respect to WTC 1 (see
Fig. B3). This was done to correctly account for the different core orientation in WTC 2 versus WTC 1
with respect to the point of aircraft impact. Since the CONTAM model considers the interior office space
as one large open zone, the core can be effectively rotated by placing damage openings in the proper
locations, rather than physically rotating the location of the core zones on the CONTAM SketchPad.
Northerly Wind (328
o
)
North (exterior wall)
South (exterior wall)
E
a
s
t
(
e
x
t
e
r
i
o
r
w
a
l
l
)
W
e
s
t
(
e
x
t
e
r
i
o
r
w
a
l
l
)
N
o
r
t
h
(
c
o
r
e
w
a
l
l
)
West (core wall)
S
o
u
t
h
(
c
o
r
e
w
a
l
l
)
East (core wall)
Core Area
Tenant Area
Damaged WTC 1 CONTAMW Model Damaged WTC 2 CONTAMW Model
Northerly Wind (328
o
)
North (exterior wall)
South (exterior wall)
E
a
s
t
(
e
x
t
e
r
i
o
r
w
a
l
l
)
W
e
s
t
(
e
x
t
e
r
i
o
r
w
a
l
l
)
E
a
s
t
(
c
o
r
e
w
a
l
l
)
North (core wall)
W
e
s
t
(
c
o
r
e
w
a
l
l
)
South (core wall)
Core Area
Tenant Area
Figure B3. CONTAM SketchPad, floor 15, exterior and core wall directional orientations
for damaged CONTAM models.
Table B18 shows the estimated impacted damage to the exterior wall of WTC 2. The impact damage
was added in CONTAM to the external walls as two-way flow openings. The estimated impacted damage
to the south exterior wall of WTC 2 is based on Fig. 227 in the FEMA report (McAllister 2002).
Impacted damage to the north and east exterior walls of WTC 2 are based on video data showing post-
impact window damage (Custer 2002).
Table B18. Estimated impact damage to exterior walls, WTC 2.
North South East West
(ft
2
) (ft
2
) (ft
2
) (ft
2
)
78 0 360 0 0
79 0 550 0 0
80 24 472 690 0
81 119 824 690 0
82 143 415 690 0
83 95 241 0 0
84 0 97 0 0
Floor
Appendix B
158 NIST NCSTAR 1-4D, WTC Investigation
Table B19 shows the estimated impact damage to the concrete floor slabs and core walls of WTC 2. The
impact damage was added in CONTAM to the floors in the tenant zone and to the core walls as a two-
way flow opening. The estimated impact damage is based on NIST preliminary impact estimates shown
in Appendix A. The floor slab damage is based on the solid line up to the core are in Figs. A7 through
A12. The core wall damage is based on the dotted lines across the core area walls in Figs. A7 through
A12.
Table B19. Estimated impact damage, floor slabs
and core walls, WTC 2.
(ft
2
) (ft
2
)
78 830 0
79 2,210 240
80 650 265
81 910 265
82 1,685 265
83 600 0
84 240 0
Floor
Floor Slab
Damage
Core Wall
Damage
Table B20 shows the estimated impact damage to the stair and elevator shafts in WTC 2. The impact
damage was added in CONTAM to the associated shaft walls as a two-way flow opening. The estimated
impact damage is based on the NIST impacts estimates shown in Appendix A. The stair and elevator
damage is based on the dotted lines in Figs. 71 and 72.
Table B20. Estimated impact damage, stair and elevators, WTC 2.
Express Elevator
2 3 6-7 81-83 84-86 87-89 90-92 93-95
(ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
) (ft
2
)
78 288 0 0 0 0 0 0 0
79 120 120 300 432 432 432 264 84
80 0 0 468 144 120 84 0 0
81 0 0 468 144 120 84 0 0
Stairs Local Elevators
Floor
CONTAM Model Construction
NIST NCSTAR 1-4D, WTC Investigation 159
B.3.3 REFERENCES
Custer, R. L. P. 2002. World Trade Center Expert Report: Documentation and analysis of fire spread in
World Trade Center events of September 11, 2001. Arup Fire, New York, NY, August.
HAI (Hughes Associates, Inc.) and DCE (Dillon Consulting Engineers). 1996. Smoke Management
Evaluation Study of The World Trade Center Complex Excluding the Concourse and Plaza. Prepared
for the PANYNJ (Port Authority of New York and New J ersey). New York, NY.
Klote, J . H., and J . A. Milke. 1992. Design of Smoke Management Systems. American Society of Heating,
Refrigerating and Air-Conditioning Engineers (ASHRAE). Atlanta, GA.
Klote, J . H., and J . A. Milke. 2002. Principles of Smoke Management. American Society of Heating
Refrigeration and Air-Conditioning Engineers (ASHRAE). Atlanta, GA.
LEA, F. C. 1942. Hydraulics for Engineers and Engineering Students: Sixth edition. London, UK.
McAllister, T., ed. 2002. World Trade Center Building Performance Study: Data Collection, Preliminary
Observations, and Recommendations. FEMA 403. Federal Emergency Management Agency.
Washington, DC, May.
UL (Underwriters Laboratories). 2003b. Standard for Safety Leakage Rated Dampers for Use in Smoke
Control Systems. UL Standard 555S, Northbrook, IL.
Walton, G. N. (1997). CONTAM96 User Manual, NISIR 6065, National Institute of Standards and
Technology.
Appendix B
160 NIST NCSTAR 1-4D, WTC Investigation
This page intentionally left blank.
NIST NCSTAR 1-4D, WTC Investigation 161
Appendix C
CONTAM MODEL CALIBRATION
As part of the 1996 Hughes Associates Incorporated (HAI)/Dillon Consulting Engineers (DCE) study,
measurements were made of actual airflows and pressures at various locations in World Trade Center
(WTC) 1 for different ventilation modes. At the time, a limited model calibration was performed in order
to perform the analysis in the 1996 study.
To complete the modeling summarized in this report, many changes were made to the model (summarized
in Appendix B) to add building features not present in the 1996 model and to add additional leakage path
types to provide additional degrees of freedom in performing the simulations. As such, it was desired to
revisit the calibration of the model in order to attempt to provide better agreement between the model and
the measured data.
This appendix discusses the efforts that went into the recalibration of the model, and the relative
agreement between the actual WTC pressure measurements (see Tables B2 through B5) and pressures
simulated using the 2004 WTC CONTAM model. Calibration efforts were conducted for the four
different ventilation modes (Modes 1 through 4), discussed in Appendix B.1.4.
C.1 PRE-CALIBRATION EFFORTS
The 2004 WTC CONTAM model was constructed using the 1996 WTC CONTAM model as a starting
point. Changes were made to the model to add building features that were identified as missing from the
1996 model and to reflect lessons learned over the intervening years with regard to the use of CONTAM
to model the performance of smoke management systems in high-rise buildings. This knowledge was
applied to the calibration process and construction of the 2004 WTC model.
The first step in the 2004 WTC CONTAM model construction involved setting all of the corrected door
leakage areas used in 1996 model back to the values measured during the 1996 WTC site visit. Measured
door leakage areas are given in Appendix B, Table B11. Since the door leakage values (core, elevator,
and stair doors) reflect actual measured values, these leakages remained fixed during the calibration
process.
The leakage areas of several key construction elements were not included as individual elements during
the 1996 study and therefore needed to be estimated. The non-measured construction elements included:
Exterior walls
Stair walls
Elevator walls
Floors
Appendix C
162 NIST NCSTAR 1-4D, WTC Investigation
Leakage area estimates for these construction elements were kept within the bounds of tightness typically
observed in commercial buildings (Klote and Milke 2002). Typical leakage areas of walls and floors of
commercial buildings are listed in Table C1.
Table C1. Typical leakage areas of walls and floors in commercial buildings.
Construction Element Tightness Leakage Area Ratio
a
Tight 0.01 in.
2
/ft
2
Average 0.03 in.
2
/ft
2
Loose 0.06 in.
2
/ft
2
Exterior Building Walls
(included construction cracks and cracks
around windows and doors)
Very Loose 0.19 in.
2
/ft
2
Tight 0.02 in.
2
/ft
Average 0.19 in.
2
/ft
Stairway Walls
(includes construction cracks but not cracks
around windows or doors)
Loose 0.61 in.
2
/ft
Tight 0.31 in.
2
/ft
Average 1.45 in.
2
/ft
Elevator Shaft Walls
(includes construction cracks but not cracks
around doors
Loose 3.11 in.
2
/ft
Tight 0.00095 in.
2
/ft
2
Average 0.00749 in.
2
/ft
2
Floors
(includes construction cracks and gaps around
penetrations)
Loose 0.02448 in.
2
/ft
2
a.
For the stairway and elevator shaft walls, the leakage area ratio is the area of the leakage through the
wall per linear ft based on a 12 ft slab-to-slab height.
Measurements made by Tamura and Shaw (Klote and Milke 2002) show stairway walls to be
approximately one order of magnitude tighter than elevator shaft walls (see Table C1 for details). In the
1996 WTC CONTAM model, elevator shaft walls and stairway walls were combined together as one
building construction element for simplicity. In the 2004 WTC CONTAM model, elevator shafts and
stairway walls were modeled as two separate construction elements (airflow elements), to account for
their considerable difference in tightness. This allowed the modeler to independently vary the leakage
areas for the elevator shaft walls and stairway walls.
C.2 CALIBRATION EFFORTS
Initial calibration efforts involved matching simulated data with measured data under Mode 2, where all
fans are set to summer normal. The heating, ventilating, and air conditioning (HVAC) system within
CONTAM was set up to depict ventilation Mode 2, where all fans were running at 85 percent of design
capacity. Using an iterative process, leakage areas for the non-measured construction elements were
varied until simulated results were in reasonable agreement with measured values. Under Mode 2,
simulated results showed the best agreement with measured values when:
Exterior walls were set to tight
Elevator shaft wall were set to tight
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 163
Stairway walls were set to average
Floors were set to tight
This tightness configuration for the non-measured construction elements also produced a reasonable
agreement under Mode 1.
The iterative process also involved varying outside temperature and wind velocities. CONTAM
simulations, running under Mode 2, showed that outside temperature differences of up to 10 F and wind
velocity variation of up to 20 mph produced relatively small changes in simulated pressures, as illustrated
in Figs. C21 through C30. Mode 1 was calibrated using an outside temperature of 63 F and Modes 2
through 4 were calibrated using an outside temperature of 70 F. The temperatures used were based on
measurements made during the site survey (see Sec. B.1.4) and the sequence in which the four different
HVAC modes were conducted, as outside temperature would have risen during the day from the early
morning to late morning period in which the pressure tests were performed.
The next calibration step involved testing the tightness configuration that produced a reasonable
agreement under Modes 1 and 2, for Modes 3 and 4. When the CONTAM model was set up for Mode 3
(core pressurization for the full height of the tower), using the tightness configuration used to match
Modes 1 and 2, the model simulations produced very high positive pressures across the tenant doors.
Simulated pressures across the tenant door were greater than 1.0 in. H
2
O. These high pressures across the
tenant doors would have been apparent during the site survey, because the doors would have been pinned
shut due to the substantial overpressure. In actuality, under the core pressurization configuration (Modes 3
and 4), positive pressures were measured between 0.01 in. H
2
O and 0.05 in. H
2
O (see Tables B5 through
B8).
At this point in the calibration process the results showed that the 2004 CONTAM model agreed
reasonably well with Modes 1 and 2, but not with Modes 3 and 4. Because this study involved analyzing
the effectiveness of various smoke management approaches, including the core pressurization mode,
Modes 3 and 4 were used to calibrate the final version(s) of the 2004 WTC CONTAM model. Principally,
the positive pressures measured across the tenant doors and elevator 50 during core pressurization
(Modes 3 and 4) were used to calibrate the final version(s) of the 2004 WTC CONTAM model.
C.2.1 Architectural and Mechanical Building Features Added
The initial simulation results for Modes 3 and 4 warranted further investigation into the buildings
architecture and HVAC system. Investigation of the buildings ventilation system led to modeling shaft
dampers and mechanical equipment room (MER) outdoor air louvers located in WTC 1. This
modification dramatically reduced pressures across the tenant doors in Modes 3 and 4, by adding
additional pressure relief from the core and interior space to the exterior of the building. Investigation of
the buildings architecture lead to modeling the wall transfer grills located over many of the tenant doors,
providing an additional airflow path between the core and interior space. Leakage around the perimeter of
the tenant space due to ventilation penetrations for the high-pressure air was also modeled in CONTAM.
As stated in Appendix B, perimeter gaps of 1/8 in. and 1/4 in. correlate to total perimeter leakages of 8 ft
2
and 16 ft
2
, respectively, given the approximate 800 ft perimeter of the building.
Appendix C
164 NIST NCSTAR 1-4D, WTC Investigation
After the architectural and mechanical building features described above were added to the 2004 WTC
model, leakage values for the non-measured construction elements were set based on the construction
characteristics of each element. The construction characteristics of each element are as follows:
FloorsFloor construction typically consisted of 4 in. of lightweight concrete on 1 1/2 in.,
22 gauge non-composite steel deck (McAllister 2002). In the core area, slab thickness was
5 in. Based on this construction, the floor leakage was set and held at tight
(0.00095 in.
2
/ft
2
).
Exterior WallsConstruction of the buildings exterior walls consisted of closely spaced
built-up box columns. A total of 59 perimeter columns were present along each flat face of
the building. The columns were built up by welding four plates together to form an
approximately 14 in. spare section, spaced at 3-4 in. on center. This column spacing resulted
in relatively small window sizes, approximately 2 ft wide. Based on this construction, leakage
areas for the exterior wall were set between average (0.03 in.
2
/ft
2
) and tight (0.01 in.
2
/ft
2
).
Elevator Shaft WallsElevator shaft walls were constructed of taped drywall. Using drywall
greatly reduced the natural weight load of the building. Based on this construction, leakage
areas were set between average (1.45 in.
2
/ft) and loose (3.11 in.
2
/ft).
Stairway WallsStairway walls were constructed of taped drywall. Using drywall greatly
reduced the natural weigh load of the building. Stairway walls contained various types of
penetrations, adding leakage. Based on this construction, leakage areas were set between
loose (0.6 in.
2
/ft) and very loose (1.0 in.
2
/ft).
C.2.2 Model Configurations
Three different model configurations were used in order to represent the uncertainty involved with this
type of modeling. These three model configurations attempt to bracket the estimated leakage rates for the
various building construction elements, as well as the efficiency of the HVAC system. The model
configurations include:
Configuration AThis model configuration is based on best estimates of building variables.
The HVAC efficiency is set at 85 percent based on WTC building engineer estimates (HAI
and DCE 1996). Door leakage areas are set to measured leakage values (Appendix B,
Table B11). Non-measured building construction elements are set to expected leakage
values based on the construction type of each element, as discussed in Sec. C.2.1. A complete
summary of leakage areas for this configuration is given in Appendix B,
Table B12.
Configuration BSame leakage areas as configuration A. This model configuration reduces
the HVAC efficiency from 85 to 60 percent, and results in an overall reduction of pressure
within the building.
Configuration CThis model configuration increases the leakage areas (i.e., loosens) of
configuration A. The HVAC efficiency is held at 85 percent. A complete summary of leakage
areas for this configuration is given in Appendix B, Table B12.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 165
C.3 CALIBRATION RESULTS
A summary of the simulated pressures verses the measured pressures, under HVAC Modes 1 through 4,
are shown in Table C1 (Sec. C.3.1). The simulation results for model configurations AC under HVAC
Modes 1 through 4 are shown in Tables C2 through C12 (Sec. C.3.2). Graphical comparisons between
simulated and measured pressures, for HVAC Modes 1 through 4, are shown in Figs. C1 through C20
(Sec. C.3.3). Temperature effects are shown in Figs. C21 through C25 (Sec. C.3.4). Wind effects are
shown in Figs. C26 through C30 (Sec. C.3.5).
C.3.1 Summary of Results
The values listed in Table C1 are the minimum and maximum pressures under HVAC Modes 1 through
4, observed versus predicted, for the full height of the building, across the tenant doors or within a given
group of shafts. Shaft types included the following:
Stairs: Stair 1, 2, and 3
Elevator 50: Elevator 50 only
Shuttle and freight elevators: Elevator 5, 6, 17, 49, 12, 34, 67, 89, 1011, 1213, 1415,
1617, 1819, 2021, and 2223
Local elevators: Elevator 2426, 3032, 3638, 3941, 5153, 5456, 5759, 6062, 6365,
6668, 6971, 7274, 7577, 7880, 8183, 8486, 8789, 9092, and 9395
Table C1. Comparison of measured and simulated pressures for
HVAC Modes 1 through 4.
Min Max Min Max Min Max Min Max Min Max
Measured -0.04 0.01 -0.01 0.01 -0.05 0.04 -0.02 0.02 -0.02 0.01
Configuration A/B -0.04 0.02 -0.02 0.05 -0.04 0.05 -0.04 0.03 0.00 0.00
Configuration C -0.03 0.02 -0.01 0.03 -0.03 0.04 -0.03 0.02 0.00 0.00
Measured -0.25 0.02 -0.24 0.02 -0.07 0.01 -0.20 0.12 -0.07 0.06
Configuration A -0.01 0.01 -0.01 0.00 -0.01 0.00 0.00 0.01 -0.01 0.00
Configuration B -0.01 0.01 -0.01 0.00 0.00 0.00 0.00 0.01 0.00 0.00
Configuration C -0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Measured -0.03 0.10 -0.02 0.05 0.05 0.07 -0.01 0.04 0.00 0.03
Configuration A 0.00 0.05 -0.01 0.05 0.01 0.02 -0.01 0.04 0.01 0.02
Configuration B 0.00 0.03 0.00 0.02 0.00 0.01 0.00 0.01 0.00 0.01
Configuration C 0.00 0.01 0.00 0.02 0.00 0.01 0.00 0.01 0.00 0.01
Measured -0.26 0.08 -0.32 0.08 -0.25 0.02 -0.19 0.11 -0.08 0.05
Configuration A -0.05 0.14 -0.12 0.15 -0.02 0.11 -0.07 0.15 -0.01 0.05
Configuration B -0.04 0.08 -0.07 0.08 -0.02 0.05 -0.04 0.08 0.00 0.03
Configuration C -0.02 0.07 -0.05 0.07 0.00 0.05 -0.04 0.08 0.00 0.02
Tenant Doors
Mode 2
Mode 3
Mode 4
Local Elev. HVAC
Configuration
Data Description
Mode 1
Elevator 50 Stairs Shuttle Elev.
Appendix C
166 NIST NCSTAR 1-4D, WTC Investigation
For Mode 1, the 2004 WTC CONTAM model agreed fairly well with measured pressures. Maximum
stairs pressure agreed by a factor of 1.7, elevator 50 by 2.5, shuttle and freight elevators by 1.2, local
elevators by 1.4, and tenant doors by 2.5. Tenant doors agreement factors are high because simulated
pressures were very close to zero.
For Mode 2, the 2004 WTC CONTAM model agreed poorly with measured pressures. Simulated
agreement factors were between 3.1 and 85.7. These agreement factors are high because simulated
pressures were close to zero. Simulated pressure in Mode 2 is close to zero for two main reasons:
1. Non-measured building construction elements were loosened to better match Modes 3 and 4.
2. The ventilation system modeled assumed a relative air balance between the supply and return
air.
Supply and return air rates were modeled based on mechanical drawings, as shown in Tables B4a
through B4e, and not on actual measured ventilation rates. Tenant changes on individual floors to the
ventilation system were not known nor modeled. As shown Tables B6 and B8, relatively high negative
pressures were measured on floors 29 and 92 under Modes 2 and 4. Mode 4 is similar to Mode 2 on
floors 29 and 92 in that the ventilation is set to summer normal. The pressures measured on floor 29 and
92 are over an order of magnitude greater than the other measured floors. These abnormally high
pressures may be attributed to tenant changes to the ventilation system on these floors (or neighboring
floors). Tenant changes to the ventilation system on an individual floor would not only alter the pressures
on that floor, but could alter the entire network of pressures within the WTC complex. These abnormally
high pressures may also be attributed to unique architectures on these floor or neighboring floors, as
discussed in Sec. C.2.3. It should be noted that in many cases, the floors directly above/below floors 29
and 92 had fairly good agreement with measured data for these floors, indicating that pressures were in
the correct range for that area of the building.
For Mode 3, the 2004 WTC CONTAM model agreed reasonably well with measured pressures.
Maximum stair pressures agreed by a factor of 2.1, elevator 50 by 1.0, shuttle and freight elevators by 3.2,
local elevators by 1.0, and tenant doors by 1.7. As previously discussed, calibration of the 2004 WTC
CONTAM model depended heavily on the positive pressures measured during the site survey across the
tenant doors and elevator 50, under Modes 3 and 4 (core pressurization). These principal measured
pressures are highlighted in Table C1. As shown in Table C1, simulated pressures agree very well with
the principal measured pressures across the tenant doors and elevator 50.
For Mode 4, the 2004 WTC CONTAM model agreed reasonably well with measured pressures.
Maximum stair pressures agreed by a factor of 1.1, elevator 50 by 1.0, shuttle and freight elevators by 2.5,
local elevators by 1.3, and tenant doors by 1.1.
C.3.2 Simulated Results for Modes 1 through 4
The simulation results for model configurations A through C under HVAC Modes 1 through 4 are shown
in Tables C2 through C12. The measured data is summarized in Tables B5 through B8 (Appendix B)
as well as the calibration measurements taken during the 1996 HAI/DCE study. The figures discussed in
Sec. C.3.3 provide a comparison between the two sets of data.
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
6
7
C
O
N
T
A
M
M
o
d
e
l
C
a
l
i
b
r
a
t
i
o
n
Table C2. CONTAM simulated pressure differentials using configuration AMode 1 (all fans off).
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23
1 -0.023 -0.018
29 0.013 0.010 0.000 -0.017 0.000
44 -0.001 0.001 0.000 0.010 -0.042 -0.002 -0.010 -0.014 -0.002
56 -0.023 -0.027 0.000 -0.014 0.012 0.012 -0.042 0.004 0.013 -0.002
60 0.004 0.000 0.022 0.014 0.003
78 0.000 -0.001 -0.023 0.023 0.050 0.037 0.039 0.040 0.041 0.049 -0.001
83 -0.036 -0.003 -0.021 0.022 0.000
92 -0.003 0.019 0.010 0.040 0.003
104 0.004 0.020 0.022 0.046 0.004
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 -0.007 -0.008 -0.023
29 0.008 0.004
44 0.000 0.000 -0.006 -0.007 -0.022 -0.026 -0.029 -0.026
56 -0.004 -0.005
60 0.023 0.023
78 -0.001 -0.002 -0.008 -0.013 -0.030 -0.025 -0.042 -0.042
83 0.000 -0.001
92 0.029 0.024
104 0.008 0.007
Tenant
Doors
Floor
Stair Shuttle & Freight Elevators
III-C Floor
Local Elevator Banks
II-D III-A I-C II-A II-B III-B III-D II-C
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate air flowing out of a space.
1
6
8
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
A
p
p
e
n
d
i
x
C
Table C3. CONTAM simulated pressure differentials using configuration CMode 1 (all fans off).
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23
1 -0.013 -0.010 -0.016
29 0.006 0.004 0.000 -0.011 0.000
44 0.000 0.000 -0.001 0.008 -0.029 -0.003 -0.007 -0.010 0.008 0.008 -0.029 0.002 0.009 -0.001
56 -0.013 -0.017 -0.003 -0.010 -0.001
60 0.003 -0.011 0.013 0.006 0.002
78 0.000 -0.001 -0.018 0.001 0.013 0.038 0.030 0.032 0.034 0.032 0.038 -0.001
83 -0.025 -0.002 -0.015 0.002 0.013 0.000
92 -0.003 0.010 0.007 0.014 0.023 0.002
104 0.003 0.010 0.017 0.020 0.027 0.003
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 -0.005 -0.005 -0.012
29 0.005 0.003
44 0.000 0.000 -0.003 -0.004 -0.011 -0.015 -0.016 -0.015
56 -0.002 -0.002
60 0.016 0.015
78 -0.002 -0.004 -0.007 -0.012 -0.020 -0.013 -0.029 -0.030
83 0.001 -0.001
92 0.021 0.016
104 0.011 0.010
Floor
Stair Shuttle & Freight Elevators
III-B Floor
Local Elevator Banks
I-C II-A
Tenant
Doors
II-B III-D II-C II-D III-A III-C
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate air flowing out of a space.
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
6
9
C
O
N
T
A
M
M
o
d
e
l
C
a
l
i
b
r
a
t
i
o
n
Table C4. CONTAM simulated pressure differentials using configuration AMode 2 (all fans set to summer normal).
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23
1 0.001 0.002 0.002
29 -0.003 -0.003 -0.002 -0.002 -0.001
44 -0.005 -0.002 -0.004 -0.008 -0.004 -0.005 -0.003 -0.005 -0.005 -0.006 -0.004 -0.005 -0.005 -0.006
56 -0.003 -0.004 -0.005 -0.003 -0.003 -0.006 -0.001
60 -0.002 -0.004 -0.005 -0.003 -0.003 -0.006 -0.001
78 0.000 -0.002 0.006 -0.002 -0.005 -0.003 -0.004 -0.005 -0.005 -0.004 -0.004 -0.004 -0.005
83 -0.009 -0.002 -0.001 -0.006 -0.007 -0.001
92 -0.012 -0.005 -0.003 -0.008 -0.008 -0.002
104 -0.010 -0.002 0.004 -0.004 -0.004 -0.002
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 -0.001 0.001 0.001
29 -0.002 -0.003
44 0.001 0.001 0.001 0.000 0.000 -0.001 -0.004 -0.005
56 0.000
60 0.000 0.000
78 0.006 0.006 0.009 0.009 0.009 0.007 -0.002 -0.002
83 -0.001 -0.002
92 -0.001 -0.001
104 -0.004 -0.004
Tenant
Doors
III-C III-D III-B
Floor
Stair Shuttle & Freight Elevators
Floor
Local Elevator Banks
I-C II-A II-B II-C II-D III-A
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate air flowing out of a space.
1
7
0
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
A
p
p
e
n
d
i
x
C
Table C5. CONTAM simulated pressure differentials using configuration BMode 2 (all fans set to summer normal).
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23
1 -0.001 0.000 0.000
29 -0.002 -0.002 -0.002 -0.002 -0.001
44 -0.005 -0.001 -0.002 -0.005 -0.002 -0.003 -0.002 -0.003 -0.002 -0.002 -0.002 -0.002 -0.002 -0.002
56 -0.002 -0.004 -0.003 -0.002 -0.003 -0.004 -0.001
60 -0.002 -0.004 -0.003 -0.002 -0.003 -0.004 -0.001
78 0.000 -0.002 0.005 0.000 -0.001 -0.001 0.000 -0.001 -0.001 0.000 0.000 0.000 -0.004
83 -0.008 -0.001 -0.001 -0.004 -0.005 -0.001
92 -0.010 -0.003 -0.002 -0.005 -0.005 -0.001
104 -0.008 -0.002 0.003 -0.001 -0.002 -0.001
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 -0.001 -0.001 0.000
29 -0.001 -0.002
44 0.001 0.001 0.001 0.001 0.001 0.000 -0.002 0.000
56 0.000
60 0.000 0.000
78 0.005 0.005 0.007 0.006 0.007 0.007 -0.002 -0.002
83 -0.001 -0.001
92 0.000 -0.001
104 -0.004 -0.004
III-B Floor
Local Elevator Banks
I-C II-A II-B II-C II-D III-A III-C III-D
Floor
Stair Shuttle & Freight Elevators Tenant
Doors
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space. Positive values indicate air flowing out of a space.
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
7
1
C
O
N
T
A
M
M
o
d
e
l
C
a
l
i
b
r
a
t
i
o
n
Table C6. CONTAM simulated pressure differentials using configuration CMode 2 (all fans set to summer normal).
1 2 3 5 6 17 49 50 1-2 3-4 6-7 8-9 10-11 12-13 14-15 16-17 18-19 20-21 22-23
1 0.000 0.000 0.000
29 0.001 0.001 -0.001 -0.001 0.000
44 -0.003 -0.001 -0.002 -0.002 -0.002 -0.001 -0.001 -0.002 -0.002 -0.002 -0.002 -0.002 -0.002 -0.001
56 -0.001 -0.002 -0.002 -0.001 -0.001 -0.002 0.000
60 -0.001 -0.002 -0.002 -0.001 -0.001 -0.002 0.000
78 0.000 -0.001 0.003 -0.001 -0.001 -0.001 -0.001 -0.001 -0.001 -0.001 -0.001 -0.001 -0.002
83 -0.004 -0.001 -0.001 -0.002 -0.002 -0.001
92 -0.005 -0.002 -0.001 -0.003 -0.003 -0.001
104 -0.006 -0.003 0.002 -0.001 -0.001 0.000
I-A I-B
24-26 30-32 36-38 39-41 51-53 54-56 57-59 60-62 63-65 66-68 69-71 72-74 75-77 78-80 81-83 84-86 87-89 90-92 93-95 96-98
1 -0.001 0.000 0.000
29 -0.001 -0.001
44 0.000 0.000 0.000 0.000 0.000 -0.001 -0.001 -0.002
56 0.000
60 0.000 0.000
78 0.003 0.003 0.004 0.004 0.004 0.003 0.000 0.001
83 -0.001 -0.001
92 0.000 0.000
104 -0.001 -0.001
III-B Floor
Local Elevator Banks
I-C II-A II-B II-C II-D III-A III-C III-D
Floor
Stair Shuttle & Freight Elevators Tenant
Doors
Note: Listed differential pressure values have units of in. H
2
O. Negative (-)values indicate air flowing into a space. Positive values indicate air flowing out of a space.
1
7
2
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
A
p
p
e
n
d
i
x
C
Table C7. CONTAM simulated pressure differentials using configuration AMode 3
(all fans set to provide core pressurization).
1 2 3 17 49 50
1
29 0.048 0.039 0.005 -0.014 0.018
44
56
60 0.000 -0.004 0.019 0.022 0.012 0.020 0.009
78
83
92 0.006 0.002 0.004 0.050 0.013
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.040 0.032
44
56
60 0.012 0.012
78
83
92 -0.008 -0.004
104
Floor
Stair Shuttle & Freight Elevators Tenant
Doors
Floor
Local elevator Banks
I-C II-B III-C
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing
into a space. Positive values indicate air flowing out of a space.
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
7
3
C
O
N
T
A
M
M
o
d
e
l
C
a
l
i
b
r
a
t
i
o
n
Table C8. CONTAM simulated pressure differentials using configuration BMode 3
(all fans set to provide core pressurization).
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space.
Positive values indicate air flowing out of a space.
1 2 3 17 49 50
1
29 0.025 0.021 0.001 -0.011 0.009
44
56
60 -0.001 -0.004 0.008 0.009 0.004 0.007 0.004
78
83
92 0.001 0.000 0.001 0.025 0.006
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.021 0.016
44
56
60 0.006 0.006
78
83
92 -0.004 -0.002
104
Floor
Local Elevator Banks
I-C II-B III-C
Floor
Stair Shuttle & Freight Elevators Tenant
Doors
1
7
4
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
A
p
p
e
n
d
i
x
C
Table C9. CONTAM simulated pressure differentials using configuration CMode 3
(all fans set to provide core pressurization).
1 2 3 17 49 50
1
29 0.013 0.010 0.002 -0.005 0.006
44
56
60 -0.001 -0.003 0.004 0.009 0.001 -0.004 0.003
78
83
92 0.003 0.002 0.002 0.018 0.006
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.014 0.013
44
56
60 0.003 0.003
78
83
92 -0.002 0.000
104
Floor
Local elevator Banks
I-C II-B III-C
Floor
Stair Shuttle & Freight Elevators Tenant
Doors
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space.
Positive values indicate air flowing out of a space.
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
7
5
C
O
N
T
A
M
M
o
d
e
l
C
a
l
i
b
r
a
t
i
o
n
Table C10. CONTAM simulated pressure differentials using configuration AMode 4
(floors 59 through 91 set to provide core pressurization, all other floors set on 100 % outside air).
1 2 3 5 6 17 49 50
1
29 0.003 0.001 -0.049 -0.001 -0.121 -0.007
44
56
60 0.035 0.027 0.139 0.105 0.091 0.146 0.050
78
83
92 -0.004 -0.035 -0.011 -0.024 0.096 -0.003
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.001 -0.002
44
56
60 0.145 0.142
78
83
92 -0.071 -0.048
104
Floor I-C
Shuttle & Freight Elevators
Local elevator Banks
II-B III-C
Floor
Stair Tenant
Doors
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space.
Positive values indicate air flowing out of a space.
1
7
6
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
A
p
p
e
n
d
i
x
C
Table C11. CONTAM simulated pressure differentials using configuration BMode 4
(floors 59 through 91 set to provide core pressurization, all other floors set on 100 % outside air).
1 2 3 5 6 17 49 50
1
29 0.006 0.002 -0.036 -0.001 -0.075 -0.004
44
56
60 0.019 0.013 0.083 0.054 0.047 0.075 0.027
78
83
92 -0.017 -0.033 -0.012 -0.023 0.044 -0.004
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.001 -0.002
44
56
60 0.083 0.082
78
83
92 -0.043 -0.030
104
Floor
Stair Tenant
Doors
Floor I-C
Local elevator Banks
II-B III-C
Shuttle & Freight Elevators
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space.
Positive values indicate air flowing out of a space.
N
I
S
T
N
C
S
T
A
R
1
-
4
D
,
W
T
C
I
n
v
e
s
t
i
g
a
t
i
o
n
1
7
7
C
O
N
T
A
M
M
o
d
e
l
C
a
l
i
b
r
a
t
i
o
n
Table C12. CONTAM simulated pressure differentials using configuration CMode 4
(floors 59 through 91 set to provide core pressurization, all other floors set on 100 % outside air).
1 2 3 5 6 17 49 50
1
29 0.003 0.002 -0.020 -0.001 -0.051 -0.003
44
56
60 0.021 0.021 0.070 0.049 0.036 0.068 0.025
78
83
92 -0.005 -0.022 -0.008 -0.004 0.037 -0.002
104
36-38 39-41 57-59 60-62 81-83 84-86
1
29 0.001 -0.002
44
56
60 0.080 0.078
78
83
92 -0.044 -0.028
104
Floor
Stair Shuttle & Freight Elevators Tenant
Doors
Floor I-C
Local elevator Banks
II-B III-C
Note: Listed differential pressure values have units of in. H
2
O. Negative (-) values indicate air flowing into a space.
Positive values indicate air flowing out of a space.
Appendix C
178 NIST NCSTAR 1-4D, WTC Investigation
C.3.3 Comparison of ResultsMeasured versus Predicted
Graphical comparisons between measured and predicted pressure differentials, for HVAC Modes 1
through 4, are shown in Figs. C1 through C20. When considering these figures it would appear that in
some cases the agreement between measured and predicted values appears quite good and in some cases
measured and predicted values vary substantially. Since measured pressures were often in the range of
0.01 in. H
2
O to 0.2 in. H
2
O, the overall agreement between measured and predicted values was good, off
by a factor of 12 rather than an order of magnitude. Differences between measured and predicted values
are due to uncertainties and simplifications in the model, as discussed in Sec. C.3.4.
-0.05
-0.04
-0.03
-0.02
-0.01
0
0.01
0.02
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration C
Figure C1. Mode 1, all fans off, stair 1.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 179
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration C
Figure C2. Mode 1, all fans off, stair 2.
-0.03
-0.025
-0.02
-0.015
-0.01
-0.005
0
0.005
0.01
0.015
0.02
0.025
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration C
Figure C3. Mode 1, all fans off, stair 3.
Appendix C
180 NIST NCSTAR 1-4D, WTC Investigation
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration C
Figure C4. Mode 1, all fans off, elevator 50.
-0.02
-0.015
-0.01
-0.005
0
0.005
0.01
0.015
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration C
Figure C5. Mode 1, all fans off, tenant doors.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 181
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C6. Mode 2, all fans set to summer normal, stair 1.
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C7. Mode 2, all fans set to summer normal, stair 2.
Appendix C
182 NIST NCSTAR 1-4D, WTC Investigation
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C8. Mode 2, all fans set to summer normal, stair 3.
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C9. Mode 2, all fans set to summer normal, elevator 50.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 183
-0.1
-0.08
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0.08
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C10. Mode 2, all fans set to summer normal, tenant doors.
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C11. Mode 3, fan set for core pressurization
for the full height of tower, stair 1.
Appendix C
184 NIST NCSTAR 1-4D, WTC Investigation
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.05
0.06
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C12. Mode 3, fan set for core pressurization
for the full height of tower, stair 2.
-0.02
0
0.02
0.04
0.06
0.08
0.1
0.12
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C13. Mode 3, fan set for core pressurization
for the full height of tower, stair 3.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 185
-0.03
-0.02
-0.01
0
0.01
0.02
0.03
0.04
0.05
0.06
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C14. Mode 3, fan set for core pressurization
for the full height of tower, elevator 50.
0
0.005
0.01
0.015
0.02
0.025
0.03
0.035
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C15. Mode 3, fan set for core pressurization
for the full height of tower, tenant doors.
Appendix C
186 NIST NCSTAR 1-4D, WTC Investigation
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C16. Mode 4, floors 59 through 91 set for core pressurization,
all other floors 100 percent recirculation, stair 1.
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C17. Mode 4, floors 59 through 91 set for core pressurization,
all other floors 100 percent recirculation, stair 2.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 187
-0.4
-0.3
-0.2
-0.1
0
0.1
0.2
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C18. Mode 4, floors 59 through 91 set for core pressurization,
all other floors 100 percent recirculation, stair 3.
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0.1
0.15
0.2
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C19. Mode 4, floors 59 through 91 set for core pressurization,
all other floors 100 percent recirculation, elevator 50.
Appendix C
188 NIST NCSTAR 1-4D, WTC Investigation
-0.1
-0.08
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0.08
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A Configuration B Configuration C
Figure C20. Mode 4, floors 59 through 91 set for core pressurization,
all other floors 100 percent recirculation, tenant doors.
C.3.4 Affect of Temperature Variation
As discussed in Appendix B (Sec. B.1.4), testing of Modes 1 through 4 conducted in support of the
1996 HAI/DCE study took place on J une 30, 1996, between 1:00 a.m. and 3:00 p.m. The National
Oceanic and Atmospheric Administration reported a temperature range of 60 F to 75 F on this date.
Since the 1996 study does not document when during the stated time period each test was performed, it is
unclear what the temperature was during each test. This had the potential to impact predicted pressure
differentials due to the impact of stack effect. As seen in Figs. C21 through C25, predicted pressure
differentials varied only slightly over a wide range of assumed temperatures. Therefore, Mode 1 was
calibrated using an outside temperature of 60 F since it was assumed that this test would have logically
been run first, during the early morning period. Modes 2 through 4 were calculated assuming an outside
air temperature of 70 F.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 189
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 70F Configuration A - 60F
Figure C21. Mode 2, all fans set to summer normal,
configuration A at 60 F and 70 F, stair 1.
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 70F Configuration A - 60F
Figure C22. Mode 2, all fans set to summer normal,
configuration A at 60 F and 70 F, stair 2.
Appendix C
190 NIST NCSTAR 1-4D, WTC Investigation
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 70F Configuration A - 60F
Figure C23. Mode 2, all fans set to summer normal,
configuration A at 60 F and 70 F, stair 3.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 191
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 70F Configuration A - 60F
Figure C24. Mode 2, all fans set to summer normal,
configuration A at 60 F and 70 F, elevator 50.
-0.1
-0.08
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0.08
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 70F Configuration A - 60F
Figure C25. Mode 2, all fans set to summer normal,
configuration A at 60 F and 70 F, tenant doors.
Appendix C
192 NIST NCSTAR 1-4D, WTC Investigation
C.3.3 Affects of Wind Velocity
As discussed in Appendix B (Sec. B.1.4), the 1996 HAI/DCE reported that during the testing performed
using Modes 1 through 4, wind was fairly constant at 2 mph. In Figs. C26 through C30, wind was
varied between 0 mph and 20 mph to test whether variation in wind conditions would affect the
calibration of the model. As seen in these figures, wind had no notable effect on the predicted pressure
differentials, and was therefore not considered in the calibration of the model.
-0.2
-0.18
-0.16
-0.14
-0.12
-0.1
-0.08
-0.06
-0.04
-0.02
0
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 0 mph Configuration A - 20 mph
Figure C26. Mode 2, all fans set to summer normal,
configuration A, 020 mph wind, stair 1.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 193
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 0 mph Configuration A - 20 mph
Figure C27. Mode 2, all fans set to summer normal,
configuration A, 020 mph wind, stair 2.
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 0 mph Configuration A -20 mph
Figure C28. Mode 2, all fans set to summer normal,
configuration A, 020 mph wind, stair 3.
Appendix C
194 NIST NCSTAR 1-4D, WTC Investigation
-0.3
-0.25
-0.2
-0.15
-0.1
-0.05
0
0.05
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 0 mph Configuration A - 20 mph
Figure C29. Mode 2, all fans set to summer normal,
configuration A, 020 mph wind, elevator 50.
-0.1
-0.08
-0.06
-0.04
-0.02
0
0.02
0.04
0.06
0.08
0 20 40 60 80 100
Floor
P
r
e
s
s
u
r
e
D
i
f
f
e
r
e
n
t
i
a
l
(
i
n
.
H
2
0
)
Measured Configuration A - 0 mph Configuration A - 20 mph
Figure C30. Mode 2, all fans set to summer normal,
configuration A, 020 mph wind, elevator 50.
CONTAM Model Calibration
NIST NCSTAR 1-4D, WTC Investigation 195
C.4 FACTORS AFFECTING CALIBRATION
There are a number of factors or combinations of factors which may help explain disagreement between
the measured pressures and simulated pressure including:
Pressure gauge accuracy and user error
Actual outdoor temperature as a function of height was not known. The outdoor temperature
in CONTAM is assumed to be uniform across the exterior face of the building
Architectural changes to individual floors that may substantially alter pressure differentials at
individual points. Potential changes include
Number of tenant doors located between core and interior space
Number and size of wall transfer grills
Configuration of tenant doors (i.e., open, partially open, closed)
Presences of convenience stairs
Configuration of stair doors (i.e., open, partially open, closed)
Configuration of stair transition doors (i.e., open, partially open, closed).
Configuration of elevator doors (i.e., open, partially open, closed)
CONTAM simulations assumed all stair doors, stair transition doors, tenant doors, and
elevator doors over the entire height of the building, were completely closed. CONTAM
simulations assumed the number of tenant doors remained consistent on each floor.
CONTAM simulations also assumed the number and size of wall transfer grills remained
consistent on each floor. Additional factors that may have contributed to variations
between simulated and measured pressure are discussed in Section C.2.3.
Tenant changes on individual floors to architecture and/or ventilation system.
Actual supply and return rates on each floor were not measured. Simulations show that the
rate of air being supplied and returned, and more importantly the ratio between them, will
have a dramatic effects build pressures.
Instructions given to building engineers to set various HVAC modes may have been
misinterpreted. For example, when a purge command was given, building engineers may
have aligned equipment differently than what was assumed. Small fans not thought to impact
the sequence may have been inadvertently left on, affecting overall pressure differentials.
C.5 CALIBRATION SUMMARY
During the calibration of the 2004 CONTAM model there were too many independent variables (leakage
types, wind/weather conditions, HVAC system alignments, design versus measured supply/return airflow,
Appendix C
196 NIST NCSTAR 1-4D, WTC Investigation
unknown architectural changes to individual floor plans) to allow for a single model configuration that
matched all of the data for the various ventilation modes tested.
It was therefore determined that in order to capture the uncertainty involved in the type of modeling being
performed, and the potential range of results given this uncertainty that the three separate building
configurations discussed in this appendix (configurations A, B, and C) would be carried forth in the
modeling performed for this report.
C.6 REFERENCES
HAI (Hughes Associates, Inc.) and DCE (Dillon Consulting Engineers). 1996. Smoke Management
Evaluation Study of The World Trade Center Complex Excluding the Concourse and Plaza. Prepared
for the PANYNJ (Port Authority of New York and New J ersey). New York, NY.
Klote, J . H., and J . A. Milke. 2002. Principles of Smoke Management. American Society of Heating
Refrigeration and Air-Conditioning Engineers (ASHRAE). Atlanta, GA.
McAllister, T., ed. 2002. World Trade Center Building Performance Study: Data Collection, Preliminary
Observations, and Recommendations. FEMA 403. Federal Emergency Management Agency.
Washington, DC, May.