/

KtH^ OF

\

^"fEAU O^ ^<^

O

Q

/

IMBS TECHNICAL NOTE 991

U.S. DEPARTMENT OF COMMERCE/National Bureau of Standards

C0MPF2 -A Program for

Calculating Post-Flashover

Fire Temperatures

53

991
NATIONAL BUREAU OF STANDARDS

The National Bureau of Standards' was established by an act of Congress March 3, 1901. The

Bureau's overall goal is to strengthen and advance the Nation's science and technology and

facilitate their effective application for public benefit. To this end, the Bureau conducts

research and provides: (1) a basis for the Nation's physical measurement system, (2) scientific and technological services for industry and government, (3) a technical basis for equity in trade, and (4) technical services to promote public safety. The Bureau's technical work is performed by the National Measurement Laboratory, the National Engineering Laboratory,

and the Institute for Computer Sciences and Technology.

THE NATIONAL MEASUREMENT LABORATORY provides the national system of

physical and chemical and materials measurement; coordinates the system with measurement

systems of other nations and furnishes essential services leading to accurate and uniform

physical and chemical measurement throughout the Nation's scientific community, industry,

and commerce; conducts materials research leading to improved methods of measurement,

standards, and data on the properties of materials needed by industry, commerce, educational

institutions, and Government; provides advisory and research services to other Government

Agencies; develops, produces, and distributes Standard Reference Materials; and provides

calibration services. The Laboratory consists of the following centers:

Absolute Physical Quantities^ — Radiation Research — Thermodynamics and

Molecular Science — Analytical Chemistry — Materials Science. THE NATIONAL ENGINEERING LABORATORY provides technology and technical

services to users in the public and private sectors to address national needs and to solve national problems in the public interest; conducts research in engineering and applied science

in support of objectives in these efforts; builds and maintains competence in the necessary

disciplines required to carry out this research and technical service; develops engineering data and measurement capabilities; provides engineering measurement traceability services; develops test methods and proposes engineering standards and code changes; develops and

proposes new engineering practices; and develops and improves mechanisms to transfer

results of its research to the utlimate user. The Laboratory consists of the following centers:

Applied Mathematics — Electronics and Electrical Engineering^ — Mechanical

Engineering and Process Technology^ — Building Technology — Fire Research —

Consumer Product Technology — Field Methods.

THE INSTITUTE FOR COMPUTER SCIENCES AND TECHNOLOGY conducts

research and provides scientific and technical services to aid Federal Agencies in the selection, acquisition, application, and use of computer technology to improve effectiveness and

economy in Government operations in accordance with Public Law 89-306 (40 U.S.C. 759), relevant Executive Orders, and other directives; carries out this mission by managing the Federal Information Processing Standards Program, developing Federal ADP standards

guidelines, and managing Federal participation in ADP voluntary standardization activities; provides scientific and technological advisory services and assistance to Federal Agencies; and

provides the technical foundation for computer-related policies of the Federal Government.

The Institute consists of the following divisions:

Systems and Software — Computer Systems Engineering — Information Technology.

'Headquarters and Laboratories at Gaithersburg, Maryland, unless otherwise noted; mailing address Washington, D.C, 20234. 'Some divisions within the center are located at Boulder, Colorado, 80303.

^1

The National Bureau of Standards was reorganized, effective April 9, 1978.
C0MPF2--A Program for

Calculating Post-Flashover

Fire Temperatures

?!a*!8n«I BnrMa of $tan4ar<»

JUN 2 ^ ym

u

Vytenis Babrauskas

Center for Fire Research

National Engineering Laboratory

National Bureau of Standards

Washington, D.C. 20234

Kt»»f OP c.

^^ftSAV ov Cr

U.S. DEPARTMENT OF COMMERCE, Juanita M. Kreps, Secretary

Jordan J. Baruch, Assistant Secretary for Science and Technology

NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Director

Issued June 1979
National Bureau of Standards Technical Note 991

Nat. Bur. Stand. (U.S.), Tech. Note 991, 76 pages (June 1979) CODEN: NBTNAE

U.S. GOVERNMENT PRINTING OFFICE

WASHINGTON: 1979

For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, D.C. 20402

Stock No, 003-003-02080-0 Price $3.50

(Add 25 percent additional for other than U.S. mailing).
TABLE OF CONTENTS

Page

LIST OF TABLES iv

LIST OF FIGURES ^^ NOMENCLATURE V

ABSTRACT 1

1. INTRODUCTION , 1

2. HISTORY OF DEVELOPMENT 1

3. THEORY 2

4. PYROLYSIS RATES 5

4.1 Liquid or Thermoplastic Pools 5 4.2 Solid Fuels 6

5. DETAILS OF SUBROUTINES 7 5.1 C0MPF2 7 5.2 CRIB 7

5 . 3 DEQNS 7

5.4 ECHOID 8

5.5 ICONDS 8 5.6 INC 8 5.7 OUTPUT 8

5.8 PFLFIX • 8

5.9 POOL 8

5.10 PP 9

5.11 PVTFIX 9

5.12 RPFIX 9

5.13 STFLOW 9

5.14 TLU 9

5.15 TRIDGF 9

6. AGREEMENT WITH EXPERIMENT 9

7. INPUT INSTRUCTIONS 9

7.1 Deck Set-up 9

7.2 Namelist VARS 10

7.3 Modes of Program Operation 10

8. FILES USED H

9

.

IMPLEMENTATION 12

10, LIST OF VARIABLES 12

11, ACKNOWLEDGMENTS 12

12 , REFERENCES 12

APPENDIX A — SA!4PLE PROBLEMS 20

APPENDIX B — PROGRAM LISTING 38

111
LIST OF TABLES

Page

Table 1. Variables specified in the input Namelist VARS 14

Table 2. List of variables 16

LIST OF FIGURES

Page

Figure 1. Flow chart for main program C0MPF2 19

IV
NOMENCLATURE

D

g

h

h

h

cAh

c Ah

E km

M

oq

Q

r

T

V

area (m^

)

combustion efficiency (-)

discharge coefficient (-)

heat capacity (J kg K ) smallest fuel dimension (m)

vj molecular weight (kg mol ) X thickness dimension (m)

e emissivity (-)— 3 p density (kg m ) -2 -k

a Stefan-Boltzmann constant (W m K ) -1 gravitational acceleration (m s )

^jT. • ^ ,VT -2 ^^-1, Subscripts convective coefficient (W m K ) enthalpy (J)

combustion enthalpy (J)-1 calorific value ( J kg ) -1 total heat of pyrolysis (J kg ) thermal conductivity (W m K ) mass (kg)

initial fuel mass (kg)

heat (J)

heat flow (W)

stoichiometric air/fuel ratio (-)

temperature (K)

regression velocity (m s )

air air

b vaporization

ep excess pyrolysate

f hot gases; pool

o ambient

p pyrolysis

r window radiation

v ventilation, window

w walls, including ceiling

Superscripts

time rate

' ' ' per volume

C0MPF2~A PROGRAM FOR CALCULATING

POST-FLASHOVER FIRE TEMPERATURES

Vytenis Babrauskas

C0MPF2 is a computer program for calculating the

characteristics of a post-f lashover fire in a single

building compartment, based on fire-induced ventilation

through a single door or window. It is intended both

for performing design calculations and for the analysis

of experimental burn data. Wood, thermoplastic, and

liquid fuels can be treated. In addition to the capability of performing calculations for compartments with

completely determined properties, routines are included

for calculating fire behavior by an innovative variable

abstraction method. A comprehensive output format is

provided which gives gas temperatures, heat flow terms,

and flow variables. The documentation includes input

instructions, sample problems, and a listing of the

program. The program is written in Fortran and constitutes an improved version of an earlier program, COMPF

.

Key words: Computer programs--f ire protection; fire

protection; fire resistance; fire tests; fire walls;

safety engineering--fires.

1. INTRODUCTION

With increasing efforts [1-4] towards rational methods of providing

fire endurance for structural building components, it becomes highly desirable for both the designer and the researcher to have available computer

programs for calculating expected fire temperatures and heat transfer through

the building components. A fire is not considered as becoming a threat to a structure and its fire barriers until it reaches the flashover stage. Flashover of a room is defined as that fire stage when the bulk of the room volume

becomes involved in flames. Operationally, this roughly coincides with flames

coming out the door or window, or an upper gas space temperature of around

600° C, or a radiant heat flux at floor level of about 20 kW/m^ . For the

purpose of designing for fire endurance, then, only post-f lashover fires are considered. The present report describes a computer program for calculating

the expected temperatures, heat and mass flows and other variables in

post-f lashover building fires. Different routines are incorporated for

producing design time-temperature curves and for permitting comparative

theoretical curves to be generated based on experimental mass loss rates.

2. HISTORY OF DEVELOPMENT

The first computer program for calculating post-f lashover fire temperatures was developed by Kawagoe [5] , in conjunction with his pioneering studies

leading to a theoretical room fire model. This model was an adaptation of an earlier graphical technique. The main limitations of both the computer program

Numbers in brackets refer to the literature references listed at the end of

this paper.

aai
and the theoretical model was the restriction to ventilation-limited fires.

Fuel-limited fires could only be expressed in terms of an empirical temperature change rate. Magnusson and Thelandersson [6] studied heat release rates

in more detail and produced a model. An unpublished computer program was used to implement that model. The normalized shape of the fire time-temperature

curves was an input variable in this model; the shape was based on sets of

typical measurements. Based on Magnusson and Thelandersson' s theory, Fedock

[7] published a similar computer program with emphasis on prestressed concrete

structures. The first program to provide for theoretically based calculations

of both ventilation-limited and fuel-limited burning was written by Tsuchiya

[8] . It was restricted to fires starting in ventilation control and to fuel

consisting of sparsely-packed wood sticks.

The predecessor to the present program, COMPF [9] , was issued in 1975

and incorporated several new advances, including the ability to treat entirely

fuel-limited fires, to allow for temperature-dependent wall properties, to

permit the optional use of numerical input fuel weight loss rates, and to

perform certain variable abstraction ("pessimization") calculations as an aid

to design. (These techniques enable an input variable to be treated non- deterministically. ) Program C0MPF2 is intended to replace program COMPF and

differs from it in the following main ways:

1. A subroutine has been added to allow treatment of fires where

thermoplastic or liquid fuel exists in the form of a pool on the

floor. The routine implements the theory discussed in reference [10]

and outlined in section 4.1; examples of calculations are also given

and discussed in that reference.

2. The deterministic wood fuel burning model has been extended to include

the possibility of densely-packed cribs.

3. Both pool fire and densely-packed crib options have been incorporated

into the pessimization routines.

'4. In addition to performing transient calculations, the program can now also treat steady-state solutions, for both lossy and adiabatic

walls.

5. The program is now in S.I. units throughout.

6. Certain corrections and improvements have been incorporated in the

calculation routines. The method for the iterative solution of the

heat balance equation has especially been improved.

3 . THEORY

The post-f lashover compartment fire theory has been given in some detail

in reference [11], thus, only a brief summary will be given here. The main

assumptions are:

• The compartment represents a well-stirred reactor, i.e., spatial

temperature variations in the hot fire gases are ignored.

• The model is quasi-steady. Time variations in fuel release rate and

in conduction losses are fully included. However, time rate of change

terms in gas phase mass and energy balance are dropped.

• Air supply and gas outflow is through a single window in a vertical wall

and is the result of fire-induced convection.

The thermal discontinuity away from the window region is at a level below

the bottom of the window. The volume below the discontinuity is occupied
by cold incoming air. In a flashed-over fire this discontinuity is close

to the floor. Its exact location below the window bottom is immaterial [12]

Burning is limited by rates of air or fuel supply rather than by gas

phase chemical kinetics.

Walls (including the ceiling) are modeled as portions of a homogeneous

solid of finite thickness. Temperature-dependent material properties

are allowed for.

The heat balance equation is

^c -"^f ^^T.

- h

298 -Q-Q-Q =0 w ^r ep

(1)

where h denotes enthalpy and the definition of symbols is given in the

Nomenclature section. The subscripts on the enthalpy terms denote the temperature at which they are evaluated. The window radiation loss is, simply

1+ u A a {T_ - T ) (2)

The wall loss term has a radiative and a convective component,

-1 k k Q = A [o -^, =--i-7 =- (T^ - T ) + h (T^ - T )] (3) w w 1/e^ + 1/e -If w f w ^ r w

The convective coefficient h above is not well known since the exact flow

conditions at the wall and ceiling surfaces in a post-f lashover fire are not

known in detail. The convective fraction is much less than the radiative

fraction, permitting a rather simplified treatment. For turbulent-free con- vection flow over flat plates the value for h should depend [13] on

(T_ - T ) ^^3^ p^ value of

f w

h = 5.0 (T. - T ) 1 w

1 / 3 (4)

was selected as being in reasonable agreement with data.

The analysis in [11] shows that for compartments greater than about 2 m

on a side, a flame emissivity of e^=0.9 may be used.

of

The enthalpy evolved from combustion, h , must be evaluated as the lesser

mpep A h b (5)

or

m . Ah b air c p, r(6)
When equation (5) is limiting, the combustion is known as "fuel-limited",

while if equation (6) is smaller, combustion is "ventilation-limited". Here

bp represents the maximum combustion efficiency and is a largely unknown

number. Since it represents, effectively, the "unmixedness" of the combustion,

there may be a scale effect, with smaller compartment spaces showing less

complete mixing. Experimental data can generally be correlated within a range

of about 0.7 < b < 0.9.

p -

The value of m . is obtained from the Bernoulli equation at the window

T . air and is

air

2

3 do 2g1-W^T /W T^

f o' o f [1+(W T./W.T [l+(m /m . )]2)l/3]3

o f f o p air

A /h v V.

(7)

The discharge coefficient has been determined by

to be 0.68 for normal-shaped windows. This value does

where the window takes up almost an entire wall. For

patterns have not been studied, but data can be correl

about one-half its normal value. The molecular weight

is not exactly known since the composition of the gase

unburned fuel gases (excess pyrolysates) is generally

their molecular weight has been assumed equal to that

contribution of carbon monoxide and other minor combus

ignored. The main dependence of m^ij- is on the window

reasonable values of temperature the whole expression

equal to

Prahl and Emmons [14]

not hold in cases

such windows the flow

ated by taking C, at

of the -products, W_,

s, especially the

unknown. For simplicity

of nitrogen. The

tion products is also

parameter A^ ^^v- •^°^ becomes approximately

air ( 0.45 to 0.50 ) A /h v V

(8)

But this approximation has not been employed here.

Ventilation through multiple openings has not been provided for in this

program. An approach for treating such problems is given in [11]

.

The heat of combustion. Ah , is taken as the net value since the hot gas

outflow is above 100° c- The stoichiometric ratio, r, is a constant for a pure

material; a tabulation of values is given in [15].

The rate of pyrolysis, Ap , is one of , the hardest quantities to determine.

A discussion of available values is given in the next section.

Vdr The outflow mass rate, mf , is by mass conservation the sum of A^irThe enthalpy of the outflow products, hm and h " 29 1 and

is evaluated on the

assumption that the combusted fuel goes to CO2 and H„0. No account is

taken of CO for two reasons: because the effect on a mass basis would be

very small, and because it was considered advisable not to introduce any

reaction kinetics. Also, only elements C, H, 0, and N have been considered

for the fuel composition.

The excess pyrolysate term Qgp/ is the heat required to vaporize the excess pyrolysates. Note that with the conventional definition of heat of

combustion, the loss for vaporization of combusted pyrolysates is already

included in A h . c The second major equation to be solved is for heat conduction through

the wall.
3T ^ 3T , , ,

The wall is initially at ambient temperature, T , and is subjected to

boundary conditions at the fire side of: °

w 4 It _k __ = h [T^ - T^(o)] + ea [T^ - T^(o)] (10)

and on the unexposed side (x = L)

,

ci T -^ ^ = ^ f^w(^) - ^o^ ^ ^^ f^w (L) - T^ (11)

For the fire side the convective coefficient has been given above. For

the unexposed side a value of

h = 1.87 [T - T (L) ] 1^3 (12) o w ^ '

was taken.

4. PYROLYSIS RATES

4.1 Liquid or Thermoplastic Pools

There is currently only one fuel arrangement where the pyrolysis rate may

adequately be predicted from theory. It consists of a pool of thermoplastic

or liquid fuel on the floor. The fuel is pyrolyzed solely by radiant flux

and "sees" the compartment with a view factor of 1.0 and itself with a view

factor of zero. In addition, the fuel must pyrolyze at a known surface

temperature, T, , and with a known heat of pyrolysis. Ah . Then:

m = A. ^^ (^f - ^b) (13)

p "f Ah

P

Tewarson and Pion [16] have measured heats of pyrolysis for numerous thermoplastic materials.

The above simple model is fully adequate for steady-state solutions.

At the start of the fire, however, the radiation feedback is small from the

hot gas volume but may be larger from the local plume above the pool itself.

Thus, a plume term should be added in to model the starting transient. Very

limited experimental data by Burgess [17] and J^y Modak [18] can be used to

derive an empirical relationship for the plume pyrolysis rate as:

m = A 0.0014 ^^c (kg/s) (14)

Ah

P

This relationship does not take into account differences in flame

emissivities for various materials; as a result, it only provides a crude

measure. In the present application, however, the contribution of this term

is minor; therefore, an approximate expression is adequate. Also, as the room radiation increases, the effect of plume radiation on pyrolyzing the

5
fuel decreases. For a radiatively black room, at high temperature, the plume

term should properly be negligible. This interaction is crudely modeled by

multiplying the plume term by a proportionality factor before adding to the

far-field term. The proportionality factor, X/ has been set equal to

X = 1.0

4 4 T - T

f b (15)

IVOO'* - Tj^

With x>0.

4.'2 Solid Fuels

Empirical data are available for the mass loss rates of wood planks in

flashed-over fires. Because of the nature of wood combustion, these rates are not especially sensitive to room radiation and can be specified [11] using

a regression velocity of 7-15 ym/s. This relationship is adequate to describe

the burning of large, isolated wood panels. For pieces thin in two or three

dimensions, yet still widely spaced, the following expression is suitable:

1-1/^

(16)

Here M is the original mass, m is the mass at a given time, and F is a constant equal to 2 for cylinders or rectangular sticks and equal to 3 for

m

P .. F . m

' M C ^ M

o

spheres or cubes. C is given by

c =

°^^p

with D being the smallest fuel dimension and v the regression velocity. For

thin fuels v is approximately

_6 _0 .

6

v - 1.7 X 10 D (m/s) (17)

The final arrangement for wood fuel for which data are available is a crib, or a regular stacked array. From the data of Nilsson [19] and Yamashika

[20] , a set of simplified relationships has been evolved for the three crib

burning regimes.

Fuel Surface Control

% = ^ ^p ( S ) ^ ^o (^^5/^) (^^) ^ ^ o

_6 _0.

6

V = 1.7 X 10 D

P
Crib Porosity Control

*p = 4.4 X 10-' (Vh ) % (,9)

c D

\

= ratio of stick clear spacing to crib height

Room Ventilation Control

m = P

0.12 A^ /h^ (20)

In calculations, each of the three rates above are determined and the

lowest rate taken as governing.

5. DETAILS OF SUBROUTINES

The program routines are written in Fortran language. A complete

listing is given in appendix B. The following are brief descriptions of the

operation of each subroutine.

5.1 C0MPF2

C0MPF2 is the main program. It calls most of the calculational routines.

A flow chart of C0MPF2 is given in figure 1. The program starts with the

initialization of certain constants and default values. The input title and

namelist are then read in. If tabular data are specified, subroutine INC is

called. ICONDS is then called in to set initial starting values. The input

data are echoed in ECHOID. After that, the appropriate computational routine

is called in. If no iteration failure has occurred the program then loops

back to the start and goes to the next problem. In case of iteration failure,

the program returns to the same problem, this time printing out additional

intermediate calculation values. This intermediate output can also be forced

to appear by specifying KTRACE=1.

5.2 CRIB

Subroutine CRIB calculates the burning of wood crib fires. A trial gas

temperature value is assumed for the first time step. This value is preset,

but may be overridden by specifying a value of TINPT. The flow quantities

are computed, then the wall losses are determined by calling DESOLV. The heat

balance is then determined. If the normalized residue is greater than 0.002,

the iteration continues. The new temperature is normally determined by the

Newton method. If divergence results, a scanning technique is used initially

and a splitting of differences once a bounded oscillation results. After

successful convergence a new wall temperature profile is established by calling

RSTA. The calculation then proceeds to the next time step. Computation is

terminated at the end of time MTIME, or when gas temperature drops to 353 K, or if errors or convergence failure is detected.

5.3 DEQNS

Subroutine DEQNS computes wall heat conduction using the Crank-Nicolson

method [21]. DEQNS has two entry points: DESOLV and RSTA. The radiation

boundary condition is linearized; updating every iteration rather than every

time step ensures minimal error. An additional within-loop iteration is also

used. DEQNS calls TRIDGF to solve the equation matrix.
5.4 ECHOID

Subroutine ECHOID echoes the input data. The complete data set is given

for each run, rather than just the changed values. Care has been taken to

give physical meaning for the variables printed.

5.5 ICONDS

Subroutine ICONDS initializes starting values and does some preliminary

calculations on the input data. It also makes a few checks on the validity

of the input data. The user, however, is cautioned that this checking is

very rudimentary and in case of error exit or iteration failure the input

data must be carefully examined.

5.6 INC

Subroutine INC is called in when tabular input data are to be read.

5.7 OUTPUT

Subroutine OUTPUT is the primary output routine. It writes at each

time step a large number of variables to output files (logical units) 2 and

3. The temperatures, burning rates, and other primary variables are put on file 2, while the heat balance values and the mass fractions are v/ritten on file 3. OUTPUT also converts temperatures from Kelvin to degrees Celsius

before printing them out.

5.8 PFLFIX

PFLFIX is a pessimization design routine. Fuel pyrolysis rate is

calculated according to governing equations, but the ventilation is pessimized

by instantaneously adjusting the window width to give the highest possible

temperatures. Wood stick or wood crib fuel is assumed unless PLFUEL=T,

in which case a pool fire is used. The window width is not allowed to exceed

a maximum, as set by AWDOW/HWDOW. Calculations stop when the fuel, as specified by FLOAD, is exhausted, since the window width would be undefined

beyond that point. Calculational procedures are similar to those in CRIB.

5 .

9

POOL

POOL is a pool fire burning routine. Computational details are similar

to those as in CRIB. The pyrolysis rate is based on equations 13, 14, and 15.

Three modes of subroutine operation are possible. If STOICH=T, the steadystate temperatures and pool area are determined for stoichiometric burning.

If EISCAN=T, the steady-state solution is found for a given pool area greater

than stoichiometric. The pool area is specified by use of the parameter EITA,

defined as [10]

A /h V V

n = A /h V V

(21)

stoich

For constant window size, this becomes simply a ratio of pool areas. No

solutions are possible for n^l- Finally, a transient calculation can be made,

which proceeds similarly as in the other transient calculations. The user must make sure that the pool size given is sufficiently large so that n^l.
5.10 PP

Subroutine PP is a plotting routine. Details are not given since plotting

routines are dependent on the hardware used.

5.11 PVTFIX

PVTFIX is a pessimization routine, and is effectively the inverse of

PFLFIX. In PVTFIX a fixed ventilation opening is specified. The fuel release

rate is instantaneously varied to always result in the highest possible burning

temperature. Temperatures drop after the fuel load is consumed. Computational

details are similar to those in PFLFIX.

5.12 RPFIX

For comparison of measured data against numerical predictions a routine

is needed which can accept m rates as an input tabular function of time.

RPFIX provides for this type of checking calculation. The case of measured

combustion rate input (as provided, for instance, by oxygen depletion measurements in the window outflow) can also be treated by dividing the measured rate

by A h (net) and setting b = 1.0.

5.13 STFLOW

Subroutine STFLOW is a wall heat conduction routine. It is similar to

DEQNS, except that only the steady-state temperatures are determined.

5.14 TLU

Function TLU is a tabular data interpolating function used in several

subroutines. If the independent variable entered is smaller than the smallest

data point or larger than the largest data point, the output is set equal to

the smallest, or largest dependent value, respectively.

5.15 TRIDGF

Subroutine TRIDGF uses a Gauss elimination procedure to solve a set of

tri-diagonal matrix equations.

6 . AGREEMENT WITH EXPERIMENT

A comparison of numerical predictions with experimental results has been

given in [22] for the program COMPF. Similar agreement should hold for C0MPF2,

since C0MPF2 is improved mainly in operational features, especially increased

versatility, while retaining the same theoretical model as in COMPF. For pool

fires useable full-scale experimental data are not available.

7. INPUT INSTRUCTIONS

7.1 Deck Set-up

The input is assigned to file 1. Each problem run consists of two or three card groups, as follows:

1. Title card (20A4). One card only. Card must be present. The

identifying information from the title card is printed at the head

of the output.

2. Namelist card(s). One or more cards. Details are given in the next

section.
3. Tabular input (optional). This input group is contained only for

the first run and for those ensuing runs where NEWPRP=T. If no tabular input is present, then blank cards must not be inserted. If

tabular input is present, it is arranged as follows:

First card: NCN, NCP, NEM, NR, NQG (1013) . These are the number of

points for the wall thermal conductivity, wall heat capacity, wall

emissivity, mass pyrolysis rate, and wall heat generation rate,

respectively. The number of points may be 0, 1, or greater than 1.

If N=0, then the previous run value is unchanged. If N=l, then it

is assumed the value is a constant, independent of temperature or time. If N>1, then an array is inputted.

Ensuing card(s): These are in the format (8F10.0) and arranged in

pairs (independent, dependent). For wall thermal properties, temperature is the independent variable, while for mass pyrolysis rate it is

time. The order is: CNDA, CPW, EMSA, RPX, QGEN . First all the points

(if any) for CNDA are read in, four pairs per card. Then a new card

is started even if the last card is part-full, and CPW is read in.

The process is continued for EMSA, RPX, QGEN. No blank cards may be

inserted. If N=l for an array, then the constant value is entered

in columns 11-20.

After cards for one run are finished, the cards for the next one are stacked, again with no blank cards.

7.2 Name list VARS

For all non-tabular data, the namelist format was adopted. This

undeservedly obscure Fortran feature is highly advantageous for the present

application. Its features include:

• Semi free- format input

• Variables may be in any order

• Unneeded variable values need not be specified

• Variables needed, but not specified in current run are automatically

set equal to the prior given value.

The namelist card(s) must contain the following information: the first

card must start with $VARS in columns 2-6, then a space, then the desired

values, separated by commas. Input may be continued on continuation cards,

each of which must have columns 1-2 blank; The stream is terminated by a $

after the last variable.

The user is cautioned to check the input carefully, since namelist

format provides for only rudimentary error messages. The namelist VARS

values are written to file 5 when read in. In normal operation file 5 can

be rewound or discarded. If error failures occur, however, the VARS listing

on file 5 may be useful in determining input errors. Table 1 lists all the variables inputted in namelist VARS.

7.3 Modes of Program Operation

Time

Three possibilities are available: complete time-temperature curve calculation, calculation of steady-state temperature for a given wall, or the

calculation of a steady-state temperature for adiabatic walls. To select

adiabatic walls, set ADIA=TRUE. To select steady-state solution for real walls,

10
set STEADY=TRUE. To obtain complete time-temperature curve, set ADIA and

STEADY both FALSE. Note that for some fuel pyrolysis conditions below not

all three possibilities are available.

Fuel Pyrolysis

The following modes of operations are available:

1) Pool fire

Set PLFUEL=TRUE.

a. Time- temperature curve, for given ventilation and pool area. Specify SIZE, Set STOICH and EISCAN both FALSE.

b. Burning conditions at steady state for stoichiometric pool size,

that is, determine values for EITA=1. Set STOICH=TRUE. Do not

input SIZE. Do not set EISCAN=TRUE.

c. Burning conditions for any other EITA. Set EISCAN=TRUE. Specify

EITA. This option must be preceded by the stoichiometric problem

(option lb, above) . SIZE input is not used; if given, the value

is disregarded.

2) Wood crib fire. This is the default option. Set FLSPEC, PLFUEL,

RPSPEC, and VTSPEC all FALSE.

a. Simple stick burning. Must specify a value for REGRES greater

than zero.

b. Nilsson's crib formulas for crib burning in three possible

regimes. Specify REGRES=0 . (default). Also specify SH.

3) Checking option when tabular input pyrolysis rates are given.

Set RPSPEC=TRUE. Also must set NEWPRP=TRUE and give an appropriate

array of RPX.

4) Pessimization over ventilation. Set FLSPEC=--TRUE . Window width is

automatically adjusted, but is no greater than determined by the

