Transcription

1 BUOYANCY-DRIVEN NATURAL VENTILATION OP ENCLOSED SPACES M. R. Marshall* and P. L. Stewart-Darling* A simple mathematical model for the buoyancy driven ventilation of an enclosed space, using a two-pipe ventilation system, is described. Experimental data which support the predictions of this model are presented and practical means of achieving recommended ventilation levels are discussed. (Key words: ventilation, buoyancy, natural gas, confined spaces). INTRODUCTION The accidental release of a flammable gas or vapour from process plant housed in a building could result in the formation of potentially hazardous mixtures unless means of preventing such a situation developing are available. An effective measure is to provide sufficient ventilation air so that any credible accidental leakage can be dispersed safely. Normally, the ventilation air will be provided by natural ventilation and in most cases this will derive from the pressure differentials created across a building by the effect of wind. However, in some circumstances, natural ventilation will be dependent on buoyancy forces derived from the difference in densities of the atmosphere within and outside a building. Por example, at low wind speeds (i.e. below about 1ms -1 ) provided that ventilator openings are suitably positioned, the natural ventilation of a building will be provided mainly by buoyancy forces rather than wind effects. In addition, for the particular case of a below ground installation the provision of adequate natural ventilation will depend almost entirely on the effective use of available buoyancy forces. Normally the density difference that promotes air movement due to buoyancy will result from a difference in temperature between the atmospheres within and outside a building but the *British Gas Corporation, Midlands Research Station, Wharf Lane, Solihull, West Midlands B91 2JW. 7

