Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

On The Design of Thermal Breathing

Download as pdf or txt
Download as pdf or txt
You are on page 1of 6
At a glance
Powered by AI
The paper presents a simulation procedure and design criteria for sizing thermal breathing devices for liquid storage tanks that considers non-uniform tank temperatures and vapour condensation.

The analysis assumes an initial hot, sunny condition that changes to a sudden drop in temperature, solar radiation and heavy rain at time t=0, causing the tank enclosures' temperatures to decrease stepwise.

The simulation procedure considers a multi-zone representation of the tank to account for non-uniform temperatures, and includes the role of vapour condensation on the tank's dynamic behavior.

0957 5820

Institution of Chemical Engineers


Open access under CC BY-NC-ND license.
Trans IChemE, Vol 77, Part B, Novem ber 1999

ON THE DESIGN OF THERMAL BREATHING


DEVICES FOR LIQUID STORAGE TANKS
P. SALATINO*, G. VOLPICELLI* and P. VOLPE**
*Dipartimento di Ingegneria Chimica, Universita degli Studi di Napoli Federico II, Napoli, Italy
**Kuwait Raf nazione e Chimica S.p.A., Napoli, Italy

T
he dynam ic behaviour of an atmospheric liquid storage tank upon sudden change of
meteorological conditions is analysed with the aim of establishing criteria for the safe
and effective design of thermal breathing devices. A simulation model is presented
based on a lumped-parameter multiple-zone representation of the tank. Simpli ed design
criteria based on reasonable and conservative approxim ations of the governing equations are
derived. Application of the simulation procedure and of design criteria is exempli ed. Results
are discussed and critically compared with prescriptions of existing codes.
Keywords: storage tank; thermal breathing; dynamic simulation.

INTRODUCTION of the tank enclosures. The role of vapour condensation


Venting of atmospheric liquid storage tanks is normally on the dynamic behaviour of the tank is also critically
required in order to maintain pressure inside the tank at the considered. A simpli ed design criterion is obtained on
atmospheric level upon occurrence of one of the following the basis of approximations of the governing equations.
events: Application of the analysis is exempli ed and results
compared with prescriptions of design criteria provided
feeding/withdrawal of liquid streams to/from the tank; by the API code 1 , by the Naumann and PTB Formulas and
changes in the density of gases/vapours in the tank by the draft EN standard 5 . Sensitivity of the procedure to
associated with changes in meteorological conditions (either selected variables is presented and discussed.
temperature or pressure);
changes in density of gases/vapours in the tank due to
DEVELOPMENT OF A SIMULATION PROCEDURE
overheating associated with exposure to accidental res.
The Reference Scenario
General design criteria for venting of storage tanks
are provided by the ANSI/API Standard 2000-1992 .
1 A schematic representation of the tank is given in
Prescriptions cover either minimum thermal venting Figure 1. The temperature TG of bulk gases in the tank
capacities or emergency venting to be provided in case of has been assumed uniform throughout the vessel at any
accidental re exposure. The prescribed thermal inbreathing time. The temperature TL of the liquid in the tank has been
capacity is evaluated (Table 2 in Reference 1) on the basis of assumed constant over the time and averaged over the
a maximum estimated rate of heat loss (20 Btu per hour and period considered.
per square foot of tank surface) upon sudden decrease of the The following scenario has been assumed as a reference:
external temperature (as might happen, for instance, after a
summer rainstorm). Maximum thermal outbreathing capa- Phase 1
city is rated at 60% of the maximum inbreathing capacity. The weather is initially hot and sunny. The tank is
Alternative design criteria for venting devices are provided exposed to solar radiation and to average ambient
by the Naumann Formulas, developed by ESSO, and by the temperature TA ,1 .
German PTB-TRbF Formulas. Hoechst Formulas, reported
2
by Sigel et al. , take into account the heat capacity of the
tank metal enclosures when heat transfer rate between
the tank enclosures and the environm ent is nite. Fullarton
and co-workers 3,4 embodied consideration of vapour con-
densation in the assessment of venting requirements. A draft
5
European normative has recently been developed .
All the cited methods are based on the assumption that
the enclosures of the gas in the tank are at uniform tem-
perature. In the present paper a design procedure for thermal
venting devices for atmospheric liquid storage tanks is
presented which takes into account temperature non-
uniformity by the recourse to a multiple-zone description Figure 1. Schematic representation of the storage tank.