inputted value of AWDOW/HWDOW. Program stops when fuel is exhausted.

a. Simple stick burning. Must specify a value for REGRES greater

than zero and set PLFUEL=FALSE.

b. Nilsson's crib formulas for crib burning in three possible

regimes. Set PFLUEL=FALSE and REGRES=0. Also specify SH

.

c. Pool burning. Set PLFUEL=TRUE.

5) Pessimization over fuel pyrolysis rate. Set VTSPEC=TRUE. Fuel

pyrolysis rate is automatically adjusted for pessimal burning

conditions

.

8. FILES USED

The Fortran file logical units must be declared as follows:

File 1 -- Input

File 2 -- Output (echoed input and main calculated variables)

File 3 -- Output (heat balance and mass fractions)

File 4 -- Output (intermediate tracing output - used only if KTRACE=1)

File 5 -- Output (listing of namelist VARS contents)

.

11
File 5 can be arranged to be rewound after each problem so that it will

contain data only in case of error failure.

9. IMPLEMENTATION

Program C0MPF2 has been successfully implemented on a UNIVAC 110 8 computer.

The predecessor program, COMPF, was run on a CDC 6400 computer. The program

uses, as much as possible, only standard Fortran expressions. Minor unavoidable

implementation differences exist, however, in commands associated with file

usage.

10. LIST OF VARIABLES

Table 2 gives a list of all the major problem variables.

11. ACKNOWLEDGMENTS

Ulf Wickstrom (Lund Institute of Technology) assisted in .program

development; Richard Peacock (NBS) helped implement the program.

12. REFERENCES

[1] Pettersson, O., Magnusson, S-E., and Thor, J., Fire Engineering Design

of Steel Structures {Publication 50), StSlbygnadsinstitutet, Stockholm

(1976) .

[2] [Law, M. ] , Design Guide for Fire Safety of Bare Exterior Structural

Steel, prepared for American Iron and Steel Institute by Ove Arup and

Partners, London (1977).

[3] Centre Scientifique et Technique du Batiment, The Behavior of Concrete

Structures in Fire -- A Method for Prediction by Calculation, Nat. Bur.

Stand.' (U.S.), Tech. Note 710-10 (1978), translation from French.

[4] Centre Technique Industriel de la Construction M^tallique, Methode du

Provision par le Calcul du Comportement au Feu des Structures en Acier,

Construction Metallique, No. 4 (1976).

[5] Kawagoe, K. , Estimation of Fire Temperature - Time Curve in Rooms

(Research Paper 29), Building Research Institute, Tokyo (1967). \

[6] Magnusson, S-E., and Thelandersson, S., Time-Temperature Curves of

Complete Process of Fire Development, Acta Polytechnica Scandinavica

(Civil Engineering and Building Construction Series) , No. 65 (1970)

.

[7] Fedock, J. J., Combustion Gas Temperatures in Prestressed Concrete

Apartment Fires, M.S. Thesis, University of Colorado, Department of

Civil and Environmental Engineering (1973).

[8] Tsuchiya, Y., Computer Program for the Behavior of Fire in an Enclosure

(CP 32) , National Research Council, Division of Building Research,

Ottawa (1972) .

[9] Babrauskas, V., COMPF, A Program for Calculating Post-Flashover Fire

Temperatures (Report UCB FRG 7 5-2) , Fire Research Group, University

of California, Berkeley (1975) .

[10] Babrauskas, V., and Wickstrom, U. G., Thermoplastic Pool Compartment

Fires, Combustion and Flame, 34_' 195-202 (1979) .

12
[11] Babrauskas, V., and Williamson, R. B., Fire and Materials, 2_, 39-53

(1978) .

[12] Rockett, J. A., Corabustion Science and Technology, 1_2, 165-175 (1976).

[13] McAdams, W. M. , Heat Transmission, McGraw-Hill Book Co., New York (1954).

[14] Prahl, J., and Emmons, H. W., Combustion and Flame, 25_, 369-385 (1975).

[15] Throne, J. L., and Griskey, R. G., Modern Plastics, 49_, 96-100 (Nov. 1972]

[16] Tewarson, A., and Pion, R. F., A Laboratory-Scale Test Method for the

Measurement of Flammability Parameters (Serial 22524), Factory Mutual

Research Corporation, Norwood, Mass. (1977).

[17] Burgess, D. S., Strasser, A., and Grumer, J., Fire Research Abstracts

and Reviews, 3, 177-192 (1961) .

[18] Modak, A., and Croce, P. A., Plastic Pool Fires (Serial 22361-3),

Factory Mutual Research Corporation, Norwood, Mass. (1976).

[19] Nilsson, L., The Effect of Porosity and Air Flow on the Rate of

Combustion of Fire in an Enclosed Space, Bulletin 18, Division of

Structural Mechanics and Concrete Construction, Lund Institute of

Technology, Lund (1971)

.

[20] Yamashika, S., and Kurimoto, H., Report of the Fire Institute of Japan,

No. 41, 8-15 (1976) .

[21] Carnahan, B. , Luther, H. A,, and Wilkes, J. 0., Applied Numerical

Methods, John Wiley i- Sons, Inc., Somerset, N. J. (1969).

[22] Babrauskas, v.. Fire Endurance in Buildings, Ph.D. dissertation.

University of California, Berkeley (1976).

13

BBBIDnDBm BBBH^mi
Variable

Table 1. Variables specified in the input Namelist VARS

Default Values

First

run Following

run Information

ADIA FALSE pv

AFLOOR none pv

AWALL none pv

AWDOW none pv

BPF none pv

CD 0.68 pv

CFLPC 44.4 pv

CPPYR(2) CPN2 pv

CVGROS none pv

DENSW none pv

DHP none pv

DTIME none pv

EF 0.9 pv

EI SCAN FALSE FALSE

EITA none pv

FLOAD none pv

FLSPEC FALSE FALSE

HFLPC 5.4 pv

HWDOW none pv

I RUN 1 sequential

IX 10 pv

KTRACE

MTIME 360. pv

MWPYR 28.97 pv

NEWPLT FALSE FALSE

NEWPRP TRUE FALSE

NFLPC 0. pv

OFLPC 0. pv

PLFUEL FALSE FALSE

PLOT FALSE pv

PNCH FALSE pv

PRNT none pv

REGRES 0. pv

RPSPEC FALSE FALSE

SH pv

SHAPE 2. pv

SIZE none pv

STEADY FALSE FALSE

STOICH FALSE FALSE

if true, walls area adiabatic and only

steady-state solution is sought

area of floor (m^)

gross area of walls and ceiling (m^) area of window (m^)

maximum fraction of pyrolyzed fuel burned

to be ^1.0

discharge coefficient

percent, by weight, of carbon in fuel

heat capacity of pyrolysis gases (J/kg-K)

upper calorific value for dry fuel (J/kg)

wall density (kg/in-^)

total heat of gasification for fuel (J/kg)

increment of time step (s)

gas emissivity, assumed gray

if true, solve steady-state problem in

POOL for a given EITA

normalized air-fuel parameter for pool

burning

fuel load (kg/m^ floor area)

if true, pessimize ventilation for a specified pyrolysis rate

percent of hydrogen, by weight, in fuel

window height (m) run problem number

number of wall slices, to be _< 10

print intermediate output if =1

maximum fire time (s)

molecular weight of pyrolysis gases (g/gmole)

if true, start new plot frame

if true, new data arrays will be given

percent of nitrogen, by weight, in fuel

percent of oxygen, by weight, in fuel

if true, fuel is a pool fire

if true, plot time-temperature curve

if true, punch time-temperature curve interval at which results are to be

printed (s)

rate of fuel regression (m/s)

if true, use tabular input fuel pyrolysis

ratio of clear spacing between sticks/crib

height for crib

shape factor in pyrolysis equation for

wood sticks

for cribs: smallest dimension of

stick (m)

for pools: pool area (m^)

if true, only steady-state solution is to

be sought

if true, EITA=1 solution is sought in

POOL

14
Table 1 TBOILC 0. pv

THICKW none pv

TINPT 0. 0. VTSPEC FALSE FALSE

WFLPC 0. pv

Note:

(continued)

fuel vaporization temperature for pools (C)

wall thickness (m)

optional input iteration gas temperature

(K)

if true, pessimize pyrolysis rate for a specified ventilation

percent of water, by weight, in fuel

pv = previous value

15

noonm
Table 2. List of variables

AD I

A

AFLOOR

AWALL

AWALLN

AWDOW

BPF

BWDOW

BWORST

C

CD

CFLDC

CND

CNDA

CNG

CNL

CNV

CPA

CPCO

CPC02

CPH2

CPH20

CPN2

CP02

CPPYR

CPW

CVGROS

CVNET

DENF

DENSA

DENSW

DENU

DERIVl

DERIV2

DIF

DTGAS

DTIME

DX

EF

EI SCAN

EITA

EMS

EMSA

FC

FLOAD

FLREM

FLSPEC

FUELPC

Fl

F2

G

IIHCP

HF

HFLPC

HU

HIN

HRATIO

HWDOW

ILINE

true if walls are adiabatic

area of floor (m^)

gross area of walls and ceiling (m^)

AWALL minus window area

area of window (m^)

maximum fraction of pyrolyzed fuel burned

width of window (m)

window width (.LE. BWDOW) which maximizes gas temperatures (m)

moles of carbon in fuel (mole/kg fuel)

discharge coefficient

percent of carbon, by weight, in fuel

conductivity of a given wall slice (W/m-K)

conductivity of the wall, as a function of temperature (W/m-K)

average conductivity, next to higher numbered slice

average conductivity, next to lower numbered slice

numerical factor in heat transfer coefficient

heat capacity of ambient air (J/kg-K)

heat capacity of CO, as a function of temperature (J/kg-K)

heat capacity of CO2/ as a function of temperature (J/kg-K)

heat capacity of H2, as a function of temperature (J/kg-K)

heat capacity of H2O, as a function of temperature (J/kg-K)

heat capacity of N2, as a function of temperature (J/kg-K)

heat capacity of O2/ as a function of temperature (J/kg-K)

heat capacity of pyrolysis gases, as a function of temperature

(J/kg-K)

wall heat capacity, as a function of temperature (J/kg-K)

upper calorific value for dry fuel (J/kg)

lower calorific value for moist fuel (J/kg)

Biot Number/2--f ire side

ambient air density (kg/m^)

wall density (kg/m^)

Biot Number/2--unexposed side

current derivative of heat balance remainder (W/K)

previous derivative of heat balance remainder (W/K)

temperature error in iteration (K)

increment in gas temperature (K)

time increment (s)

wall thickness increment (m)

effective flame grey body emissivity

true if seeking contant EITA^^l solution

dimensionless air/fuel parameter for pool burning

computed wall emissivity for parallel plane problem

wall emissivity, as a function of temperature

true if in fuel control

fuel load (kg/m^ of floor area) mass of fuel remaining at a given time (kg)

true if fuel pyrolysis rate is fixed and ventilation pessimized

percent of original fuel supply still remaining

current heat balance error (W)

previous heat balance error (W)

acceleration of gravity (m/s^)

moles of hydrogen in fuel (mole/kg fuel)

variable in solving differential equation

effective heat transfer coefficient, fire side (W/m^-K)

percent of hydrogen, by weight, in fuel

effective heat transfer coefficient, unexposed side (W/m^-K)

height of neutral plane (m)

fractional height of neutral plan above window bottom

height of window (r") line number

16
Table 2, (continued)

IPG

I RUN

IX

IXC

IXL

J

JM

JPRINT

K

KD

KITER

KNTRL

KTRACE

MTIME

MWIN

MWOUT

MWPYR

N

NCND

NCPW

NEMS

NEWPLT

NEWPRP

NFLPC

NQGEN

NRP

OFLPC

OPENF

PLFUEL

PLOT

PNCH

PRNT

QCONW

QFIRE

QFLOW

QFUEL

QGEN

QRADO

QFLOW

QFUEL

QGEN

QRADO

QRADW

QWLSUM

R

RO

RC

REGRES

RMA

RMF

RP

RPSPEC

RPX

SCAN

SH

SHAPE

SIGMA

SIZE

SIZEl

page number

run number

number of wall slices

number of middle slice

number of penultimate slice

number of current time step

maximum number of time steps

output to be printed every JPRINT time steps

number of trial iterations at any given time step

number of iterations to converge differential equation

equals for normal operation, equals 1 for convergence failure

parameter indicating exit status

print intermediate tracing output if KTRACE-1

maximum time for fire simulation (s)

molecular weight of ambient air (g/g-mole)

molecular weight of exhaust gases (g/g-mole)

molecular weight of pyrolysis gases (g/g-mole)

moles of nitrogen in fuel (mole/kg fuel)

number of points in CNDA table

number of points in CPW table

number of points in EMSA table

true if start new plot frame (not overlay previous one)

true if read in new set of tabular data

percent of nitrogen, by weight, in fuel

number of points in QGEN table

number of points in RPX table

moles of oxygen in fuel (mole/kg fuel)

percent of oxygen by weight, in fuel

opening factor ration (m ) true if pool fire configuration

true if plot time-temperature curve

true if punch time-temperature curve number of times per second output is to be printed

heat transfered to walls by convection (W)

heat generated by combustion (W)

net flow enthalpy (exhaust minus inflow) (W)

heat lost in heating up unburned fuel fraction (W)

wall heat generation, as a function of temp (W/m^)

heat radiated out the window (W) new flow enthalpy (exhaust minus inflow) (W)

heat lost in heating up unburned fuel fraction (W)

wall heat generation, as a function of temp (W/m^)

heat radiated out the window (W)

heat transfered to walls by radiation (W)

total heat removed from compartment and passing into the walls (J)

stoichiometric air/fuel mass ratio

stoichiometric oxygen/fuel mass ratio

rate of burning (kg/s)

rate of fuel surface regression (m/s) mass inflow rate of air (kg/s) mass outflow rate of hot gases (kg/s)

rate of pyrolysis (kg/s)

true if rate of pyrolysis is prescribed as input

rate of pyrolysis, as a function of time (kg/s)

true if search for solution by scanning temperatures

ratio of clear spacing between sticks to crib height

contant indicating shape of fuel sticks

Stef an-Boltzmann constant (w/m^-K**)

thickness of crib sticks (m) area of pool (m^)

pool area for EITA=1 condition (m^)

17

"TM"'"""TM" HHff ri^^fl
Table 2. (continued)

STEADY true if only steady-state calculation to be made

STOICH true if pool fire and EITA=1

TAMB ambient temperature (K)

TGAS gas temperature (K)

TGASl previous value of TGAS (K)

TGAS2 previous value of TGASl (K)

TGASC gas temperature (C)

TGASN closest gas temperature, lower than true (K)

TGASP closest gas temperature, higher than true (K)

TGOLD value of TGAS from prior time step (K)

THICKW wall thickness (m)

TINPT input trial starting gas temperature (K)

TITLE title of this run

TSF wall surface temperature, fire side (K)

TSU wall surface temperature, unexposed side (K)

TTIME total time (s)

Tl old wall temperature profile (K)

T2 new wall temperature profile (K)

T2C wall temperature profile (C)

VAVGIN average inflow velocity (m/s)

VTSPEC true if ventilation is fixed and pyrolysis rate pessimized

W moles of water in fuel (mole/kg fuel)

WA format constant

WB format constant

WFLPC percent of water, by weight, in fuel

WTFUEL initial total mass of fuel (kg)

YC02 mass fraction of CO2 in outflow

YH20 mass fraction of H2O in outflow

YN2 mass fraction of N2 in outflow

Y02 mass fraction of O2 in outflow

YPYR mass fraction of pyrolysates in outflow

18

HP,
r START J

RESET COUNTERS

^1 READ l INPUT

PLTRST 4 rINC

ICONDS

PFLFIX

DOUT

KTRACE=1

KrTER=1

Figure 1. Flow chart for main program C0MPF2

19

HflflBflH nsasHira Kaoi^siaa^Ki
APPENDIX A — SAMPLE PROBLEMS

Given below is a set of ten concatenated input problems. Each problem

is intended to test out a subroutine or other feature of the program. The

output for the problems is given following the input.

TEST PROGRAM FOR POOL FIRE. STEADY STATE* EITA=1.0

SVARS AFLOOR*20« . AWALL=eO« . AWDGIK=4# . BPF = .7 .CD=0 . 68. CFLPC=85 . 7 . CVGR0S=46«5E6.DENSW=790« .DHP=2.4E6.EITA=1 .O.FLOAD=20e .HFLPC=14o3.

HI*DOW=1.5.0FLPC=O.O.P1_FUEL=T.STOICH=T.TBOILC=390«.

THICK W=0« 038. WFLPC=0«0$

001000001

0.17

0*5

TEST PROGRAM FOR POOL FIRE. STEADY STATE. EITA=0«Oi

SVARS EISCAN=T.EITA=0«01 .PLFUEL=T$

TEST PROGRAM WITH DELIBERATE ERROR TO CHECK KTRACE OPERATION

SVARS £ISCAN=T. PLFUEL=T.TBOILC=2000»

$

TEST PROGRAM FOR POOL FIRE. TRANSIENT CASE. SIZE=7,5 M2

SVARS OTIME=60« .MTIME=3600«.NEWPRP=T. PLFUEL=T. PRNT=60. . SI ZE=7. 5.

TBOILC=390«S

000001

840*

TEST PROGRAM FOR WOOD CRIB FIRE. REGRES SPECIFIED

SVARS CFLPC=44«4.CVGROS=i8«8E6.FLOAD=:S0«0«HFLPC=5«4.OFLPC=38.2.

REGRES-1«5E-5.SHAPE=2«0.SIZE=0.05.WFLPC=12.0S

TEST PROGRAM 1 FOR WOOD CRIB FIRE. NILSSON'S FORMULAS

SVARS REGRES=0*0.SH=0«10S

TEST PROGRAM 2 FOR WOOD CRIB FIRE. NILSSON'S FORMULAS

SVARS FLOAD=20«.SH=0«20S

TEST PROGRAM FOR PVTFIX ROUTINE. VARIABLE WALL PROPERTIES

SVARS NEWPRP=T.VTSPEC=TS

004006

£73* 0*2i 372* 0«2i 373« 0«&6

273. 1090. 372» 1090. 373. 47300.

3e4« 5000* 413. 5000* 4&4« 840*

TEST PROGRAM FOR PFLFIX ROUTINE. POOL OPTION

SVARS AWDOW=10«.CFLPC=e5e7.CVGROS=46«5E6.FLSPEC=T.HFLPC=14.3.NEWPRP=T.

OFLPC=0«0.PLFUEL=T.SIZE=5.0.WFLPC=0.0$

OOiOOl

0*17

640*

TEST PROGRAM FOR RPFIX ROUTINE

SV^.RS MTIME=1903« . N£WPRP=T . RPSPEC=TS

000000000003

0*0 0«i2 S20« 0*12 i2ie 0.25

a 073. 0«26

363* 47300

i073« 840*

20
7.D

a: o in N H in I u <

>. a s

-^

, OJ I/) o o t- in z u u <

> a £

ai M I/) o u s

< a I

a. J

UJ