2 density difference can also be caused by the release of a buoyant gas within a building. A simple mathematical model of buoyancy driven ventilation, based on a two-pipe ventilation system, has been developed in which the pressure differential to cause air movement is derived from the release of a buoyant gas within an enclosure. A design formula based on this simple model is used to determine the ventilator area requirement of underground pits and similar installations operated by British Gas. Use of this formula should ensure that the levels of ventilation achieved are adequate to disperse 3afely any credible natural gas leakages that occur in these installations. This paper presents the results of an experimental study of buoyancy driven ventilation (caused by the release of a buoyant gas) and the application of the data to the design of an effective natural ventilation system is discussed. MODEL OP BUOYANCY-DRIVEN VENTILATION A relationship has been derived between ventilator area, gas leakage rate and steady state gas concentration on the basis that the driving force for air movement is provided by the presence of a layer of buoyant gas-air mixture within an enclosure. The situation considered was that of an underground pit having a two-pipe ventilation system, the air inlet pipe terminating close to the floor of the pit and the outlet pipe being located at high level. (Figure 1). In deriving the relationship, it was assumed that gas released in the enclosure would produce a layer of gas-air mixture of uniform concentration extending downwards from the top of the enclosure to the level of the point of leakage (1). On this basis, a layer depth of 0-4d (where d is the depth of the enclosure) was taken to be representative. The flow of mixture out of the enclosure is given by (2): Q M = C d A eff (2AP//>m)* where, P =V Gh o and (D 1 / A Jff 1 / A i 2 n + 1 /Ag ut Optimum ventilation air flow is achieved when k± ^out and for this condition, A eff = A in A/2 = A out / /2 8

4 EXPERIMENTAL In order to test the validity of the simple model of buoyancy driven ventilation, experiments using releases of natural gas were conducted in several vessels, which differed in both volume and shape. All of these were fitted with a two-pipe ventilation system as shown diagramatically in Figure 1. To simulate practical leakage scenarios, the natural gas was released into the test enclosures as a pressurised jet, via small nozzles ranging in diameter from 0.6mm to 4-Omm. The influence on steady state gas concentrations of the gas delivery pressure and the leak orientation, as well as the actual gas release rate, were investigated. The majority of the experiments were conducted with a vessel having dimensions of 1.23m x 1.54m x 0.92m deep and an internal volume of 1.7m 3. Using this enclosure, the effect on the steady state gas concentration of changes in pipe diameter and pipe height were investigated, for natural gas release rates between and ~ 3 m 3 s- 1. The influence of vessel volume and shape on the steady state gas concentration was investigated by conducting additional experiments in two other enclosures: a duct-like vessel 0.61m x 1.83m x 0.61m deep, of internal volume 0.68m 3, and a cubical enclosure of 3ide length 2.7m and an internal volume of approximately 20m 3. In practice, it is likely that cowls would be fitted bo the vent pipes to prevent ingress of water into an enclosure. Consequently, the effect of these on steady state gas concentrations was investigated. During an experiment the build up of gas concentration within the test enclosure was continuously monitored by means of a multi-stream rapid gas chromatograph. Sampling points were located specifically under the base of the inlet and outlet pipes, others being located at various positions throughout the enclosure volume. RESULTS _AND _DISCUSSION A typical curve of the build-up of gas concentration with time is shown in Figure 2. This illustrates clearly that a buoyancy driven flow has been established i.e. fresh air is sampled at the base of the inlet pipe and gas-air mixture is leaving the enclosure via the outlet pipe. The curve also indicates that the mixture accumulation within the enclosure is of essentially uniform concentration. Data from all the tests performed showed that the time taken to reach steady state conditions was inversely proportional the magnitudes of both the gas leakage rate and the available ventilator area (hence, the air ventilation rate). The build-up of gas concentration illustrated in Figure 2 is representative of experiments in which the flow is established 10

5 correctly. However, in some cases flow reversal occured; the effect of this on the gas concentration is described in a later section. M^SSi of Pipe Diameter All the tests to investigate the effect of pipe diameter on the gas build-up were carried out in the 1.7m 3 vessel, using an horizontal leak orientation in the middle of the vessel and a gas delivery pressure of 0-35bar. In any given test, the pipes were of the same diameter with heights of 4m and 3m for the outlet and inlet pipes respectively. Figure 3 summarises the variation of the measured steady state gas concentrations with gas leakage rate for each of the pipe diameters used. The trend observed is as expected from the predictions of the 'simple model' i.e. the larger the pipe diamter, the lower the final gas concentration for a given gas leakage rate. Effect of Pipe Height To determine the effect of both the relative and actual heights of the vent pipes, tests were carried out with inlet and outlet pipe heights of 1m and 2m and also 1m and 4m respectively, in addition to those in which the corresponding pipe heights were 3m and 4m. With a difference in vent pipe height of 1m (i.e. a height combination of 3m and 4m or 1m and 2m), the results obtained for the three diameters of pipe used were similar: for an outlet pipe height of 4m, the measured steady state gas concentrations were about 20$ - 30$ higher than predicted, whereas for an outlet pipe height of 2m the measured gas concentrations were approximately 20$ - 30$ lower than predicted by the simple model. Data for the 0.154m diameter pipes is presented in Figure 4- Tests carried out with a larger difference in the pipe heights (e.g. an inlet pipe height of 1m and an outlet pipe height of 4m) resulted in lower steady state gas concentrations than predicted. It is considered that the differences between measured and calculated gas concentrations may be due to frictional losses caused by flow along the vent pipes. This factor can be included in the mathematical model. However, in practice its effect on the calculation of the required ventilator area is not significant for the small leakage rates typical of real situations (i.e. up to about lo -4 ra 3 s -1 ). Effect of Vessel Geometry The influence of vessel geometry is included in the simple model only indirectly, through defining the depth of the layer of gas-air mixture as being equal to 0.4d. The effect on steady state gas concentrations of variations in vessel volume and shape is illustrated in Figure 5- It is clear that, for a given gas release rate, the measured steady state gas concentration decreases as the vessel volume is increased. This agrees with the model predictions but there is 11

6 some evidence to suggest that the shape of the enclosure also has an influence. Thus, for the 20m 3 cubical vessel the measured and predicted gas concentrations were almost identical whereas concentrations measured in the smaller, non-cubical vessels were higher than expected and, in fact, increased as the ratio of Lmax: Lmin for the vessel increased. However, this apparent effect of vessel shape could also be a consequence of the initial momentum of the gas jet release, which would become less significant as the volume of the enclosure increased. Effect of Leak Source Characteristics In deriving the simple model the gas leak source characteristics were not specified i.e. no attempt was made to include such parameters as gas delivery pressure, leak orientation or the position of the gas release within an enclosure. Consequently, the influence of these factors on steady state gas concentrations was investigated experimentally. However, the data indicates that, in practice, the influence of leak source characteristics is relatively minor. Thus, measured gas concentrations were independent of the gas supply pressure over the range investigated i.e to 0-35 bar (Figure 6) and, whilst the leak orientation did affect the steady state gas concentration (Figure 7). the variation was limited to a maximum of about ±15$ of the mean, decreasing to about ±10$ at higher gas leakage rates. In the 1.7m 3 vessel used for the experiments, with a vertical leak orientation (upwards or downwards) the jet momentum would tend to induce an essentially upwards flow in the enclosure and, hence, assist the outward flow of gases. However, a horizontal gas release could not and this, combined with the fact that the horizontal gas releases were always directed away from the outlet pipe, might explain the higher gas concentrations obtained with this leak orientation. In a larger volume enclosure, it is anticipated that the influence of leak orientation on gas concentration would be less. Although the effect was relatively small, the influence on steady state gas concentrations of the position of the gas release within an enclosure was in general as predicted by the simple model. That is, the higher the level of the leak (and hence the shallower the layer of gas-air mixture formed) the greater the magnitude of the steady state gas concentration achieved. m Jl Reversal During some experiments, usually with high gas release rates (greater than m 3 s _1 ), flow reversal(s) occurred - particularly under gusty wind conditions. In some cases, a steady flow pattern was never established whilst in others a reverse flow was established and maintained. i.e. fresh air was drawn into the outlet pipe and mixture flowed out of the test vessel via the (shorter) inlet pipe. The effect of this flow reversal is to reduce the effective height of the outlet pipe and hence lead to an increase in the steady state gas concentration due to the reduced flow of ventilation air. 12

7 Table 1 compares some measured gas concentrations obtained under 'correct' and reversed flow conditions. TABLE 1. EFFECT OF FLOW REVERSAL ON GAS CONCENTRATION Gas Leakage Rate Measured Gas Concentration {% gas-in-air) (m 3 s -1 ) Correct Flow Reversed Flow "" * 1.2. TO" To overcome the problem, the use of different types of cowls was investigated. These included the standard 'chinamans hat' type, a modified version of this design which has an inner inverted cone (claimed by the manufacturers to prevent back pressure in a vent pipe) and a rotating, French designed 'Aspiromatic' cowl (which was installed only on the outlet pipe). All of the cowls were successful in preventing flow reversal and, in addition, did not affect the ventilation air flow to any significant extent. The 'Aspiromatic' cowl (installed only on the oulet pipe) was particularly effective: because this rotated freely in the wind, it tended to induce gas-air mixture up the outlet pipe, hence assisting to establish flow in the correct direction. CONCLUSIONS 1. Given the assumptions inherent in the derivation of the simple mathematical model of buoyancy driven ventilation, the model predictions and experimental data were in reasonable agreement. 2. Ventilator areas calculated using this model will provide levels of ventilation adequate to disperse 3afely gas leakages of the magnitude anticipated in governor pits and similar below ground installations. (i.e. leakages up to about 10-4 m 3 s- 1 ). 3- The use of cowls on vent pipes will normally prevent flow reversal. The most effective arrangement would appear to be a 'Chinamans hat' type of cowl on the inlet pipe and an 'Aspiromatic' cowl on the outlet pipe. 13

University Physics (Prof. David Flory) Chapt_15 Thursday, November 15, 2007 Page 1 Name: Date: 1. The principle of fluid pressure that is used in hydraulic brakes or lifts is that: A) pressure is the same