354
THERMAL BREATHING DEVICES FOR LIQUID STORAGE TANKS 355

Phase 2 coef cients between the roof/shell and ambient air and
At time t = 0 weather conditions change. The combined the ambient temperature TA,1 are known.
effect of sudden drop of solar radiative ux incident on The energy balance on the gas within the tank, under
the tank roof, decrease of ambient temperature to TA,2 pseudo-steady state conditions, reads:
and strong downpour are assumed to determine stepwise
decrease of the enclosures temperature. The temperature
0 = AR hRG (TG,1 2 TR,1 ) + AS hSG(TG,1 2 TS,1 )
of gases and vapours in the tank starts decreasing, and so + AL hLG (TG,1 2 TL ) (4 )
does the pressure. As pressure reaches a threshold vacuum
level (that depends on inbreathing valves characteristic where AR, AS and AL are the areas of the surfaces enclosing
curve), inbreathing starts and ambient air is admitted into the gas, i.e. the roof/shell internal surface and the gas-liquid
the tank. interface. hLG is the heat transfer coef cient between the
gas and the bulk of the liquid.
Solution of equation (4) with respect to TG yields:

=S
Evaluation of Variables Before Change of Ai hiG Ti
Weather Conditions TG (5 )
S Ai hiG
The energy balance over a surface element of the tank
enclosure (assuming that temperature is uniform across the where summation is extended over all the surfaces enclosing
metal sheet) yields: the gas.
dT
re ce s = ks$ 2
T + qe + hA (TA 2 T ) + hG (TG 2 T)
Evaluation of Variables After Change of
dt
(1 ) Weather Conditions
In equation (1) T is the element temperature, $ 2 T is its At t = 0 it is assumed that sudden change of weather
Laplacian, re ce and k the heat capacity per unit volum e conditions associated with: (a) vanishing of radiative
and the thermal conductivity of the metal sheet, respec- ux absorbed by the roof/shell; (b) decrease of ambient
tively, s its thickness, q is the incident radiative ux, e is the temperature and, possibly, (c) downpour, occurs. Under
surface emissivity, hA and hG the heat transfer coef cients such conditions it is likely that roof/shell temperatures
of the element with the ambient and the gases in the tank, might approach, over a relatively limited time-scale, the
respectively. The latter embody contributions from both wet bulb temperature corresponding to the new ambient
radiation and convection (either free or forced). values of temperature and humidity. Due to the dif culty
Order of magnitude calculations based on typical of assessing the actual heat transfer rate between the gas
parameters of storage tanks suggest that the rst term at enclosures and the environm ent, it is (conservatively)
RHS of equation (1) is much smaller than the others, i.e., assumed that stepwise drop of roof/shell temperatures to
conduction along the metal enclosures can usually be the values TR,2, TS,2 occurs upon change of weather
neglected. This implies that spatial variation of T is related conditions. Correspondin gly and in the absence of inbreath-
only to change of parameters qe , hA and hG. Provided that ing the temperature of the gas/vapours in the tank decreases
these parameters can be assumed locally uniform , the towards an ultimate value TG ,2 that can be calculated from
enclosure can be lumped into multiple zones at piecewise equations (2), (3) and (5) referring to the new conditions.
uniform temperature. This approach was followed in the In the absence of inbreathing and neglecting conden-
present study, and the enclosures lumped into the roof sation of vapours at temperatures below the dew point, the
and the shell, each characterized by uniform values of maximum vacuum level in the tank would be:

( )
parameters TR, qR e , hRA, hRG and TS, qS e , hSA, hSG, respec- TG,2
tively. The liquid free-surface, at temperature TL, completes | p2 2 p1 | = p1 12 (6 )
TG,1
the enclosure of the gas in the tank.
Analysis of equation (1) further indicates that the term at The actual vacuum level established in the tank under
LHS is typically much smaller than the others as far as the combined effect of gas cooling and of inbreathing has
heating (or cooling) rates associated with daily excursion the value given by equation (6) as the upper limit. The actual
of ambient parameters are considered. This implies that, maximum vacuum level in the tank can be calculated by
as far as the dynamic behaviour of the tank over the time- following the dynam ic behaviour of the system.
scale of a day is concerned, the pseudo-steady state
approximation can be invoked.
In the light of the above approximations, the pseudo- Simulation of the Dynamics of the Tank upon
steady state roof and shell temperatures can be obtained Sudden Change of Boundary Conditions
by solution of the following energy balances: The energy balance on the gas in the tank, under transient
qR e = hRA (TR,1 2 TA,1 ) + hRG (TR,1 2 TG,1 ) (2 ) conditions, modi es into:

qS e = hSA (TS,1 2 TA,1 ) + hSG (TS,1 2 d TG


TG,1 ) (3 ) nc =S Ai hiG (TG 2 Ti,2 ) + nc(T
d A,2 2 TG )
dt
The last term in both equations (2) and (3), representing
dp
the heat ux from the roof/shell to the gas inside the tank, + V dt (7 )
is generally negligible. Equations (2) and (3) in the
unknowns TR,1 and TS,1 are solved once the radiative where n is the number of moles of gas/vapour in the tank,
heat uxes absorbed by the roof/shell, the heat transfer V the volume occupied by gas/vapours, c the gas molar

Trans IChemE, Vol 77, Part B, November 1999


356 SALATINO et al.

