10.08.2013 Views

Code Manual for CONTAIN 2.0 - Federation of American Scientists

Code Manual for CONTAIN 2.0 - Federation of American Scientists

Code Manual for CONTAIN 2.0 - Federation of American Scientists

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

where N is the number <strong>of</strong> the surfaces in the enclosure, Fj~is the user-input view factor from surface<br />

j to surface k (the fraction <strong>of</strong> radiation from surface j that is directed into the solid angle subtended ~<br />

by surface k), e~J~is the emissivity <strong>of</strong> gas mixture, &,Jis the surface emissivity, and B~= o(T~)4is<br />

the Planck blackbody radiation <strong>for</strong> the gas. The modeling <strong>of</strong> e~~~,which is a function <strong>of</strong> the userinput<br />

beam length Lj~between surfaces j and k, is discussed in Section 10.3.3.<br />

The dry surface ernissivities, view factors Fj~,and characteristic beam lengths Lj~between surfaces<br />

are provided by the user as “emsvt,” “vufacn,” and “beaml,” respectively, in the RAD-HEAT input<br />

block described in Section 14.3.1.5. When the surface is dry, the “emsvt” value is used <strong>for</strong> &,,j.<br />

Whenever a water fdm is present on the surface, 0.94, the emissivity <strong>of</strong> water, is used instead <strong>of</strong> the<br />

dry surface emissivity. Because the surface and gas temperatures <strong>of</strong> the enclosure are known, the<br />

radiosity <strong>of</strong> each surface can be obtained from Equation (10-60) through use <strong>of</strong> a standard linear<br />

equation solver. The flux q~~,is then obtained from Equation (10-58).<br />

If the user has not input the beam length ~~ between surfaces <strong>for</strong> the net enclosure model, a default<br />

value given by<br />

3.6V;<br />

L= (10-61)<br />

5A.<br />

i=l<br />

5,1<br />

is used, where ~ is the initial cell gas volume, and Asi is the area <strong>of</strong> surface i.<br />

For the view factor between surfaces, the user should only supply the view factor in one direction<br />

(say, from surface i to j). The code will then calculate the view factor in the opposite direction (from<br />

surface j to i) according to the following equation:<br />

A, i Fij = As,j Fji (10-62)<br />

where ~J is the area <strong>of</strong> surface j, Fijis the view factor from surface i to j, and Fjiis the view factor<br />

from surface j to i.<br />

10.3.2 Direct Radiative Heat Transfer Model<br />

The direct radiative heat transfer model in <strong>CONTAIN</strong> calculates the radiative heat transfer between<br />

the atmosphere and each heat transfer structure, between the atmosphere and the uppermost lower<br />

cell layer, and between the uppermost lower cell layer and each structure inner surface. The name<br />

<strong>of</strong> this model arises from the fact that secondruy reflections involving three or more enclosure<br />

elements are treated in an approximate fashion, in contrast to the net enclosure radiation model. A<br />

further limitation is that direct radiation between structure surfaces is not modeled. Unlike the net<br />

enclosure radiation model, the direct model permits radiative heat transfer from the atmosphere to<br />

both the inner and outer surfaces <strong>of</strong> a structure if both lie within a ceil. The direct model <strong>for</strong> gas to<br />

surface heat transfer is invoked by speci~ing either GASWAL or GEOBL in the RAD-HEAT input<br />

block. The modeling <strong>of</strong> direct radiative exchange between the lower cell and each heat transfer ~<br />

O 10-38 6/30/97

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!