Manual M000 rev0 03/00 THE WAY THE VENTURI AND ORIFICES WORK CHAPTER All industrial combustion systems are made up of 3 main parts: ) The mixer which mixes fuel gas with combustion air in the correct ratio

LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/ This experiment will introduce you to the kinetic properties of low-pressure gases. You will make observations on the rates with which selected

Increase in Evaporation Caused by Running Spa Jets swhim.com Nomenclature A pipe cross-section area, m D water inlet diameter of the venturi tube nozzle, mm diameter of small end of the throat of the venturi

Exercise 3-2 Orifice Plates EXERCISE OBJECTIVE In this exercise, you will study how differential pressure flowmeters operate. You will describe the relationship between the flow rate and the pressure drop

CHAPTER 7 : SMOKE METERS AND THEIR INSTALLATIONS 1 Scope : This Chapter covers the requirements of smoke meters and their installation on engines for full load and free acceleration tests, mentioned in

PNEUMATIC CONVEYOR OVERVIEW conveying is the transport of dry granular or powdered material though a pipeline using air (or an inert gas) as the carrier. It is used extensively in industry where powders

POP Safety Valve POP Safety Valve INTRODUCTION The effects of exceeding safe pressure levels in an unprotected pressure vessel or system, can have catastrophic effects on both plant and personnel. Safety

Pressure Variation with Depth Pressure in a static fluid does not change in the horizontal direction as the horizontal forces balance each other out. However, pressure in a static fluid does change with

Chapter 10.4 Learning Goals Define buoyancy. Explain the relationship between density and buoyancy. Discuss applications of Archimedes principle. 10.4 Buoyancy is a force Buoyancy is a measure of the upward