speci c heat at constant pressure, p the absolute pressure differentiation of the equation of state:
in the tank. Summation in equation (7) is again extended d p RTG nd Rn T
over all the i-th surfaces enclosing the gas. The material = + V d tG = pQ + p d TG
(16 )
balance on gas in the tank yields: dt V d V TG d t
dn where Q is the volum etric ow rate of the inbreathed air.
= nd (8 ) In order to prevent pressure from decreasing, the inbreath-
dt ing valves must provide at least the air ow rate required
where nd , the molar ow rate of gas entering the tank via the to ensure that dp/d t = 0, that is:
inbreathing valves, depends on the vacuum level inside
V d TG
the tank according to characteristic curves of the inbreathing Q$ 2 (17 )
valves: TG d t
nd = n(
d p 2 pA ) (9 ) At the beginning of the cool-down process the inbreathing
required is maximum. Particularization of equation (17) at
Pressure in the tank is related to the number of moles of t = 0 yields:
gas by the ideal gas equation of state:
nRT Q$
V 1
(TG,1 2 TG,2 ) =
S Ai hiG R
(TG,1 2
p= (10 ) TG,1 t c p
TG,2 )
V
Simulation of the transient behaviour requires integration (18 )
of the two ordinary differential equations (7) and (8), with
initial conditions:
p1 V APPLICATION TO DESIGN/VERIFICATION:
for t = 0, TG = TG,1 and n = AN EXAMPLE
RTG,1
Application of the above procedure is exempli ed in the
following. A cone-roof tank characterized by the following
geometrical parameters is considered: diameter: 70 m;
DEVELOPMENT OF A SIMPLIFIED
height: 15 m; liquid storage capacity: 56,000 m 3; roof
DESIGN CRITERION 3
tangent: 0.167; total tank volume: 63,000 m . The tank is
A useful approximation of equation (7) can be developed equipped with four breather valves whose characteristic
when considering that the inbreathing valves have to curve is reported in Figure 2. It is assumed that at the time of
provide enough air for tank pressure not to change appre- the sudden change of weather conditions the tank is almost
ciably. If the pressure derivative related term is dropped completely empty.
in equation (7), and if the further simpli cation is made The evaluation of input parameters to the simulation
of neglecting the term associated with the enthalpy of procedure is detailed in the Appendix. It is: qR =
incoming air, equation (7) becomes: 800 Wm2 2 ; e = 0 .35; TA,1 = 309 K; TL = 298 K. Further-
more it is assumed that the radiative heat ux incident
d TG
nc = S Ai hiG (TG 2 Ti,2 ) on the shell is negligible and that shell temperature can be
dt taken equal to the ambient temperature (TS,1 = TA,1 =
= TG S 309 K). hRA in equation (2) sums up contributions from
Ai hiG 2 S Ai hiG Ti,2 (11 )
convection, either free or forced, and radiation. The radia-
or: tive contribution accounts for about 5 W m 2 K 1. Free
d TG convection regime is assumed and the convective contri-
t = TG 2 TG,2 (12 ) bution is 5 W m K . Accordingly hRA = 10 W m 2 K2 .
2 1 2 1
dt
Equation (2) yields the equilibrium roof temperature
where TG,2 has already been de ned and t is a time constant TR,1 = 337 K (64C) before change of weather conditions.
Assuming: AR = 3900 m ; AS = 3200 m ; AL = 3800 m ;
2 2 2
given by:
nc pV c
t = = (13 )
S Ai hiG RTG S A i hG
The solution of equation (12) is straightforward and
yields:
t
TG,2 )[ 1 2
TG = TG,1 2 (TG,1 2
( t) ]
exp 2 (14 )

From this equation the meaning of the characteristic time


t can be appreciated: it is the time-scale over which gas
temperature decreases from the value TG ,1 to approach the
ultimate (equilibrium ) gas temperature TG ,2 .
The initial rate of change of gas temperature is:

d TG 1
=2 (TG,1 2 TG,2 ) (15 )
dt t
The rate of change of pressure in the tank is, from Figure 2. Characteristic curve of inbreathing valves.

Trans IChemE, Vol 77, Part B, November 1999


THERMAL BREATHING DEVICES FOR LIQUID STORAGE TANKS 357

hRG = hSG = hLG = 3 W m2 2 K2 1 , equation (5) yields


TG,1 = 315 K (42C).
Immediately after the change of weather conditions both
ambient and roof/shell temperatures change. It is assumed
TR,2 = TS,2 = TA,2 = 293 K. From equation (5), and assum-
ing hRG = hSG = 5 W m2 K2 ; hLG = 3 W m2 K2 , it is
2 1 2 1

TG,2 = 294 K (21C).

Dynamic Simulation
Results of dynam ic simulations are presented in Figures 3
to 6. They report the pressure in the tank (Figure 3), the
volum etric ow rate of inbreathed air (Figure 4), the tem-
perature of gas in the tank (Figure 5) and the total heat ux
from the roof/shell to the ambient (Figure 6) as functions
Figure 4. Results of numerical simulation: inbreathed ow rate vs time.
of time. The equilibration of temperature in the tank with
the new ambient temperature takes place over a time
interval of about one hour. Table 1 (column 1) reports the unchanged. Results are compared with those of the base
main parameters of the calculations, namely the maxi- case in Table 1 (column 3). It can be noted that maximum
mum inbreathing rate, the maximum vacuum level and the underpressure rises to 48 mm w.g., while maximum load
maximum heat ux from the tank to the atmosphere after on inbreathing valves is in the order of 13,000 m 3 hr 1.
change of weather conditions. A maximum vacuum level The analysis crucially relies on the knowledge of gas-
of about 40 mm w.g. is reached shortly after the beginning roof/shell and of gas-liquid heat transfer coef cients, whose
of cool-down. Correspondin gly, maximum ow rate of determination is at present the most critical aspect of the
inbreathed air is 9600 m 3 hr2 1 . Maximum heat ux from whole procedure. Values adopted in the calculations were
the tank to the atmosphere is in the order of 110 W m 2. This determined on the basis of literature correlations for free-
gure, which corresponds to about 35 BTU/hr sqft, is far convective heat transfer from at surfaces6 . The validity
larger than the value (20 BTU/hr sqft) that provides the of these equations, however, was extrapolated well beyond
basis for the ANSI/API 2000 standard 1. the ranges of Nu and Gr numbers that provide their
Computations were repeated considering the case in experimental background, for lack of better suited correla-
which the tank is half lled. Correspondingly, AS = 1600 m
2
tions. In order to check the sensitivity of the simulation
and the total volum e occupied by gas/vapours would be procedure to the values of the heat transfer coef cient hRG,
V = 35 , 000 m 3 . Results, summarized in Table 1 (column 2), additional computations were carried out by taking twice
indicate that maximum vacuum level and inbreathing the value used in the base case calculation. Results of this
rate would be about 37 mm w.g. and 8400 m 3 hr 1 , respec- sensitivity analysis are reported in Table 1 (column 4). It
tively. This situation turns out to be, therefore, less critical is observed that maximum vacuum level increases to about
than that occurring in the case of an empty tank. 48 mm w.g.. Correspondingly inbreathing demand is
Results reported so far are relative to a computed initial 3
raised to more than 13,000 m hr .
1

gas temperature of TG,1 = 42C, with gas temperature


drop after change of weather conditions of 21C. The API
design directives [Appendix A of Reference 1] indicate The In uence of Vapour Condensation
that temperature drops of roof plates as large as 60F (33C) Computations have been developed so far neglecting
can be observed upon change of weather conditions. In vapour condensation as roof/shell temperatures are brought
order to check the sensitivity of the model to the gas tem- below the dew point of the gas/vapour mixture. Fullarton
perature drop, further calculations have been carried out
assuming TG,1 = 50C, while leaving the other parameters

Figure 5. Results of numerical simulation: temperature of gas/vapours in


Figure 3. Results of numerical simulation: vacuum level vs time. the tank vs time.

Trans IChemE, Vol 77, Part B, November 1999


358 SALATINO et al.

Simpli ed Design Criterion


The design criterion expressed by equation (18) has been
applied to the example. Using temperatures and heat trans-
fer coef cients reported above, the characteristic time for
temperature decay is t = 25 min. Correspondingly the
maximum volumetric ow rate of inbreathed air required
is Q = 10 , 000 m hr2 and is reported in Table 1 (column
3 1

1). The gure obtained with the approxim ated criterion


equation (18) compares very satisfactorily (within a 4%
error) with that obtained after numerical simulation of the
tank dynamics, and it lies on the safe side. Similar compu-
tations have been carried out with reference to the other
cases considered in Table 1. It is suggested that the design
criterion expressed by equation (18), with t given by
Figure 6. Results of numerical simulation: heat ux from the roof/shell to equation (13) and TG1 , TG,2 by equation (5), be used to
the atmosphere vs time. evaluate the maximum load on inbreathing valves.
When developing a conservative design criterion the
3,4 most critical case, i.e. that corresponding to the empty tank,
and co-workers addressed the in uence of vapour
condensation on the establishment of vacuum in storage is to be considered. Equation (18) can be further simpli ed
tanks. They recognized the existence of two mutually if an average value is assumed for the heat transfer
coef cients hiG . Accordingly, equation (18) yields:
opposing phenom ena, namely the release of the heat of
condensation at the wall, which decreases the wall cooling Q = K A (19 )
1
velocity, and the additional volum e ow rate associated where the constant K takes the value 1 m hr if it is assumed
with the condensation ux. The rst effect is ruled out in the hiG > 4 W m 2 2 K2 1 and TG1 2 TG2 > 25C. The area A of
present study as the metal enclosure temperature is assumed surfaces enclosing gas/vapours can be simply expressed
to change stepwise upon change of weather conditions. in this case in terms of the volume V of the tank (supposed
Assessment of the maximum extent of vapour condensation cylindrical) and of the aspect ratio r between the tank
is based on the evaluation of the molar fraction y of height H and its diameter D. It is:
condensables in the gas phase prior to change of weather 1 + 2r
conditions. It is worth noting that the value of y should be Q = K V 2 /3 3 W W (20 )
upperly limited by saturation at the lowest temperature r 1
of the enclosures (that would also be the upper limit of the With the above assumptions K takes the value 2.6 m hr .
gaseous mixture dew point) and not at the gas tempera- Application of equation (20) to the design/veri cation
ture TG (this has been a common assumption in the of inbreathing devices for the tank considered in the case
2 4
literature ). The maximum amount of condensables study is exempli ed in Table 2. Its results are compared
present in the gas phase prior to change of weather with those obtained from the application of equation (18)
conditions should be evaluated accordingly. With reference as well as from direct numerical simulation. Results
to the case in hand and assuming water vapour as the only obtained from the application of the ANSI/API Code, of
condensable specie, the maximum dew point of the mixture the PTB-TrbF and EN Formulas (both based on expressions
would be TL = 25C and vapour condensation as gas is
m
of the type Q 3 V , m being 0.71 and 0.70, respectively)
cooled to TG,2 = 21C would account for no more than 10% and of the Naumann Formulas (based on an expression
of the type Q 3 D (1 + 4 r)) are also reported for com-
2
of the total inbreathing demand. By no means should this
nding be considered a general conclusion. parison. It can be noted that the ANSI/API prescriptions
In addition to the effects discussed above, some enhance- signi cantly underestimate the maximum inbreathing
ment of heat transfer between gas and the enclosures demand with respect to results of the present analysis. On
might be expected as a consequence of condensation: it the other hand, both the PTB-TrbF and the Naumann
7
has been observed that heat transfer can be dramatically Formulas appear to give conservative prescriptions. Appli-
augmented by occurrence of phase change even in the cation of the draft EN code yields values of the inbreathing
presence of only small amounts of condensables. According demand that can be either smaller or larger than those of
to equation (18), this would lead to larger inbreathing the present analysis depending on the value of the constant
requirements. C adopted in the computation.

Table 1. Results of numerical simulation for the case study.

1 2 3 4
base case half- lled tank TG,1 = 50C hRG = 10 W m 2 2
K2
1

Maximum inbreathing ow rate , m h 2


1 3 1
9600 8400 13300 13100
(10,000) (8650) (13,800) (14,500)
Maximum vacuum level, mm w.g. 2 40 2 37 2 48 2 48
Maximum heat ux, W m2
2
111 116 150 172
1
The inbreathing ow rate calculated according to equation (18) is given in parentheses.

Trans IChemE, Vol 77, Part B, November 1999


THERMAL BREATHING DEVICES FOR LIQUID STORAGE TANKS 359

Table 2. Maximum thermal inbreathing ow rate for the case study (empty Heat transfer coef cient, free convection, horizontal
tank) calculated according to different design criteria. plate, cold-facing up or hot-facing down: 60% of the
Thermal inbreathing ow rate, value calculated with the previous equation.
Design criterion m hr 2
3 1
Radiative heat transfer coef cient from the roof to the
1
atmosphere: h = 4 seT 3R .
ANSI/API Code 4100
Naumann Formula 12,800
PTB-TRbF Formula
5
11,300 NOMENCLATURE
pr EN 265001 6300 10,500
Present work dynamic simulation 9600 A area
Present work equation (18) 10,000 c molar speci c heat at constant pressure
Present work equation (20) 9900 ce speci c heat of wall material
D tank diameter
h heat transfer coef cient
H tank height
k thermal conductivity of wall material
K proportionality constant
CONCLUSIONS n number of moles of gas in the tank
Detailed simulation of the dynamics of an atmospheric nd molar rate of air inbreathing
liquid storage tank upon sudden change of weather condi- p pressure
q radiative ux
tions has been carried out by means of a simple lumped- Q volumetric ow rate of inbreathed air
parameter model. Simpli ed criteria for the design of r tank aspect ratio (= H/D)
inbreathing valves have been developed, based on reason- R gas constant
able and conservative approximations. s tank wall thickness
t time
The dynamics of a cone-roof liquid storage tank subject
T temperature
to a sudden change of weather conditions has been simu- V gas volume in the tank
lated. Maximum inbreathing capacities and vacuum levels y molar fraction of condensables in the gas
calculated by direct numerical solution of the governing e emissivity
equations compare well with indications coming from the re density of wall material
t characteristic time for gas temperature decay
simpli ed design criteria. The simpli ed criterion provides
a reasonably good, slightly conservative tool for the design Subscripts
of inbreathing devices, easier to use than the direct A ambient
numerical simulation. Results of the analysis are extremely G gas
L liquid
sensitive to the values of the heat transfer coef cients
LG between liquid and gas
assumed in the computations. Their evaluation is at present R roof
the most critical point of the evaluation procedure. RA between roof and ambient
Analysis of results relative to the case study indicates RG between roof and gas
that rates of heat transfer from the tank to the environment S shell
SA between shell and ambient
well in excess of 20 BTU/sqft hr (i.e., the value assumed SG between shell and gas
by the ANSI/API Standard 1) might establish after sudden 1,2 before, after change of weather conditions
change in weather conditions. Inbreathing requirements
calculated as per the Naumann and PTB Formulas are
slightly larger than those predicted with the design criterion REFERENCES
5
developed in this paper. The draft EN prescription might 1. Venting Atmospheric and Low-Pressure Storage Tanks, API Standard
not lie on the safe side. Application of the ANSI/API 2000, 4th edition, 1992 (American Petroleum Institute, Washington
D.C.).
code appears to be not conservative enough, and consistent 2. Sigel, R., Kuxdorf, B, Meiss, R. and Schwarz, H., 1983, Chem Ing Tech,
oversizing of the thermal breathing devices in respect to 55: 377.
its prescriptions should be considered as a safer design 3. Fullarton, D., 1986, Chem Ing Tech, 58: 38.
measure. 4. Fullarton, D., Evripidis, J. and Schlunder, E.U., 1987, Chem Eng
Process, 22: 137.
5. Final Draft of Annex K Venting Systems for pr EN 265001 Part 1,
APPENDIX 1998.
6. McAdams, W. H., 1954, Heat Transmission, 3rd edition (McGraw Hill,
Evaluation of Heat Fluxes and New York).
Heat Transfer Coef cients 7. Kern, D. Q., 1965, Process Heat Transfer (McGraw Hill, New York).
8. Guida al Controllo Energetico della Progettazione, 1985 (Consiglio
Total radiant heat ux over a horizontal surface: Data Nazionale delle Ricerche, Rome).
published by C.N.R. . The gure corresponds to 42 latitude
8
9. Kreith, F., 1973, Principles of Heat Transfer (Donnelley, New York).
N, summer, noon.
Emissivity of the roof surface: 0.35 is emissivity to solar
radiation, 0.9 the emissivity for heat ux leaving the surface ADDRESS
at TR = 50C 9. In either case the roof is assumed to be Correspondence concerning this paper should be addressed to
painted in white. Professor P. Salatino, Dipartimento di Ingegneria Chimica, Universita
degli Studi di Napoli Federico II, P. le Tecchio 80-80125 Napoli, Italy.
Heat transfer coef cient, free convection, horizontal
plate, hot-facing up or cold-facing down: Nu = 0 .14 The manuscript was received 12 March 1999 and accepted for
(Gr Pr)1 /3 6 publication after revision 21 September 1999.

Trans IChemE, Vol 77, Part B, November 1999

You might also like