J (-

HI U 3 a

(M K in N

z u < N

> a E •

J 1- O -< <t • J Z O <t J u 11 1 > o

< (M

JC Q liJ II J Z Of II II 4. < < a. \- < o

z >- UJ II X UJ h- 4 (T (j o q: u

X. t- <<-<<<

Ji u. a

UJ in

Q. «

O Q

«

(M O

S + N O

O < ^ t

nj

II UJ in o

«t > o

UJ _i o> a o n

< o:

3fe2

> a

cr a UJ

UJ J Q. a < s

X UJ Q t- i- < z O UJ 13

_l z J "

UJ < -J

q: h- -< « o o

U. t- ID

I- t- z z

UJ UJ

I 13 13

i y N i:

: V o> N

-> • T O

I 00 N

lO CM ~ •

I o •

I + o

lO II •^

I (0 OH I • in -

I P) UJ o

4 h- UJ < + z

in a. > in 3

a «

V in *

Ol s o

• 13 • u Si X

UJ

in N

U V O

a 13 •4 ^ •

in fO a s Ifl a. 13 "<

a 0>n

U. 3 «

1 •

o

o

in Ol

I- t- z z

lU UJ u u

q: a. a. a cc UJ UJ UJ UJ UJ z

a a a a a o

zo N n o o o 1-

»-« • « • • a in H in * 3 « 10 — 03 in £ O II II II II II O a 10 D O

z z Z N * O Ul UJ t • U. O Z 13 Z 13 * ri o

o

J 10 q: UJ q; a

3 < >

U. U I

.. J - O 13 z a . o > » UJ >- a N -

I- <M <

<n K - in II 3 Z « >- m a • J > 2 3 OH

o CD «» q: •- u z >. >

3 II a •<

u. in

u. in li. in o < o -

K 13 s: < I- Z Ul

UJ X m q

1 (3 « « UJ » in I- £

J UJ > u <

< 3 J < J

3 o a u.

»- ct q: u. u < > w

< J a 2 < 3 3 13

q: u u. s I UJ UJ o « >

S J X UJ o o Q. < cr

_l 2 U S U

£ 13

in II J CO J UJ >

< Z 1- « ^ "•

1 u in

1 « z

1 X UJ

1 \- Q

-J « J <

< *M o J Q

3 • Q < < N a (D a t- in o

*4 u ^~ Q a. w* II < z

UJ (- UJ X N in X M <t in J U • Q. 0. • (^ J O 3 £ U) « UJ E o

UJ < 3 a -1 1- K 13 '^ U. UJX

u. u

a. o 1- in<

Ul • a o

£ o 1- « in UJ

1- zo

Iu

UJ£

•^ a

t21

n ttiiuMJiimMiiHiiiiinummi "TM'TM°'°TM Hil KHHinH
3 Z

u.

u

CJ I- 1/1 I u < > Q. I o H m o o < > a z

II z u < > a. £ z

Ml (0V a o M ^ a> in

Q.

V • 15 Z o t • O « (\J (M

tu o z z Z

oc 1- (D

ta4 < (\j o o cri IL J z in o C7> « » CO 00 J h- o ^ * • « Z o •* • O 111 II a > « o

II

aOu II • o z q: II II u.

IL < <

1- <

a: u. O UJ z >• ai II X UJ 1- o 4 a o u a o u a t- < < < < o m u. UJ lU U. UJ o z a a: I «

0! o 3 < o a a UJ W s z < o q; o -" I

t- 1 J o Q z u

in 1 J o z UJ fl lU 1 < _i IH a. « H a u. 3 a o I ^ N ^ S 01 V

-) • -> o

(9 -0 CM «« f >£ o o • V o + o -> +o

•0 II o II 0 in • in 4

CVJ o • n UJ o z + •0 » t- u ^ N o * < z ^ o •* in a 1

:»: OJ t- II

>-

m< 3

CD o zN o I II o 15 o S • u u z a t- in CJi 1 o M J o V * UJ 1 O 1

01

II 1/) UJQ

UJ3

UJ3

u. tV)

a

1-

IM

1-<

in II 1 in

N

II tn ,

>-o > 3

CD z

ex

"*

J

V

UJ

in < > o IZ Z 3 o 1- tUJ UJ J o 1- 1- t- t- 1- Q o ID «^ a M CJ IZ

a: o n z z z z z u Z > > z UJ ro II 1- < a. UJ UJ UJ UJ UJ U. 3 II a a z

cr > II

u

u

u <Ju o IL in

N o s >

UJ

a

a UJ

aa (Xa q: o U.

m u. in

O IZ Z o 1-

a o a: UJ UJ UJ UJ UJ z o < o a ^ k^ Q o U. 3 a. a a a 0. o K IS z <t <0 >

IT

J o 1- z < K z UJ

NJn •

Q. u. < o N

n o o o 1- UJ I in o < o o l-

>- (T k4 • • t • t in I 19 ti^ 1.^ UJ z • 0-

U

Q IX

a UJ

K in * 3 CO t- z a

N 3 < UJ J 0. f^ 00 CD J UJ >• u < UJ C3 o a < z V) £ < X J < J I II Z

-1 I UJ o II II II II IIO 3 o a. u. IIK o a i- H a. CO n U K cr IZ u. < in II u

J < z z z z N <t

u < V CO UJ

J in

UJ o UJ o o lU UJ • u. < -1 0. z < a

J UJ > -J 3 J z o z o z o t n o 3 3 15 < < z H <

u.

J h^ o o UJ o X •-* crU U. Z 1 3 ^ z 1 UJ < J J CD q: o a m 1- UJ Ul o M > -1 1 V in a 1 a. 1- •H UJ X o > H t- II II < 3 J X UJ UJ 1 z UJ 1 1^ o Q 3 < >- X ^ < o UJ O a Q. < q: 3 1 I UJ X 1 u. 1- m u. u I o Z 3 a. a I J z u z o U. t- Q K

J

U CM

in in

J

a

Q.N 00 4 CM

CZ CS •n

in >>ZO COo

in rt

J Q. o

a.

J

UJ s. UJ 4

IZ UJ • IZ S

^4 fN

M UJ

u.

J

u CM u. z d J J 4 4 a J

O o X t o 4 4 CM a o a m IZ 1-U

in

u in

u

HQ Q. • M 4 Z CM

t- »- UJ 3

in in I < 4 J J

a.

u • a o a * in o 3 z z in CM z o •0

IZ UJ 4 in IZJ

K CM

UJX

1-

O UJX

ILOain

22
TEST PROGRAM MITH DECISERATE ERROR TO CHECK KTRACE OPERATION GEOMETRY AND VENTILATION

C0MPF2 VERSION 1.1 - RUN NO. 3

WALL SURFACE AREA = 80.0 M2 FLOOR AREA = 20.00 M2 WINDOW HEIGHT = 1 . 50 M AREA = 4.00 M2 OPENING FACTOR = 4.899 M2 .

S

DISCHARGE COEFF.= .68

-FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 20.0 KG/M2

TOTAL ENTHALPY OF PYROLYSIS= 2.40*06 J/KG BOILING TEMPERATURE=2000. OEG C

PERCENT 8Y WEIGHT

PERCENT

PERCENT PERCENT

PERCENT

FUEL COMPOSITION

CARBON = 85. HYDROGEN = 14. OXYGEN = NITROGEN = . WATER =

R = 14.78

RO- 3.43

HEAT OF COMBUSTION OF DRY FUEL = LOWER ACTUAL HEAT OF COMBUSTION = MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( .1I27»TGAS 1010.) J/KG- MAXIMUM FRACTION OF PYROLYSATES BURNED = .70 GREY-GAS FLAME EMISSIVITY = .900 FUEL AREA= 174.07 M2

46.50+06 J/KG 43.36+06 J/KG

WALL THERMAL PROPERTIES- THICKNESS

DENSITY = .038 M 790. KG/M3

THERMAL CONDUCTIVITY EMISSIVITY = .50

.170 W/M-K

THERMOPLASTIC POOL FIRE PAGE NO. 1 RUN NO.

TEMP

GAS.C WALL TEMPS RC KG/S

EXC .PYR. KG/S

FUEL AIR IN PCT KG/S

N.P. VELOCITY MOL.WT FUEL

M/S CNTRL

ITERATION FAILURE—

TGAS= 1890.00

TGAS.LT.TBOIL TGAS= 1890.0 GO TO NEXT CASE

RUN NO. 3 TGASl TGAS2 1800.00 .00 1810.00 1800.00 1900.00 1810.00

.43*08

.40+08 .10+08 2.43+08

2.40+08

DERI VI K KO KH 1.35+05 1 3 -2.93+05 2 3 -3. 33+05 3 3

1730.86 1791.43

1740.49 1801.51 1827.13 1892.16

4.028+05 -1. 2.891+04 -1. 2.405+03 -9,

009+08

008+08

183+07

4.052+05-59.938

4.083+0 5-59.068 4.363+05-50.558

.009

.001 .000

23
TEST PROGRAM FOR POOL FIRE, TRANSIENT CASE. SIZE=7.5 M2

GEOMETRY AND VENTILATION

WALL SURFACE AREA = 80.0 M2

FLOOR AREA = 20.00 M2

WINDOW HEIGHT = 1.50 M AREA = 4.00 M2

OPENING FACTOR = 4.899 M2.5

DISCHARGE COEFF.= .68 — FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 20.0 KG/M2

TOTAL ENTHALPY OF PYROLYSIS= 2.40+06 J/KG

BOILING TEMPERATURE= 390. DEG C

FUEL COMPOSITION

CARBON = 85.7 PERCENT BY WEIGHT

HYDROGEN = 14.3 PERCENT

OXYGEN = .0 PERCENT

NITROGEN = .0 PERCENT

WATER = .0 PERCENT

R = 14.78

R0= 3.43

HEAT OF COMBUSTION OF DRY FUEL = 46.50+06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 43.36+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( .1127*TGAS + 1010.) J/KG-K

MAXIMUM FRACTION OF PYROLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

FUEL AREA= 7.50 M2

—-WALL THERMAL PROPERTIES r- THICKNESS = .038 M DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY = .170 W/M-K

HEAT CAPACITY = 840. J/KG-K,

EMISSIVITY = .50

I

24
THERMOPLASTIC POOL FIRE

33

34

TIME

S

60.

120.

180. 240.

300. 360.

420.

480. 540.

600. 660.

720.

780.

84 0. 900. 960.

1020. 1080^

1 140. 1200. 1260.

1 320.

1 380. 1440. 1500.

1560.

1620. 1680.

1740.

1 800.

1 86 0. 1920. 1980. 2040. 2100.

TEMP GASiC

739.

826.

852.

870.

882. 891.

898.

903.

907.

91 1. 914. 917.

919.

921. 923.

925. 926.

358.

255.

211.

184.

164. 149.

138. 129.

121.

1 14. 108.

103. 99.

95.

82.

80.

WALL TEMPS RP RC EXC .PVR. FUEL AIR n

C KG/S KG/S KG/S PCT KG/S

513. 26. 25. • 306 .108 . 198 95 .4 2.29

712. 34. 25. .36 2 . 106 . 256 90 .0 2.23

763. 53. 25. .381 .105 .277 84 .3 2.21 798. 83. 26. .396 . 104 .292 78 .3 2.20

820. 117. 27. .406 .104 .302 72 .2 2.19

836. 151. 30. .413 .103 .310 66 .0 2.18

847. 183. 34. .419 .103 .316 59 .7 2.18

857. 213. 39. .424 .103 .321 53 • 4 2.17

86 4. 240. 46. .428 . 103 .325 47 .0 2.17

871 . 265. 53. .431 .103 .328 40 .5 2.17

876. 288. 61. .434 .102 .331 34 .0 2.16

881. 309. 69. .436 .102 .334 27 .5 2.16

885. 329. 77. .4 39 .102 .336 20 .9 2.16

888. 347. 85. • 440 .102 .338 14 .3 2.16

892. 364. 93. .442 .102 • 340 7 .6 2.16

394. 380. 101. .444 .102 .342 1 .0 2.15

897. 394. 108. .445 .102 .34 3 .0 2.15

534. 407. 115. .000 • 000 .000 .0 2.47

409. 414. 122. .000 • 000 .000 .0 2.39

350. 412. 128. .000 • 000 .000 .0 2.32

310. 402. 133. .000 • 000 .000 .0 2.26

280. 387. 137. .000 • 000 .000 .0 2.21

257. 371. 141. .000 • 000 .000 .0 2.16

238. 355. 143. .000 .000 .000 . 2.12

222. 339. 145. • 000 .000 .000 .0 2.09

209. 324. 145. .003 .000 .000 • 2.05

196. 310. 145. .000 .000 .000 .0 2.02

188. 296. 144. .000 .000 .000 .0 1 .99

179. 284. 142. .000 .000 .000 .0 1.97

171. 272. 140. .000 .000 .000 .0 1.95

164. 261. 137. .000 .000 .000 .0 1.92

157. 251. 135. .000 .000 .000 .0 1 .90 151. 24 1. 132. .000 .000 .000 .0 1.88

146. 232. 129. .000 .000 .000 .0 1.86

141. 223. 126. .000 .000 .000 .0 1.85

136. 215. 123. .000 ,000 .000 .0 1.83

60.

120.

180. 240. 300. 360.

420.

480.

540. 600.

660. 720.

780.

840. 900.

960. 1020.

1080.

1140. 1200. 1260.

1320.

1380^ 1440^

ISOO.

1560.

1620. 1680. 1740.

ISOO.

1860. 1920.

1980. 2040. 2100.

THERMOPLASTIC POOL FIRE

HEAT BALANCE

GAS FLO* WND RAD WALL C> PCT PCT PCT 44.886 5.322 1 1.092 51.859 7.174 4.559

54.089 7.950 3.297 55.726 8.546 2.523 56.787 8.943 2.091 57.S86 9.250 1 .787 58.196 9.486 1.574 58.693 9.681 1.41 1 59.095 9.840 1 .283 59.440 9.978 1 .180 59.737 10.097 1.094 59.996 10.201 1.022

6 0.2 24 10.294 .96 60.427 10.376 .907 60.610 10.451 .860 60.7 76 10.518 .818 60.927 10.580 .781 96.257 3.743 -40.770

97.357 2.643 -52.654

97.717 2.283 -58.391 97.914 2.086 -62.064 98.036 1 .964 -64.780

98.120 1.880 -67.042 98.180 1.820 -68.677

98.223 1.777 -69.805

98.258 1 .742 -71 .093 98.285 1.715 -71 .919 98.306 1.694 -72.779

98.324 1.676 -73.334 98.338 1.662 -73.959

98.349 1.651 -74.362 98.359 1.641 -74.862

98.367 1.633 -75.260

98.374 1.626 -75.620

98.379 1.621 -75.777

98.384 1.616 -76.079

WALL RAO

PCT

28.91 22.952

20.091

17.707

16.097

14.79 4 13.785

12.954

12.252

1 I .659

1 1.144 10.691

10.29 9.929

9.604 9.308

9.037

-59.079 -47.486

-41.795 -37.859

-35.129

-33.129

-31 .50 5 -30.128

-29.103 -28.159

-27.407

-26.712 -26.146

-25.603

-25.162 -24.756

-24.394 -24.028 -23.735

4.694*06

4.582+06

4,543+06 4.514+06

4.496+06

4. 481+06 4.470+06 4.461+06

4. 453+06

4. 447+06

4. 442+06 4,437+06 4.433+06 4.429+06

4.426+06

4.423+06

4.4 20+06 0.000

0.000

0. 000

0. 000 0.000

0. 000 0.000

0. 00 0.000

0, 000 0,000

0.000

0. 000

0. 000 0.000

0. 000 3.000

0, 000 0.000

a-WALL SUM

1 .127+08

1 .883 + 08 2.520+08 3.068+08

3,559+08

4.005+08

4.417+08 4.801+08

5.163+08

5.505+08 5.832+08 6.143+08

6.443+08 6.731+08

7.008+08

7.277+08 7.537+08

7.537+08

7.537+08 7.537+08

7.537+08

7.537+08

7.537+08 7.537+08

7.537+08

7.537+08

7.637+08

7.537+08

7.537+08

7.537+08 7.537+08

7.537+08

7.537+08

7.537+08

7.537+08

7.537+08

N.P. VELOCITY MOL. WT FUEL

M/S CNTRL

.38 1.86 29.87 F .37 1.87 29.85 F .37 1.87 29.85 F .37 1.87 29.84 F .36 1.87 29.84 F .36 1.87 29.83 F .36 1.87 29.83 F .36 1.87 29.83 F .36 1.87 29.83 F

• 36 1.87 29.83 F

• 36 1.87 29.82 F

• 36 1.87 29.82 F .36 1.87 29.82 F .36 1.87 29.82 F .36 1.87 29,82 F .36 1.87 29.82 F .36 1.87 29.82 F .43 1.78 26.90 T .45 1.66 26.95 T .45 1.59 26.94 T .46 1.54 26.90 T .46 1.49 26.87 T .46 1.45 26.84 T .47 1 .42 26.80 T .47 1.39 26.78 T .47 1,36 26.75 T .47 1,34 26.73 T .47 1 .31 26.70 T .47 1.30 26.67 T .47 1.28 26.64 T .48 1.26 26,62 T .48 1.24 26,6 T .48 1.23 26,58 T .48 1.22 26,56 T .48 1 .20 26.54 T .48 1.19 26.52 T

PAGE NO. 1 RUN NO.

y02 rN2 YC02 YH20 YPYR

PCT PCT PCT PCT PCT MASS MASS MASS MASS MASS

.063 .679 . 131 .054 .076 .058 .663 . 128 .052 .099

.358 .657 . 127 .052 . 107 .357 .652 . 126 .052 .113 .057 .650 .126 .051 .116 .057 .647 . 125 ,051 . 1 19 .057 .646 . 125 ,051 . 122 .057 .644 . 125 .051 .124

.0 56 .643 .124 .051 .125

.0 56 .642 . 124 .051 .126

.0 56 • 64 1 .124 .051 .128

.056 .64 1 . 124 .051 .129

.056 • 640 . 124 .051 .130

.056 .639 ,124 .051 .130

.0 56 .639 • 123 .051 .131

.0 56 .638 . 123 .051 . 132

.0 56 .638 . 123 ,050 . 132 .230 .770 .000 .000 .330

.230 .770 .000 .000 .000

.230 .770 .000 .000 .300

.230 .770 .033 .300 .030

.230 .770 .000 .000 .000

.230 .770 .000 .003 .000

.230 .770 .030 .330 .000

.230 .770 .000 .000 .000

.230 .770 ,000 .000 .000

.230 .770 .000 .003 .000

.230 .773 • 300 .000 .000

.230 .770 .300 .000 .300 .230 • 770 .000 .000 .000

.230 .770 .000 .000 .000

.230 .770 .000 .000 .000

.230 .770 .030 .000 .000

.230 .770 ,000 .000 .000

.230 .770 .000 .000 .000

.230 .770 .000 .000 .000

25

MMMIKMIMJI
TEST PROGRAM FJDR WOOD CRIB FIRE, REGRES SPECIFIED

GEOMETRY AND VENTILATION

WALL SURFACE AREA = 80.0 M2

FLOOR AREA = 20,00 M2

WINDOW HEIGHT = 1.50 M

AREA = 4.00 M2

OPENING FACTOR = 4.399 M2.5

DISCHARGE COEFF.= .68

FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 10.0 KG/M2

FUEL COMPOSITION

CARBON = 44.4 PERCENT BY WEIGHT

HYDROGEN = 5.4 PERCENT

OXYGEN = 38.2 PERCENT

NITROGEN = .0 PERCENT

WATER = 12.0 PERCENT

R = 5.32

R0= 1.23

HEAT OF COMBUSTION OF DRY FUEL = 18.80+06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 15.07+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( .1127«TGAS + 1010.) J/KG-K

MAXIMUM FRACTION OF PYROLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

RATE OF REGRESSION = 15.00-06 M/S

FUEL DIMENSION = .050 M SHAPE FACTOR = 2.00

WALL THERMAL PROPERTIES

THICKNESS = .038 M

DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY = .170 W/M-K

HEAT CAPACITY = 840. J/KG-K

EMISSIVITY = .50

26
CRIB FIRE

TIME TEMP WALL TEMPS RP RC EXC.PYR. FUEL AIR IN

S GAS.C C KG/S KG/S KG/S PCT KG/S

i 0. 510. 296. 26. 25. .240 .168 .072 92 8 2,35

2 60. 560. 422. 30. 25. .231 .162 .069 85 9 2,35

3 120. 572. 461. 41. 25. .222 .156 • 067 79.2 2.36

4 180. 574. 482. 58. 25. .214 . 150 .064 72 8 2.36

5 240. 571. 490. 78. 26. .205 . 143 .061 66 6 2.37

6 30 0. 563. 492. 98. 28. .1 96 .137 .059 60.8 2.37

7 360. 553. 489. 117. 30. .187 .131 .056 55 2 2.38

a 420. 54 0. 483. 134. 33. .178 .125 .053 49.3 2.38

9 480. 526. 474. 149. 37. .169 .119 .051 44 7 2.39

10 540. 510. 463. 162. 42. .161 .1 12 .048 39.9 2.39

1 1 600. 4 94. 451. 174. 46. .152 .106 .045 35 4 2.40

12 660. 476. 437. 184. 51. .143 . 100 .043 31 1 2.40

13 720. 457. 423. 192. 56. .134 .0 94 .040 27.1 2.40

14 780. 438. 408. 199. 61. .125 .0 87 .037 23.3 2,40

15 840. 417. 392. 204. 66. .116 .081 .035 19.3 2.40

16 900. 396. 375. 209. 70. .107 .0 75 .032 16 6 2.40

17 960. 374. 358. 212. 74. .098 .069 .029 13 7 2.39

18 1020. 352. 34 1. 214. 78. .0 89 .062 .027 11 2.38

19 1080. 328. 324. 216. 81. .080 .056 .024 8.6 2.37

20 1 140. 304. 307. 216. 84. .071 .049 .021 6.5 2,35

21 1200. 279. 289. 216. 87. .061 .043 .018 4 7 2.33

22 1260. 253. 27 1, 214. 89. .052 .036 .016 3.1 2.29

23 1320. 226. 252. 213. 91. .042 .030 .013 1.9 2.25

24 1380. 199. 233. 210. 93. .033 .023 .010 .9 2.19

25 1440. 169. 213. 207. 94. .022 .016 .007 2 2,10

26 1500. 134. 190. 203. 95. .011 .008 .003 1 .96 27 1560. 100. 166. 198. 96. .000 .000 .000 .0 1,74

28 1620. 92. 152. 193. 96. .000 .000 . 000 1,68

29 1680. 87. 143. 137. 96. .000 .000 .000 1.63

30 1740. 83. 135. 181. 96. .000 .000 .000 1.59

31 1800. 79. 129. 175. 95. .000 .000 .000 1.55

CRIB FIRE

TIME HEAT BALANCE Q--WALL GAS FLOW WNO RAD WALL CNV WALL RAO Q-FIRE SUM

PCT PCT PCT PCT w J

0. 55.544 3.303 19 249 21 . 906 2.531*06 6.250+07

60. 63 693 4.410 1 1 .122 20. 80 2.4 38'f06 1 ,092+08

120. 6 7.4 96 4.843 a .615 19. 040 2.3454-06 1.481+08

ISO. 70.479 5.1 13 7 081 17. 334 2.252+06 1 811+08

240. 72.815 5.238 6 .119 15. 82 8 2. 159+06 2.095+08

300. 74 815 5.278 5 .434 14. 47 3 2. 066+06 2.342+08

360. 76 584 5.253 4.919 13. 244 1.973+06 2 557+08

420. 78.203 5.181 4 507 12. 109 1. 880+06 2.745+08

480. 79.721 S.074 4 .160 1 1 . 046 1. 786+06 2 .908+08

540. 81.173 4.941 3 .852 10. 036 1.693+06 3 049+08

600. 82 590 4.787 3 .564 9. 062 1. 599+06 3 .170+08

660. 83 .993 4.61 9 3 .280 8. 112 1.505+06 3 .273+08

720. 85 ,405 4.440 2 .987 7. 173 1 .4 11 + 06 3 ,359+08

780. 86.848 4.254 2.673 6. 232 1.317+06 3.429+08

840. 88 347 4.064 2 .324 5. 281 1.222+06 3 485 + 08 900. 89 926 3.872 1 .924 4. 299 1 . 127+06 3 ,527+08

960. 91 618 3.681 1.461 3. 268 1. 032+06 3 ,556+08

1020. 93 446 3.492 .919 2. 150 9.368+05 3 ,573+08

1080. 95.478 3.308 .318 . 902 8.407+05 3 580+08

1 140. 96 899 3.101 - .220 -. 632 7.440+05 3.580+08

1200. 97.133 2.867 -1 .288 -2. 188 6. 465+05 3 580+08

1260. 97.353 2.647 -3 .035 -3. 814 5.479+05 3.580+08

1320. 97 .559 2.441 -5 .736 -5. 614 4.476+05 3 580+08

1380. 97 .750 2.250 -10 .04 1 -7. 776 3.446+05 3.580+08

1440. 97.922 2.078 -17 571 -10. 734 2.366+05 3 580+08

1500. 98.072 1.928 -34 .79 1 -16. 210 1. 134+05 3.580+08

1560. 98 .139 1.861 -73 .76 -26. 340 0.000 3.580+08

1620. 98 .135 1.865 -74 .639 -25. 412 0.000 3 580+08

1680. 98 .127 1.873 -75 .239 -24. 810 0.000 3. 560+08

1740. 98 .1 16 1.884 -75 .698 -24. 354 0.000 3.580+08

1800. 98 .104 1.896 -76 .063 -23. 994 0.0 00 3 580+08

PAGE NO. 1 RUN NO.

N.P. VELOCITY MOL.WT FUEL

M/S CNTRL

.41 1,80 29,90 T .40 1,82 29,86 T .40 1,83 29,83 T .40 1,83 29,79 T .40 1.83 29.76 T .40 1.83 29.72 T .41 1.83 29.69 T .41 1,82 29,65 T .41 1.82 29.62 T .41 1,81 29.58 T .41 1,81 29.55 T .42 1.80 29.51 T .42 1,79 29.48 T .42 1.78 29.44 T .42 1.76 29.41 T .43 1.75 29.37 T ,43 1.73 29.33 T ,43 1.71 29,30 T ,44 1.69 29,26 T ,44 1.66 29.23 T ,45 1.63 29.19 T ,45 1,59 29.15 T ,45 1,54 29.12 T ,46 1.48 29.08 T ,47 1,40 29.03 T ,47 1.29 28.98 T .48 1.13 23.92 T .48 1.08 28.92 T .48 1.05 28.92 T .49 1.02 28.92 T .49 .99 28.92 T

PAGE NO. 1 RUN NO. yo2 »N2 YC02 yH20 VPYR

PCT PCT PCT PCT PCT MASS MASS MASS MASS MASS

.129 • 699 , 106 .043 .028

.132 ,701 . 102 .041 .027

.136 ,704 .09 8 .040 .026 .139 ,706 .09 4 .038 .025 .143 .709 .091 .037 .024 .147 ,71 1 .087 .035 .023

.150 ,714 ,083 .034 .022

.154 ,716 . 079 ,032 .021 .158 .719 .075 ,030 ,020

.161 .722 ,072 .029 .019 .165 .724 ,068 ,027 ,018

,169 .727 ,064 .026 ,017

, 172 .729 ,06 .024 .016

.176 ,732 ,056 .023 .015

.180 ,735 ,052 .021 .014

. 183 ,737 ,049 .020 .013

.187 .740 ,045 .018 .012

.191 .74 2 .04 1 .017 .01 1 ,194 .74 5 ,037 .015 .010

. 198 .748 .033 .013 .009

.202 .75 .029 .012 .008

.206 .753 ,025 .010 .007

.210 .756 ,021 ,009 .006

.214 ,759 ,017 ,007 .004

.218 .762 .012 .005 .003

.224 .766 .006 .003 .002 .230 .770 ,000 ,000 .000

.230 .770 ,000 .000 .000

.230 ,770 .000 .000 .000 .230 ,770 . 000 .000 .000

.230 ,770 .000 .000 .000

27

l>u.lllu»MIMmnnBBi!cTMDTM sflraB
TEST PROGRAM 1 FOR WOOD CRIB FIRE. NILSSON'S FORMULAS

GEOMETRY AND VENTILATION

WALL SURFACE AREA = 80.0 M2

FLOOR AREA = 20.00 M2

WINDOW HEIGHT = 1.50 M AREA = 4.00 M2

OPENING FACTOR = 4.899 M2.5

DISCHARGE COEFF.= .68

FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 10.0 KG/M2

FUEL COMPOSITION

CARBON = 44.4 PERCENT BY WEIGHT

HYDROGEN = 5.4 PERCENT

OXYGEN = 38.2 PERCENT

NITROGEN = .0 PERCENT

WATER = 12.0 PERCENT

R = 5.32

R0= 1.23

HEAT OF COMBUSTION OF DRY FUEL = 18.80-I-06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 15.07+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( .1127*TGAS + 1010.) J/KG-K

MAXIMUM FRACTION OF PYROLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

RATE OF REGRESSION = 00.00 M/S

FUEL DIMENSION = .050 M SHAPE FACTOR = 2.00

CRIB SPACING/HEIGHT RATIO= .100

WALL THERMAL PROPERTIES

THICKNESS = .038 M DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY = .170 W/M-K

HEAT CAPACITY = 840. J/KG-K

EMISSIVITY = .50

28

BftM
CRIB FIRE PAGE NO. 1 RUN NO.

TIME TEMP KACL TEMPS RP RC EXC .PYR. FUEL AIR IN N.P. VELOCITY MOL. WT FUEL

S GAS.C C KG/S K.G/S KG/S PCT KG/S M/S CNTRL

1 0. 408. 217. 25. 25. .176 . 123 .053 94 .7 2.36 .42 1.74 29.65 T

2 60. 452. 315. 28. 25. .176 . 123 .053 89 .4 2.38 .42 1.78 29.65 T

3 120. 472. 357. 36. 25. . 176 .123 .053 84 .2 2.38 .41 1.79 29.65 T 4 180. 487. 386. 49. 25. .176 .123 .053 78 .9 2.38 .41 1.80 29.65 T

S 240. 498. 408. 64. 26. .176 .123 .053 73 .6 2.38 .41 1.80 29.65 T 6 300. 506. 424. 80. 27. .176 .123 .053 68 3 2.38 .41 1 .81 29.65 T

7 360. 513. 438. 96. 29. .176 . 123 .053 63 2.38 .41 1.81 29.65 T 8 420. 519. 449. 111. 31 . .176 . 123 .053 57 .8 2.39 .41 1.81 29.65 T

9 480. 524. 458. 126. 34. .176 .123 .053 52 .5 2.39 .41 1.32 29.65 T 10 540. 524. 463. 139. 38. .173 . 121 .052 47 .3 2.39 .41 1.82 29.64 T

1 1 600. 512. 458. 152. 42. .165 .115 .049 42 3 2.39 .41 1.81 29.60 T 12 660. 497. 449. 164. 46. .156 .109 .047 37 7 2.40 .41 1.81 29.56 T 13 720. 481 . 438. 174. 51 . .147 .103 .044 33 3 2.40 .42 1.80 29.53 T ! 780. 464. 425. 183. 55. .138 .0 97 .041 29. 1 2.40 .42 1.79 29.49 T 15 840. 445. 411. 191. 60. .129 .090 .039 25 2 2.40 .42 1.78 29.46 T 16 900. 425. 396. 198. 64. .120 .084 .036 21 6 2.40 .42 1.77 29.42 T 17 960. 405. 381 . 203. 68. .111 .078 .033 18 3 2.40 .43 1.76 29.39 T 18 1020. 384. 364. 207. 72. .102 .072 .031 15 2 2.39 .43 1 .74 29.35 T 19 1080. 362. 348. 210. 76. .093 .065 .023 12.4 2.39 .43 1.72 29.32 T 20 1 140. 339. 33 0. 212. 79. .084 .059 .025 9.9 2.38 .44 1.70 29.28 T 21 1200. 315. 313. 213. 82. .075 .053 .023 7.6 2.36 .44 1.67 29.25 T 22 1260. 291. 296. 213. 85. .066 .046 .020 5.6 2.34 .44 1.64 29.21 T 23 1320. 266. 279. 213. 87. .057 .040 .017 3.9 2.31 .45 1.61 29.17 T 24 1380. 240. 26 1 . 211. 90. .047 .033 .014 2.5 2.27 .45 1.57 29. 14 T 25 1440. 213. 242. 209. 91. .038 .027 .011 1.4 2.22 .46 1.51 29.10 T 26 1500. 185. 222. 206. 93. .028 .020 .008 5 2.15 .46 1.45 29.06 T 27 1560. 154. 20 1. 203. 94. .017 .012 .005 2.04 .47 1.36 29.01 T 28 1620. 109. 175. 199. 94. .002 .001 .001 1 .80 .48 1.17 28.93 T 29 1680. 95. 157. 194. 95. .000 .000 .000 1.70 .48 1.10 28.92 T 30 1740. 89. 146. 189. 95. .000 .000 .000 .0 1.65 .43 1.06 28.92 T 31 1800. 64. 138. 183. 95. .000 .000 .000 1.60 .48 1.03 23.92 T 32 1860. 81 . 131. 177. 95. .000 .000 .000 .0 1 .56 .49 1.00 28.92 T 33 1920. 77. 125. 171. 94. .000 .000 .000 .0 1 .53 .49 .98 23.92 T CRIB FIRE PAGE NO. 1 RUN NO. TIME HEAT BALANCE 0--WALL V02 YN2 YC02 YH20 Y?YR GAS FLOW UNO RAO WALU CNV WALL RAD 0-FIRE SUM PCT PCT PCT PCT PCT PCT PCT PCT PCT W J MASS MASS MASS MASS MASS

a. 57.683 2.530 22 .491 17.294 1.856«^06 4 .431+07 .154 .717 .079 .032 .021 60. 64.928 3.282 14.508 17.280 l.856<-06 7.971+07 .155 .717 .079 .032 .021

120. 68.196 3.673 11 .511 16.622 1.856+06 .1 10+08 .155 .717 .078 .032 .021 180. 70.573 3.978 9 .568 15.88 1 1 . 8564-06 1 394+08 . 1 55 .717 .078 .032 .021 240. 72.319 4.215 8 266 15.200 1.856+06 .655+ 08 .155 .717 .078 .032 .021 300. 73.7 06 4.41 1 7 306 14.578 1 . 856*06 1 899+08 .155 .717 .078 .032 .021 360. 74.838 4.576 6 .571 14.020 1.8S6+06 2 128+08 . 155 .717 .078 .032 .021 420. 75.787 4,71 7 5.985 13.515 1 .856+06 2 .345+08 .155 .717 .078 .032 .021 430. 76.597 4.341 5 .508 13.059 1.856+06 2 .552+08 .155 .717 .073 .032 .021 540. 77.663 4.912 5.042 12.382 1.829+06 2.743+08 . 156 .718 .077 .031 .020

600. 79 542 4.86 2 4 .485 11.111 1. 736+06 2 906+08 . 160 .720 .073 .030 .019

660. 81 169 4.755 4 .077 10.00 1 1.643+06 3 045+08 .163 .723 .070 .023 .018

720. 82.701 4.618 3.722 8.962 1.550+ 06 3 .162+08 .167 .726 .066 .027 .017

780. 84.190 4.461 3 .388 7.965 1.456+06 3 .262+08 .171 .728 .06 2 .025 .016

840. 85 673 4.291 3 .052 6.989 1.362+06 3 344+08 .174 .731 .058 .023 .015

900. 87 180 4.113 2 .697 6.017 1.268+06 3 410+08 .178 .733 .054 .022 .014 960. 88.743 3.930 2.307 5.036 1. 174+06 3 462+08 .181 .736 .051 .020 .013 1020. 90.393 3.745 1 .864 4.022 1.080+06 3 500+08 .185 .738 .047 .019 .012 1080. 92.159 3.560 1 .351 2.951 9.850+05 3.525+08 .189 .74 1 .043 .017 .01 1

1 140. 94.085 3.378 .759 1.783 8.897+05 3.539+08 .193 .744 .039 .016 .010

1200. 96.227 3.200 .135 .453 7.938+05 3 542+08 .196 .746 .035 .014 .009 1260. 97 ,021 2.979 - .524 -1.158 6.973+05 3 542+08 .200 .749 .03 1 .013 .008

1320. 97 247 2.753 -1 .857 -2.747 5.998+05 3.542+08 .204 .752 .027 .011 .007 1380. 97 458 2.542 -3 97 -4.448 5. 009+05 3 542+08 .208 .754 .023 .009 .006

1440. 97.655 2.345 -7.264 -6.383 4.001+05 3 542+08 .212 .75 7 .019 .008 • 005

1500. 97 836 2.164 -12.695 -8.823 2.958+05 3 542+08 .216 .760 .015 .006 .004

1560. 97.997 2.003 -23 .094 -12.532 1.840+05 3 542+08 .221 .763 .010 .004 .003

1620. 98.133 1.867 -62 659 -23.888 2.191+04 3.542+08 .229 .769 .00 1 .001 .000 1680. 98 .137 1.863 -74 .329 -25.735 0.000 3.542+08 .230 .770 .000 .000 .000

1740. 98.130 1.870 -75 024 -25.025 0.0 00 3 542+08 .230 .770 .000 .000 .000

1800. 98 .120 1.880 -75.537 -24.512 0.000 3. 542+08 .230 .770 .000 .000 .000 1860. 98 .108 1.892 -75 .938 -24.116 0.000 3 542+08 .230 .770 .000 .000 .000 1920. 98.095 1.905 -76 262 -23.798 0.000 3.542+08 .230 .770 .000 .000 .000

29

HHBHl
TEST PROGRAM 2 FOR WOOD CRIB FIRE, NIUSSON'S FORMULAS

GEOMETRY ANO VENTILATION

WALL SURFACE AREA =: 80.0 M2

FLOOR AREA = 20.00 M2

WINDOW HEIGHT = 1.50 M AREA =^ 4.00 M2

OPENING FACTOR = 4.899 M2 .

5

DISCHARGE COEFF.= .68

FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 20.0 KG/M2

FUEL COMPOSITION

CARBON = 44.4 PERCENT BY WEIGHT

HYDROGEN = 5.4 PERCENT

OXYGEN = 38.2 PERCENT

NITROGEN = .0 PERCENT

WATER = 12.0 PERCENT

R = 5.32

R0= 1.23

HEAT OF COMBUSTION OF DRY FUEL = 18.80+06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 15.07+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( .1127*TGAS 1010.) J/KG-K

MAXIMUM FRACTION OF PYROLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

RATE OF REGRESSION = 00.00 M/S

FUEL DIMENSION = .050 M SHAPE FACTOR = 2.00

CRIB SPACING/HEIGHT RATIO= .200

WALL THERMAL PROPERTIES—

THICKNESS = .038 M DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY = .170 W/M-K

HEAT CAPACITY = 840. J/KG-K

EMISSIVITY = .50

30
CRIB FIRE PAGE NO. 1 RUN NO^

TIME TEMP WALL TEMPS RP RC EXC.PYR. FUEL AIR IN N.P. VELOCITY MOL.WT FUEu

S GAS.C C K6/S KG/S KG/S PCT KG/S M/S CNTRL

J 0. 739. 514. 26. 25. .479 .2 88 . 191 92 .8 2.19 .37 .83 30.60 F

2 60. 858. 745. 34. 25. .461 .2 86 .176 35 .9 2.17 .36 .85 30.60 F

3 120. 902. 815. 54. 2S. .444 .236 . 158 79 .2 2.17 .36 .86 30.60 F

4 160. 936. 366. 35. 26. .426 .286 .141 72 .8 2.17 .36 .87 30.59 F 5 240. 960. 900. 122. 27. .409 .286 .123 66 .7 2.17 .36 .87 30.59 T 6 300. 955. 90S. 159. 30. .391 .274 .117 60 .8 2.18 .36 .88 30.52 T

7 360. 942. 397. 194. 34. .374 .261 .112 55 .2 2.20 .36 .88 30.45 T

8 420. 924. 884. 227. 40. .355 .249 .107 49 .9 2.21 .37 .88 30.38 T

9 480. 903. 366. 255. 47. .338 .237 . 101 44 .8 2.23 .37 .83 30.31 T to 540. 879. 846. 281 . 55. .321 .224 .096 40 .0 2.25 .37 .38 30.24 T

11 600. 853. 82 2. 302. 63. .303 .212 .091 35 .5 2.26 .37 .88 30.17 T 12 660. 824. 797. 32 1. 72. .285 .200 .086 31 .2 2.23 .38 .83 30.09 T 13 720. 794. 769. 336. 80. .267 .187 .080 27 .2 2.30 .38 .33 30.02 T 1* 780. 762. 740. 349. 89. .250 .175 .075 23 .4 2.31 .33 .38 29.95 T IS 640. 728. 709. 360. 97. .232 . 162 .070 20 .0 2.33 .39 .87 29.87 T 16 900. 693. 678. 368. 104. .214 . 150 .064 16 .8 2.35 .39 .87 29.80 T 17 960. 656. 644. 374. 111. .196 .137 .059 13 .8 2.37 .40 .86 29.73 T IS 1020. 617. 610. 378. 117. .178 .125 .053 11 .1 2.38 .40 .85 29.65 T 19 1080. 576. S75. 380. 122. .160 .1 12 .048 8 .7 2.39 .40 .84 29.58 T 20 1 140. 533. 540. 380. 127. .142 .0 99 .042 6 .6 2.41 .41 .83 29.51 T 21 1200. 489. 504. 379. 131 . .123 .086 .037 4 .8 2.41 .42 .81 29.43 T 22 1260. 442. 467. 376. 134. .105 .073 .031 3 .2 2.42 .42 .78 29.36 T 23 1320. 393. 429. 372. 137. .036 .060 .026 1 .9 2.41 .43 .75 29.28 T 2« 1380. 342. 390. 366. 139. .066 .045 .020 .9 2.39 .44 .71 29.21 T 25 144 0. 286. 350. 359. 141 . .0 46 .032 .014 .2 2.35 .44 • 64 29.12 T 26 1500. 223. 307. 350. 142. .023 .016 .0 07 .0 2.25 .46 .54 29.03 T 27 1560. 155. 261. 341. 143. .000 .000 .000 .0 2.06 .47 .37 28.92 T 28 1620. 140. 236. 330. 143. .000 .000 .000 .0 1 .99 .47 .31 28.92 T 29 1680. 129. 219. 318. 143. .000 .000 .000 .0 1.94 .47 .27 28.92 T 30 1740. 121. 20S. 306. 142. .000 .000 .000 .0 1 .69 .48 .24 28.92 T 31 1800. lis. 194. 294. 141 . .000 .000 .000 .0 1.85 .48 .21 28.92 T 32 186 0. 109. 184. 282. 139. .000 .000 .000 .0 1.31 .48 .18 28.92 T 33 1920. 105. 175. 270. 137. .000 .000 .000 .0 1.78 .43 .15 28.92 T 34 1980. 101 . 168. 260. 135. .000 .000 .000 .0 1.75 .48 .13 28.92 T 35 2040. 97. 161. 250. 133. .000 .000 .000 .0 1.72 .48 .11 23.92 T 36 2100. 94. 155. 240. 130. .000 .000 .000 .0 1.69 .48 .09 28.92 T 37 2160. 91. 149. 231. 128. .000 .000 .000 .0 1.66 .48 .07 28.92 T 38 2220. 88. 144. 223. 125. .000 .000 .000 .0 1.64 .43 .05 28.92 T 39 2280. 85. 139. 214. 122. .000 .000 .000 .0 1.61 .48 • 04 28.92 T 40 2340. 83. 135. 207. 119. .000 .000 .000 .0 1.59 .49 • 02 28.92 T 41 2400. 81 . 131. 199. 117. .000 .000 .000 .0 1 .56 .49 .00 28.92 T 42 2460. 79. 127. 192. 114. .000 .000 .000 .0 1.54 .49 .99 28.92 T

CRIB FIRE PAGE NO. 1 RUN NO. TIME HEAT BALANCE Q -WALL Y02 YN2 YC02 YH20 YPYR

GAS FLOII UNO RAD *ALL CNV WALL RAO Q-FIRE SUM PCT PCT PCT PCT PCT PCT PCT PCT PCT w J MASS MASS MASS MASS MASS

0. SI 239 5 .446 1 1 .986 31 .325 4.342+06 1 .128+08 .055 .632 . 176 .0 74 .071

60. 59.890 a .578 4 .838 26.713 4.305+06 1 .943+08 .056 .635 . 177 • 074 .067

120. 62.785 10 .011 3 .335 23.867 4.301 + 06 2 .647+03 .056 .639 . 178 .0 73 .061

lao. 64 837 11 .207 2 .533 21 .406 4.303+06 3 .255+08 .057 .644 . 179 .073 .054

240. 66 133 12 .113 2 .058 19.70 6 4. 310+06 3 .828+08 .057 .648 . 180 .073 .047

300. 68.439 12 .473 1 .715 1 7 . 38 3 4. 125+06 4 .300+08 .064 .653 .173 .070 .046

360. 70.315 12 .500 1 .524 15.674 3.939+06 4 ,707+08 .071 .658 . 166 .067 .044

420. 72.040 12 .361 1 .389 14.222 3.754+06 5 ,058+08 .078 .663 . 158 .064 .042

480. 73.682 12 .102 1 .287 12.939 3.568+06 5 ,363+08 .036 .669 . ISO .061 .040

540. 75 291 1 1 .7S4 1 .203 1 1.759 3.331+06 5 ,626+03 .093 .674 . 142 .058 .037

600. 76.901 11 .339 1 .123 10.638 3. 195+06 5 ,351+08 .101 .679 . 135 .054 .035

660. 78 537 10 .370 1 .056 9.543 3.0 03+06 6 043+08 . 103 .684 . 127 .051 .033

720. 80. 222 10 .360 .978 8.444 2.821+05 6 202+08 .1 16 .690 . 119 .048 .031 780. 81 9 79 9 .318 .888 7.317 2.633+06 6.332+08 .123 .695 .111 .045 .029

840. 83 833 9 .251 .778 6.136 2.445+06 6.433+08 .131 .700 • 103 .042 .027 900. 85 814 8 .668 .64 4.374 2. 256+06 6 508+03 .1 39 .706 .09 5 .033 .025

960. 87. 96 a .074 .464 3.49 6 2.067+06 6 557+08 . 146 .711 .037 .035 .023

1020. 90. 318 7 .476 .243 1.956 1.877+06 6 582+08 . 154 .716 . 079 .032 .021 1080. 92. 946 6 .878 .010 .165 1.686+06 6 583+08 .162 .722 .071 .029 .019

1140. 93. 846 6 .154 - .314 -1.899 1.494+06 6 583+08 .169 .727 .063 .026 .017

1200. 94, 558 5 .442 -1 .020 -4.049 1. 300+06 6 583+08 .177 .733 .055 .022 .015

1260. 95. 238 4 .762 -2 .239 -6.376 1. 103+06 6, 533+08 .185 .738 .047 .019 .012

1320. 95. 883 4 .117 -4 .264 -8.925 9.040+05 6 583+08 .192 .744 .039 .016 .010

1380. 96. 485 3 .515 -7 .700 -1 1.364 6.995+05 6.583+08 .200 .749 .031 .012 .008

1440. 97. 049 2 .951 -14 .093 -15.683 4.853+05 6 583+03 .209 • 755 .022 .009 .006

1500. 97.574 2 .426 -28 .573 -21 .936 2. 439+05 6 583+03 .219 .762 .012 .005 .003

1560. 97. 982 2 .018 -66 .264 -3 3.90 9 0.000 6 583+08 .230 .770 .000 .000 .000

1620. 98. 049 .951 -68 .413 -31 .701 0.000 6 583+03 .230 .770 .000 .000 .000

1680. 93 085 .915 -69 .814 -30.274 0.000 6. 583+08 .230 .770 .000 .000 .000

1740. 98. 107 .893 -70 .393 -29.183 0.000 6.583+08 .230 .770 .000 .000 .000

1800. 98. 122 .878 -71 .769 -28.312 0.000 6 583+08 .230 .770 .000 .000 .000

1860. 98 131 .869 -72 .495 -27.590 0.000 6.533+03 .230 .770 .000 .000 .000

1920. 98. 136 .864 -73 .116 -26.981 0.000 6 583+08 .230 .770 .000 .000 .000

1980. 98. 138 .362 -73 .665 -26.462 0.000 6 583+03 .230 .770 .000 .000 .000

2040. 98. 138 .362 -73 .953 -25.954 0.000 6 583+03 .230 .770 .000 • 000 .000

2 100. 98. 136 .864 -74 .573 -25.606 0.000 6 583+08 .230 .770 .000 .000 .000 2160. 98. 133 .867 -74 .763 -25.201 0.000 6 583+08 .230 .770 .000 .000 .000 2220. 98 128 .372 -75 . 1 14 -24.888 0.000 6. 583+08 .230 .770 .000 .000 .000

2280. 98. 123 .877 -75 .413 -24.604 0.000 6. 583+08 .230 .770 .00 .000 .000 2340. 98. 1 16 .884 -75 .679 -24.347 0. 000 6 583+08 .230 .770 .000 .000 .000

2400. 98. 109 .891 -75 .769 -24.079 0.000 6, 583+08 .230 .770 .000 .000 .000 24-0. 98.100 .900 -75 .979 -23.867 0.000 6, 583+08 .230 .770 .000 • 000 .000

31

.,.i,l.l—^Li.Hn.,..-^^^. mj^gmig HHI
TEST PROGRAM FOR PVTFIX ROUTINE. VARIABLE WALL PROPERTIES

GEOMETRY AND VENTILATION

WALL SURFACE AREA = 80.0 M2

FLOOR AREA = 20.00 M2

WINDOW HEIGHT = 1.50 M AREA = 4.00 M2

OPENING FACTOR = 4.899 M2.5

DISCHARGE COEFF.= .68

FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 20.0 KG/M2

FUEL COMPOSITION

CARBON = 44.4 PERCENT BY WEIGHT

HYDROGEN = 5.4 PERCENT

OXYGEN = 38.2 PERCENT

NITROGEN = .0 PERCENT

WATER = 12.0 PERCENT

R = 5.32

R0= 1.23

HEAT OF COMBUSTION OF DRY FUEL = 18.80+06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 15,07+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = < .1127*TGAS + 1010.) J/KG-K

MAXIMUM FRACTION OF PYROLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

WALL THERMAL PROPERTIES

THICKNESS = .038 M DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY ARRAY (W/M-K)

TEMPERATURE CONDUCTIVITY

273.0 .210

372.0 .210

373.0 .160

1073.0 .260

HEAT CAPACITY ARRAY (J/KG-K)

TEMPERATURE HEAT CAPACITY

273.0 1090.

372.0 1090.

373.0 47300.

383.0 47300.

384.0 5000.

413.0 5000.

414.0 840.

1073.0 840.

EMISSIVITY = .50

32
VENTILATION SPECIFIED. FUEL PYROLYSIS ADJUSTED FOR WORST CONDITIONS PAGE NO.

TIME TEMP

S GAS.C

1 0. 720.

2 60. 853.

3 120. 867.

4 180. 913.

5 240. 933.

6 300. 947.

7 360. 959.

8 420. 969.

9 480. 978.

10 540. 986.

11 600. 992.

12 660. 998.

13 720. 1002.

14 780. 1007.

15 840. 1010.

16 900. 1013.

17 960. 1015.

18 1020. 426.

19 1080. 296.

20 1 140. 239.

21 1200. 203.

22 1260. 178. 23 1320. 160. 24 1380. 145. 25 1440. 134. 26 1500. 124. 27 1560. 1 15. 28 1620. 108. 29 1680. 102. 30 1740. 96.

31 1800. 92.

32 1860. 87.

33 1920. 84.

34 1980. 80. 35 2040. 77.

WALL TEMPS RP RC EXC.PYR. FUEL AIR IN

C KG/S KG/S KG/S PCT KG/S

453. 26. 25. .419 .293 . 126 93 .7 2.23

729. 32. 25. .414 .290 .124 87 .5 2.20

787. 48. 25. .412 • 289 . 124 81 .3 2.19

329. 75. 26. .411 .2 38 .123 75 .2 2.19

860. 96. 27. .410 .287 . 123 69 . 2.18

881 . 131. 29. .409 .287 . 123 62 .9 2.18

898. 146. 32. .409 .2 86 . 123 56 .7 2.17

914. 214. 36. .408 .286 .122 50 .6 2.17

927. 252. 41. .408 .285 . 122 44 .5 2.17

938. 278. 47. .407 .285 . 122 38 . 4 2.17

947. 312. S3. .407 .285 .122 32 .3 2. 17 955. 343. 57. .407 .285 .122 26 .2 2. 16 961. 363. 61 . .407 .285 .122 20 . 1 2.16

967. 379. 64. .406 .284 . 122 14 .0 2.16

972. 392. 67. .406 .284 .122 7 .9 2.16

976. 403. 69. .406 .284 . 122 1 .8 2.16

980. 416. 72. .406 .284 . 122 .0 2.16

606. 428. 73. .000 .000 .000 .0 2.47

459. 433. 74. .000 .000 .000 .0 2.39

385. 419. 76. .000 .000 .000 .0 2.30

335. 392. 77. .000 .000 .000 . 2.22

298. 364. 78. .000 .000 .000 .0 2.15

269. 337. 79. .000 .000 .000 .0 2.08

245. 312. 80. .000 .000 .000 .0 2.02

226. 290. 81. .000 .000 .000 .0 1.96

209. 271. 81. .000 .000 .000 .0 1.91 194. 255. 82. .000 .000 .000 .0 1.85

181. 240. 82. .000 .000 .000 .0 1.80

170. 227. 82. .000 .000 .000 .0 1.76

160. 216. 82. .000 .000 .000 .0 1.71

152. 205. 83. .000 .0 00 .000 .0 1.67

144. 196. 83. .000 .000 .000 .0 1 .63 137. 187. 83. .000 .0 00 .000 .0 1.60

131. 179. 83. .000 .000 .000 .0 1.56

125. 172. 83. .000 .000 .000 .0 1.53

VENTILATION SPECIFIED, FUEL PYROLYSIS ADJUSTED FOR WORST CONDITIONS

60.

120.

180. 240. 300.

360. 420.

480.

540. 600.

660. 720.

780. 840.

900.

960.

1020. 1080.

1 140. 1200.

1260. 1320.

1380.

1440.

1 500. 1560.

1620.

1680. 1740.

1800.

1860. 1920. 1980.

040.

GAS FLOW

PCT 48.339

58.154

60.694 62.608

64.103 65.155

66.052 66.825

67.492

68.080

68.577

68.996

69.354

69.673 69.940

70.160 70.342

95.334 96.914

9 7.433

97.701 97.857

97.958

98.025

98.071

96.101

98.121 98.133

98.138

98. 136 98.135

98.126

98.1 19 98.108

98.095

HEAT BALANCE

WNO RAO WALL CNV PCT PCT 4.946

8.318 9.415

10.308

11.047

1 1.590 12.068

12.493

12.868 13.204 13.494

13.741 13.956

14. 149 14.313

14.448

14.561 4.666 3.086

2.567

2.299

2.143 2.042

1.975 1.929

1 .899 1.879

1 .867 1.862

1 .862 1.865

1.872

1.881 1.892 1.905

14.745 5.410

4.060

3.244

2.695

2.344

2.081

1.862 1.686

1 .545 1.429

1.335

1.254 1.187

1.132

1.088

1 .052 -34.635 -47.715

-54.773 -59.342

-62.819

-65.406

-67.489 -69.232

-70.589

-71.719 -72.665

-73.461 -74.133 -74.703

-75. 188 -75.603

-75.960

-76.269

WALL RAD

PCT 32.022

28.134 25.768

23.808

22. 137

2 0.86 19.804

18.80 7 17.929 17.190

16.52 9 15.963 15.443

15.000

14.625

14.315

14.056 -65.520

-52.475

-45.402

-40.481 -37.109

-34.51

2

-32.461 -30.849

-29.463 -28.347

-27,396

-26.598

-25.925

-25.355

-24,871 -24,459

-24.108

-23.809

4.419+06

4.3634-06

4. 347+06 4.334+06

4.324+06

4.3 17+06 4.310+06

4. 305+06

4. 301+06

4. 296+06 4.293+06

4. 290+06 4.288+06

4. 285+06 4.284+06

4.282+06 4.281+06

0.000 0.000

0.000

0. 000 0.000

0.000 0.000

0.000

0.000 0.000

0.000

0.000

0. 000 0.000

0.000

0. 000

0.0 00 0.000

0-WALL

SUM

1.240+08

2.1 18+08 2.896+08

3.599+08

4.244+08 4.845+08

5,41 1+08 5.945+08

6.451+0 8 6.934+08

7.396+08

7.841+08

6.271+08

8.667+08

9.092+08 9.486+08

9.676+08

9.876 + 08 9.876+08

9.676+08

9.876 + 08 9.876+06

9.876+08 9.876+08

9.876+08

9.876+08

9.876+08 9.876+08

9.676+06

9.876+08 9.876+08

9.876 + 08 9.876+08

9.876 + 08 9.S76+08

N.P. VELOCITY MOL.WT FUEL

M/S CNTRL

,38 1.84 30.59 F .37 1.86 30.59 F .37 1.87 30.59 F ,36 1.67 30.59 F .36 1.67 30.59 F .36 1.67 30.59 F .36 1.87 30.59 F

• 36 1.88 30.59 F .36 1.88 30.59 F .36 1.88 30.59 F .36 1.88 30.59 F .36 1.88 30.59 F .36 1.88 30.59 F .36 1.88 30.59 F ,36 1.88 30.59 F .36 1.88 30.59 F .36 1.88 30.59 F .43 1.79 28.92 T .45 1.67 26.92 T .45 1.58 28.92 T .46 1.50 28.92 T .47 1 .44 28,92 T .47 1.38 28.92 T .47 1.33 28.92 T .47 1.29 28,92 T .46 1.25 28.92 T .48 1.21 28.92 T .46 1.17 28.92 T .48 1 .14 28.92 T .48 1.11 28.92 T .48 1 .08 28.92 T .48 1.05 28.92 T .46 1.03 28.92 T .49 1.00 28.92 T .49 .96 26.92 T

PAGE NO. 1 RUN NO.

Y02 YN2 YC02 YH20 YPYR PCT PCT PCT PCT PCT MASS MASS MASS MASS MASS

.057 .648 .180 .073 .047 .057 .646 . 180 .073 .047

.057 .64 8 . 180 .073 .047 ,057 .648 . 180 .073 .047

,057 .648 . 180 .073 .047

.057 .648 . 180 .073 ,047 .057 .648 , 180 .073 .047

.057 .648 . 180 .073 .047 .057 .648 .180 .073 .047 .057 .648 . 180 .073 .047 .057 .646 . 160 .073 .047 .057 .648 . 180 .073 .047 .057 .648 . 180 .073 .047

.057 .64 8 . 180 .073 .047

.057 .648 , 180 .073 .047 .057 .64 8 . 180 .0 73 .047 .057 ,64 8 . 180 .073 ,047 .230 .770 .000 .000 ,000

.230 .770 .000 .000 ,000

.230 .770 .000 .000 .000

,230 .770 .000 .000 .000

,230 ,770 ,000 .000 .000

,230 .770 .000 .000 ,000

,230 .770 .000 .000 ,900 .230 .770 ,000 ,000 .000

.230 .770 .000 .000 .000

.230 .770 ,000 .000 .000

.230 .770 .000 .000 .000

.230 .770 .000 .000 .000

.230 ,770 .000 .000 .000

.230 .770 .000 .000 .000

.230 .770 .000 .000 ,000 .230 .770 .000 .000 ,000 .230 .7 70 .000 .000 .000

.230 .770 .000 .000 .000

3^

''''TMTMTMTM- nn
TEST PROGRAM FOR PFLFIX ROUTINE. POOL OPTION

-GEOMETRY AND VENTILATION

WALL SURFACE AREA = 80.0 M2

FLOOR AREA = 20.00 M2

WINDOW HEIGHT = 1.50 M AREA = 10.00 M2

OPENING FACTOR = 12.247 M2.5

DISCHARGE COEFF.= .68

FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 20.0 KG/M2

TOTAL ENTHALPY OF PYROLYSIS= 2.40+06 J/KG

BOILING TEMPERATURE= 390. DEG C

FUEL COMPOSITION

CARBON = 85.7 PERCENT BY WEIGHT

HYDROGEN = 14.3 PERCENT

OXYGEN = .0 PERCENT

NITROGEN = .0 PERCENT

WATER = .0 PERCENT

R = 14.78

R0= 3.43

HEAT OF COMBUSTION OF DRY FUEL = 46.50+06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 43.36+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( •1127*TGAS + 1010.) J/KG-K

MAXIMUM FRACTION OF PYROLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

FUEL AREA= 5.00 M2

WALL THERMAL PROPERTIES

THICKNESS = .038 M DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY = .170 W/M-K

HEAT CAPACITY = 840. J/KG-K

EMISSIVITY = .50

34

I

If
O I O l- Z Q 3 X

-I

-Ja

Ui t- 3 Z

ILU « {ft O V o s J

ai>

<^t^\^\C\OvO'OOiO*0^v0^vO

U.U.U.U.U.U.U.U.U.U.I1.U.U.

a <n > H o • a o < > a. £ zz3a o in

CM t- in I u < > a. X • (M (n Q o 1- in Z u o < > Q. £

LUO<

a. in

CM H- in z u < > a £

c^•<4>o^-o>CJloo-t-«cJCM<M 0OC500OOOO0OOOO

iniflinininininininininininin OOOOOOOOOOOOOO — cMoceoNi\0>oioiflioin«»-* CM-«-<00000000000

az

V)zo

1^HM z O M (fl z V o (X o u - ^ < -Maa3 J K

UJ

u

K 3 a a u.

u.o

111 •

ta

V) > in D a V

-1 • o o <j a < X

UJ zo

•4 in H W N < a o J ^ MHZ

lU in > a s a o • s^ Q

UJ

»rf u.u

UI in

a. Q. U) z

UJ «% H J

U o J o J

a. <X

U>M«>Joa> u

a. Q. • £ in J UJ <

UJ 1- tf 3

U.

UI£

-.

w

H

^r)HC^aO(SSN'0>o>Olnu)U)

N««ioininuiu)inif)i/)ioinin • •##•••••••••• inininininininininiflininmin

cMr)r)OvO-«>oonif.(7>cM-*o

soninooo(oinN(7>cM*

o>ooiDf>-'Oioin<»ncM

— •«

iOoinoo)<»o>nvoo>'<nin>o — o>ocM<Mi^n<t-**u)ininin — -CMCMCMCMCMCMCMCMCMCMCMCM

CMCMCMNCMCMCMCMCMCMCMCMCMCM

co^-•ln^oc^r)^oc^l^^oooc7' (DinNCO(>0>00""-"-"-<""-' m^it^^^inininininininin

ininiOiO<oc\j(0>oinin'Of>.NN (MCMCMCMCMl»)r)*lflv0N(0»O

c^o^nln-^<tloc^-«<5N<M•lt -» lr>^'^>o-•u)<^cMlOC^-<*^o _j «^[MCMCMnr)n*<»* o

>ocMo><fcoeoin-<'Oo*'00>'H

ioinNo-«cMt<)<t'*inioininvo (DOO-««-'H« — -<•<-<-•<

*4\00^<^vOO^>OCOOCMrOlO

inCM^invOiONNNNIOIDlOflO

OOOOOOOOOOOOOO <OCMIO^O\OCMOO<tO'OCMCD M•4CMnn'<t^u)'04)^-^>

•<cM^<^l/>'0^a>c^o•"CMlo*

(M 1- in o u <

en > a z zo

N4

1-^azo J u J< £ H SOT in 1 in

a. o J o < u u

UJ

J

a

a

J

in <

o^ UI

-1

u o z o < a J

Q < <

CO o:H

inu

^4

KQa

in < z

>• UJ 3 J r oa>a o

-1

J

H

UJ U.U 3 a

U. <n

ni^CMCMCMCMNCMCMCMCMCMCMCM OOOOOOOOOOOOOO

cDCOQOCoasoococoooaooocoioio OOOOOOOOOOOOOO + + + + + + + 'f + + + + + + ^-cM<^^0ln0l•«ortO5l-<c^n coocon**(ocMNo»o*»oroN^ NN*— N.CMN-HinO>nNO* • ••••••••••••• •"CMM'**inin>o>oioNNioa)

OOOOOOOOOOOOOO + + + + + + + + + + + + + SNDO^Nin^nCMCM-lOOO

o-<eoocMln*lJ>*<ocMl^J^^o< N'*NO«CMNCM<0U)O'0inN

-<CM^nC^O>CM^OCMCOlnCMOOO^O

cMcMoi0Nr»>04)ininu)u)<«'4'

CM - •<

OOO^rt»4\0»HvOlOvON«-«ID»0

to<tr<>ir)«o>i/)-«iD'0<tr)-<o

-.c^^ln<^^ol«)Ioc^JcMcMCMcMf^J

|»^n>oNO>*c^r»>o^N*<^cMlo cMl/)lO<^*o«lDlnol^<t<»cM

l«)Nlflma)CM<t>0C0O-<CMnl<t'

O00OiOO.<-<»^

— CMCMCMCMCM rt— -<CMCMCMCMIMCMCMtM(MCMCM

ln(^•*l«)c^^o^•l')ocMNcM•<lO noOinCMNONO'NCMiOIDIO

cMo-<-"N-<ins(?i

— n*in>o

C\IIS0>OO<4>4M>4CMCMCMCMCM

OOOOOOOOOOOOOO lOCMCD^OvOCMOJ-^OvOCMCO

'^^c\lmn<t'tlt)•o•ot'^^

35

u KHUUUUIUIUrai IHHia '^'"•"'""'^ I^IHHHHliHiM
TEST PROGRAM FOR RPFIX ROUTINE

GEOMETRy AND VENTILATION

WALL SURFACE AREA - 80.0 M2

FLOOR AREA = 20.00 M2

WINDOW HEIGHT = 1.50 M AREA = 10.00 M2

OPENING FACTOR = 12.247 M2.5

DISCHARGE COEFF.= .68

FUEL LOAD PROPERTIES

FIRE LOAD PER FLOOR AREA = 20.0 KG/M2

FUEL COMPOSITION

CARBON = 85.7 PERCENT BY WEIGHT

HYDROGEN = 14.3 PERCENT

OXYGEN = .0 PERCENT

NITROGEN = .0 PERCENT

WATER = .0 PERCENT

R = 14.78

R0= 3.43

HEAT OF COMBUSTION OF DRY FUEL = 46.50+06 J/KG

LOWER ACTUAL HEAT OF COMBUSTION = 43.36+06 J/KG

MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = 28.97

CP OF PYROLYSIS GAS = ( .1127*TGAS + 1010.) J/KG-K

MAXIMUM FRACTION OF PYRCLYSATES BURNED = .70

GREY-GAS FLAME EMISSIVITY = .900

RATE OF PYROLYSIS (KG/S)

TIME RP

0. . 120

120. .120

121 . .250

WALL THERMAL PROPERTIES

THICKNESS = .036 M DENSITY = 790. KG/M3

THERMAL CONDUCTIVITY = .170 W/M-K

HEAT CAPACITY = 840. J/KG-K

EMISSIVITY = .50

36

i
INPUTTED VALUES OF RP ARE USEO

TIME TEMP ALU TEMPS RP RC EXC .PVR. FUEL AIR IN N.P. VELOCITY MOL. WT FUEL

S GAS.C C KG/S KG/S KG/S PCT KG/S H/S CNTRL

1 0. 432. 234. 26. 25. .120 .084 .036 98.2 6.11 .43 .79 29,23 T

2 60. 457. 323. 29. 25. .120 .0 84 .036 96..4 6.13 .42 .80 29.23 T 3 120. 467. 355. 37. 25. .120 .084 .036 94.6 6.13 .42 .81 29.23 T 4 180. 805. 666. 51 . 25. .250 .175 ,075 90..9 5.98 .39 ,91 29.56 T 5 240. 843. 756. 71. 26. .250 .175 .075 87.. 1 5.95 .39 .92 29,56 T 6 300. 856. 784. 97. 27. .250 .175 .075 83..4 5.94 .39 .92 29.56 T

7 360. 866. 804. 128. 29. .250 .175 .075 79.,6 5.94 .39 .92 29.56 T 8 420. 872. 817. 160. 32. .250 .175 .075 75.9 5.93 .38 .92 29.57 T 9 480. 877. 827. 190. 37. .250 .175 . 075 72..1 5.93 .38 .92 29.57 T 10 54 0, 881 . 835- 218. 42. .250 .175 .075 68.4 5.93 .38 .92 29.57 T 11 600. 884. 84 1. 244. 49. .250 .175 .075 64.6 S.92 .38 .92 29.57 T 12 660. 887. 846. 267. 56. .250 .175 .075 60.9 5.92 .38 .92 29.57 T 13 720. 889. 851. 289. 64. .250 .175 .075 57.. 1 5.92 .38 .92 29.57 T 14 780. 891. 855. 309. 72. .250 .175 .075 53.,4 5.92 .38 .92 29.57 T 15 840. 893. 858. 327. 80. .250 .175 .075 49..6 5.92 .38 .92 29.57 T 16 900. 894. 86 1. 344. 88. .250 .175 .075 45..9 5.92 .38 .92 29.57 T 17 960. 896. 86 4. 360. 95. .250 .175 .075 42., 1 5.92 .38 .92 29.57 T 18 1020. 897. 366. 375. 102. .250 . 175 .075 38.4 5.91 .38 .92 29.57 T 19 1080. 898. 868. 389. 109. .250 . 175 .075 34..6 5.91 .38 .92 29.57 T 20 1 140. 899. 870. 401. 116. .250 ,175 .075 30.,9 5.91 .38 .92 29.57 T 21 1200. 900. 872. 413. 122. .250 .175 .075 27.. 1 5.91 .38 .92 29.57 T 22 1260. 901 . 874. 425. 127. .250 .175 .075 23.,4 5.91 .38 .93 29.57 T 23 1320. 902. 876. 435. 133. .250 .175 .075 19.6 5.91 .38 .93 29.57 T 2* 1380, 903. 877. 445. 138. .250 .175 .075 15.,9 5.91 .38 .93 29.57 T 25 1440, 904. 878. 454. 142. .250 . 175 .075 12.. 1 5.91 .38 .93 29.57 T 26 1500. 904. 380. 463. 147. .250 ,175 . 075 8..4 5.91 .38 .93 29.57 T 27 1560, 905. 88 1. 471 . 151. .250 .175 .075 4,.6 5,91 .38 .93 29.57 T 28 1620, 90S. 882. 478. 154. .250 . 175 .075 ..9 5.91 .38 .93 29.57 T 29 1680, 906. 883. 485. 158. .250 .175 .075 .,0 5.91 .38 .93 29.57 T 30 1740. 206. 476. 491 . 161 . .000 .000 .000 ,0 5.58 .46 .51 28.92 T 31 1800. 148. 357. 491. 164. .000 .000 .000 ,,0 5.08 .47 .34 28.92 T 32 1860, 126. 307. 481 . 167. .000 ,000 .000 .0 4.80 .48 .26 28.92 T 33 1920, 1 12. 27 2. 464. 169. .000 .000 .000 •,0 4.58 .48 .19 28.92 T

INPUTTED VALUES OF RP ARE USED PAGE NQ. 1 RUN NO.

TIME HEAT BALANCE Q-WALL Y02 YN2 YC02 YH20 YPYR GAS FLOW WND RAO WALL CNV WALL RAD 0-FIRE SUM PCT PCT PCT PCT PCT

PCT PCT PCT PCT w J MASS MASS MASS MASS MASS

0. 75,.919 3.719 1 1 .044 9. 317 3.642»-06 4.450+07 . 179 .755 ,04 2 .017 .006

60. 80..960 4.296 6 .596 8. 147 3.642+06 7.672+07 . 179 .755 .042 .017 .006

120. 82,.877 4.533 5 ,171 7. 419 3.642*06 1 .042 + 08 . 1 80 .755 • 04 2 .017 .006

180. 72 .523 10.016 3 ,312 14. 154 7.588+06 1.837+08 .124 ,739 .088 .036 .012

240. 76-.015 1 1.532 1 .779 10. 659 7.588+06 2.404+08 , 124 ,739 . 089 .036 .012

300. 77 .201 12.086 1 .379 9. 318 7.588+06 2,891+08 .124 ,739 .089 .036 .012 360. 78 .063 12.502 1 . 122 8. 308 7.588+06 3.320+08 .124 .739 .089 .036 .012

420. 78,.650 12.792 .963 7. 606 7.588+06 3.710+08 . 124 .739 .089 .036 ,012

480. 79 .090 13.013 .848 7. 049 7.588+06 4.070+08 .124 .739 .089 .036 .012

540. 79 .441 13.191 .762 6. 606 7.588+06 4.405+08 .124 .739 .089 .036 ,012

600. 79 .729 13.339 .694 6. 239 7.588+06 4.721+08 , 123 .739 .089 .036 .012

660. 79 .970 13.464 .639 5. 927 7. 588+06 5.020+08 .123 -739 .089 .036 .012

720. 80,.176 13,571 .594 5. 659 7.588+06 5.305+08 .123 .739 .089 .037 .012 780. 80 .355 13.665 .556 5. 424 7.588+06 5.577+08 .123 .73 9 .089 .037 .012

840. 80 .513 13.748 .523 5. 217 7.588+06 5.838+08 .123 .739 .089 .037 .012

900. 80 .653 13.822 .494 5. 031 7.588+06 6.090+08 .123 .739 .089 .037 ,012 960. 80,.778 13.389 ,469 4. 86 4 7.588+06 6.332+08 ,123 .739 .089 .037 .012 1020. 80,.892 13,949 ,447 4, 71 2 7.588+06 6.567+08 ,123 .739 .089 .037 .012 1080. SO .996 14.005 .427 4. 573 7.588+06 6.795+08 . 123 .739 .089 .037 .012

1140. 81,.090 14.055 .409 4. 445 7.588+06 7.016+08 .123 .739 -089 .037 .012 1200. 81 .178 14. 102 .39 3 4. 32 8 7.588+06 7.231+08 .123 .739 .089 .037 .012

1260. 81,.258 14. 146 .378 4. 219 7.588+06 7.440+08 ,123 .739 .089 .037 .012

1320. 81 .333 14. 186 .364 4. 11 7 7.588+06 7.644+08 .123 .739 .089 .037 .012

1380. 81 .402 14.224 .352 4. 023 7. 588+06 7.643+08 , 123 .739 .089 .0J7 .012

1440. 81 .466 14.258 .340 3. 935 7.588+06 8.038+08 .1 23 .739 . 089 .037 .012

1500. 81,.526 14.291 .330 3. 853 7.588+06 8.228+08 .123 .739 .089 .037 .01 2 1560. 81 .582 14.322 ,320 3. 776 7.588+06 8.415+08 .123 .739 .089 .037 .01 2 1620. 81 .635 14.350 .311 3. 704 7. 588+06 8.598+08 .123 .739 .089 .037 .012 1680. 81 .683 14.377 .303 3. 637 7.588+06 8.777+08 .123 .739 .089 .037 .012

1740. 97 .678 2.322 -55 .416 -44. 732 0.000 8.777+08 ,230 .770 .000 .000 .000

1800. 98,.014 1.986 -64 .641 -35. 282 0.000 8.777+08 ,230 .770 .000 .000 .000 1860. 98,.094 1 .906 -68 .298 -31 . 506 0.000 8.777+08 .230 .770 .000 .000 .000 1920. 98 .127 1.873 -70 .817 -29. 06 6 0. 000 8.777+08 .230 .770 .000 ,000 .000

37

mammmaBom iaBB]
APPENDIX B — PROGRAM LISTING

38

.
cC COMPF2 MAIN PROGRAM

C CONSTITUTES REVISION OF PROGRAM COMPF

C COMPF2 VERSION 1.0 PROGRAMMED 1 MARCH 1978 BY V. BABRAUSKAS

C VERSION 1,1 MINOR REVISIONS 18 AUG 1978

C

COMMON /CNSTS/ AW ALLNt BWOOW, DENS A. G , GASCNT . KTRACE ,MTI ME

COMMON /CP/ CPA.CPC0(2) tCPC02(2} tCPH2C2) .CPH20(2) .CPN2(2)

t

1 CP02(2) .CPPYR(2)

COMMON /FUEL/ C.CFLPC. CVGROS. CVNET t H, HFLPC. MWPYR . Nt NFLPC. Ot

1 OFLPC.Rt ROtREGRES.SH,SHAPE,SIZE,W,WFLPC, WTFUEL

COMMON /GP/ AWDOWiBPFtCDtCNV.DTIME,EMS( 2) .HWDOWtIX.IXC» IXLt

1 J» JM.JPtJPRI NT tK,KD,KH, KIT ER,KNTRL,MWlN,MWOUT«RC»RP .SIGMA

COMMON /LOGIC/ FC.FLSPEC . KRI T. NE WPLT . NEWPRP.PLFUEL. PLOT . PNCH»

1 RPSPEC. VTSPEC

COMMON /PLAST/ TBOILC .DHP , STOICH, SI ZE 1 . EITA . EISCAN

COMMON /PRBLM/ ADI A. AFLOOR. AW ALL . DENS W. FLO ADf IRUN. OPENF,

1 PRNTiSTEADY.THICKW

COMMON /TEMP/ OENF, DENU» T AMB»TGAS,TI NPT «T 1 ( 20 ),»T 2( 20 ) . TSF ,TSU

COMMON /THERML/ CNDA( 2. 1 ) . CPW ( 2. 1 0) » DX. EF, EMSA ( 2 . 10 )

.

1 NCNOtNCPWtNEMStNQGEN«NRP.QGEN( 2t 10) «RPX( 2 f 50

)

COMMON /TITLE/ TITLE(14)

COMMON /CPLOT/ BUFX (500) . BUFY( 500 ), SCALX .SCALY . SPECS ( 30)

INTEGER TITLE

LOGICAL AOIA.EISCAN.FCtFLSPECiKRIT.NEWPLT.NEWPRP. PLFUEL.

1 PL OT.PNCH, RPSPEC. STEADY, STOICH, VTSPEC

REAL MWIN.MWOUT.MWPYR.MTIME.N.NFLPC

C

DATA ADIA.CD.CFLPC»CNV»CPA»DENSA.EF.G

1 /.FALSE. (0.6 8. 44.

4

.5.0. 1005..1.18.0.9.9.8/*

2 HFLPC.IRUN.IX.MTIME.MWPYR.NEWPRP

3 /5. 4.0. 10.360. .28.97, .TRUE./.

4 NFLPCOFLPC. PLOT. PNCH.REGRES.SH. SHAPE

5 /O.O, 38.2, .FALSE., .FALSE. ,0.0, 0.0. 2./.

6 SIGMA, SIZEI, STEADY. TINPT.WFLPC

7 /5.6697E-8.-10.. .FALSE.. 0. 0. 12.0/

C HEAT CAPACITIES ARE GIVEN IN THE FORM CP( 1) TEMP+CPC 2)

DATA CPCO /0.1185, 1018./.

1 CPC02 /0.2114. 931./.

2 CPH20 /0.3E49. 1814./,

3 CPH2 /0.68e2, 13966./,

4 CP02 /0.3704. 931./.

5 CPN2 /0.1127, 1010./.

6 CPPYR /0.1127. 1010./

C PROPERTIES OF PYROLYSIS GASES ARE ASSUMED SAME AS FOR NITROGEN

C

NA MEL I ST /VARS/ ADI A. AFLOOR, AWALL . AWDQW , BPF .CD, CFLPC , CPPYR , C VGROS

.

1 OENSW.DHP.DTIME.EF.EISCAN. EITA.FLOAD.FLSPEC.HFLPCHWOOW.

2 I RUN. I X,KTRACE,MT I ME, MWPYR.NEWPLT.NEWPRP. NFLPCOFLPC, PLFUEL.

3 PLOT, PNCH,PRNT,REGRES,RPSP EC, SH, SHAPE, SIZE, STEADY, STOICH,

4 TINPT.THICKW,TBOILC, VTSPEC, WFLPC, SCALX , SCALY

10 READ ( 1,900, END=150) TITLE

900 FORMAT (13A6,A2)

WRITE(2,910) TITLE

910 FORMAT( IHl, 13A6,A2)

EISCAN=. FALSE.

STOICH=. FALSE.

39

"''^"''"''''""llli I IHMMiri f l —
STEADY=. FALSE.

KITER=

TINPT= 0.

KNTRL= 1 KTRACE=

NEWPLT= .FALSE.

RPSPEC= .FALSE.

FLSPEC= .FALSE.

VTSPEC= .FALSE.

PLFUEL= .FALSE.

IRUN= IRUN+1 lljj 20 READ (l.VARS)

WRITE (SsVARS)

IF(ADIA.OR.EISCAN.OR.STOICH) STEAOY= .TRUE.

IF (PNCH) PUNCH 900. TITLE

IF (NEWPLT. AND. PLOT) CALL PLTRST

IF (NEWPRP) CALL INC

NEWPRP= .FALSE.

30 CALL ICONOS

IF (KNTRL.EQ.2) GOTO 10

IF ( KT RACE. NE. 1 ) GOTO 50

40 WRITE CAt901} IRUN

901 FORMAT {•! RUN N0.«.I4//

1 • TGASl TGAS2 Fl F2 DERIVl K KO • 2 'KH J T2(l) TSF QFIRE QFLOW QRAOW • 3 • RP RC •/)

50 IF (KITER.EQ.I) GO TO 60

CALL ECHOIO

60 IF (RPSPECJ GOTO 70

IF (FLSPEC) GOTO 80

IF (VTSPECJ GOTO 90

IF CPLFUEL. AND.. NOT. FLSPEC) GOTO 100

CALL CRIB

GOTO 110

70 CALL RPFIX

GOTO 110

80 CALL PFLFIX

GOT 110

90 CALL PVTFIX \ GOT O 110

100 CALL POOL

110 GO TO ( 120* 10.130«120)«KNTRL

C KNTRL= 1 INITIAL VALUE

C KNTRL= 2 INPUT ERROR DETECTED, PROCEED TO NEW RUN

C KNTRL= 3 ITERATION FAILURE—PRINT OUT STEPS EVEN IF KTRACE= O.

C KNTRL= 4 SIMULATION TIME LIMIT EXCEEDED

120 CONTINUE

C INSERT HERE ANY REWIND COMMAND TO BE DONE IF NO ERROR

IF (PLOT-OR.PNCH) CALL DOUT

GO TO 10

130 IF (KTRACE.EQ.l ) GOTO 10

140 KTRACE= 1 KITER= 1 WRITE (2<.d03) TGAS

903 FORMAT ( 'O—ITERATI ON FAILURE— •/• TGAS=« .F16 .2 ) GO TO 30

150 CONTINUE

C ENDFILE 5 STOP

END

40
SUBROUTINE CRIB

CRIB FIRE ROUTINE

EQUATIONS FOLLOW NILSSON'S DATA FOR WOOD CRIBS.

OTHER FUEL CRIBS CAN BE TREATED IF PYROLYSIS CONSTANTS

ARE KNOWN.

COMMON /CNSTS/ AWALLNt BW DOW, DENSA, G . GASCNT . KTRACE ,MTI ME

COMMON /CP/ CPA,CPC0(2),CPC02(2) tCPH2(2) tCPH20(2) fCPN2(2)

t

1 CP02(2).CPPYR(2J

COMMON /FUEL/ C t CFLPC tCVGROS.CVNET t H, HFLPC . MWPYR t N. NFLPC, O,

1 OFLPC.R.R0tREGRES,SH»SHAPE,SIZE, W.WFLPC. WTFUEL

COMMON /GP/ AWDOWtBPF,CD,CNV,OTlME,EMS( 2) .HWDOWtlX.lXC. iXLt

1 Jt JMtJPtJPRI NT* KtKDtKH.K ITER, KNTRL,MWlN«MWOUT*RC«RP« SIGMA

COMMON /LOGIC/ FC,FLSPEC,KRIT,NEWPLT,NEWPRP,PLFUELtPLOT.PNCH,

1 RPSPECVTSPEC

COMMON /PLAST/ TBOILCDHP , STOI GH, SI ZE 1 ,EITA ,EISCAN

COMMON /PRBLM/ ADI A, AFLOOR, AWALL , DENSW.FLOAD.I RUN .OPENF

,

1 PRNT, STEADY, THICKW

COMMON /OS/ QC0NWtQFlRE«QFLOW*QRADO,QRADW,QWLSUM

COMMON /TEMP/ DENF, DENU»T AMB, TGAS.TINPT ,T 1 ( 20 > , T2( 20 ) »TSF ,T SU

COMMON /THERML/ CNDAC 2, I ) , CPW( 2, 1 0) ,DX . EF, EMSA( 2 . 1 0)

.

1 NCND,NCPW,NEMS,NQGEN,NRP.QGEN(2, 10) tRPX(2.50)

COMMON /WOUT/ BWORST,FLREM, HRATIO,RMA,RMF,TTIME»VAVGlN.

1 WA, WB,YC02. YH20, YN2, Y02,YPYR

LOGICAL ADIA,EISCAN,FC,FLSPEC,KRIT,NEWPRP,PLFUEL,

1 PLOT. PNCH,RPSPEC, SCAN, STEADY, STOI CH.VTSP EC

REAL MWIN,MW0UT,MWPYR,MTIME,N,NFLPC

IF (STEADY) GOTO 190

FC= .FALSE.

SCAN= .FALSE.

QRAOW=0.

QCONW= 0.

F2 = 0.

F1=0.

DTGAS=10.

CALL HEAONG

START TIME LOOP

DO 170 J=1,JM

KH=

0ERIV1= 1. TGAS2= 0.

TGAS1= 0.

TGASP= 2000.

TGASN= TAMB

20 CONTINUE

K=

30 CONTINUE

IF (FLREM.GT.O.) GOTO 40

RC= 0.

RP= 0.

GO TO 50

40 IF (REGRES.LE.0.0) GOTO 45

USE THIS FORMULA IF INPUT REGRES IS SPECIFIED

RP= REGRES*2.*SHAPE/SIZE*FLREM*«( I .- 1 ./SHAPE )*WTFUEL**( 1. /SHAPE)

GO TO 50

45 CONTINUE

41

'^^niHH!aUHIlBIIUIMIUMMIIMH«l .,MI.I-.l»l«Li,..IHMI..- i^HHHIld
<; FUEL SURFACE CONTKOL

C ASSUME CRIB STICK DENSITY RHOCR= 500 KG/M**3

RHOCR= 500.

REGREN= l.24E-3/RHOCR*SIZE**-0.6

RP1= RE6REN*2.*SHAPE/SIZE*FLREM**( l.-l ./SHAPE )*WTFUEL** ( 1 ./SHAPE)

C; CRIB POROSITY CONTROL

RP2= 0.22«WTFUEL/(RHOCR*SIZE)*SH

C ROOM VENTILATION CONTROL

RP3= 0.120*AWOOW*SQRT(HWDOW)

RP= AMINI (RPl,RP2.RP3)

50 RMF= RMA-I-RP YC02= 3.66667*CFLPC*RC/100./RMF

YH20= ( WFL PC*RP+9 • 0* HFL PC*RC ) / 1 . /RMF

Y02= (0.23*RMA-R0«RC)/RMF

YN2= 0.77*RMA/RMF +NFLPC*RP/1 00 ./RMF

YPYR= (RP-RC)/RMF

IF(YPYR.LT..0> YPYR= 0.

MWOUT= 44.«YC02+1 8.*YH20+28.*YN2+32 .*Y02+MWPYR*YPYR

HRATIO= l./( l.+( (TGAS/TAMB>*(MWlN/MWOUT)*( I . +RP/RMA)**2}

1 **0. 3333333333)

C NOTE HIN IS TAKEN AS POSITIVE

HIN= HWOOW* HRATIO

ZW= I .-MWOUT*TAMB/MWI N/TGAS

IF(ZW)195i55.55

55 VAVGIN= 0.666667*SQRT(2.*G*HIN*ZW)

RMA= C0*VAVG1N*H1N*BW00W*DENSA

RMF= RMA+RP

IF (RMA/R-RP) 60.60.65

60 RC= BPF*RMA/R

GO TO 70

65 RC= BPF*RP

FC= .TRUE.

70 CONTINUE

QFLOW= RMF*(YCO2*(TGAS*(0.5*CPCO2( i ) *TGAS+CPC02{ 2 ) )-TAMB*(0.5*

1 CPC02( 1 )*TAMB+CPC02(2)) ) +YH20*( TGAS*( . 5*CPH20( 1 ) *T6AS+

2 CPH20(2) )-'TAMB*(0.5*CPH20{ 1 )*TAMB+CPH20(2) )) +Y02*(TGAS*(

3 0.5*CP02< 1 )*TGAS+CP02( 2) )-T AMB* ( .5*CP02{ 1 )*T AMB+CP02 ( 2) )

)

4 +YN2*(TGAS*(0.5*CPN2(1 )*TG AS+CPN2( 2 ) )-TAMB*( 0. 5*CPN2( 1 )*

5 TAMB+CPN2(2) ) ) +YPYR* ( T6AS*( .5*CPPYR( 1 ) *TGAS+CPPYR( 2)

)

6 -TAMB*(0.5*CPPYR(1 )*TAMB+CPPYR( 2 ) )) ) QFIRE= RC*CVNET

IF (ADIA) GOTO 90

CALL DESOLV

QRADW= AWALLN*EMS( I )*SI GMA*(T6AS**4.-TSF**4 .

)

QCONW= AWALLN*(TGAS-=TSF)*CNV*( (TGAS-TSF) *(T GAS-TSF )) **0. 1 666666 7 90 CONTINUE

QRADO= AWDOW*SIGMA*(TGAS**4.-TAMB**4.)

K= K + l F3=F2

F2 = Fl

Fl= QFIRE-=QFL0W-QRADO-QRADW-QC0NW

TGAS3=TGAS2

TGAS2=T6AS1

TGASi=TGAS

IF (F1.LT.0..ANO.TGAS.LT.TGASP) TGASP=TGAS

IF (Fl.GT.O. .AND.TGAS.GT.TGASN) TGASN=TGAS

DERIV2= DERIVI

IF (TGAS1.EQ.TGAS2) GOTO 130

42
DER1VI=(FI-F2)/{TGAS1-TGAS2)

IF (KTRACE.GT.O) WRITE (4*99) TGAS 1 »TGAS2.F 1 ,F2 .DERI VI ,K .KD

.

1 KH. J.T2( 1 ).TSF»QFIRE,QFLOW,QRAOW,RP,RC

99 F0RMAT(2F9.2.3( I PE9. 2) ,313. 15 , 2( 0PF9.2) . 3 ( IPEIO .3) ,2{0PF7,3))

IF {.NOT. SCAN) GOTO 95

IF (F1/F2.GE.0.0) GOTO 93

SCAN= .FALSE.

GOTO 100

93 TGAS= TGAS-DTGAS

IF (TGAS.LT.TAMB) GOTO 200

GOTO 120

95 IF (DERIVl .LT..0.ANO.ABS(F2) .GT..0001)GOTO 100

IF(OERIV2.LT-.0.AND.J.GT.2) GOTO 100

TGAS= TGASl+DTGAS/5.

GOTO 120

100 DIF= ABS(F1/QFL0W)

IF (OIF.LT.0.002.AND.ABS(TGAS2-TGAS1) .LT.2. ) GOTO 130

TGAS = (F1*TGAS2-F2*TGAS1 )/(Fl-=F2)

IF (K.GT.10.ANO.FI.LT.0..AND.TGAS.GT.TGASP) GOTO 105

IF (K.GT.lO.AND.Fl .GT.O.. AND.TGAS.lt .TGASN) GOTO 105

IF (K.EQ.l.AND.KH.EQ.O) TGAS= TGASl +10.

IF (TGAS.GT.2000.) GOTO 110

IF (TGAS.LT.(TAMB+30.)) GOTO 110

GOTO 120

105 TGAS= (TGASN+TGASP)/2.

GOTO 120

110 SCAN= .TRUE.

TGAS= 1900.

120 CONTINUE

IF (K-200) 30t30»200

130 CONTINUE

CALL RSTA

FLREM= FLREM-RP*DTIME

IFCFLREM.LT.O.) FLREM=0.

IF (QCONW.GT.O. ) QWLSUM= QWLSUM+( QRADW+QCONW )*DTI ME

IF (TTIME .GE. MTIME) GO TO 210

IF (TGAS.LE.353..ANO.J.GE.10) GO TO 210

IF ( J.EQ.l ) GO TO 150

IF ( JP.LT. JPRINT) GO TO 160

JP=

150 CALL OUTPUT

160 JP= JP+1

TTIME= TTIME+DTIME

170 CONTINUE

END TIME STEP DO-LOOP

180 CONTINUE

185 CALL OUTPUT

RETURN

ERROR IN INPUT

190 CONTINUE

KNTRL= 2 WRITE (2.910)

910 FORMAT (///• CRIB ROUTINE DOES NOT ACCEPT STEADY-STATE CAS£*)

RETURN

SQUARE ROOT ERROR

195 CONTINUE

IF(KTRACE.EQ.l) WRITE(2.930) TGAS.RC. RP. YPYR . ZW .RMA.MWOUT

930 FORMAT (/• TGAS= • . F5.0. • RC=«.E10.4.' RP=».E10.4.

43

Bsa^muB mmsmmsmaamaa^
I • YPYR=« ,E10.A.« ZW=«,F6.4.' RMA=« ,E1 0.4, • MWOUT=« , F6.1

)

C FAIL TO CONVERGE. ERROR EXIT

20 CONTINUE

KNTRL=3

RETURN

C FIRE IS OVER (TRANSIENT CASE)

210 CONTINUE

CALL OUTPUT

RETURN

END

44

HP
SUBROUTINE DEQNS

DIFFERENTIAL EQUATION SOLVER BASED ON CRANK-NICOLSON METHOD.

COMMON /GP/ AWDOW,BPFtCDiCNV,DTlME.EMS(2}»HWDOWtIX»IXC.IXLf

1 J. JM.JPt JPRINTtKtKD,KH,KlTERtKNTRLiMWlNtMWOUT.RCtRP»SIGMA

COMMON /PRBLM/ AOIAt AFLOOR. AWALL t DENSWtFLO AD»I RUN, OPENF,

1 PRNTtSTEADY.THlCKW

COMMON /TEMP/ OENF » DENU.TAMB.TGAS ,TI NPT tT I ( 20) ,T2( 20 ) ,TSF .TSU

COMMON /THERML/ CNDA (2» 1 > . CPW( 2, 1 0) , DX, EF, EMSA( 2 . 10)

.

1 NCNO.NCPWtNEMStNQGEN*NRP,QGEN(2* 10) *RPX( 2 « 50

)

DIMENSION A (20) • B( 20)

t

C(20 ) • CND( 20) , D( 20 ) • HCP(20

)

ENTRY RSTA

ENTER HERE WHEN READY FOR NEW TIME STEP (FINISHED ITERATING)

DO 10 I^lflX

TKI )= T2{I ) 10 CONTINUE

TGOLD= TGAS

IF (J.EQ.O) TGOLD=TAMB

RETURN

ENTRY DESOLV

SOLVE DIFFERENTIAL EQUATION

KD= 1 DX1= DX

DO 20 1=1. IX

CND(1)= TLU(CNDA.NCN0.T1(I) ) 20 HCP(I)= DEN£W«OX/OTlME*TLU(CPW,NCPW«Tl(I))

EMS(1)= l./( 1 ./TLU(EMSA,NEMS.TSF) +1 ./EF -1.)

EMS(2)= TLU(EMSA,NEMS.TSU)

DO 50 I=2tlXL

CNL= l./(DX/CND(I-l)+DX/CND( I)

)

IF (1.EQ.2 ) CN_= l./(DXl/CND( 1)+0X/CND(2))

CNG= 1 ./(DX/CND(I )+DX/CND(I+l)

)

A(I)= -CNL

B(I)= HCP( I )+CNL+CNG

C(I )= -CNG

0( I )= (HCPd )-CNL-CNG)*Tl (I)+CNL*Tl( I-l )+CNG*Tl ( I + l ) 50 IF (NQGEN.GT.O) D(I)= D( I ) +DX*TLU C QGEN.NQGEN.Tl ( I )

)

CNG= l./(DXl/CND( 1)+DX/CND( 2)

)

C( 1 )= -CNG

CNL= 1./(0X/CND(IXL)+DX/CND(IX)

)

A(IX)= -CNL

ENTER HERE WHEN KD.GT.l SINCE PRIOR EXPRESSIONS DO NOT CHANGE

30 TSFOLD= TSF

ZRF= TGAS*CTGAS*(TGAS+TSF)+TSF«TSF)+TSF*TSF«TSF

ZCF= CNV*( (TGAS-TSF)*(TGAS-TSF))**0. 1666666 667

HF= ZCF+SIGMA*EMS( 1)*ZRF

DENF= HF»OXl/2./CND( 1)

ZF= HF/2./(DENF+l.)

45

iHurai
B(l)= HCP( 1 )+ZF+CNG

D( 1)= (HCP( n-ZF-CNG)«Tl( 1 )+ZF*( TGAS-I-TGOLD) +CNG*T 1 C 2>

IF (NQGEN.GT.O) 0(l)= D( I ) +DX*TLU( QGEN. NQGEN .Tl ( 1 )

)

ZRU= TAMB«{TAMB*(TAMB+TSU)+TSU*TSU)+TSU*TSU*TSU

ZCU= 1 .87*( (TAMB-TSU)«(TAMB-TSU)>**0. 1666666667

HU= ZCU+SIGMA*EMS{2)*ZRU

DENU= HU*DX/2./CND(IX}

ZU= HU/2«/(DeNU+l«)

B(IX)= HCPCIX)+ZU+CNL

D(IXJ= {HCP(IX)-ZU-CNL)*T1C1X)+ZU*2.*TAMB+CNL«T1(IXL)

IF (NQGEN.GT.O) D{IX)= D( I X > + DX*TLU( QGEN . NQGEN.T 1 ( I X)

)

CALL TRIOGF ( A* B» C. 0tT2

t

IX)

TSF= (DENF«TGAS+T2< 1))/(DENF+1 .)

TSU= (0ENU*TAMB+T2(IX) )/(DENU+l •> KD= KD+1

IF {ABS(TSF-TSF0LD).LT.4) RETURN

IF (KD.LE.6) GO TO 30

WRITE (2»100) TSF.TSFOLD

100 FORMAT CO FAIL TO CONVERGE D.E. TSF=«,F7.2. • TSFOLO=« F7.2)

IF (XD.LE.aO) GO TO 30

RETURN

END

46

II
SUBROUTINE ECHOlD

SUBROUTINE TO ECHO INPUT DATA

COMMON /CP/ CPA,CPC0<2)tCPC02(2) .CPH2(2) .CPH20(2> tCPN2(2)t

CP02(2} tCPPYR(2J

COMMON /FUEL/ C t CFLPC t C VGROS. CVNET tH .HFLPC . MWPYR.N* NFLPC* O.

OFLPC.R.RO.REGREStSH, SHAPEtSlZE.W.WFLPC.WTFUEL

COMMON /GP/ AWDOW.BPFtCD,CNV,DTlME.EMS( 2).HWD0W,lX.lXCt iXLt

J. JM, JPt JPRINTtK,KOtKH,KlTER,KNTRLiMWIN»MWOUT«RC.RP«SIGMA

COMMON /LOGIC/ FC.FLSPEC t KRI Tt NE WPLT t NEWPRP ,PLFUEL» PLOT » PNCH

.

RPSPEC.VTSPEC

COMMON /PLAST/ T BOlLC.OHPt STOlCHt SI ZEl t EIT A t EISCAN

COMMON /PRBLM/ AOIA. AFLOOR tA WALL . DENSW tFLO AOt I RUNt OPENFt

PRNT.STEADYtTHICKW

COMMON /THERML/ CNDA( 2t 1 ) t CPW( 2t I 0) f DX . EF» EMSA( 2 1 1 0) i NCN0*NCPW«NEMStNQGENtNRP.QGEN(2t 1 ) • RPX( 2 « 50

)

LOGICAL AOIA. EI SCAN. FC.FLSPECt KRl T« NEWPRPiPLFUELt

PL OT.PNCH.RPSPEC. STEADY, ST OICH.VT SPEC

REAL MWIN. MWOUT.MWPYRtN.NFLPC

WRITE (2.90) IRUN

WRITE (2*91) AWALL* AFLOORt HWOOW* AWOOW. OPENF. CO

WRITE (2.92J FLOAD

IF (PLFUEL) WRITE (2,93) DHP,TBOILC

WRITE (2,94) CFLPC,HFLPC,OFLPC,NFLPC,WFLPC,R,R0

WRITE (2,95) CVGROS,CVNET»MWPYR,CPPYR.BPF,EF

IF (.NOT. (PLFUEL. OR. VTSPEC. OR. RPSPEO) WRITE (2,96) REGRES.SIZE,

I SHAPE

IF (PLFUEL. AND. .NOT. STOICH) WRITE (2,97) SIZE

IF ( .NOT. (PLFUEL. OR. VTSPEC. OR. RPSPEO. AND. SH.GT. 0.0) WRITE (2,908)

1 SH

IF (RPSPEC.AND.NRP.EQ.l) WRITE (2,913) RPX(2,1)

IF (RPSPEC.AND.NRP.NE.l) WRITE (2,91A) ( ( RPX ( I , J ) , 1= 1 , 2 ) , J= 1 , NRP)

IF(ADIA) GOTO 200

WRITE (2,98) THICKW,DENSW

IF (NCND.EQ.l) WRITE (2,900) CNDA(2,1)

IF (NCND.GT.l) WRITE ( 2,90 1 ) ( (CNDA( I , J) . I >1 , 2 ) , J=l , NCNO

)

IF (NCPW.EQ.l) WRITE (2,904) CPW(2,l)

IF (NCPW.GT.l) WRITE ( 2, 905) ( (CPW( I , J ) , I =1 . 2 ) . J-1 , NCPW

)

IF (NEMS.EQ.l) WRITE (2,902) EMSA(2,1)

IF (NEMS.NE.l) WRITE ( 2, 903) ( ( EMSA( I , J ) , 1 = 1 ,2 ) . J=l , NEMS

)

IF (NQGEN.EQ.l) WRITE (2,906) QGEN(2,1)

IF (NQGEN.GT.l) WRITE ( 2, 907) ( ( QGEN( I , J) , 1 = 1 ,2) . J = l . NQGEN)

200 IF(AOIA) WRITE(2,909)

RETURN

90 FORMAT (1H+, T86,«C0MPF2 VERSION 1.1 - RUN N0.»,I4)

91 FORMAT ('0—• GEOMETRY AND VENTILATION----*//

1 • WALL SURFACE AREA = • , F8 . 1 , • M2«/

2 • FLOOR AREA = •F8.2.« M2 •

/

3 • WINDOW HEIGHT = •,F8.2,« M«/

4 • AREA = •,F8.2.« M2«/

5 • OPENING FACTOR = •,F7.3, • M2.5«/

6 • DISCHARGE COEFF.= •,F4.2/)

92 FORMAT ('0 -FUEL LOAD PROPERTIES —*//

1 • FIRE LOAD PER FLOOR AREA =«,F6.1,« KG/M2M

913 FORMAT ('ORATE OF PYROLYSIS =',F7.2,' KG/S«)

914 FORMAT ( 'ORATE OF PYROLYSIS (KG/S)'/' TIME RP'

47

i^H HHHIHHii
I 50(/3XtF5.0.F9.3)

)

93 FORMAT (• TOTAL ENTHALPY OF PYROLYSIS= • 1 .lPE10.2t« J/KG«/« BOILING TEMPERATURE= • , 0PF5 . . • DEG C«)

9A FORMAT COFUEL COMPOSITION'/

PERCENT BY WEIGHT'/

PERCENT'/

PERCENT'/

PERCENT'

/

PERCENT'

/

1 • CARBON = '.F4.1

2 • HYDROGEN = ' . FA .

I

3 • OXYGEN = '.F4,i.

4 ' NITROGEN = ' .FA.

I

5 ' WATER = ' ,FA.l

6 • R = '.FS.a/

7 • R0= '.F5.2)

95 FORMAT ( • HEAT OF COMBUSTION OF DRY FUEL = '.2PE10.2,' J/KG'/

1 • LOWER ACTUAL HEAT OF COMBUSTION = ',E10.2,' J/KG'/

2 ' MOLECULAR WEIGHT OF UNBURNT PYROLYSATES = •.0PF6.2/

3 ' CP OF PYROLYSIS GAS = ( • .F6. A. ' TGAS + '.F6.0.

4 •) J/KG-K'/

5 • MAXIMUM FRACTION OF PYROLYSATES BURNED = •.F5.2/

6 ' GREY-GAS FLAME EMISSIVITY = •F5.3)

96 FORMAT (• RATE OF REGRESSION = ',2PE9.2,' M/S'/

1 • FUEL DIMENSION = '.0PF5.3.' M'/

2 ' SHAPE FACTOR = '•F4.2 /)

97 FORMAT (• FUEL AREA= • . F 1 . 2» ' M2'

)

908 FORMAT (' CRIB SPACING/HEIGHT RAT IO= ' . F6.3)

98 FORMAT ('0 --WALL THERMAL PROPERTIES— —•//

1 • THICKNESS = '.F5.3,' M'/

2 ' DENSITY = '.F6.0»' KG/M3')

909 FORMAT CO——WALL THERMAL PROPERTIES '//

I ' ADIABATIC WALL'//J

900 FORMAT C'OTHERMAL CONDUCTIVITY = '.F7,3.' W/M-K'}

901 FORMAT COTHERMAL CONDUCTIVITY ARRAY (W/M-K)'/

1 ' TEMPERATURE CONDUCTIVITY'. 1 ( /3X . F7 . 1 .4X. F 1 0.3 >

)

902 FORMAT COEMISSIVITY = '.FA. 2)

903 FORMAT COEMISSIVITY ARRAY'/' TEMPERATURE EMISSIVITY'

1 10(/3X.F7.1,AX,F10.3))

90A FORMAT COHEAT CAPACITY = '.F7.0t' J/KG-K'J

905 FORMAT ('OHEAT CAPACITY ARRAY (J/KG-K)'/

1 ' TEMPERATURE HEAT CAPACITY'. I 0( /3X, F7. 1 . AX . FIO . )

)

906 FORMAT COWALL HEAT GENERATED = ',F9.3.' W/M3')

907 FORMAT ('OWALL HEAT GENERATED ARRAY (W/M3)'/

1 ' TEMPERATURE QGEN'. 10 ( /3X.F7. 1 . 4X.F i 0.3 )

)

END

48

k
SUBROUTINE ICONDS

SET INITIAL CONDITIONS AND CONSTANTS

COMMON /CNSTS/ AWALLN, BWDOW, DENSA , G , GASCNT

COMMON /CP/ CPA,CPC0(2) .CPC02(2) .CPH2(2

1 CP02(2l .CPPYR(2)

COMMON /FUEL/ C .CFLPC , CVGR0S» CVNET.H.HFLP

1 OFLPCR.RO.REGRES.SH.SHAPE.SIZE.Wt WFLPC.

COMMON /GP/ AWDOW,BPF,CO,CNV,DTIME,EMS(

1 Jt JM, JP. JPRINT.K.KD.KH.KITER.KNTRL.MWIN.

COMMON /LOGIC/ FC . FLSPEC. KRI T, NE WPLT. NEWPR

1 RPSPEC. VTSPEC

COMMON /PLAST/ TBOILC .DHP , STOICH, SI ZE 1 , EI TA

COMMON /PRBLM/ AOI A, AFLOOR. AWALL , DENSW .FLO

1 PRNT, STEADY. THICKW

COMMON /QS/ QC0NW»QFIRE.QFLOW.QRAD0.QRA

COMMON /TEMP/ DENF» DENU. TAMB. TGAS. T INPT ,T

COMMON /THERML/ CNDA ( 2. 1 ) . CPW( 2. 1 ) . DX, EF.

1 NCN0«NCPW.NEMS«NQ6EN,NRP.QGEN(2. 10) «RPX(

COMMON /WOUT/ BWORST. FLREM, HRAT I O. RMA, RMF

1 VKA,WB.YC02.YH20fYN2. Y02.YPYR

LOGICAL ADIA,EISCAN,FC,FLSPEC.KRIT.NEWPRP,P

1 PLOT, PNCH, RPSPEC. STEADY. STOICH

REAL MWIN.MWOUT.MWPYR.MTI ME, N, NFL PC

FC= .FALSE.

KNTRL= 1 AWALLN= AWALL-AWDOW

BWOOW= AWDOW/HWDOW

WTFUEL= FLOAD*AFLOOR

OPENF= AWOOW*SORT(HWOOW)

MW1N= 28.97

MWOUT= MWIN

TAMB= 298.

TGAS= 1800.

IF (TINPT.GT.O.

)

IXL= IX-1

IXC= IX/2

TSF= TAMB

TSU= TAMB

DENF= 0.

DENU= 0.

WA= 6H

WB= 5H

IF < .NOT. FLSPEC)

WA= 6HWIN00W

WB= 5HWI0TH

10 JP=

IF ( .NOT. STEADY)

IF (STEADY) GOTO

IF (DTIME.GT. 0.00001) GOTO 15

WRITE (2.95)

KNTRL= 2 RETURN

95 FORMAT (///• FOR NON-STEADY PROBLEMS MUST S

1 'GREATER THAN ZERO')

15 IF (MTIME.GT.DTIME) GOTO 20

.KTRACE.MTIME

) .CPH20(2),CPN2(2)

.

C.MWPYR.N.NFLPC. O. WTFUEL

2) .HWDOW, IX. IXC. I XL.

MWOUT.RC.RP, SIGMA

P. PL FUEL. PLOT, PNCH.

,EISCAN

AD.IRUN.OPENF,

DW.QWLSUM

1(20). T2(20) .TSF.TSU

EMSA(2. 10)

.

2.50)

.TTIME.VAVGIN,

LFUEL.

TGAS= TINPT

GO TO 10

JPRINT= PRNT/DTIME + 20

(l.O-l.E-6)

PECIFY DTIME '.

49

""TM"TM""*TMTM HI ^gg^m
MRITE (2.92)

92 FORMAT ( /// • FOR NON-STEADY PROBLEMS MUST SPECIFY MTIME ••

1 'GREATER THAN DTIME*

)

KNTRL= 2

RETURN

20 CONTINUE

IF < .NOT, STEADY) JM= MTIME/DTIME +2

IF (STEADY) JM= 2

DX= THICKW/IX

FLREM= WTFUEL ^ TTIME= 0.

EMS(1)= 1 ./( 1 ./TLU(EMSA,NEMS, TAMB) +1./EF-1.)

QWLSUM= 0.

DO 60 1=1, IX

60 T2(I)= TAMB

J=

CALL RSTA

TOTAL= CFLPC+HFLPC+OFLPC+NFLPC+WFLPC

IF (TOTAL.LT. 101.1 .AND. TOTAL. GT. 98. 9) GOTO 70

CHECK FOR ERRORS IN FUEL COMPOSITION

KNTRL= 2

WRITE (2.90)

90 FORMAT (///• SUM OF FUEL COMPOSITION INPUT IS INCORRECT')

RETURN

70 C= CFLPC*( 10./12.)

H= HFLPC+10.

0= OFLPC*( 10./16.)

W= WFLPC*( lO./ia.

)

N= NFLPC*( 10./14.

R0= (C+H/4.-0/2. )*32./1000.

R= RO/0.232

CVNET= CVGROS*( 1 .-WFLPC/ 1 00 . )-( WFLPC+9. 0*HFLPC) /I 00.*2440 .E+3

LATENT HEAT OF H20 EVAPORATION= 2440E+3 J/KG AT 25 C

RMA= 0. 16*AWDOW*OENSA*SQRT(G*HWDOW)

RMF= RMA

RC= BPF*RMA/R

Y02= 0.10

IF (STOICH) EITA= 1.

IF (EISCAN) SIZE= SIZEl/EITA ^ RETURN

END

50

I^P
SUBROUTINE INC

ROUTINE TO INPUT ALL CONSTANTS NOT GIVEN IN NAMELIST -VARSCOMMON /THERML/ CNDA( 2t I ) t CPW( 2 . 10) .DX EFt EMSA( 2 . 1 0)

.

1 NCN0,NCPW.NEMStNQGEN,NRP«QGEN(2. 10) .RPX(2t50)

NQGEN=

READ (lt91) NCNtNCPt NEMtNR«NQG

IF (NCN.EQ.O) GO TO 10

NCND= NCN

READ (1,90) (CNDA( l.I),CNDA{2.I ),I=1 ,NCND)

10 IF (NCP.EQ.O) GO TO 20

NCPW= NCP

READ (1,90) (CPW( 1,1 ),CPW(2,I ),I=1,NCPW)

20 IF (NEM.EQ.O) GO TO 30

NEMS= NEM

READ (1,90) (EMSA( 1,I),EMSA(2,I ),I = 1,NEMS)

30 IF (NR.EQ.O) GO TO 40

NRP= NR

READ (1,90) (RPX( 1,I),RPX(2,I ),I = 1,NRP)

NQGEN= NQG

40 IF (NQG.EQ.O) GO TO 50

READ (1,90) (QGEN( 1,1) ,QGEN(2*I ),I=1 ,NQGEN)

50 RETURN

90 FORMAT (8F10.0)

91 FORMAT (1013)

END

51

"'°'"'°"*"TM'" maS!^a!mm,«Mm,w,l«,m«mm,,,,«,«^
SUBROUTINE OUTPUT

PRINTS OUTPUT DATA

COMMON /FUEL/ C t CFLPC. C VGROS iCVNET . Ht HFLPC t MWPYRt N. NFLPC. Ot

1 OFLPC.R,R0tREGRES»SH,SHAPE.SIZEtW,WFLPC.WTFUEL

COMMON /GP/ AWDOW.BPFtCOtCNV,DTlMEtEMS( 2) .HWDOWtlX.IXC. IXL.

1 Jt JMtJPtJPRI NT .RtKOtKH, KIT ER»KNTRL.MWlN,MWOUT.RC.RP» SIGMA

COMMON /LOGIC/ FCt FLSPEC. KRITiNE WPLT» NEWPRP iPLFUEL.PLOT t PNCHt

1 RPSPEC.VTSPEC

COMMON /PLAST/ TBOILC tOHP , STOI CH, SI 2E I .EIT A . EISCAN

COMMON /PRBLM/ ADI A, AFLOOR. AWALL . DENSW .FLO AO,I RUN» OPENF

I PRNT.STEADYtTHICKW

COMMON /OS/ QCONW.QFlREtQFLOW.QRAOOtQRADWtQWLSUM

COMMON /TEMP/ DENF , DENUtTAMB.TGAStTlNPT .T I ( 20) .T2( 20 J »TSF .TSU

COMMON /WOUT/ BWORST t FLREMt HRAT 10. RMA.RMF «TTI ME , VAVGI N

,

1 WA,WB.YC02.YH20.YN2.V02.YPYR

LOGICAL AOIA.EISCAN.FC.FLSPEC.KRIT.NEWPRP.PLFUEL.

1 PLOT. PNCH.RPSPEC. STEADY, ST O I CH.VT SPEC

REAL MWIN.MWOUT.MWPYR.N.NFLPC

LOGICAL DATPRT

DIMENSION T2C(3)

DATA DATPRT /.FALSE./

IF (KITER.EQ.I) RETURN

IF (PLOT.OR.PNCH) CALL DSTO (TTIME.TGAS)

IF (ILINE.LE.47) GO TO 50

DATPRT= .TRUE.

GO TO 300

50 TGASC= TGAS-273,

T2C(1)= TSF-273.

T2C(2)= T2(lXC)-273.

T2C{3)= TSU-273.

FUELPC= FLREM/WTFUEL*100.

QNORM= QFIRE/100.

IF (QCONW.LT.Oi QNORM= ( QFLOW + QRADO) /1 00.

ZFLOW= QFLOW/QNORM

ZRADO= QRADO/QNORM

ZCONW= QCONW/QNORM

ZRADW= QRADW/QNQRM

EXCE£S= RP-RC

ILINE = ILINE+1 \ WRITE (2.90) J. TTIME»TGASC.T2C.RP,RC. EXCESS. FUELPC.RMA.HRATIO,

1 VAVGIN, MWOUT.FC

IF (FLSPEC) WRITE (2.91) BWORST

WRITE (3.92) TTIME.ZFLOW.ZRADO.ZCONW.ZRADW.QFIRE.QWLSUM.

1 Y02. YN2. YC02.YH2O.YPYR

IF (STOICH) WRITE (2.901) SIZEl

RETURN

ENTRY HEADNG

START NEW PAGE

IF (KITER.EQ.I) RETURN

IPG = 1 300 CONTINUE

ILINE = IF { .NOT.RPSPEC) GOTO 315

WRITE (2.94) IPG.IRUN

52
315

325

335

(3«9A) IPG*IRUN

400

(. NOT.FL SPEC) GOTO 325

(.NOT.PLFUEL) WRITE (2

(.NOT.PLFUEL) WRITE (3

(PLFUEL) WRITE (2,905J

(PLFUEL) WRITE (3,905)

TO 400

(.NOT.VTSPEC) GOTO 335

WRITE (2<96) IPG.IRUN

WRITE (3t96) IPG.IRUN

GO TO 400

IF (.NOT.PLFUEL) GO TO

WRITE (2»97) IPG.IRUN

WRITE

GO TO

IF

IF

IFIF

IF

GO

IF

,95) IPG.IRUN

,95) IPG.IRUN

IPG.IRUN

IPG.IRUN

345

345

400

WRITE

GO TO

WRITE

WRITE

IPG = WRITE

WRITE

IPG.IRUN

IPG.IRUN

IPG.IRUN

WA, WB

90

(3.97)

400

(2,96)

(3,98)

lPG+1

(2,99)

(3,900)

IF ( .NOT.DATPRT) RETURN

OATPRT= .FALSE.

GO TO 50

FORMAT (IH ,S4,T6,F8.0,T15,F7.0,T25.3Fe.0.T44.F6.3.T51,F6.3,

I T6 0,F6.3,T69,F5.l .T79,F4. 2 , T87 ,F7. 2 ,T96 . F9.2 .T 1 09. F6 .2 . 5X , L I

)

91 FORMAT (1H+, T124, Fe.2)

92 FORMAT (IH , F7. 0, 4F1 1 .3, I PE 15 . 3, E 1 5. 3, IX , OP ( 5F8.3 )

)

94 FORMAT (1H1,T14,« INPUTTED VALUES OF RP ARE USED*

1 TIOO.'PAGE NO. •,I3,Tll5.'RUN NO.*, 14//)

95 FORMAT (1H1,T14, • FUEL PYROLYSIS (CRIB) SPECIFIED, VENTILATION •

1 • ADJUSTED FOR WORST CONDITIONS'. TIOO.'PAGE NO. • , 13 .T 1 1 5.

2 'RUN NO. ',14//)

96 FORMAT (1H1,T14.' VENTILATION SPECIFIED, FUEL PYROLYSIS ADJUSTED'

1 • FOR WORST CONDITIONS' .TIOO.'PAGE NO. ' . I3.T 1 1 5

, 'RUN NO. '.14/

2 /)

97 FORMAT (1HI,T14,' THERMOPLASTIC POOL FI RE • , T 100, ' PAGE NO.'

1 I3,T115,'RUN NO.', 14//)

98 FORMAT ( IH I ,T 14, 'CRIB Fl RE ' .T 1 00 .' PAGE NO. '.13.

I T115,'RUN NO. ',14/)

99 FORMAT ( IH0,T10, 'TIME' .T17. 'TEMP' .T30. 'WALL TEMPS ' .T47 .' RP'

,

1 T54, 'RC ,T59, 'EXC.PYR. ' .T70. 'FUEL' .T78. 'AIR I N ' .T90 .

' N.P .

'

2 T98, 'VELOCITY' ,T 1 09. ' MOL . WT ' .Tl 18, 'FUEL' .T124, A6/

3 T11,'S',T17, 'GAS.C .T34 . ' C ' .T46. • KG/S ' .T53.'KG/S' .T61 . 4 'KG/S' .T71 'PCT' .T79.'KG/S' .TlOl . 'M/S' .T108.

5 Tl IS.'CNTRL' ,T124,A5/)

900 FORMAT (T5. 'TIME' .T24. 'HEAT B ALANCE • ,T75, ' Q- WALL • , T88 , • Y02' ,T96»

2 •YN2',T10A, 'YC02' ,T 1 12, ' YH 20 ' ,T 120 ,

' YPYR •

/

2 T13,'GAS FLOW' ,T25, 'WND RAD' ,T36, ' WALL CNV ,T47 ,

' WALL RAD',

3 T6 0, 'Q-FIRE',T76, 'SUM' .T88.'PCT' ,T96, 'PCT' .T I 4, 'PCT' .Tl 12.

4 'PCT' .T 120. 'PCT'/

5 T16, 'PCT' ,T27, 'PCT' .T39. 'PCT' .T49, 'PCT' .T62.' W' ,T77,

6 ' J' ,T8a, 'MASS' ,T96. 'MASS' .T104. 'MASS' ,T112.'MASS'.

7 T120, 'MASS'/)

901 FORMAT ( // ' STOICHIOMETRIC FUEL SIZE= '.F8.3.' M2 '

)

905 FORMAT (1H1,T14, ' FUEL PYROLYSIS (POOL) SPECIFIED, VENTILATION'

1 • ADJUSTED FOR WORST CONDITIONS', TIOC'PAGE NO. ' , 13 .T 1 1 5

,

2 'RUN NO. • ,14//)

END

53

mma msjms^sass^s^^assa
SUBROUTINE PFLFlX

PESSIMIZATION ROUTINE

FIXED FUEL PYROLYSISt WINDOW WIDTH VARIED FOR WORST

BURNING CONDITIONS.

EITHER CRIB OR POOL FUELS ACCEPTED

FOR POOL FUELS t MUST ALSO SET PLFUEL=.TRUE.

COMMON /CNSTS/ AW ALLN t BW DOW» DENSA, G . GASCNT t KTRACE ,MTI ME

COMMON /CP/ CPA,CPC0(2}tCPCO2(2) tCPH2(2) tCPH20(2) .CPN2(2).

1 CP02(2) tCPPY«(2)

COMMON /FUEL/ C . CFLPC tCVGROS t CVNET t H iHFLPC . MWPYR t Nt NFLPCt O.

1 OFLPC.R. RO.REGRESiSH.SHAPEfSIZEtW.WFLPC. WTFUEL

COMMON /GP/ AWDOW.BPF,CD.CNV,DTIMEtEM£( 2) »HWDOW,IX.IXC. IXL.

1 Jt JMtJPtJPRI NT tKtKOiKHt KIT ERtKNTRLtMWiNtM WOUTtRCtRP • SIGMA

COMMON /LOGIC/ FC.FLSPECtKRIT t NEWPLT. NEWPRP . PLFUELtPLOT. PNCH

t

I RPSPEC.VTSPEC

COMMON /PLAST/ T BOILCi DHP. STOICH, SI ZEl t EITA , EISCAN

COMMON /PRBLM/ ADI A, AFLOOR, A WALL , DENS W. FLOAD . I RUN» OPENF.

1 PRNT.STEAOY.THICKW

COMMON /QS/ QCONW,QFIRE,QFLOW»QRADO,QRADW,QWLSUM

COMMON /TEMP/ DENF, DENUt T AMBt TGASt T INPT «T 1 ( 20 ) »T2{ 20 ) .TSF ,TSU

COMMON /THERML/ CNDA< 2 1 I ) i CPW( 2. I ) . DX EF. EMSA( 2 . 10)

.

1 NCNDf NCPWtNEMS«NQGENtNRP*QGEN(2» 10) «RPX(2«50}

COMMON /WOUT/ BWORSTt FLREMt HRATI Ot RMA.RMF .TTI ME . V AVGI Nt

1 WA»WB.YC02fYH20.YN2» Y02tYPYR

LOGICAL ADIA.EISCAN.FCtFLSPEC.KRITtNEWPRP.PLFUEL,

1 PLOT.PNCHtRPSPEC. SCANt STEADY. STOI CH, VTSPEC

REAL MWIN.MWOUT. MWPYR. MTI ME.Ni NFLPC

IF (STEADY) GOTO 190

FC= .FALSE.

SCAN= .FALSE.

QRADW=0.

QCONW= 0.

F2=0.

F1=0.

DTGAS=10.

TB0IL=TB0ILC+273.

CALL HEADNG , START TIME LOOP

DO 170 J=1.JM

KH=

DERIV1= 1. TGAS2= 0.

TGASl= 0.

TGASP= 2000.

TGASN= TAME

20 CONTINUE

K=

30 CONTINUE

KR=

RMA= 0.666667*CD*0.5*AWDOW*DENSA«SQRT(G*HWDOW*C 1 .-TAMB/TGAS )

)

IF (FLREM) 220,220.35

AS SOON AS FUEL IS EXHAUSTED PROGRAM MUST STOP.

SINCE WINDOW SIZE WOULD NOT BE WELL DEFINED.

35 IF { .NOT.PLFUELI GOTO 40

54
RP= SIZE*EF*SI GMA*(TGAS**4.-TB0IL**4,>/0HP

PLUME= S1ZE*0.0014.*CVNET/DHP

PROP= 1 .->(TGAS**4.-TBOIL**4.}/(1700.**4.=TBOIL«*4.)

IF CPROP.LT.O.) PROP= 0.

RP= RP+PROP*PLUME

GO TO 50

40 IF (REGRES.LE.0.0) GOTO 45

C USE THIS FORMULA IF INPUT REGRES IS SPECIFIED

RP= REGRES*2.*SHAPE/S1ZE*FLREM**{ 1 .-1 ./SHAPE )*WTFUEL**( 1. /SHAPE)

GO TO 50

45 CONTINUE

C FUEL SURFACE CONTROL

C ASSUME CRIB STICK DENSITY RHOCR= 500 KG/M**3

RHOCR= 500.

REGREN= 1.24E-3/RHOCR*SI2E**-0.6

RPl= RE6REN*2.*SHAPE/SIZE*FLREM*«( 1 .- I ./SHAPE)* WT FUEL** ( 1 ./SHAPE)

C CRIB POROSITY CONTROL

RP2= 0.22*WTFUEL/{RHOCR*SIZE)*SH

RP= AMINl (RP1.RP2)

50 RMF= RMA+ RP

YC02= 3.66e67*CFLPC*RC/100./RMF

YH20= (WFLPC*RP+9.0*HFLPC*RC)/100./RMF

Y02= (0.23*RMA-R0*RC)/RMF

YN2= 0.77*RMA/RMF +NFLPC*RP/l 00 ./RMF

YPYR= (RP-RC)/RMF

IF(YPYR.LT..O) YPYR= 0.

MW0UT= 44.*YC02+18.*YH20+28.*YN2+32.*Y02+MWPYR*YPYR

HRATIO= l./( l. + ( (TGAS/TAMB)*( MWIN/MWOUT ) ( 1 . +RP/RMA )**2

)

1 0.3333333333)

C NOTE HIN IS TAKEN AS POSITIVE

HIN= HWDOW* HRATIO

ZW=l .-MWOUT^TAMB/MWIN/TGAS

IF(ZW)195t55.55

55 VAVGIN= 0.ee6667^SQRT(2.^G+HlN^ZW)

RMA= CD^VAVGIN^HlN^BWDOW»DENSA

IF (RP-RMA/R) 60.60t65

60 RC= RP«BPF

BWORST= BWDOW^RC^R/RMA

RMA= RC*R/BPF

GO TO 70

65 RC= BPF*RMA/R

BWORST= BWDOW

C RECALCULATE Y- VALUES SINCE RP.RC HAVE BEEN CHANGED

70 KR= KR+1

IF (KR-3) 50.75*75

75 CONTINUE

QFLOW= RMF^(YCO2+(TGAS^(0.5«CPCO2( I ) TGAS+CPC02( 2) )-TAMB*(0 .5*

1 CPC02( 1 )«TAMB+CPC02(2) ) ) +YH20+(T GAS*( . 5*CPH20( 1 ) *TGAS+

2 CPH20(2) )-TAMB*{0.5*CPH2O( 1)*TAMB+CPH20(2) )) +Y02*(TGAS*(

3 0.5*CP02( 1 )*TGAS+ CP02(2) )-T AMB* ( .5*CP02 ( I )*TAMB+CP02 ( 2) )

)

4 +YN2*(TGAS*<0.5*CPN2( 1 ) *TGAS+CPN2 ( 2 ) )-T AMB*( 0. 5*CPN2 ( I )*

5 TAMB+CPN2(2)) ) +YPYR* (T GAS *( 0. S^CPPYR ( I ) *TGAS+CPPYR{ 2)

)

6 -TAMB^(0.5*CPPYR( I ) TA MB+CPPYR ( 2) ) ))

QFIRE= RC^CVNET

IF (ADIA) GOTO 90

CALL DESOLV

QRADW= AWALLN*EMS(

1

)*SI GMA^( TGAS^^4.-TSF^*4 .

)

QCONW= AWALLN*(TGAS-TSF)*CNV*( { TGAS-TSF )*( T GAS-TSF ) )**0 . 1 66 66667

55
90 CONTINUE

QRAOO= HWOOW*BWORST*SIGMA*(TGAS4i4>4.-TAMB*4>4 • ) K= K-l-l F3=F2

F2= F1

Fl= QFIRE-QFLOW-QRADO-QRADW-QCONW

TGAS3=TGAS2

TGAS2=TGAS1

TGAS1=TGAS

IF (Fl.LT.O. .AND.TGAS.lt .TGASP) TGASP=TGAS

IF (Fl .GT«0..AND.TGAS.6T.TGASN> TGASN=TGAS

0ERIV2= DERIVI

IF (TGAS1.EQ.TGAS2) GOTO 130

DERI V1=(F1-F2)/(T GAS X-TGAS2}

IF (KTRACE.GT.O) WRITE (4t99) TGASl tTGAS2 t Fl ,F2 .DERI V 1 ,K. KD

,

1 KH,JtT2(l) tTSFfQFIRE,QFLOW,QRADW,RP,RC

99 FORMAT (2F9.2t3( IPE9. 2) . 31 3 1 15 t 2( 0PF9. 2> t3( IPE10.3> .2C0PF7.3) ) IF (.NOT. SCAN) GOTO 95

IF {F1/F2.GE.0.0) GOTO 93

SCAN= .FALSE.

GOTO 100

93 TGAS= TGAS-DTGAS

IF ( (PLFUEL.AND.(TGAS.LT.TBOIL)) .OR.(TGAS.LT.TAMB) ) GOTO 200

GOTO 120

95 IF (DERIV2. LT..O. AND. ABS(F2) .GT.. 0001 >GOTO 100

IF(OERIV2.LT..O.ANO.J.GT.2) GOTO 100

TGAS= TGASl + OTGAS

GOTO 120

100 DIF= ABSCFl/QFLOW)

IF (DIF.LT.0.002.AND.ABS(TGAS2-TGAS1).LT.2. ) GOTO 130

TGAS=(F1*TGAS2-F2*TGAS1)/(F1-F2)

IF (K. GT. 10. AND. Fl.LT.O.. AND. TGAS.GT. TGASP) GOTO 105

IF (K. GT. 10. AND. Fl .GT.O.. AND.TGAS.lt .TGASN) GOTO 105

IF (TGAS.GT. TGASP. OR.TGAS.lt. TGASN) TGAS=(T GASP+TGASN) /2.

IF (K.EQ.l.AND.KH.EQ.O) TGAS= TGASl+10.

IF (TGAS.GT. 2000.) GOTO 110

IF (TGAS.LT.(TAMB+30.) ) GOTO 110

IF (PLFUEL.AND.TGAS.LT.TBOIL) GOTO 110

GOTO 120

105 TGAS= (TGASN+TGASP)/2.

GOTO 120

110 SCANS .TRUE.

TGAS= 1900.

120 CONTINUE

IF (K-200) 30.30.200

130 CONTINUE

CALL RSTA

FLREM= FLREM-RP*OTlME

IF(FLREM.LT.O) FLREM=0.

IF (QCONW.GT.O. ) QWLSUM= QWLSUM+( QRADW+QCON W )*DTI ME

IF (TTIME .GE. MTIME) GO TO 210

IF (TGAS.LE.353..AND.J.GE.10) GO TO 210

IF ( J.EQ.l ) GO TO 150

IF < JP.LT.JPRINT) GO TO 160

JP=

150 CALL OUTPUT

160 JP= JP+1

TTIME= TTIME+DTIME

56
170 CONTINUE

END TIME STEP OO-LOOP

PET URN

ERROR IN INPUT

190 CONTINUE

KNTRL= 2 WRITE (2.910)

910 FORMAT (///• PFLFIX ROUTINE DOES NOT ACCEPT STEADY-STATE CASE')

195 CONTINUE

1F(KTRACE.EQ. 1) WRITE(2«930) TGAS .RC.RP. YPYR . ZW .RMA, MWOUT

930 FORMAT (/• TGAS=' . F5 .0. • RC=«.E10.4.« RP=«.E10.4.

I • YPYR=',E10.*. • ZW=«.F6.4.« RMA= • , El 0.4. • MWOUT = « t F6 . 1 ) FAIL TO CONVERGE. ERROR EXIT

200 CONTINUE

KNTRL=3

RETURN

FIRE IS OVER (TRANSIENT CASE)

210 CONTINUE

CALL OUTPUT

220 CONTINUE

RETURN

END

57

mma^mas^ammsBo^
SUBROUTINE POOL

C

C PLASTIC FUEL (POOL FIRE) ROUTINE

C FUEL 'SEES* ONLY COMPARTMENT AND NOT ITSELF

C

COMMON /CNSTS/ AWALLN , BWDOW , DENSA , G . GASCNT , KTRACE , MTI ME

COMMON /CP/ CPA,CPC0(2J,CPC02(2) .CPH2<2) »CPH20(2) .CPN2(2) t 1 CP02(2> tCPPYR(2J

COMMON /FUEL/ C .CFLPC» CVGROS* CVNET tH ,HFLPC . MWPYR .N, NFLPC. O.

1 OFLPC.RtR0tREGRES»SHtSHAPEtSlZE.WtWFLPC» WTFUEL

COMMON /GP/ AWDOW»BPF,CDtCNVt DTlMEtEMS( 2) .HWDOW.IX. IXC. IXL»

1 J« JMtJPtJPRI NT tK«KO«KH,K ITER (KNTRL«MWIN«MW0UT*RC*RP. SIGMA

COMMON /LOGIC/ FC.FLSPECf KRI T« NE WPLT. NEWPRP »PLFUEL. PLOT . PNCH

.

I RPSPEC.VTSPEC

COMMON /PLAST/ TBOILC. DHP, STOI CH, SI ZEl .EIT A , EISCAN

COMMON /PRBLM/ ADI A, AFLOOR. AWALL t OENSW, FLO AD,IRUN .OPENF

,

1 PRNT.STEADYiTHICKW

COMMON /QS/ QCONtt,QFlRE»QFLOW.QRADO»QRADW.QWLSUM

COMMON /TEMP/ DENF. DENU.T AM B.TGAS. TINPT .T 1 ( 20 ) . T2( 20 ) tTSF ,T SU

COMMON /THERML/ CNDA( 2 1 10 ) . CPW( 2. I ) .DX . EF, EMSA( 2 . I 0)

1 NCN0»NCPWtNEMS«NQGENtNRPtQGEN(2t 10) «RPX( 2*50)

COMMON /WOUT/ BWORST»FLREMt HRATI O t RMA.RMF ,TTI ME . VAVGI N,

1 WAtWBt YC02t YH20.YN2t Y02tYPYR

LOGICAL AOIA, EI SCANt FCt FLSPEC » KRI T. NEWPRP. PLFUEL

.

1 PLOT .PNCH.RPSPECf SCAN, STEADY, STOICH,VT SPEC

REAL MWIN, MWOUT, MWPYR, MTI ME, N, NFLPC

C Y- DENOTE MASS FRACTIONS OF OUTFLOW

C RMA= MASS INFLOW RATE

C RMF= MASS OUTFLOW RATE

C R= STOICHIOMETRIC AIR/FUEL MASS RATIO

C R0= STOICHIOMETRIC OXYGEN/FUEL RATIO

SCAN= .FALSE.

C INITIALIZES STOICHIOMETRIC CASE

IF (STOICH) EITA =1.

QRAOW=0.

QCONW= 0.

F2=0.

F1=0.

DTGAS = 10.

TBOIL=TBOlLC+273.

FC= .FALSE.

CALL HEAONG

TF2= TGAS

C START TIME STEP

DO 170 J=1.JM

KH=

DER1V1= I. TGAS2= 0.

TGAS1= 0-

TGASP= 2000.

TGASN= TAMB

20 CONTINUE

K=

30 CONTINUE

IF (STOICH) GOTO 34

32 RP= SIZE*EF«SIGMA*(TGAS**4.-TB0IL**A.)/DHP

IF (STEADY) GOTO 33

I

58
PLUME= SIZE*0,0014*CVNET/DHP

PROP= 1 .-(TGAS**A.-TB0IL**4.)/( I700.**4.-TBOIL«*4 .

)

IF (PROP.LT.O.) PROP= 0.

RP= RP+PROP*PLUME

33 CONTINUE

RMF= RMA+RP

1F(FLREM.LE.0)RP=0.

34 YC02= 3.66667*CFLPC«RC/100./RMF

YH20= (WFLPC*RP+9.0*HFLPC*RC)/100./RMF

Y02= C0.23*RMA=R0«RC)/RMF

YN2= 0.77*RMA/RMF +NFLPC*RP/1 00 ./RMF

YPYR= (RP-RC)/RMF

IF(YPYR.LT. .0) YPYR= 0.

QFUEL= (RP-RC)*DHP

MWOUT= 44.«YC02+18.*YH20+28.*YN2+32«*Y02+MWPYR*YPYR

HRAT2 0= l./( l.+( (TGAS/TAMB)*(MWlN/MWOUT)*( 1 .+RP/RMA) 2)

1 **0 .3333333333)

NOTE HlN IS TAKEN AS POSITIVE

HIN= HWDOW* HRATIO

ZW= 1 .-MWOUT*T AM 8/ MW I N/T GAS

IF<ZW)195t35t35

35 VAVGIN= 0.ee6667*SQRT(2.«G«HIN*ZW)

RMA= CO*VAVGlN*HIN*BWDOW*DENSA

RMF= RMA+RP

RC= BPF*RMA/R

IF (STOICH) RP= RC/BPF

IF < (EISCAN. AND. ElTA.LT.l.). OR. STOICH) GOTO 37

36 IF (RC.GT.RP*BPF) GOTO 40

37 FC= .FALSE.

GO TO 45

40 FC= .TRUE.

RC= RP*BPF

RECALCULATE VALUES IF IN FUEL CONTROL REGIME.

YC02= 3.66667*CFLPC*RC/100./RMF

YH20= (WFLPC*RP+9.0*HFLPC*RC)/100./RMF

Y02= {0.23*RMA=R0*RC)/RMF

YN2= 0.77*RMA/RMF +NFLPC*RP/l 00 ./RMF

YPYR= (RP-RC)/RMF

45 CONTINUE

QFLOW= RMF»(YC02*(TGAS*(0.5*CPC02(1 )*TGAS+CPC02( 2 ) )-TAMB*{0.5*

1 CPC02( 1 )*TAMB+CPC02<2) ) ) +YH20*( TGAS* ( . 5«CPH20 ( 1 ) *T6AS+

2 CPH20(2) )-TAMB*( 0.5*CPH2O( l)*TAMB+CPH20(2))) +Y02*(TGAS*(

3 0.5*CPO2{l )*TGAS + CP02<2) )-T AMB*( 0.5*CP02 ( I )«TAMB+CP02( 2) )

)

4 +YN2*(TGAS*(0.5*CPN2( I ) TGAS+CPN2 ( 2) )-TAMB*( 0.5*CPN2( 1 )*

5 TAMB+CPN2(2) ) ) +YP YR* (T GAS* ( . 5*CPPYR ( 1 ) TGAS+CPPYR ( 2)

)

6 -TAMB«(0 .5*CPPYR( I )«TAMB+CPPYR ( 2) ) )

>

QFIRE= RC*CVNET

IF (AOIA) GOTO 90

IF ( .NOT. STEADY) CALL DESOLV

IF (STEADY) CALL STFLOW

QRADW= AWALLN*EMS( 1 )*SI GMA*(T GAS**4.-=TSF**4 . ) QCONW= AWALLN*(TGAS-TSF)*CNV*( ( TGAS-TSF) *(TGAS-TSF ) )**0 . 16666667

90 CONTINUE

QRADO= AWOOW*SI GMA*(TGAS**4.-TAMB**4 .

)

K= K+1

F3= F2

F2=F1

Fl= QFIRE-QFLOW-QFUEL-QRAOO-QRADW-QCONW

59
TGAS3=TGAS2

TGAS2=TGAS1

TGAS1=TGAS

IF tFl.LT.0..ANO.TGAS.LT.TGASP) TGA£P=TGAS

IF (Fl.GT.O. .AND.TGAS.GT.TGASN) TGASN=TGAS

DERIV2= DERIVl

IF (TGAS1.EQ.TGAS2) GOTO 130

DER1VI=(F1-F2)/(TGASI-TGAS2)

IF (KTRACE.GT.O) WRITE (4,99) T GAS! .TGAS2. F 1 ,F2.DERI V I , K .KD

,

I KH.JtT2( 1 J tTSF,QFIRE,QFLOW,QRADW.RP,RC

99 F0RMAT(2F9.2.3( 1PE9.2),3I3.15.2( 0PF9.2).3( lPEiO.3) .2(0PF7.3) ) IF (.NOT. SCAN) GOTO 95

IF (F1/F2.GE.0. 0) GOTO 93

SCAN= .FALSE.

GOTO 100

93 TGAS= TGAS-DTGAS

IF (TGAS.LT.TAMB) GOTO 200

IF (TGAS.LT.TBOIL.ANO.(FLREM.GT.O.O) ) GOTO 190

GOTO 120

95 IF (DERIVl .LT..0. AND. ABS(F2) .GT..0001)GOTO 100

IF(DERIV2.LT..0.AND.J.GT.2) GOTO 100

TGAS= TGASl+OTGAS

GOTO 120

100 DIF= ABS(Fl/QFLOW)

IF (DIF.LT.0.002.ANO.ABS(TGAS2=TGASl ).LT.2. ) GOTO 130

TGA£=(FI*TGAS2-F2*TGAS1)/(F1-F2)

IF (K.GT.10.AND.F1.LT.0..AND.TGAS.GT.TGASP) GOTO 105

IF (K. GT. 10. AND. Fl.GT.O. .ANO.TGAS.lt .TGASN) GOTO 105

IF (K.EQ.l .AND.KH.EQ.O) TGAS= TGASl+10.

IF (TGAS.GT.2000. ) GOTO 110

IF (TGAS.LT.(TAMB+30.) ) GOTO 110

IF (TGAS.LT.TBOIL.AND.(FLREM.GT.O.O)) GOTO 110

GOTO 120

105 TGAS= (TGASN+TGASP)/2.

GOTO 120

110 SCAN= .TRUE.

TGAS= 1900.

120 CONTINUE

IF (STEADY. AND. .NOT. ADIA) CALL STFLOW

IF (K-200) 30.30.200

130 IF (STEADY) GOTO 180

CALL RSTA

FLREM= FLREM-RP*DTIME

IF(FLREM.LT.O) FLREM=0.

IF (QCONW.GT.O. ) QWLSUM= QWLSUM+( QRADW+QCONW )*DTI ME

IF (TTIME .GE. MTIME) GO TO 210

IF (TGAS.LE.353. .AND.J.GE. 10) GO TO 210

IF (J.EQ.l ) GO TO 150

IF (JP.LT.JPRINT) GO TO 160

JP=

150 CALL OUTPUT

160 JP= JP+1

TTIME= TTIME+DTIME

170 CONTINUE

END TIME STEP DO-LOOP

180 CONTINUE

IF( .NOT.STOICH) GOTO 185

FIND STOICHIOMETRIC FUEL SIZE

60

I
SIZE1= RP/(EF*SIGMA*(TGAS**4.-TB0IL**4.)/DHP)

185 CALL OUTPUT

NORMAL EXIT WHEN STEADY. EQ.T

RETURN

ERROR EXIT

190 CONTINUE

IF (KTRACE.EQ.l) WRITE (2t910> TGAS . 910 FORMAT!///' TGAS.LT.TBOIL TGAS=« t F8. 1

.

• GOTO NEXT CASE'///)

GOTO 200

SQUARE ROOT ERROR

195 CONTINUE

IF(KTRACE.EQ.l ) WRITE(2t930) TGAS t RC.RP. YPYR t ZW ,RMA, MWOUT

930 FORMAT (/' TGAS=' t F5 .0 .' RC=',E10.4.' RP='.E10.4»

1 ' YPYR=',E10.4. ' ZW='.F6.4t' RMA=' .El 0.4 t • MWOUT = ' . F6. I

)

FAIL TO CONVERGE. ERROR EXIT

200 CONTINUE

KNTRL=3

RETURN

FIRE IS OVER (TRANSIENT CASE)

210 CONTINUE

CALL OUTPUT

RETURN

END

61
SUBROUTINE PP

PLOTTING SUBROUTINE

THIS ROUTINE IS LEFT BLANK SINCE IT IS MACHINE-DEPENDENT

ENTRY PLTRST

THIS ENTRY SETS UP THE INIALIZATION OF PLOTTING

ENTRY DSTO

THIS ENTRY IS CALLED EACH TIME TO STORE A DATA POINT

ENTRY DOUT

THIS IS THE LAST ENTRY FOR A GIVEN RUN

RETURN

END

I

62
SUBROUTINE PVTFIX

C

C PESSIMIZATION ROUTINE

C FIXED VENTILATIONt WORST POSSIBLE FUEL PYROLYSIS RATE.

C

COMMON /CNSTS/ AWALLN, BW DOW, DENS A, G . GASCNT , KTRACE , MTI ME

COMMON /CP/ CPA,CPC0(2J.CPC02(2) .CPH2( 2) .CPH20(2>,CPN2(2)

.

1 CP02(2) ,CPPYR{2)

COMMON /FUEL/ C, CFLPCt CVGROS, CVNET t H,HFLPC. MWPYR , N. NFLPC, O.

1 OFLPC,R.R0tREGRES.SH,SHAPE,SIZE,W,WFLPC» WTFUEL

COMMON /GP/ AWDOW,BPF,CD»CNV. DTIMEtEMS(2>,HWD0W»IXtIXCtIXL.

1 J.JMtJPt JPRI NT t K. KD » KH » KIT ER,KNTRL,MWlN.MWOUT»RCtRP, SIGMA

COMMON /LOGIC/ FC, FLSPEC, KRITt NEWPLT, NEWPRP .PLFUEL.PLOT .PNCH

,

1 RPSPECt VTSPEC

COMMON /PLAST/ T BOlLC, OHP , STOlCH, SI ZEl , EITA , EI SCAN

COMMON /PRBLM/ ADI A , AFLOOR , A WALL , DENS W»FLOAD, I RUN, OPENF

,

1 PRNT.STEADYtTHICKW

COMMON /QS/ QCONW,QFIRE,QFLOW, QRADO»QRADW,QWLSUM

COMMON /TEMP/ DENF, DENUtT AMB.TGAS, TINPT ,T 1 ( 20i .T2( 20 ) ,TSF ,T SU

COMMON /THERML/ CNDA( 2, I ) , CPW( 2, 1 0) , DX , EF, EMSA{ 2 , 1 0)

,

1 NCND,NCPW,NEMS,NQGEN,NRP,QGEN(2, 10) ,RPX(2,50)

COMMON /WOUT/ BWORSTtFLREM,HRATI0,RMAtRMF,TTlME,VAVGlN,

1 WA,WB,YC02,YH20, YN2, Y02,YPYR

LOGICAL ADIA,EISCAN,FC,FLSPEC,KRIT,NEWPRP,PLFUEL,

1 PLOT,PNCH,RPSPEC,SCAN,STEADY,STOICH,VTSPEC

REAL MWIN,MW0UT»MWPYR,MTIME,N,NFLPC

IF (STEADY) GOTO 190

FC= .FALSE,

SCAN= .FALSE.

QRADW=0.

QCONW= O.

F2=0.

F1=0.

DTGAS=10.

CALL HEADNG

C START TIME LOOP

DO 170 J=1,JM

KH=

DERIV1= 1. TGAS2= 0.

TGAS1= 0.

TGASP= 2000.

T6ASN= TAMB

20 CONTINUE

K=

30 CONTINUE

IF (FLREM.GT.O. ) GOTO 32

RC= 0.

RP= 0.

FC= .TRUE.

32 RMF= RMA+RP

YC02= 3.6e667*CFLPC«RC/100./RMF

YH20= (WFLPC*RP+9.0«HFLPC«RC>/100./RMF

Y02= (0.23«RMA-R0«RC)/RMF

YN2= 0.77«RMA/RMF +NFLPC*RP/1 00 ./RMF

YPYR= (RP-RCJ/RMF

1F{YPYR.LT..0) YPYR= 0.

63
MWOUT= 44.*YC02+18.*YH20+28.*YN2+32.*Y02+MWPYR*YPYR

HRATIO= l./( I .+( (TGAS/TAMB)*(MWlN/MWOUT)*( I . RP/RMA)**2

)

1 **0. 3333333333)

NOTE HIN IS TAKEN AS POSITIVE

HIN= HWDOW* HRATIO

ZW=1 .-MW0UT*TAMB/MWIN/TGAS

IF( ZW)195t35.35

35 VAVGIN= 0.666667*SQRT(2.*G*HIN*ZW)

RMA= CD*VAVGlN*HlN«BWOOW*DENSA

RMF= RMA+RP

IF (.NOT.FC) RC= BPF*RMA/R

IF (.NOT.FC) RP= RC/BPF

45 CONTINUE

QFLOW= RMF«(YC02*(TGAS«(0.5«CPC02< 1 )*TGAS+CPC02 ( 2 ) )-TAMB*(0 .5*

1 CPC02(1 J*TAMB+CPC02(2} ) ) +YH20*( T GAS*( . 5*CPH20 ( 1 ) *TGAS+

2 CPH20(2i )-TAMB*(0«5*CPH2O( l)*TAMB+CPH20(2J )) +Y02*(TGAS*(

3 0.5*CPO2(l )*T6AS+CP02(2) >-T AMB*( .5*CP02 ( 1 )*TAMB+CP02( 2) )

)

4 +YN2«(TGAS«(0.5*CPN2( I } *TGAS+CPN2( 2) )-T AMB*( . 5«CPN2 ( 1 )

5 TAM8+CPN2(2}> ) +YPYR* (T GAS* ( . 5*CPPYR ( 1 > TGAS+CPPYR( 2}

)

6 -TAMB*{0.5*CPPYR( I > TAMB+CPPYRt 2) ) ))

QFIRE= RC*CVNET

IF (ADIA) GOTO 90

CALL DESOLV

QRADW= AWA'_LN*EMS( 1 )*SI GMA*(T 6AS**4 .-TSF**4 . ) QCONW= AWALLN*(TGASTMTSF>*CNV*( (T6AS-TSF ) (T GAS-TSF > }**0 . I 6666667

90 CONTINUE

QRAOO= AWDOW*SI GMA* (TGAS**4 .-TAMB**4 .

)

K= K+1

F3 = F2

F2=F1

Fl= QFIRE-QFLOW-QRADO-QRADW-QCONW

TGAS3=TGAS2

TGAS2=TGASl

TGAS1=TGAS

IF (F1.LT.0..AND.TGAS.LT.TGASPJ TGASP=TGAS

IF (Fl.GT.O. .AND.TGAS.GT.TGASN) TGASN=TGAS

DERIV2= OERIVl

IF {TGAS1.EQ.TGAS2J GOTO 130 i DERIV1=(FI-F2)/(TGAS1-TGAS2)

IF (KTRACE.GT.O) WRITE (4f99) TGASl .TGAS2. F 1 .F2 .DERI V 1 , K , KD

,

t KHt J.T2( I )»TSF,QFIREtQFLOW.QRADW.RP.RC

99 FORMAT (2F9. 2 t3{ 1PE9. 2) ,313 . 15 . 2( 0PF9 .2 ) .3< IPEIO.3) ,2(0PF7,3)

}

IF (.NOT. SCAN) GOTO 95

IF (F1/F2.GE.0.0} GOTO 93

SCAN= .FALSE.

GOTO 100

93 TGAS= TGAS-OTGAS

IF (TGAS.LT.TAMB) GOTO 200

GOTO 120

95 IF (OERIVl. LT..O. AND. ABS(F2) .GT-.OOODGOTO 100

IF(DERIV2.LT..0.AND.J.GT.2} GOTO 100

TGAS= TGASl+OTGAS

GOTO 120

100 DIF= ABS(F1/QFL0W>

IF {DIF.LT.0.002.AND.ABS(TGAS2-TGAS1) .LT.2. ) GOTO 130

TGAS=(F1*TGAS2-F2*TGAS1 )/(Fl-F2>

IF (K.GT.lO.AND.Fl .LT.O..AND.TGAS.GT.TGASP) GOTO 105

IF (K.GT.lO.AND.Fl .GT.O-. ANO.TGAS.lt .TGASN) GOTO 105

64
IF (K.EQ.l.ANO.KH.EQ.O) T GAS= TGASl +10.

IF (TGAS.GT.2000. ) GOTO 110

IF (TGAS.LT.(TAMB+30.J ) GOTO 110

GOTO 120

105 TGAS= (TGASN+TGASP)/2.

GOTO 120

110 SCAN= .TRUE.

TGAS= 1900.

120 CONTINUE

IF (K-200) 30t30»200

130 CONTINUE

CALL RSTA

FLREM= FLREM-RP^DTIME

IF(FLREM.LT.O) FLREM=0.

IF (QCONW.GT.O. ) QWLSUM= QWLSUM+( QRADW+QCCNW )*DTI ME

IF (TTIME .GE. MTIMEJ GO TO 210

IF (TGAS.LE.353..AND.J.GE.10) GO TO 210

IF (J.EQ.l ) GO TO 150

IF { JP.LT. JPfiINT) GO TO 160

JP=

150 CALL OUTPUT

160 JP= JP+1

TTIME= TTIME+DTIME

170 CONTINUE

END TIME STEP DO-LOOP

RETURN

ERROR IN INPUT

190 CONTINUE

KNTRL= 2 WRITE (2.910)

910 FORMAT (///• PVTFIX ROUTINE DOES NOT ACCEPT STEADY-STATE CASE'}

RETURN

SQUARE ROOT ERROR

195 CONTINUE

IF(KTRACE.EQ.l) WRITE(2.930) TGAS. RC. RP , YPYR , ZW ,RMA , MWOUT

930 FORMAT (/• TG AS = « .F5. 0» • RC = «.E10.4.' RP=«.E10.4,

1 • YPYR=« ,E10.4. • ZW='.F6.4.* RMA= • , E 10.4, • MWOUT=« ,F 6. 1

)

FAIL TO CONVERGE, ERROR EXIT

200 CONTINUE

KNTRL = 3

RETURN

FIRE IS OVER (TRANSIENT CASE)

210 CONTINUE

CALL OUTPUT

RETURN

END

65
SUBROUTINE RPFI

X

TABULAR FUEL PYROLYSIS ROUTINE

FUEL PYROLYSIS RATE IS AN INPUT VARIABLE.

COMMON /CNSTS/ A WALLN, BWOOW, DENSA .G « GASCNT » KTRACE .MTI ME

COMMON /CP/ CPA,CPC0(2).CPC02(2) .CPH2(2) .CPH20(2) tCPN2(2) •

1 CP02{2).CPPYR(2)

COMMON /FUEL/ C. CFLPCt CVGROSt CVNET . H, HFLPC .MWPYR ,N , NFLPC . O,

I OFLPCt R.RO.REGRES.SH. SHAPE t SI ZEtW.WFLPC.WTFUEL

COMMON /GP/ AWOOWtBPF.CDtCNV. DTIME,EMS<2)»HWD0W,IX,IXC,IXL,

I J, JMt JP» JPRINT«KtKD«KH.KlTER.KNTRL»MWlN. MWOUT . RC. RP t SI GMA

COMMON /LOGIC/ FC» FLSPECt KRIT t NEWPLT. NEWPRP ,PLFUEL t PLOT .PNCH

.

1 RPSPEC. VTSPEC

COMMON /PLAST/ T BOlLCt DHP, STOICH, SI ZEl , EITA , EI SCAN

COMMON /PRBLM/ ADI At AFLOOR t AWALL t DENSW.FLO AD .1 RUN, OPENFt

1 PRNT.STEADY.THICKW

COMMON /QS/ QCONW,QFIREiQFLOW.QRADO»QRAOW,QWLSUM

COMMON /TEMP/ DENF, OENU.T AMB.TGAS.T INPT . T 1 ( 20 ) .T2( 20 ) ,TSF ,TSU

COMMON /THERML/ CNDA( 2, 1 ) , CPW( 2, I 0) .DX ,EF, EMSA( 2 . 1 )

,

1 NCND.NCPW.NEMS.NQGENf NRPt QGEN(2. 10)«RPX( 2.50}

COMMON /WOUT/ BWORST»FLREM.HRATIO.RMA.RMF,TTlME. VAVGiNt

I WAt WB»YC02tYH20.YN2» Y02,YPYR

LOGICAL ADIA.EISCAN.FC.FLSPEC.KRITtNEWPRP.PLFUEL*

I PLOT. PNCH, RPSPEC. SCAN, STEADY, STOICH, VTSPEC

REAL MWI N.MWOUT, MWPYR. MTI ME. N. NFLPC

SCAN= .FALSE.

QRAOW=0.

QCONW= 0.

F2 = 0.

Fl = 0.

DTGAS=10.

CALL HEADNG

START TIME LOOP

DO 170 J=1.JM

KH=

DERIVl= 1.

TGAS2= 0. ' TGASl= 0.

TGASP= 2000.

TGASN= TAMB

20 CONTINUE

K=

30 CONTINUE

FC= .FALSE.

IF (FLREM.GT.O. ) RP= TLU ( RPX, NRP. TTI ME>

IF (FLREM.LE.O.i RP= 0.

RMF= RMA+RP

YC02= 3.66667*CFLPC«RC/100./RMF

YH20= (WFLPC*RP+9.0*HFLPC*RC)/100 ./RMF

Y02= (0.23*RMA-'R0*RC>/RMF

YN2= 0.77*RMA/RMF +NFLPC*RP/100 ./RMF

YPYR= (RP-RC)/RMF

IF(YPYR.LT. .0) YPYR= 0.

MWOUT= 44.*YC02+18.*YH20+28.*YN2+32.*Y02+MWPYR*YPYR

HRATIO= l./( l.+( (TGAS/TAMB)*(MWIN/MW0UT)*(1.+PP/RMA)**2)

66
1 «*0 .3333333333)

NOTE HIN IS TAKEN AS POSITIVE

HIN= HWDOW* HRATIO

ZW=l.-MWOUT*TAMB/MWlN/TGAS

IF{ZW>195.35.35

35 VAVGIN= 0.6e6e67*SQRT(2.*G*HIN*ZW)

RMA= CD*VAVGIN*HIN*BW00W«DENSA

RMF= KMA+fiP

IF (RMA/R-RP) 40.*0t45

40 RC= BPF*RMA/'R

GO TO 50

45 RC= BPF*RP

FC= .TRUE.

50 CONTINUE

QFLOW= RMF*(YC02*(TGAS*(0.5*CPC02( 1 ) *TGAS+CPCD2( 2 )

) -TAME* (0.5*

1 CPC02( 1 )*TAMB+CPC02(2)) ) +YH20*(T GAS*( . 5«CPH20( 1 ) TGAS+

2 CPH20(2) )-TAMB*{0.5*CPH2O{ l)*TAMB+CPH20(2i) ) +Y02*(TGAS*<

3 0.5*CPO2(l )*TGAS+CP02(2} )-T AMB*{ .5*CP02 ( 1 )*TAMB+CP02 ( 2) )

)

4 +YN2*(TGAS*( 0.5*CPN2(1 }*TGAS+CPN2 ( 2 J )"TAMB*( . 5*CPN2( 1 )*

5 TAMB+CPN2C2)) i +YPYR*( TGAS*( 0. 5*CPPYR( 1 ) TGAS+CPPYRC 2 )

)

6 =TAMB«(0.5*CPPYR( 1 ) *T AMB+CPPYR (2) ) ))

QFIRE= RC*CVNET

IF (AOIA) GOTO 90

IF ( .NOT. STEADY) CALL OESOLV

IF (STEADY) CALL STFLOW

QRADW= AWALLN*EMS( 1 )*SI GMA*(TGAS**A.-TSF**4 .

)

QCONW= AWALLN*(TGAS-TSF)*CNV*( (TGAS-TSF) ( T GAS-TSF ) )**0. 1 6666667

90 CONTINUE

QRADO= AWOOW*SI GMA*( TGAS**4 .-T AMB**4 .

)

K= K+l

F3=F2

F2 = Fl

F 1= OFIRE-QFLOW-QRADO-QRADW-QCONW

TGAS3=TGAS2

TGAS2=TGAS1

TGAS1=TGAS

IF (Fl.LT.O. .ANO.TGAS.lt .TGASP) TGASP=TGAS

IF (Fl,GT.O..AND.TGAS.GT.TGASN) T6ASN=TGAS

DERIV2= DERIVl

IF (TGASl .EQ.TGAS2) GOTO 130

DERIV1=(F1-=F2)/(TGAS1-TGAS2)

IF (KTRACE. GT.O) WRITE (4*99) TGASl .TGAS2. F I . F2 . DERI V I . K. KO

.

I KH, J,T2( 1) .TSF.QFlRE.QFLOW.QRADW.RPtRC

99 FORMAT ( 2F9.2.3( IPE9. 2) t 313. 15. 2( 0PF9.2 ) . 3 ( 1 PEl .3

)

,2(0PF7.3) ) IF (.NOT. SCAN) GOTO 95

IF (F1/F2.GE.0. 0) GOTO 93

£CAN= .FALSE.

GOTO 100

93 TGAS= TGAS-DTGAS

IF (TGAS.LT.TAMB) GOTO 200

GOTO 120

95 IF (DERIVl .LT..O. AND. ABS(F2) .GT..0001)GOTO 100

IF(DERIV2.LT. .O.AND. J.GT.2) GOTO 100

TGAS= TGASl+OTGAS

GOTO 120

100 DIF= ABS(F1/QFL0W)

IF (DIF.LT.0.002.AND.ABS(TGAS2-TGAS1) .LT.2. ) GOTO 130

TGAS=(F1*TGAS2-F2*TGASI)/(F1-F2)

67
IF (K.GT.lO.AND.Fl .LT.O..AND.TGAS.GT.TGASP) GOTO 105

IF (K.GT.lO.AND.Fl .GT.O.. AND.TGAS.lt. TGASN) GOTO 105

IF (K.EQ.l. AND.KH.EQ.O) TGAS= TGASl+lO.

IF (TGAS.GT.2000.> GOTO I 10

IF (TGAS.LT.(TAMB+30.) ) GOTO 110

GOTO 120

105 TGAS= (TGASN+TGASP)/2.

GOTO 120

110 SCAN= .TRUE.

TGAS= 1900.

120 CONTINUE

IF (STEADY. AND.. NOT. ADIA) CALL STFLOVK

IF (K-200) 30.30.200

130 CONTINUE

IF (STEADY) GOTO 180

CALL RSTA

FLREM= FLREM-RP^DTIME

IF(FLREM.LT.O. J FLREM=0.

IF (QCONW.GT.O.) QWLSUMs: QWLSUM-l-( QRADW-I-QCON W )*DTI ME

IF (TTIME .GE. MTlMEi GO TO 210

IF (TGAS.LE.353..ANO.J.GE.10) GO TO 210

IF (J.EQ.l ) GO TO 150

IF ( JP.LT. JPRINT) GO TO 160

JP=

150 CALL OUTPUT

160 JP= JP+l

TTIME= TTIME+DTIME

170 CONTINUE

END TIME STEP 00-LOOP

180 CONTINUE

185 CALL OUTPUT

NORMAL EXIT WHEN STEADY. EQ.T

RETURN , SQUARE ROOT ERROR

195 CONTINUE

IF(KTRACE.EQ.l) WRITE(2i930) TGAS.RC .RP. YPYR» ZW ,RMAt MWOUT

930 FORMAT (/• TGAS=« iF5 .0. • RC=«.E10.4.» RP=«,E10.4.

1 • YPYR=».E10.4. • ZW=«.F6.4t« RMA= • t El 0.4. • MWOUT = « . F6 . 1

)

FAIL TO CONVERGE. ERROR EXIT

200 CONTINUE

KNTRL=3

RETURN

FIRE IS OVER (TRANSIENT CASE)

210 CONTINUE

CALL OUTPUT

RETURN

END

I

68 II
SUBROUTINE STFLOW

C

C CALCULATES WALL HEAT CONDUCTION WHEN STEADY-STATE

C CONDITION ONLY IS NEEDED,

C

COMMON /GP/ AWDOWtBPF.CD,CNV,DTlME«EMS(2) .HWOOWtlX.IXCt IXLt

1 J. JM. JP. JPRINTtK.KOtKH.KITERfKNTRLtMWlN.MWOUT.RC.RP.SIGMA

COMMON /PRBLM/ AOIAt AFLOORt AWALL, DENS W.FLO AD t I RUN.OPENF»

1 PRNT.STEADY.THICKW

COMMON /TEMP/ DENF »DENU.TAMB,TGAS.TI NPT .T 1 ( 20) .T 2( 20 } .TSF t TSU

COMMON /THERML/ CNDAC2. 1 > . CPW( 2. 1 0) . DX . EF , EMSA ( 2 . 10 )

.

1 NCNDtNCPWtNEMS*NQGEN«NRPtQGEN( 2« 10) «RPX( 2 t 50

)

C BIOT= BIOT NUMBER

KD=

TSF= TGAS -30.

TSU= TAMB +30.

10 CONTINUE

TSFOLD= TSF

TSUOLD= TSU

EMS(1)= l./( l./TLU{EMSA.NEMS.TSF) +l./EF-l.i

EMS(2)= TLU(EMSA. NEWS. TSU)

TAVG= (TGAS+TAMB)/2.

CND= TLU(CNDA.NCNO»TAVG)

ZRF= TGAS*(TGAS*(TGAS-i-TSF)+TSF*TSF)+TSF*TSF*TSF

ZCF= CNV*( (TGAS-TSF)*(TGAS-TSF) J**0. 16666667

HF= ZCF+EMS( 1 )*S1GMA«ZRF

B10TF= HF*THICKW/CND

ZRU= TAMB«(TAMB*{TAMB+TSU)+TSU«TSU)+TSU*TSU*TSU

ZCU= 1.31*( (TAMB-TSU)*(TAMB-TSU))** 0.166666 67

HU= ZCU+EMS(2)*SIGMA*ZRU

BI0TU= HU*THICKW/CND

TSF= ( (BlOTF+HF/HU)*TGAS+TAMB)/( l.+BIOTF+HF/HU)

TSU= ( (BI0TU+HU/HF)*TAMB+TGAS)/(1 . +BIOTU+HU/HF)

T2(l)= TSF- {TSF-TSU)*DX/THICKW/2.

T2(IXC)= (TSF+TSU>/2.

IF ( (ABS<TSF-TSFOLD) .LT.3.) .AND. ( ABS{TSU-TSUOLD) . LT .3 . )

)

1 RETURN

KD= KD+1

TSU= <TSU+TSU0LD)/2.

IF <KD.LT.20) GOTO 10

WRITE (2.90) TSF.TSFOLD

90 FORMAT (//• UNSUCCESSFUL ITERATION IN STFLOW'/

1 • TSF=«,F15.2. • TSFOLD=« ,F15.2)

RETURN

END

69
FUNCTION TLU ( AR«AY,NUM, VALIN

)

TABULAR LOOK-UP INTERPOLATING ROUTINE

DIMENSION ARRAY(2«NUM)

ARRAY(1.I>= INDEPENDENT VARIABLE

ARRAY(2»I>= DEPENDENT VARIABLE

INTERPOLATES LINEARLY WITHIN GIVEN DOMAIN. SETS EQUAL TO

SMALLEST OR LARGEST VALUE IF OUTSIDE THE DOMAIN.

IF (NUM.NE.l) GO TO 10

TLU= ARRAY (2. 1

)

RETURN

10 IF (NUM.NE.2) GO TO 20

1= 2 GO TO 50

20 IF CVALIN.GT.ARRAY(i.l) ) GO TO 30

TLU= ARRAY (2. 1

)

RETURN

30 DO AO 1=2. NUM

IF <VALIN.LE.ARRAY(1

,

I)) GO TO 50

40 CONTINUE

TLU= ARRAY(2.NUM)

RETURN

50 TLU= ARRAY(2.I-1) + (VALIN - ARRAY( 1 , I- I ) )*

1 ((ARRAY(2tI} - ARRAY{2»I-1)) / (ARRAYd.I) - ARRAY ( 1 , 1- 1 ) ) ) RETURN

END

SUBROUTINE TRIDGF ( A, B« C* O. E, IX ) TRIDIAGONAL GAUSS ELIMINATION PROCEDURE FOR UNSYMMETRIC

MATRictS.

A=LEFT OF DIAGONAL. B=DI AGONAL. C=RIGHT OF DIAGONAL.

D= CONSTANT VECTOR. E= SOLUTION VECTOR. I X= SIZE OF MATRIX.

DIMENSION A(20).B(20).C(20) .D(20) .E(20) .CP( 20)

CP(l)= C(l )/B(l

)

Ed )= DdJ/BC I)

C(IX)= 0.

IXL= IX1

00 10 I=29lX

J= I-l

BX= B(IJ-CP( J>*A(I)

CP(I)= C(I)/BX

E(I)= (D(I>-E(J}«ACIJ)/BX

10 CONTINUE

DO 20 1=1. IXL

J= IX-I

E(J)= E( J)-E( J+1)*CP(J)

20 CONTINUE

RETURN

END

70
NBS-114A (REV. 9-76)

U.S. DEPT. OF tOMM. BIBLIOGRAPHIC DATA

SHEET

1. PUBLICATION OR REPORT NO. NBS Tech Note 991

2,(kiv*t AccesslOT No, 3, fi6Ctj}t6ftf4Acce$5^i<?nNo,

. .^v -^N VV-VS..V^^ .^.,< .^-V^ .. ^

4. TITLE AND SUBTITLE 5. Publication Date

0CMPF2—A Program for Calculating Post-Flashover F'ire

Tanperatures

June 1979

<, f^feft^'^trganization CkKle

7. AUTHOR(S)

Vytenis Babrauskas

8. Performing Organ. Report No.

9. PERFORMING ORGANIZATION NAME AND ADDRESS M. Fm}«ct/Task/Wort< Unit Ho.

NATIONAL BUREAU OF STANDARDS

DEPARTMENT OF COMMERCE

WASHINGTON, DC 20234

11. Contract/Grant No.

12. SPONSORING ORGANIZATION NAME AND COMPLETE ADDRESS rsfreet. City. St^te, ZIP)

Same as 9 above

13. Type of Report & Period Covered Final

'14. Sponsoring Agmcy Code

15. SUPPLEMENTARY NOTES

I i Document describes a computer program; SF-185, FIPS Software Summary, is attached.

Id. ABSTRACT (A 200-word or less factual summary of most significant information. If document includes a significant bibliography or

literature survey, mention it here.) COVIPF2 Is a computer prograiTi for calculatlng the characteristics of a post-flashover fire in a single building compartment, based on fire-induced

ventilation through a single door or window. It is intended both for performing

design calculations and for the analysis of experimental burn data. Wood, thermoplastic, and liquid fuels can be treated. In addition to the capability of performin;

calculations for ccmpartments with ccxnpletely determined properties, routines are included for calculating fire behavior by an innovative variable abstraction method.

A comprehensive output format is provided which gives gas tenperatures , heat flow

teiTOS, and flow variables. The documentation includes input instructions, sample

problems, and a listing of the program. The program is written in Fortran and

constitutes an improved version of an earlier program, OCMPF.

17. KEY WORDS (six to twelve entries; alphabetical order; capitalize only the first letter of the first key word unless a proper name;

separated by semicolons) Coriputer programs—fire protection; fire protection; fire resistance; fire tests;

fire walls; safety engineering—fires.

M. AVAILABILITY IXXUnlimited

For Official Distribution. Do Not Release to NTIS

fXiI Order From Sup. of Doc, U.S. Government Government Prinfino ri Office, Washington, DC

20402, SD Stock No. SN003-OO3-O2O8O-O

[J Order From National Technical Information Service (NTIS), Springfield, VA, 22161

19. SECURITY CLASS

(THIS REPORT)

UNCLASSIFIED

20. SECURITY CLASS

(THIS PAGE)

UNCLASSIFIED

21. NO. OF

PRINTED PAGES

76

22. Price$3.50

U5COMM-DC

There's

anew

look

/ y W^X? IX ... the monthly

y/ >^L ^^v / /^ magazine of the Nationy^^^^^^^tf^/^ al Bureau of Standards.

/^ ^^^^^^ X Still featured are special ar- ^^^^^k^^x tides of general interest on \^^ y / current topics such as consum- __\:^^^y / c product safety and building ^^^'^ / technology. In addition, new sections are designed to . . . PROVIDE

SCIENTISTS with illustrated discussions

of recent technical developments and

work in progress . . . iNFORM INDUSTRIAL

MANAGERS of technology transfer activities in Federal and private labs. . . DESCRIBE TO MANUFACTURERS advances in the field of voluntary and

mandatory standards. The new DIMENSIONS/NBS also

carries complete listings of upcoming conferences to be

held at NBS and reports on all the latest NBS publications,

with information on how to order. Finally, each issue carries a page of News Briefs, aimed at keeping scientist and consumalike up to date on major developments at the Nation's physical sciences and measurement laboratory.

(please delach here) SUBSCRIPTION ORDER FORM

Enter my Subscription To DIMENSIONS/NBS at S11.00. Add S2.75 for foreign mailing. No additional postage is required for mailing within the United States or its possessions. Domestic remittances

should be made either by postal money order, express money order, or check. Foreign remittances

should be made either by intefnational money order, draft on an American bank, or by UNESCO

coupons. Send Subscription to:

I NAME-FIRST, LAST I

Remittance Enclosed (Make checks payable

to Superintendent of Documents)

Charge to my Deposit Account No.

. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 COMPANY NAME OR ADDITIONAL ADDRESS LINE

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ! 1 1 1

STREET ADDRESS

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 > > 1 1 1 1 1 1

CITY J 1 1 1 1 M 1 1 1 1 1 1 1 1 1

STATE

1

ZIP CODE

1 1 1 1

MAIL ORDER FORM TO:

Superintendent of Documents

Government Printing Office Washington, D.C. 20402

PLEASE PRINT

NBS TECHNICAL PUBLICATIONS

PERIODICALS

JOURNAL OF RESEARCH—The Journal of Research

of the National Bureau of Standards reports NBS research

and development in those disciplines of the physical and

engineering sciences in which the Bureau is active. These

include physics, chemistry, engineering, mathematics, and

computer sciences. Papers cover a broad range of subjects,

j with major emphasis on measurement methodology, and

the basic technology underlying standardization. Also in- cluded from time to time are survey articles on topics closely

related to the Bureau's technical and scientific programs. As

a special service to subscribers each issue contains complete

citations to all recent NBS publications in NBS and non- NBS media. Issued six times a year. Annual subscription: domestic $17.00; foreign $21.25. Single copy, $3.00 domestic;

$3.75 foreign. Note: The Journal was formerly published in two sections:

' Section A "Physics and Chemistry" and Section B "Mathematical Sciences." DIMENSIONS/NBS

This monthly magazine is published to inform scientists, engineers, businessmen, industry, teachers, students, and

consumers of the latest advances in science and technology,

with primary emphasis on the work at NBS. The magazine

highlights and reviews such issues as energy research, fire protection, building technology, metric conversion, pollution

abatement, health and safety, and consumer product per- formance. In addition, it reports the results of Bureau pro- grams in measurement standards and techniques, properties

of matter and materials, engineering standards and services, instrumentation, and automatic data processing.

Annual subscription: Domestic, $11.00; Foreign $13.75

NONPERIODICALS

Monographs—Major contributions to the technical liter- ature on various subjects related to the Bureau's scientific and technical activities. Handbooks—Recommended codes of engineering and industrial practice (including safety codes) developed in cooperation with interested industries, professional organizations, and regulatory bodies.

Special Publications—Include proceedings of conferences

sponsored by NBS, NBS annual reports, and other special

publications appropriate to this grouping such as wall charts, pocket cards, and bibliographies.

Applied Mathematics Series—Mathematical tables, manuals, and studies of special interest to physicists, engineers,

chemists, biologists, mathematicians, computer programmers,

and others engaged in scientific and technical work.

National Standard Reference Data Series—Provides quantitative data on the physical and chemical properties of

materials, compiled from the world's literature and critically evaluated. Developed under a world-wide program coordinated by NBS. Program under authority of National

Standard Data Act (Public Law 90-396).

NOTE: At present the principal publication outlet for these data is the Journal of Physical and Chemical Reference Data (JPCRD) published quarterly for NBS by the American Chemical Society (ACS) and the American Institute of

Physics (AIP). Subscriptions, reprints, and supplements

available from ACS, 1155 Sixteenth St. N.W., Wash., D.C.

20056.

Building Science Series—Disseminates technical information developed at the Bureau on building materials, components,

systems, and whole structures. The series presents research

results, test methods, and performance criteria related to the

structural and environmental functions and the durability and safety characteristics of building elements and systems. Technical Notes—Studies or reports which are complete in themselves but restrictive in their treatment of a subject.

Analogous to monographs but not so comprehensive in scope or definitive in treatment of the subject area. Often

serve as a vehicle for final reports of work performed at NBS under the sponsorship of other government agencies.

Voluntary Product Standards—Developed under procedures

published by the Department of Commerce in Part 10, Title 15, of the Code of Federal Regulations. The purpose

of the standards is to establish nationally recognized require- ments for products, and to provide all concerned interests with a basis for common understanding of the characteristics of the products. NBS administers this program as a supple- ment to the activities of the private sector standardizing

organizations.

Consumer Information Series—Practical information, based on NBS research and experience, covering areas of interest

to the consumer. Easily understandable language and

illustrations provide useful background knowledge for shopping in today's technological marketplace.

Order above NBS publications from: Superintendent of Documents, Government Printing Office, Washington, D.C.

20402. Order following NBS publications—NBSIR's and FIPS from

the National Technical Information Services, Springfield, Va. 22161.

Federal Information Processing Standards Publications (FIPS PUB)—Publications in this series collectively consti- tute the Federal Information Processing Standards Register.

Register serves as the official source of information in the Federal Government regarding standards issued by NBS

pursuant to the Federal Property and Administrative Services Act of 1949 as amended. Public Law 89-306 (79 Stat. 1127), and as implemented by Executive Order 11717

(38 FR 12315, dated May 11, 1973) and Part 6 of Title 15 CFR (Code of Federal Regulations).

NBS Interagency Reports (NBSIR)—^A special series of

interim or final reports on work performed by NBS for outside sponsors (both government and non-government).

In general, initial distribution is handled by the sponsor;

public distribution is by the National Technical Information

Services (Springfield, Va. 22161) in paper copy or microfiche form.

BIBLIOGRAPHIC SUBSCRIPTION SERVICES

The following current-awareness and literature-survey bibli- ographies are issued periodically by the Bureau:

Cryogenic Data Center Current Awareness Service. A litera- ture survey issued biweekly. Annual subscription: Domestic, $25.00; Foreign, $30.00.

Liquified Natural Gas. A literature survey issued quarterly. Annual subscription: $20.00.

Superconducting Devices and Materials. A literature survey

issued quarterly. Annual subscription: $30.00. Send subscription orders and remittances for the preceding bibliographic

services to National Bureau of Standards, Cryogenic Data

Center (275.02) Boulder, Colorado 80302.
U.S. DEPARTMENT OF COMMERCE

National Bureau of Standards

Washington, D.C. 20234

OFFICIAL BUSINESS

Penalty for Private Use, S300

POSTAGE AND FEES PAID

U.S. DEPARTMENT OF COMMERCE COM-215

^

SPECIAL FOURTH-CLASS RATE BOOK