IIT JEE Achiever 2014 Ist Year Physics-2: Worksheet-1 Date: 2014-06-26 Hydrostatics 1. A liquid can easily change its shape but a solid cannot because (A) the density of a liquid is smaller than that of

INTRODUCTION Multiple port vessels. Hydraulic design of side-ported RO vessel arrangements (How to get the best out of the side-ported vessels) This document has the final target to show the best multiport

Applications of Bernoulli s principle Principle states that areas with faster moving fluids will experience less pressure Artery o When blood flows through narrower regions of arteries, the speed increases

Gas/Liquids Separators Part 3 Quantifying Separation Performance Mark Bothamley, John M. Campbell/PetroSkills In this third article of a three-part series, the results of selected gas/liquid separation

CSC-2000 SERIES Reset Volume Controllers MADE IN U.S.A. DESCRIPTION The CSC-2000 series are designed for use on VAV terminal units in HVAC systems. These are submaster air velocity controllers whose velocity

EDEXCEL NATIONALS UNIT 6 MECHANICAL PRINCIPLES and APPLICATIONS ASSIGNMENT No. 4 NAME: I agree to the assessment as contained in this assignment. I confirm that the work submitted is my own work. Signature

ASME BPVC.VIII.1-2015 UG-35.2 UG-36 (4) It is recognized that it is impractical to write requirements to cover the multiplicity of devices used for quick access, or to prevent negligent operation or the

Theoretical Analysis of Opposing Air Bearing Concept This concept utilizes air bearings to constrain five degrees of freedom of the optic as shown in the figure below. Three pairs of inherently compensated

Design and Safety Document for the Vacuum Windows of the NPDGamma Liquid Hydrogen Target at SNS Prepared: Checked: Approved: H. Nann W. Fox M. Snow The NPDGamma experiment is going to run at BL13 at SNS

EXPLOSION VENTING - THE PREDICTED EFFECTS OF INERTIA By Steve Cooper - Stuvex Safety Systems Limited Explosion venting is an established and well used method of primary explosion protection within industry.

QUESTION 1 A ball is thrown horizontally from a cliff with a speed of 10 ms -1 shown in the diagram at right. Neglecting the effect of air resistance and taking gravitational acceleration to be g +9.8ms

H16 Losses in Piping Systems The equipment described in this manual is manufactured and distributed by TECQUIPMENT LIMITED Suppliers of technological laboratory equipment designed for teaching. BONSALL

Lecture 6 HYDRAULIC CIRCUIT DESIGN AND ANALYSIS Example 1.14 Design a car crushing system. The crushing force required is such that a 15 cm diameter cylinder is required at a working pressure of 16.5 kg/cm.

Lesson: Ocean Circulation By Keith Meldahl Corresponding to Chapter 9: Ocean Circulation As this figure shows, there is a connection between the prevailing easterly and westerly winds (discussed in Chapter

HAZARDS OF HEATING IMMISCIBLE LIQUIDS R.B. Keey Department of Chemical Engineering, University of Canterbury, Christchurch, New Zealand. A recent incident at an oil re-refining plant in New Zealand highlighted

Thermal design of APFC Panels The life of the power capacitors and other equipments in APFC panels depend very much on the operating temperature. In panels with detuned harmonic filter reactors and thyristor

MAKING MODERN LIVING POSSIBLE Technical brochure Regulating valves REG-SA and REG-SB REG-SA and REG-SB are angleway and straightway hand regulating valves, which act as normal stop valves in closed position.

2.4. Applications of Boundary Layer Meteorology 2.4.1. Temporal Evolution & Prediction of the PBL Earlier, we saw the following figure showing the diurnal evolution of PBL. With a typical diurnal cycle,

larner-johnsonvalve presentation History! The Larner-Johnson valve was originally designed and developed in 1920 s for use on hydraelectric power plants in the U.S.A., where it pioneered the application

ractical Guide The following article was published in ASHRAE Journal, March 2003. Copyright 2003 American Society of Heating, Refrigerating and Air- Conditioning Engineers, Inc. It is presented for educational

APPENDIX A: SENSITIVITY TESTS INPUTS AND SETUP TABLES This appendix contains two tables for each case and model; the first outlines the assumptions and inputs that were kept constant for each case, and

1.26. A certain object weighs 300 N at the earth's surface. Determine the mass of the object (in kilograms) and its weight (in newtons) when located on a planet with an acceleration of gravity equal to