Vented Gas Explosions: by Rafiziana Md. Kasmani
Vented Gas Explosions: by Rafiziana Md. Kasmani
By
Rafiziana Md. Kasmani
Submittedin accordance
with the requirementsfor the degreeof Doctorof
Philosophy
University of Leeds
Schoolof Process,EnviromnentalandMaterialsEngineering
Underthe supervisionof
and
Dr.H.N.Phylaktou,Phd,B.Eng
September2008
ABSTRACT
mb
= AS. p,, = CdeA,(2p,,Pd)o"
(1)
The equationgiven in (2) is slightly different from (1) as is about 6.5 times the mass
flow of the first methodasit takesthe effect of (E-1) whereE is the expansionratio.
A critical review werecarriedout for the applicability, validity and limitation on the
venting correlationsadoptedin NFPA 68 and EuropeanStandardwith 470 literature
experimentaldata,coveringa wide rangeof valuesfor vesselvolumeand geometries,
bursting vent pressure,P,, UD ratio, maximum reducedpressure,P"d and ignition
location. The fuels involved are methane,propane,hydrogen, town gas, ethylene,
It was observedthat end ignition leads to a higher explosion severity than central
ignition in most cases,implying that central ignition is not a worst-casescenario in
gasventedexplosionsas reportedpreviously
iii
CONTENTS
Abstract
Contents
List of figures vii
List of tables xii
Nomenclature xiii
Acknowledgements, xv
Publications by author xvi
1.0 Introduction I
.......................................................................................................
1.1 Venting gasexplosions 2
......................................................................................
1.1.1 Generaloverview 2
.......................................................................................
1.2 Overviewof the project studies 8
..........................................................................
1.2.1 Simply venting gasexplosions 9
...................................................................
1.2.1.1 Self-accelerationof sphericalflames 12
....................................................
1.2.1.2 Implication of the objectivesto the presentstudy 14
..................................
1.2.2 Ventedgasexplosionwith relief pipe attached 14
. ........................................
1.2.2.1 The influenceof ductpipe diameter 17
......................................................
1.2.2.2 Implication of the objectivesto the presentstudy 17
..................................
2.5 Instrumentationanddatacollection 34
..................................................................
2.5.1 Instrumentationtechniques 34
.......................................................................
2.5.1.1 Explosion pressure history 35
....................................................................
2.5.1.2 Measurement of flame positions 35
...........................................................
2.5.1.3 Thermocouples 36
.....................................................................................
2.5.1.4 Flame speedcalculation 37
.......................................................................
2.5.2 Data acquisition 39
.......................................................................................
2.6 Identification hazardsand safety measures 39
.......................................................
2.6.1 Vessel failure 40
..........................................................................................
2.6.2 Transmission of the explosion to auxiliary equipment 40
..............................
2.6.3 Releaseof combustible gas into the Test room 41
.........................................
2.7 Operating procedures 41
.......................................................................................
2.7.1 Fuel entrainment 41
......................................................................................
2.7.2 Ignition procedure 42
....................................................................................
2.7.3 Vessel purging 42
.........................................................................................
2.8 Ancillary equipment 42
........................................................................................
2.8.1 Isolation valves and pipe works 43
................................................................
2.8.2 Mixture pressure monitoring system 43
........................................................
2.8.3 Vacuum gate valve 47
...................................................................................
2.8.4 Evacuation system 47
....................................................................................
2.8.4.1 Vacuum pump A 47
..................................................................................
2.8.4.2 Vacuum pump B 47
..................................................................................
2.8.5 Ignition system 48
........................................................... I...I.,.......................
2.8.6 Nitrogen purging system 48
..........................................................................
3.0 Introduction 50
.....................................................................................................
3.1 Parameterinvolved in empirical equations 51
.......................................................
3.2 Venting mechanism 53
. ........................................................................................
3.3 Venting theory 58
.................................................................................................
3.4 Comparisonof Bartknecht and Swift correlationsas been adoptedin NFPA 68
usingderivedmethodson the samebasis.....................................................................
71
3.5 Reviewof the publishedexperimentaldata 74
. .....................................................
V
experimental data 76
. ...................................................................................................
3.5.2 The influence of turbulent enhancementfactor, 0 80
.....................................
3.5.3 The influence of the vessel's geometry 85
.....................................................
3.5.4 Comparison between K and A, /A, term for vessel's geometries ............... 90
3.5.5 Evidence for an additional influence of vessel volume on Pred.................. 93
3.6 Concludingremarks 96
.........................................................................................
4.0 Introduction 98
.....................................................................................................
4.1 Generalfeaturesof experimentaltests 99
..............................................................
4.2 Generalexplosiondevelopment 102
.....................................................................
4.2.1 Resultsand discussionson Test vesselI 104
................................................
4.2.1.1 Maximum pressure,P.. as a function of equivalenceratio ................104
4.2.1.2 Flamespeeds 108
.....................................................................................
4.2.2 Resultsand discussionon Test vessel2 116
..................................................
4.2.2.1 Maximum pressure,Prn. 116
....................................................................
4.2.2.2 Flamespeeds 119
......................................................................................
4.3 Deflagrationto detonationin test vessel 122
.........................................................
4.4 Influenceof vent coefficient,K, volume andburstvent pressure,P, on Prnfoc.. 126
4.5 Comparisonbetweentheoryand experimentaldata 134
........................................
4.6 Concludingremarks 138
.......................................................................................
REFERENCES 221
....................................................................................
APPENDIX 230
........................................................................................
vii
LIST OF FIGURES
aluminiurn foil. From bottom left: 100g A4 paper and magazine paper 34
Figure 2.6 Typical thermocouple traces indicating the change in voltage
the rig. 37
Figure2.8 Flow diagramillustratingsafetyandoperationalprocedure 44
Figure2.9 Exampleof operatingprocedurefor carryingexplosiontests 46
Figure3.1 Pressurebehaviourversustime for a closedand ventedexplosions 53
Figure3.2 The pressure-timevariationin a ventedexplosion.Graph
o cubic 87
Figure 3.14The influenceof rectangularvesselson Predfor propane/air
0 rectangular 89
Figure 3.15The influenceof variousgeometrieson Predfor town gas/air
0 rectangular,o cubic, o sphere 89
Figure 3.16The influenceof the vessel'sShapeOnPred-
A cylinder c,sphere 90
Figure 3.17Methane/airfor Predasa function of K 91
Figure 3.18Pmdof methane/airat variousgeometriesasa function of AdAv 92
Figure 3.19Ko asa function of vesselsizefor propane/air(4%) measured
at an initial pressure of I atm (Chippett, 1984) 94
Figure 3.20 Effect Of Predon vessel volume at constant K 95
Figure4.1 Typicalpressure-timehistoriesfor centreandendignition
in TestvesselI 103
Figure4.2 Typicalpressure-timehistoriesandflame positionin Test vessel2 103
Figure 4.3 Methane/airat different equivalenceratio, 0 105
Figure4.4 Propane/airat different equivalenceratio, 0 106
Figure4.5 Ethylene/airat different equivalenceratio, 0 106
Figure4.6 Maximumrateof pressurerise,dP/dt(max) at variousequivalence
ratios for propane/air(above)andethylene/air(below) 107
Figure 4.7 Hydrogen/airat differentequivalenceratio, 4) 107
Figure 4,8 Flamespeedsof methanelairat endandcentralignition as
venting 110
Figure 4.10 Flamespeedsof propane/airat end andcentralignition as
methane/airmixture 128
Figure4.25 Hydrogen/airat stoichiometricconcentrationfor different K 129
Figure4.26 Influenceof P,,on P for Testvessel2 133
......
historiesfor simply ventedandduct vented
Figure5.1 Pressure-time
explosionfor methane/airat 0=1.06 for endignition (initially open
venting).dP/dttraceswere given for both cases 146
Figure 5.2 Pressurerecords at selected positions along the test vessel for
ratio 172
Figure5.18 Propane/airfor end andcentralignition asa functionof equivalence
ratio 174
Figure5.19Ethylene/airfor endand centreignition asa function of equivalence
ratio 175
Figure5.20 Hydrogen/airfor endand centreignition as a function of equivalence
ratio 176
Figure5.21 Flamespeedin the duct pipe for propane/airasa function of
equivalenceratio at endignition 178
xi
LIST OF TABLES
from BradleyandMitheson(BradleyandMitcheson,1978b)except
andfuels 83
data 199
NOMENCLATURE
Ka Karlovitz number
Rred dimensionless
maximumpressure(Ped/Pj)
ill, dimensionless
vent burstingpressure(P,/Pi)
Greeksymbols
turbulent enhancementfactor
PU unbumtgasdensity
Pb burntgasdensity
X/P deflagration-outflow-interaction
number
& expansibilityfactor
(D equivalenceratio
AP pressuredifferential
xvii
ACKNOWLEDGEMENTS
Thank is also expressedto R.A. Boreham, for his help and technical expertise for
completingthis project.I would like to thankmy other supervisor,Dr.Roth Phylaktou,Ms
S.K Willacy for her help and motivation and fellow postgraduatesfor guidance
throughoutthis project,
Most especiallya special thank to my parents and family for a constantsupport and
encouragementin completion of this project and thesis. Not forgetful to Nadiah, Liz,
Toby, Emma,Danial,Brian and Bev for welcomedistractionfrom study and making my
life in Leedsmoremeaningful.
Kasmani, R.M, Andrews, G.E, Phylaktou,H.N, Willacy, S.K. (2007), Vented gas
explosionin a cylindricalvesselwith a vent duct,Third EuropeanCombustionMeeting
(ECM 2007),Crete,Chania.
1.0 Introduction
Studieson vented gas and dust explosionshave shown the evolution of venting area
with pressuredepend on the nature and state of the initial explosive mixture i. e.
composition,initial pressureandtemperature,pre-ignition turbulenceand on the vessel
characteristicsi. e. dimensionand shape,position of the ignition source,location, size,
strengthand shapeof vent, presenceof obstacleswithin the vessel.An understandingof
the mechanismsby which pressureis generatedin ventedexplosionsis importantin the
design of explosion reliefs and to the investigation of incidents. Such knowledge
providesthe basicof the developmentof prediction methodsfor use in risk and hazard
studieson an actualor proposedindustrialplant.
quiescent gas-air mixtures. Yet, the different basic physical aspects of vented
deflagrationshave still not beenstudiedsufficiently. For decades,venting experiments
witb various initial conditions and modelling model using lumped parameteras a
venting parameter were done to understand the dynamics and physical process of
venting. From the intensive works done by previous researchers(Bartknecht, 1993,
Bradley andMitcheson, 1978a,Bradley and Mitcheson, 1978b,Burgoyne and Wilson,
1960, Cousins and Cottons, 1951, Cubbage and Marshall, 1974, Cubbage and
Simmonds, 1955, Cubbage and Simmonds, 1957, Donat, 1977, Zalosh, 1980),
correlations were offered to be used in order to design the required vent area;
unfortunatelythey are strictly applied within the limits of the experimentalrange
investigated.Hence,the key knowledgeaboutthe physicaland dynamic mechanismof
venting with or without the relief pipe needsto be understoodin order to reduce the
potential of explosion severity.
1.1.1 Generaloverview
per unit areaof the relief, location of ignition as well as shelf (obstacles).
arrangement
Low breakingburst pressureslead to large flame elongation,and hencelarger name
areaand result to higher expectedvesselpressureswhereasat high breaking pressure
the effectof flamedistortionis reduced.
correlated the experimental results as shown in Eq. 1.1 where a and b given in Table 1.1.
-b (1.1)
A, = aP,,
d
NFPA 68 adopted Swifts equation (Swift, 1983) to apply for design and scale-up of
it in
venting system purposes and seem that the results obtained were good agreements
data obtained by Donat (Donat, 1977) and Harris and Briscoe (Harris
with experimental
Briscoe, 1967). Again, there is a limitation in using this approach as it can be
and
i. Pred is not more than 200 mbar and only suits to low-
applied on certain conditions e.
enclosures. Meanwhile, Swift's correlation is given as A, as a function of
strength
cross-sectionalarea of the A,
enclosure, multiplied to reactivity terms which include the
burning velocity, S,, and turbulent enhancementfactor, P. The use of the vented vessel
in
surface area vent correlations was first introduced by Runes (Runes, 1972) as the
maximum possibleflame to
area and was used calculatethe maximum flame mass
buming rate.
(Nagy andVerakis, 1983) and hencethe two parametersto quantify the reactivity are
directly related.
of volume dependence
on KG was also madeby Chippet (Chippett,1984)on his work,
on propane/airmixtures.
1978a, Bradley and Mitcheson, 1978b) have compared the results obtained by lumped-
parameter model of vented explosions to the experimental data but yet, only a
quantitative agreement was observed. The work was regarded more complete since it
also offers a sound theoretical basic in predicting vent area, however the criteria they
developed have been shown inadequate for large volumes with obstacle-generated
turbulence (Swift, 1983). Further critical examinations of various existing equations and
especially of their extrapolation within and beyond their recommended validity range
have been presentedby Molkov el al (Molkov, Dobashi, Suzuki and Hirano, 2000) who
also developed a new so-called 'universal correlation' for vent area calculation. Razus
et al (Razus and Krause, 2001) reviewed the empirical correlations developed for vent
sizing based on the reduced pressure parameter. Following the authors, equation by
Bradley and Mitcheson (Bradley and Mitcheson, 1978a, Bradley and Mitcheson, 1978b)
and Molkov (Molkov, Korolchenko and Alexandrov, 1997) gave relatively close
predictions to experimental data for most common hydrocarbon-air mixture.
Modelling on venting explosion has been carried and still, disparity betweenmodels
proposedis huge. Complex computer codes requires considerableinput data and
accuratethermophysicalpropertiesas well asaccessto super-computerfacilities are of
little use in the systemdesign (Zalosh, 1995).Their utility in exploring less complex
flame laws however,is not in question.Designtechniquesshouldbe quick and easyto
use,accurate,conservativeyet realistic andprovideresultsthat accordwith experience.
manyexperimentaldonein stoichiometricconcentration
The publishedgas venting data from 1920sto date were complied and collected with
different volumes, gas reactivity, ignition position, shapes and different initial
conditions.From the reporteddata, the comparisonbetweenBartknecht's and Swift's
equation were made as those equationswere offered in NFPA 2002 and European
Standard(2007). Further, those data were analysedto be comparedto Bradley and
9
The differenceof these two methods was the use of S,, and S,, (E-1) terms in the
equations.The simple approacheswere evaluatedwith additional of compressibility
factor and turbulencefactor to simulatethe real explosions.Thesetwo methodswere
comparedwith the publishedexperimentaldatato justify the applicability, validity and
limitation with respectto the practicalapplication.The detail discussionon the effect of
0.5820
1/K, = 0.197/Pd (1.5)
0.5820
1/Kv = 0.164/Pred (1.6)
0.5820
VK, = 0.291APýd
0.5820
1/Kv = 0.220/Pd
It will be shown that these predictions are grossly in error and even with a vent pipe the
overpressures come no where near these values. When the experimental data of
Bartknecht is examined (Kasmani, Willacy, Phylaktou and Andrews, 2006, NFPA68,
2002), it is clear that the exponent Of Predis fitted to the data with I/K > 0.2 or K<5
and the be
correlationsshouldnever usedfor higher K. This is not recognizedin NFPA
68 or in theEuropeanguidanceon venting.The limitation on the correlationis statedas
Pred<2 bar andthe experimentaldatashowsthat this is effectivelya limit of K of 10. It
is found that in cubic vessels,this condition is not valid as the flame have touchedthe
wall well before2 bar overpressurewas reached.However,none of the experimental
data for Pred> 0.8 bar fits the correlation and all is lower than the correlation, This
effectively givesthe validity of the correlationat K<5 andPred< 0.8 bar andthis is the
incompressiblevent flow regime. When the nearestBartknecht's data point to the
In this work, the effect of fuel/air mixtures concentrationat different ignition position
Further, the use of vent cover i.e. lightweight relief panel, rupture membrane or
explosiondoor has proven to have effect on explosion pressuregeneration inside the
i. e. back pressure due to inertia of the vent cover (Molkov,
vessels/enclosures
12
Grigorash, Eber, Tamanini and Dobashi, 2004). This commercial vent cover can be
explosion. The example of commercial vent cover available for venting is shown in
Fig. 1.1.
The effect of inertia has been studied extensively in order to quantify the effect of
inertia to the Predduring venting (Bartknecht, 1985, Cubbage and Simmonds, 1955,
Molkov, Gfigorash, Eber and Makarov, 2004). However, in this present work, the effect
of vent cover inertia can be neglected as the material used is easily broken and not
impending the flame propagation during venting explosion.
Although the laminar burning velocity and hence the flame speed in a spherical flame at
small then having a period of near constant flame speed and at so-called ' critical flame
radius', the flame starts to rapidly accelerate through a mechanism called ' self-
acceleration' or flame cellularity. The net result is the KG is a function of the vessel size.
Harris and Wickens (Harris and Wickens, 1989) had done series of small and large scale
13
that they propagateat the laminar flow rate before startingto self-acceleratedue to the
processof flame instability throughthe formation of cellular flames that have a greater
total flame surfacearea than a smoothlaminar flame. The detail of self-acceleration
mechanism hasbeen discussedby Bradley et al (Bradley,Cresswell
and Puttock, 2001).
In cellular flame, the magnitudeof the increasein flame speeddependson the type of
gasand the equivalenceratio of the mixtures.It is generallyknown that methanehas a
smaller accelerationthan propanewhich in turn, is smaller comparedto hydrogen.
Hydrogen/air flame becomescellular very easily as been shown by Andrews and
Bradley (Andrewsand Bradley, 1973)in 50 mm diameterof closedcylindrical vessel.
They observedthat the onsetof flame cellularity occursat the mixtures leanerthan 25
% concentration.
explosionprotectiondesign.
14
The studyof the venting gasexplosionof different fuel/air mixtures and concentrations
were further examinedin terms of the effect of duct pipe attachedto the main vesselon
P ,,. Basedon the experimentalanalysisof ventingexplosionwith and without a pipe,
it is known that the severityof the explosionis likely to be 2 or 12 fold increasewith the
presenceof a duct to
with respect simply vented vessels (Bartknecht, 1981, Ferrara,
Benedetto,Salzano,Russo,2006, Ponizy and Leyer, 1999a,1999b,Siwek, 1996). The
resultsof DeGood et al (DeGood and Chartrathi, 1991) and McCann et al (McCann,
Thomas and Edwards, 1985) for various duct lengths supportedthese results. The
surprisingfeatureof theseresultsis that the length of the duct is not a major parameter
andonly short and long ductswere separatedby Bartknecht(Bartknecht,1993).This is
becausethe increasein the overpressureis dominatedby duct entry and exit pressure
losses;a duct would have to be about 100pipe diameterslong before the flow friction
losseswereequalto the combinedinlet and exit pressurelosses.The presentwork used
15
accelerationof the flame prior to the vent duct and also increasethe flow velocity,
turbulence and flame speed in the vent duct. Tests with ducted explosion vents
generallydisplay Helmholtz oscillations, that is the pocket of burned gas within the
vesselundergoes bulk motion towards andaway from the vent opening,due to the mass
of the duct contentsandthe compressibilityof the gasin the primary enclosurethat acts
as a spring(McCann,Thomasand Edwards,1985).The occurrenceof a sharppressure
peak in the vent duct is also mentionedby Kordylewski and co-worker (Kordylewski
The increasein Predwith the additionof a vent duct is so largethat vent ductscannotbe
used without increasingthe vent area and duct size to achieve a reduction in the
overpressure.However,there is insufficient designdata for gasesto enablethis to be
done effectively and the physics of the processfor gas explosionventing is not well
understood.This contrastswith the situation for dust explosions,where a substantial
experimentaldatabaseexists(Lunn, CrowhurstandHey, 1988).The aim of the present
work was to explorethe physicsof the venting processwith a vent duct attachedwith
variousduct diametersto investigatethe effect of the vent staticburst pressureand the
duct diametersize for central and end ignition for an explosionvesselat the limit of
applicabilityof compactvesselventingcorrelationswith a lengthto diameterratio, LM
of 2. It is shownthat generally,reverseflow in ductsof largerdiameterswill result in
large scalemixing in the vesselwhich in turns generateshigherpressureincreasein the
The influence of vent ducts on gas explosions was investigatedwith the aim of
determiningwhetherthe use of a vent duct of large areathan that of the vent, would
reducethe very large increases
in overpressurethat occur whenvent ducts are used.In
the presenceof a vent duct, an increaseof venting areaand duct diameterhas been
found to not alwaysresult in a decreasein the peak over pressure(Ponizy and Leyer,
1999a).
The main reasonfor the increasein the overpressurewhenlong vent ductsare attached
to vents is due to the phasein the explosionwhen the flame is in the vent pipe with
unburnt gas mixture aheadof it. The expansionof the burnt gases in the vent pipe
greatly acceleratesthe unbumt gas flow and this increasesthe vent pipe friction, inlet
and exit pressurelosses (Lunn, Crowhurstand Hey, 1988).These effects are a function
of the dynamic in
pressure the In
vent pipe. principle the dynamic in
pressures the vent
be by
pipe can reduced simply using a larger vent diameterthan that for the vent, rather
than increasingboth the vent and vent pipe sizes.For example,if the vent pipe was
twice the diameterof the vent then the vent pipe dynamicpressurewould be reducedby
a factor of 16, if the vent massflow rate remainedconstant.Some evidence that a larger
vent pipe diameter would reduce the overpressurewith no changein the vent size was
provided by Nagy (Nagy and Verakis), is
which quoted in NFPA 68 (2002). From the
experimentsperformedon dust explosionswith vent pipes (Lunn, Crowhurst and Hey,
1988),Hey (Hey, 1991) has suggestedthat the techniqueof using a larger vent duct
diameterthan the vent diameteris effective if the duct area/ventarea is about - 2-2.5
and when Predis less than 0.5barg. It is consideredthat this approachwould be a
simpler methodof designingfor safe vent pipes and the presentwork investigatedfor
gas explosionsa vent pipe that was close to twice the vent diameter,as recommended
by Hey (Hey, 1991).
relief pipe attached to the primary vessel. It is the aim of the present study to quantify
the relationship of the Pnz,,,with ignition position, equivalence ratio and vent bursting
pressureP, as well as the duct pipe diameter and length in order to highlight the most
leading factors of the rapid increase in P.. at the presence of the duct in vented gas
explosion with respect to the simply vented explosion.
CHAPTER 2
The researchhas been carried out in Room BII which has been designatedas the
'Explosion Hazards-HighPressureTest Facility'. The site is shown schematicallyin
Fig. 2.1. Safetyconsiderationshavebeentop priorities to make sure that all tests were
carried out in a safety manner.Two separatedrooms have been designed: 'Control
room' and 'Test room'. Accessbetweenthe two roomswasvia an interlockeddoor that
formed part of the safety systemcontrolling the ignition circuit. The dump-vessel,test
vesselsand instrumentationequipmentswere locatedin the Test room. Data logging
hardwarewas situatedin this areaandwas electronicallylinked to a computernetwork
in the Control room. Data collection was synchronisedwith ignition using specialised
softwareknownasFAMOS.
Test Room
Control room
pressurewith time in venting gas exýlosion dependson the nature stateof the initial
explosive mixture, i. e. composition, initial pressureand temperature,pre-ignition
turbulenceand on the vesselcharacteristicsi. e. dimensionand shape,position of the
ignition source,location, size,strengthand shapeof the vent, presenceof the obstacles
within the vessel.Further,the addition of the relief pipe to the vesselcan significantly
increasethe violenceof the explosionand the pressurein the vesselto the factor of 10
or morein which the explosionis venteddirectly to the atmosphere(Bartknecth,1981,
Kordylewski and Wach, 1988,Molkov, Baratov and Korolchenko, 1993, Ponizy and
Leyer, 1999a,Ponizy and Leyer, 1999b).The test configurationwas employedfor the
following reasons;
VT = VI +Vd
Vb=
7.5V, (2.2)
V,
P2 = PI
V2
(2.3)
where
PI =I atm.
3)
Vi = systemvolumeafter combustion (Vb + Vd) (M
Vd
P2 ='X 7.5V,+
(2.4)
V, +Vd
P=6.5V, (2.5)
SYSV, + Vj
23
This calculation was performed for each test vessel volume and values of Py, are listed
in Table 2.1.
The Test vessel I was made up from the existing 05 in diameter and Im length
cylindrical test vesselsthat have been designedto BS1560and rated at 28 bars as
illustratedin Fig. 2.2. The vesselwas drilled with tappedbosseswelded on at positions
along their length,to allow instrumentationin the form of thermocouplesand pressure
transducersto be fitted. A removableblank plate was usedto close one end of the
vessels.This was centrallydrilled and tappedto allow a sparkplug, hoseconnectedto
the barometerandhoseconnectedto vacuumpump to be fitted. At the other end of the
vessel,a 0.162in diameterwas drilled to simulatethe vent areaof 0.012 m2.A certified
hydraulic pressuretest had been carried out on the test vesselsbefore commissioning
the tests.A scaleddrawingof test vesselis shownin Fig. 2.2 andmajor dimensionsand
pressurerating arelisted in Table2.1.
The Test vessel2 was madeup from 0.162 m diameterand 0.315in length cylindrical
test vesselsthat havebeen designedto BS1560 and rated at 35.5 bars (Fig. 2.3). The
vesselwas drilled with tappedbossesweldedon at positionsalongtheir length, to allow
instrumentationin the form of thermocouplesand pressuretransducersto be fitted. A
removableblank plate was used to close one end of the vessels.This was centrally
drilled and tappedto allow a sparkplug, hoseconnectedto the barometerand vacuum
pump to be fitted. Seriesof test on Test vessel 2 were carried out to investigatethe
effect of vent areaon explosiondevelopment.The single-holeorifice plate was usedto
simulatethe objectiveon varying the vent areawhich placedat the end of the vesseljust
beforethe gatevalve andit gives,
A, = (I - BR)xA
24
where A, is the vent area, BR is the blockage ratio of the orifice plate and A is the area
of the plate.
For ventingwith relief pipes,Im pipe lengthhasbeenconnectedto the other end of the
Test vessel 1 for a constantvent areaof 0.012 m2 which beencalled a Test vessel 3
(Fig. 2.3). For Im pipe, thereweretwo different internal diameterpipe fitted, i.e. 0.162
m and 0.315 m. Vesselsupportwas in the form of moving framesfabricatedin-house
from L-sectionmild steel.Further,the test vesselwastied with the cradleof the craneto
easethe movementto facilitate quick replacementof vent cover andpipe length without
major rig dismantlement.
Designpressure(bar) 28 35.5 28
Flanges
Numberof bolts 20 12 20
25
Pipe section
Assembled test-vessel
(mbar)
26
Li
--L-i
IIi
Test vessel I
N
i1
T1j
I.
/
Test vessel 2
cl
ý-- CO
Ln
(ID C)
,t
CU
LH
it
-`0
162
Z
315
openings were included at various positiolls, with diameter matched that of the tost
30
Blanking plateswere drilled and tappedto allow connectionof pipe work to facilitate
evacuation via an external vacuum pump or the introduction of compressedair.
Pressureratedball valveswere usedto isolateconnectionancillary equipmentfrom the
dump vesselprior to ignition in the test vessel.Before explosiontesting was started,a
certified hydraulicpressuretest was carriedout by pressuringthe water filled vesselto
11.25 baro.
Shell
Length(m) 0.672
,
Shellthickness(m) 0.015
Torispherical(2:1) Dishedends
Outerdiameter(m) 2.500
Assembled structure
Design pressure(bar) 9
Flanged openings
O/DIA
O/DIA 40/3
I class300
1 1_ L class300
FIgLiFe2.4 DLIIIIJ)VCSSCI
SCI)CIIIItICCleSiglIS
33
(mm) mass/area
(kgIM2)
Important parameters studied in the general analysis ol' gas explosions are tile I'lame
operation. For this present project, the relatively high pressure ocnerated and large scale
experimental rig meant that there were S111'ety 10RIICOut the ViSLIa11SatIOII
ColiSider,111011S
technique to monitor flarric travel.
35
2.5.1.1 Explosionpressurehistory
A thermocouplejunction placedin the path of a flame travel will register its arrival as a
changein voltagepotentialacrossthe junction. This is recordedas a distinct changein
the analogueoutput of the thermocouple.Accurate measurementof the time of this
signal changerelative to that measuredfrom an adjacentprobe allowed the average
flame speeds to be calculated at the midpoint between the two probes. This
thermocouple technique has been used to monitor the flame arrival in previous
explosion studies (Gardner, 1998) and the system was validated by comparison to
photographicrecords of explosions in a closed spherical vessel. Figure 2.6 shows
examplesof thermocouplesignalsas displayedon a computerscreenfrom which flame
arrival time measurementwere taken using a movable cursor, This thermocouple
techniquewas a primary methodof flame monitoring usedin this study.
36
The thermocouple tips had an inherent thermal mass; however since all the
thermocoupleswere at the samespecification,it was assumedthat they had the same
responsetime. It shouldbe notedthat the thermocoupleswere not usedto measurethe
temperaturebut merelydetectingany changein temperature.The time delayor response
time, usuallyrefersto the time lag for full responseto the stimulusand therefore,is not
relevantat this point.
0.005
O.OCK
0.003
0.002
0.001
0
45
TIM(s)
2.5.1.3 Thermocouples
compression fittings were used to seal the units in tapped BSP bosses. The
thermocoupleswere locatedthrough the vesselwalls and so would be subject to high
dynamicloadsresulting from high velocity gas flow impact, before and after the flame
arrival. In order to limit the bending and maintain accurate positioning of the
17
thermOCOUPIeS, structures were built into test vesselsas well as along the pipes
SUPPOI-t
that were used.
The support consisted of' aI rniri dianieter connecting wire, suspended above the axial
centreline of' the vessel f'rom the thin inctal strips at the end of' the vessel. The
percentage finear blockage of each strip was less than 0.4 I/r. CirCUlar Washers were
! t-
Figure 2.7 Exposedjunction, mineral insulated type-K thermocouples Inside the rig.
By dividing the distance between two thermocouples with the time difference in flarne
During the initial stage, flame speedswere relatively low but will dramatically increase
towards the downstream of' the vessels and to the connecting pipe. UsIng tile flame
thbý
UNN ýI UbMýl
LEED6
38
Raw data
top(ignition) XO'O
t2,0 X2,0 (x2, o-x 1,0)/(t2, o-t 1.0) (X2,0-X LOY2 +X1,0
Smoothing data
tij = (to,o+tjo)/2 xij = (xo.o+ xj,o)/2 (XI., -x 0,1)/(tl,l - to,j) (x 1,1-x o,1)/2+ xo.I
X3,0)/2
tn, l " (tn-1,0 +tn, O)/2 Xn, I = (Xn-1,0+ (Xn, tn- (X + Xn-
l -X n-1,1)/(tn, l - n,I -X n-1,0/2
to get preciseand accurateflame arrival time as well as the flame speedsfor the tests
done.
The short duration of the highly transient explosions meant that accurate analysis
required high speed data collection from the various instrumentations positioned
throughout the test geometry. Processing the large number of recorded data was aided
by specialisedanalysis software.
2.6.1 Vesselfailure
The designfeaturesof the vesselin this study were listed in Table 2.1 and 2.2. All test
vesselswerepressuretestedandcertified to withstandpressureabove8 baro. This is the
maximum adiabatic overpressurethat can be generatedby any fuel-air explosion
initially at I bar. The possibletransitionto detonationcombustionwas also taken into
consideration.The actual design pressurefor each test vessel was higher than the
overpressure with detonation.
associated
This hazardwasnot specifiedto this facility but was also applicableto any University
areawhere combustiblegaseswere presentand handled.In the event of an accidental
explosionwithin the Test room, it was expectedthat the blast would be relievedto the
outside through the large window along the external wall. This large vent area was
expectedto reducethe pressurebelow level that might damagethe partitioning walls
that isolatedthe Control room.The glasswindows werecoveredwith a clear membrane
that wasspeciallydesignedto limit the glasspigmentation.
In order to create a combustible mixture of known composition, the volume of both fuel
and air introduced into the test vessel had to be known accurately. The theory of partial
pressure was applied to an accuracy of 0.1 mbar (0.01% of composition). The final pre-
ignition mixture pressure was set at ambient pressure (- I atm) for
all tests. For
example, to make up a mixture with 10% fuel/air ratio, the test vessel was evacuated to
All isolating valveswere closedand the mixture inlet line was disconnectedfrom the
test vesselandconnectedto a safetyinterlock.The dataacquisitionsoftwarewas armed
and ignition circuit integrity was indicatedby light on the ignition control panel after
the openingof the gatevalve.The 'fire' button waspressedand this activatedthe spark
plug and by this, it triggeredthe data collection system.Once the samplingtime was
over, the raw transducersdata weredisplayedusing FAMOS software.Output signals
were recordedto check whether a pressurerise had occurredin the system and the
maximumoverpressure is not exceededthe critical adiabaticpressureof the fuel/air.
2.7.3 Vesselpurging
Purging involved evacuationof both test vesseland the dump vessel using vacuum
pump B (refer to later section).Once a sufficiently low vacuumpressurewas attained,
air was allowed to enter the systemfor more than ten minutes.This procedurewas to
ensurethat the systematmospherematchedthat of the outside ambient air. The test
facility wasthenreadyfor the next testpreparation.
The preparationand ignition of fuel/air mixture inside the test vesseland its subsequent
purging involved various equipmentswith associatedpipe works and valves. As
described in the previous section (refer to section 2.7), the procedures involved
including the evacuationof the test vessel,mixture preparationand recirculation and
control of the dump vesselpressure.Isolation of all ancillary equipmentwas required
beforeignition andthe test vesselhadto be isolatedfrom the connectedpipe (if venting
with relief pipe geometrywas carriedout) anddump vesselduring mixture preparation.
43
connectingpipe, the pipe attachedto the vacuumpump and the pipe attachedto the
barometershould be disconnectedprior to ignition. Step-by-steptest procedure is
shownin Fig 2.9.
Is vacuumpressurestable?
(Lessthan0.2 mbar/minleakrate)
Yes
No
abandontest pressuremethod
pressure? I
Yes
Yes
x
Yes
I
Yes
6-1
,. I --d by
ý;
N,, . ......... .
..... ........... .... ... .. Ppll
1--l. --. 1 -,
01, - VW, V: i K V4 50",
vessel.It consistedof a light-weight aluminium body and gate with Viton seals.Shaft-
feed-throughwas employedwith pneumaticactuation.Valve operationwas controlled
by a solenoidvalvein the compressed
air supplyline pipe.
A single conventionalcar sparkplug was usedto ignite the fuel-air mixture in the test
vessels.Ignition was actuatedfrom a standardcombustionengine spark (16 J) via a
high-capacitancedischargecircuit. The circuit diagramis included the various safety
interlocks. The system ensuredthat the following operationshad to be carried out
before spark ignition could be effected.
The gate valve used to isolate the test vesselsto pipe and dump vessel during the
closedandlocked.
0 The mixture inlet line had to be disconnectedfrom the test vessel. This
effectively meantisolationof the fuel supplyprior to ignition.
Indication lights, fitted to the ignition panel at the Control room were monitored to
check whether this action had been carried out. The data acquisition system was
connectedto the ignition circuit for triggeringthe fire button for datacapture.
If the fuel-air mixture was not ignited might be due to the failure of the ignition spark
plug or electrical faulty in ignition circuit, a relatively large volume of explosive
mixture would remain inside the test vessel.Before any repairscould be made,the test
vesselneededto be purgedof the gasmixture. A nitrogenpurging systemis one of the
special case circumstancesthat neededto be followed if this case happened.This
procedurecarriedout by injecting high pressurenitrogenthat would mix with the fuel.
A standardpressuregaugeon the Barocel could be monitored so that the volume of
49
nitrogen introduced was sufficient to inert the fuel-air mixtures. At the same time,
simultaneousoperation of vacuum pump B allowed the mixture to be-safety purged and
3.0 Introduction
ST=I +,8S.
concepthas been applied by Cubbage and Simmonds (Cubbage and Simmonds, 1955,
Cubbageand Simmonds,1957),Cubbageand Marshal(Cubbageand Marshall, 1974)
and Rasbash
(Rasbash,
Drysdaleand Kemp, 1976).Bartknecht(Bartknecht,1993) used
K terms to expressthe vent flow pressureloss and vent bursting pressurein his
KO to the
correlationand value express gas/air mixturereactivity. KG is defined as
where dP/dt is the rate of pressurerise (bar/s) and V is the vesselvolume (M3). Kr, is
given as 50,100,140 and 550 bar-m/s for methane, propane, town's gas and hydrogen
in 5 litre sphericalvesseldoneby Bartknecht(Bartknecht,1993).Thesevaluescan only
be used with the Bartknecht's correlation as KG value is not to be constant but to
increasewith vesselvolumeasbeeninvestigatedby Chippet(Chippett, 1984).Thus, KG
the flame speeds.This assumption is not valid as the higher compression pressure and
temperature of the unburnt gases changes the burning velocity from that of the initial
2/3y I/Y]2/3
dP 3Su ýpm 1/y_pl/r)1/3 Po 8)
Tt (P
-R *p (2-11Y - 6) '00p 1_ý p(3-2ly-, (3.3)
01
From Eq. 3.3, onecaneasily derive the dependentof Ko (from Eq.3.2) to S,,.Kumar at
al (Kumar, Bowlesand Mintz, 1992) also derived an equationfor KG as a function of
S.. This equationalso assumesthat S,,is the constantand doesnot vary with radius as
well aspressureandtemperature.The equationstatesas,
This formulation derived by Andrew (Andrews, 2004) with regards on the constant
burning velocity through the flame travel. From Eq. 3.3 3.5, it is clear that the
-
dependantof dP/dton S,,.From this relationship,it can be seenthat KC1= constantS,,.In
dust explosion,the relationshipof K., andturbulent burning velocity, STexistedand this
can be relatedwith laminar burning velocity for dust i.e. ST= PSuas shown in Eq. 3.1.
53
P,
LD
a
Time (s)
Figure 3.1 Pressurebehaviorversustime for a closedandventedexplosions.
54
The basic principal of explosion venting provides for the rapid opening of a vent of
sufficient area to allow unburned gas and explosion products to escape,thus limiting the
pressurerise to an acceptablelevel. On the general basis, the reduced pressure attained
in the vessel before venting must be less than the mechanical strength pressure of the
containers. The vent opening must be correctly sized to allow the expanding gasesto be
vented at a rapid rate so as to limit the internal development of pressure.The acceptable
pressurerise is determined by the requirement that the vessel should not rupture and in
some cases,it should not deform. Essentially, a vented explosion can be considered to
have three phasesas shown in Figure 3.2:
A venting phase. A phase during which gas flow out of the open vent is
established.In this case,the pressurewill riseagainafter the vent cover hasbeen
removeduntil the flame front reachesit maximumareawhich denotesas P2 in
Figure 3.2. Cubbageand Simmonds(Cubbageand Simmonds,1955, Cubbage
and Simmonds, 1957) thoroughly investigated the development Of Predin
industrial ovens using town gas/air and methane/airmixtures by varying the
location of the vents and different breaking pressure of vent covers. The
significantresult from their work was the doublepeakpressureprofile attained.
Doublepeakpressure-timeprofile occurswith low breakingvent pressurewhilst
single pressurepeak profile occurs with vent cover having a high breaking
pressure(Harris, 1983).A higher rise in vesselpressureis observedfor lower
venting pressuresand this is explained by increasedflame area and a longer
period of unburned gas venting, which is less efficacious than burned gas
venting (Bradley andMitcheson,1978b,Harris and Briscoe, 1967).Effectively,
the total surface area of the flame front is greatly enhancedby the sudden
openingof a vent, and this in turn increasesthe overall combustionrate of gases
within the vessel.Thus, very strong secondpressurepeaks are often detected
following the initial pressuredrop whenthe vent opens.
Turbulence and stretching of the flame surface towards the vent opening can also
result in an increase in the combustion rate (Buckland, 1980). The turbulence may
be caused by increased gas flow velocities and by the opening of vent covers.
Burgoyne and Wilson (Burgoyne and Wilson, 1960) have demonstrated that, all
other things being equal, a vent which opens smoothly can give a lower Predthan
does a diaphragm which burst abruptly at the same opening pressure. The
very
presenceof a vent cover can increaseturbulence in an explosion,
external explosion that gives rise to the P3 is due to the fast flame propagation in the
turbulent external vented premixed gases. This pressure increases as the vent area
decreases or high K value as this will produce higher and turbulent jet velocity out
from the vessel. Usually P3 develops because the size and shape of the enclosures,
the position of ignition source as well as the stretching and movement of the flame
towards the vent as it opens, allowing burnt gas to be vented while combustion is
still taking place inside. Transient pressure may also be affected by oscillation being
set up in the system (Anthony, 1978) as represented by P4. This oscillatory type
pressure gradually damped out as the flame expands. However, the burning rate
during this phase is enhanced by the turbulence generated in the shear layer between
the outflow burnt gas and unburnt gases within the vessel (Cooper, Fairweather and
Tite, 1986). The strong interaction of the flame front, the shock/pressure wave and
physical back flow into the vessel generally enhance the rate of combustion inside
the vessel and further induced Taylor instability on the section of flame front surface
farthest away from the vent. As the flame is continued to expand, it eventually
encounters the walls of the vessel and this will suddenly decrease the flame area and
hence the rate of pressure rise. Work by Solberg et al (Solberg, Pappas and
Skramstad, 1980) confirmed the existence of third pressure peaks assembling the P4
trace in 35m 3 vessel in which the vent were initially covered. They found that a
strong flow field will be generated at the opening of the vent which the flame will
be strongly accelerated towards the vent and Taylor instability may develop at the
worst case of central ignition.
57
1ýr.
10
ci
m
E to
tA
vi
i2 4G
CD 21
-S
TIME/ s
The relative size of peaks P, and P2 are determined by the size of the vent relative to the
vessel, the magnitude of P,, the flame speed and the scale and intensity of turbulence
when the vent breaks. As the vent area of the vessel or enclosure is smaller, the pressure
required for vent cover removal is increased. It is found that the bursting vent pressure,
P, is inversely proportional to the area of the vent cover (Cubbage and Marshall, 1974).
Further, the smaller the A, i. e. high K, the smaller is the vent area available upon its
removal and the less rapid, therefore, the egress of combustion products were impeded.
Due to inertia of the vent cover, a finite time is required immediately after P, is reached
to totally remove the vent cover sufficiently far from the enclosure for the flow of gases
out of the enclosure and thus, the weight per unit area of the vent cover is another
important factor influencing Pred.The heavier vent cover is not as effective as the lighter
material in minimising the pressure developed (Cubbage and Marshall, 1974, Molkov,
Grigorash, Eber, Tamanini and Dobashi, 2004).
An increase in P, will cause P, to increase with respect to P2; while high flame speed
and high turbulence will cause P2 to increase with respect to Pl- P2 also increases with
respect of P, as the vent size becomes smaller for a given vessel (Anthony, 1978).
Further, at certain condition, there is only one peak is
observed. The merging Of P2 with
P, corresponds to a situation where dP/dt becomes
negative after the operation of the
vent which occurs when the vent is relatively large or flame speed is low. However, the
58
For the purposeof correlating the theory and experimental results, P2 will be assumed as
the dominant overpressure as P, is usually regarded as the bursting burst pressure. The
problem arises as most investigators of vented explosions do not report the pressure-
time profile but simply report the maximum overpressurethey attained in that particular
experiments. Hence, it is normally not possible to know whether this is caused by PI, P2
or P3.
correlations apply for a compact vessel which Bartknecht defined as those with LJD <2
although the draft European Standard (2007) has applied the correlations for UD < 3.
Bartknecht'scorrelation(Bartknecht,1993)is given asthe vent area,A, as a function
the vessel volume, V2/3 multiplied by a complex term that includes the mixture
reactivity KG,reducedoverpressurein the ventedexplosion,Predand static vent burst
pressure,P, The V213dependenceof overpressureon the test vessel volume is a
characteristicof spherical or compact vessel explosions,where the flame remains
mostly in sphericalshapeduring ventingprocess.If the sphericalflame propagatdsat a
constantrate irrespectiveof the vesselvolume, Predshouldbe only dependenceon Av
andV213or K. Bartknechtexpressedhis correlationas:
I
0.175(Pstat-0-1) VY3
-0.0567 (3.6)
=[0.1265logloKc,
pl 182 0172
p j
,.,
(bar] -
0
2,0 - ------------------ --- - --------------------------
\OON1, o\
1,0 - ---- ---------------
0.8 -
0,6 -
0 v=i rn,
3ý-
0-
--------------------------
ErgebnIsse:
CL Vý 2 M"
0,4
0
A V-10
0 V=30
m
m3
U, V-60 m3
%
noch kub'ischem Gesetz cl
0,2 euf 1 rri3umgerechnet 0\
['m 21
0,02 0,040,06 0.1 o,2 0,4 0,6 0,8
O,OB F
(bar)
Bartknecht
2,0 -
correlation
------------------- -----------------
line
0.11 111! II
ýII-ILI. 0,2 0.4 0.6 0,8(M"]
0,01 0.02 0,04 0.06 0.1
0,0a 1/K
Figure 3.3 Bartknechtdata for propane-airat different volumesat E= 10J,P, = 100
no other publishedventing data includedas notedby Siwek (Siwek, 1996)in his study
on publishedcorrelationsin ventedgasanddust/air-mixtures.
Figure 3.3 showsthe original data from Bartknechtfor propane/airat different volume.
From the figure, Bartknechtportrayedhis correlation basedon the experimentaldata
obtained for V= 10 m3- It suggestedthat the correlation developedby Bartknecht
(Bartknecht,1993)was purely empirical with no safetyfactor offeredfor ventingdesign
At it
purposes. the x-axis, representedthe vent area, denoted as F (m2) and Pm for y-
axis. However,when the result data is re plotted for V=2m3, representedby the (+)
sign, it is found that only one data agreeswith the line that drawn from the correlation
line. After re-analysed,it is found that F parameterat x-axis is a mistakeand it should
be correctedto I/K. To justify the alterations,the methaneexperimentaldata,denoted
as * sign from Bartknecht'swork for methane/airat V= 30 m3 were re plotted and it
showeda poor agreementwith the correlationline drawn. From the figure, it is implied
that the same vent area is required irrespective of the vessel's volume which is
obviously incorrect.Further,the resultsimplied that lower overpressureattainedfor the
sameK in the largestvolumei.e. 60 m3
NFPA 68 adoptedSwifts equationto apply for design and scale-upof venting system
purposes and it seem that the results obtained were in good agreementswith
experimentaldataobtainedby Donat (Donat, 1977)and Harris and Briscoe (Harris and
Briscoe, 1967). Again, there is a limitation in using this approach as it can be applied on
certain conditions i. e. Predis not more than 200 mbar and only suits to low-strength
enclosures.
CA,
A, = (3.8)
Pled 0.5
I= constant
(3.9)
K
The use of A, insteadof V2/3by Swift is followed by Bradley and Mitheson (Bradley
and Mitcheson, 1978a,Bradley and Mitcheson, 1978b) as Swift predicts that a larger
vent areais requiredfor cubic vesselrather than sphericalvesselat the samePredand
volumewhich is not in the caseof Bartknecht'scorrelation.This meansthat A, is the
key parameter to include the influence of the vessel geometry/shapeson the
The simplicity and easeof useof Eq. 3.9 is questioned,underestimatethe
overpressure.
complexity of large-scalecombustionphenomenaand did not take into consideration
the effect of vent openingpressure.Swift (Swift, 1983)basedhis design equationon
theoretical development that indicates the independence Of PredOf the vent opening
APred : ý/-S
versusa dimensionlessvent parameter, asfollows for initially closedventing
Le, with the ventcoveron place;
Pl, /S)-0,6993
d ý 2.43(A for P 2ý1barg (3.10)
Cd AC,, (p,,
with j (3.12)
s A$Su A
C,, is
where i the speedof sound(m/s), p, is the unburntgas density 3)
(kg/m and Pbis the
burnt gas density (kg/m3) The dimensionlessvent parameterhas been termed by
.
Molkov (Molkov, 1995, Molkov, 2001) as turbulent Bradley number, Brt where the
only different from the original equationis the addeddeflagration-outflow-interaction
number,X/g.The model suits ideally for sphericalvesselflame propagationas it is the
basisfor the equationsdeveloped.
If the Eq. 3.10 and 3.11 changedto K term as followed in Bartknecht and Swift's
equationabove,it gives,
1 0.28 (E - 1) contant
CdCuipred
, 1.43
= -ýred
1 43 (3.13)
K .
1=3.53S (E - 1) contant
,= Pud 0.5 (3.14)
Kc dC. ip, ýd0-5
venting of the unburnt gasesrather than the burnt gasesand this is taken as a 'worst
case' scenariofor a safe vent design. The pressurereduction enhancedby burnt gas
venting is supportedby Maisey (Maisey, 1965), Rasbashand Rogowski (Rasbash,
1986)andYao (Yao, 1974).
gas flow aheadof it and the presenceof openingvent alters the flow patterns.These
conditionsgenerateintenseturbulencetowardsthe flame propagationin which increase
the burning velocity and hencethe maximum overpressure.
At this particularcondition
where the postulateddiscrepancybetween the theory and experimentalresults in a
presenceof turbulenceassociatedwith the increasein burning velocity, Bradley and
Mitheson (Bradley and Mitcheson, 1978a, Bradley and Mitcheson, 1978b)
recommendedthat a five fold increasein burning velocity,abovethat of the laminar
value would benecessaryto fit the experimentaldata.
In the correlation,deflagration-outflow-interaction
numberor can be called turbulence-
factor,Z1,U and Bradley-Mitchesonnumber, Br are introduced where the so-called
Br, = Br ýu (3.17)
V3
-6; r z
65
A, c",
Br = V2/3 -
S,,(E - 1)
0.37
x (1+10V"')(l+O. 5Br) 1
0.9 (3.19)
ß11+x,
P,(bar absolute)
Where TrV=- Y is the specific heat
is the expansionratio Of Pu/Pb,
Pi ,E
ratio for unburnedmixture,c,ýjis the speedof sound (m/s) and R. is 'pi' number. As the
the
samewith othercorrelations, original equations (Eq. 3.15 and 3.16) will be changed
to K term. This givesas,
V3-
-6z S,,(E - 1) 1 constant (3.20)
TT
K Cid Pred 0.4 Pred0.417
2= 2
I= V3--6z 7- P,,
X S,,(E - 1) 7-P,, d d
(3.21)
JE-1y', u' constant
K C,,, 66
and flame acceleration induced by flame instabilities. Further, this method requires
users to determine turbulent Bradley number, Bt and thermodynamic (1,,, Ej, qj) as well
p -UA
As far as the author's concern,there are basically two different approachesof the
explosion venting flow modelling applied. In order to obtain the mathematical
expressionsfor flame flow and pressuredevelopmentduring venting explosions,it is
necessaryto establishthe appropriatethermodynamicand physics fundamentalsi.e.
conservationof mass,gas law, energybalanceetc. For a vent to give no increasein
overpressureother than that due to the pressuredifferencecreatedby the massflow of
unbumt gasesthrough the vent, the vent massflow rate is assumedto be equal to the
maximummassburning rate of the flame and this considerationshould be usedas the
design massflow through the vent. This approachin calculating the maximum vent
flow could be said to overestimatethe maximum flow as the assumption of the
maximumflameareais very conservative.
Mb
= SoA, p,, (3.23)
67
2
where S,,is burning velocity (m/s), A, is the internal surface area of the enc osure rn,
p,, is the (kg/m 3) The assumption is valid as the flame
and unburnt gas velocity .
the
approaches wall at higher Predas the flame area is closer to the wall. The maximum
vent mass flow and pressure loss always occurs before the flame exits the vent at which
during the period when the flame is expanding inside the vessel and pushing unburnt
gases ahead of the flame due to the thermal expansion of the burnt gases. Once the
flame exits the vent, the burnt gasesare free to escapeand expand outside the vent
the
whereas expansion inside the vesselceasesand hencethe flame speedsinside the
vesselwill fall to S,,.Further,the burnt gasflow throughthe vent is at a low gasdensity
and thus, will lower the pressureloss at the vent entrance.Thus, the prediction of the
maximum ventingoverpressureis reducedto that of predictingthe maximum unburnt
gasflow throughthe vent and its associatedflow pressureloss.
J'nb ý- SgAs
pa :--Su(E 1)A,
p,, (3.24)
-
Cd'FAv (3.26)
Mb= A, S, p,, = A, S,,(E - I)p,, ý-- (2p,,P,,,
d)0'5
1P a "d
I- (0.41+ 0.35 2) (3.27)
+ P,,,
K, 7(pi
d
p
redMW (3.28)
RT
where MW is the molecular weight of air (0.029 kg mol"), R is the gas constant
(82.0552 x 10-6M3bar/Kmol) and T is the temperature(K). The result of taking the
densityas a functionof Pedis to increasethe massburning rateand thus requirea larger
vent areafor the sameoverpressureor a large overpressurefor thesamevent area.
69
Again, the problem arises as temperature is also assumed constant. During the initial
portion of the explosion, the vented gases are cool. As the flame expands, the gases
flowing through the vent will soon have the high temperature of the combustion
reaction and basedon the gas law, high temperatureassociatedwith high pressure
attained.The assumptionof constanttemperatureis not valid but the computationof the
temperaturerise due to compressionis difficult without constantvolume combustion.
However, Lunn (Lunn, Crowhurst and Hey, 1988) has introduced approximate
temperatureasa function of P,ýd.
empiricalmethodto estimatethe compressed
r. 286
)/P
T= 298[(P,,,,+ Pi i (3.29)
( ) 1/2
A, 0.00243S, P,, Po
= d+ (3.30)
P,
A, CCd
-,d
( Pred
A, 0.00243S,,(E - 1) + PO
= (3.31)
A, -cCd
Pred
Noted that A, and V2/3is in relation in terms of A, =C jjV3 The constant value is given
.
in Table 3.2 below. For spherical vessel, A, = 4.84 V2/3and this will result to
( P,, ) 1/2
0.0118S,,(E - 1) d
+ PO
(3.33)
K -ECd
Pld
Equation 3.32 and 3.33 are suitable only when the flow through the
vent is subsonic.
When Pred > 1.89 Po, the flow is chocked and sonic flow is
occurred. For spherical
vessel it gives,
70
S,, P,.,, + po
dz
0.0145 1/2 (3.34)
KT
1)( Pred + PO
1 S,,(E -
0.014.5 1/2 ý
P,,
(3.35)
KT d
It seemsthat the underlying theory behind venting explosion makes other published
correlationsgiven by the previous investigatorscan be compared.It is important to
comparethe theory approachdiscussedabovewith the correlationsoffered by the vent
designstandardsi.e. NFPA 68 andEuropeanGasVenting Guidance2006.
P,
-,d=
for P., k1 barg
2.43(A/ý)-'0*69" Bradley and Mitcheson
(Bradley and Mitcheson,
P,. =12.46(A/ý)-' forP: gl barg
d 1978a, Bradley and
Mitcheson,1978b)
; r,,
d = Br, -2" for <1 Molkov (Molkov, 1995)
'T,, d
0*5 for
; rr,
d= 7- 6Br, 'Tred
>1
I ( Pred
+ PO
1/2
Method2
0.00243CS,(E - 1)
P,,
K ECd
d
I=C
-0.5817
(3.36)
K 3Pred
1=
Cl C2 E-1 Pr ed -"" (3.37)
K
= 7.98 and 0.078 for hydrogen with E=7.29. If these C1 values were compared to
Swift's C constant and C3 in Bartknecht's equation, it will give the value as shown in
Table3.2.
propane and 0.41 for hydrogen. To get agreementfor the turbulent effect or self
is
factor,0,
acceleration the ratio Of C3/CIC2
as Bartknecht'swork usedin
calculated
larger vesseland it was suspectedthat self-accelerationis already presentduring the
explosions.Thelist of predictedP is given in Table 3.3 below.
73
Cuhe
Methane 0.0097 0.062 0.164 16.9 2.6
Propane 0.011 0.071 0.200 18.2 2.8
Hydrogen 0.114 0.41 0.290 2.5 0.7
Sphere
Methane 0.0078 0.050 0.164 21.0 3.3
From the listed P value for both Method I and 2, it seemsthat Method 2 gives a good
which V is up to 1000 m3 where others are not. In case of P, > 500 mbar, the different
analysis will be formed in order to predict the suitability, availability and validity of
thoseequationson the experimentaldata studied.All data sets are distinguishedwith
shapeof the vesseland different volumeof the samevesselshapei. e. ci for cubic, o for
sphere,0 for rectangularand A for cylinder. This is done to show whether there are
independentvolume effects which is not been accountedin K and if there is any
differencesdueto the vessel'sgeometry,for the sameK andvolume
Table 3.4 KG value used for vent area calculation. All data were excerpted from Bradley
and Mitheson (Bradley and Mitcheson, 1978b) except for acetone-airdata (Molkov,
Dobashi,SuzukiandHirano,2000)
Fuel-air n-dxtures KG S, C P. YU0 7 bo
(bar-m/s) ( M/S) ( bar)1/2 A
Further, since the applicability of Swift's equationis UPto Pred 100 mbar which is
-
subsonicflow regime,the sonic flow will be appliedfor Ped> 900 mbar. The sonic flow
applied when the pressureupstreamof the vent reacheda critical pressure,Pc which
76
experimental data.
10 m v. 1m*
V. Ims
,,_*ýV. 22m3
0
ED V- 30m'
-- -- Bartknecht
rj - Swift
- Sonic flow line
AO 13
Cz A Sonicventing
.0 ........ .................... ................................................
Subsonicventing
CL
13
Methane -air
10 100
K
Figure3.4Methane/airexperimentaldataasa function of K.
00 V-1m,
13 V- 20
ED (D V-10m,
10 30M3
00V .
60M3
E3 E) V-
Bartknecht
Swift
..........Sonic flow line
Ea
(D O(D M0
(1) CO
00
ie
C13 Sonicventing
.......................... =din .... m..................................................
a.
/
//te
/'
Subsonicventing
/
/', e
Propane-air
0.1
10 100
K
Figure3.5Propane/airexperimentaldataasa function of K.
V=1m 3
10
Bartknecht
Swift
Sonic flow line
Sonicventing
C51 - I
....................................................
.0:..................................
*? Subsonicventing
CL
Hydrogen-air
10 100
K
10 * V. Irr?
V-1fe
14 Vc loff?
Badknechl
21 Sonicventing
ca I-
..............................
.. .........................
LP Subsonicventing
Towngas-air
10 100
K
Figure 3.7Town's gasexperimentaldataasa function of K
From Fig.3.8 and 3.9, turbulent enhancementfactor, 0 =1 and Cd = 0.61 are used in
order to correlate the same Cd constantusing by other comparativeequations.For
CH4/adr
10
A
* V--'D. 8 rTO(DragosaNAC)
....... Molkov
'0
Method 1
40 x
Method 2
0.1
1 10 100 1000
0 Wl 0 m3(Donad)
0 \k. 25 m3(8romma)
z V-10 m3(Barfinedt)
0 WM m3(Dork-d)
0 \A.40 m3ftomma)
V660 m3(Bardoectl)
x+
V660 m3(Donat)
f0
/" *1 V-70
0 0 m3(Bromma)
+0* V68lm3(Howards&KArabis)
00
\ýL-" m3ftomma)
+00
0 V--)M. 8 m3(P. F.Thome)
0:
4 Bartimedi
....... rvblkov
0.1
- MAIOd 1
0.1 1 10 100 off
Mtiod 2
Ka
From those figures above, it is clearly seen that Bradley and Mitcheson's equation has a
adequate and safe protection in venting design can be applied for Bradley and
82
valuesin practicalapplication.
For smoothly opening vents, Yao (Yao, 1974) recommendedto use 0=3 and for
bursting diaphragms,P=4 while Pasmanet al (Pasman,Groothuizen and Gooijer,
1974)gaveP =3 for small vent area( Av :50.1 rn2) and P =1 for the larger openings(A
> 0.1 2).
m Swift and Epstein (Swift and Epstein,1987,Swift, 1983)and Swift (Swift,
1983)indicatedthat, exceptfor the mostextremecases,a turbulentenhancementfactor
of 5 is sufficient to envelop the effect of vent opening pressureand used to give
deflagrationcharacteristicvalue, C in calculatingthe desiredvent areaof the protected
Further, Bradley and Mitcheson (Bradley and Mitcheson,
equipments/containments.
1978a,BradleyandMitcheson, 1978b)alsoindicatedthat in orderto get the theory and
experimentto be in agreement,a turbulent factor of 5 would have to be assumedas
beingcreatedby the vent. This recommendation applieson the initially coveredvent in
their correlation,known assafevent designfor ventedexplosion.
( ) 112
1=0.00243CPS,, P,, Po
d+ (3.38)
K ECd P,,
d
I=0.00243CgS,, 1/2
(E 1) Pred + PO
(3.39)
K -Cd
83
Figure 3.10,3.11 and 3.12 showed the exampleson how the possible P had been
processedin order to fit with the experimentaldata. The summary of the average
turbulentfactor,P for different vessel'sgeometryandfuels given by Table 3.5,
shape gas
P P P P P 0 P P P P P P
MI M2 ml M2 MI M2 MI M2 M1 M2 MI M2
Cubic 20 4 20 3 -20 5 - - - -
Sphere - - -20 4 - - 20 3 - -
Rectangular 20 4 3 >20 1 10 1 10 2 I- 7 2
-20 -
Cylinder 20 3 15 13 1- 1- 1- 1- 1-
j
The turbulentfactor, 0 tabulatedin Table 3.5 seemsto agreewith most of the previous
but slightly lower that beengivenby Swift (Swift, 1983)andBradley and
experimenters
Mitcheson(Bradley and Mitcheson, 1978a,Bradley and Mitcheson, 1978b).However,
if the valuesin Table 3.5 is comparedto valuescalculatedusing Bartknecht's constant
in Table 3.3, it can be said that both values seemsto have a good agreement.This
implies that the vesselvolumeis the importantrole in determiningthe Predin ventedgas
lo Tomygas
KM 10 T--19-
V-0,23 niG
0V-0.23 M3
am
.01 V-1.0 rTo
0V-1.0 rra
V-1.47 rTG a
V-1.47 ffl3
V- 28nG am
0V-28 n13
TLYOAM at 2
of 20
TL#tLAerit
0.1 4t
---- TurbUýt of 3
TabLdemd 10
CF6/air 10 C3Wair
10 mi
a V. 077 ,
I NE a
V-a77
.,, a
V-1.0 V-1.0
cr. , ..0
x V--20 x V-20
x x++ Ix
I ool x V-10LO V. lao
/ + V-2ZO V-220
a x a/ -x
V. Oao , V-MO
+ /+ a
+ V-60LO W8M
/, ' o
0.1 RrOLdert d 2D Tb, J«t of 2
+,, /, /, o
TLMLMM 09 3
+1
---- TLgbLA&ld 10 ----
0 TLIrtLAertCt 7 aa TLrbLffl of 4
I
1- .
.
. TurbJert d5 / " TUrtLiErt Gf 5
0.01 - 0.01
0.1 1 10 100 a1 1 10 100
However, it seemsthat Method 1 with high 0 do not give satisfactoryresult to fit the
experimental data. In this it
case, can be argued on the application of S,, (E-1) and Su
terms for demonstratingthe unbumt gasaheadof the flame characteristic.This different
approachto the maximum vent flow rate of unburnt gasesis a crucial in vent prediction
assumptionas the difference is betweenSg Sj- S,,
or and where usually a factor of 7 or 8
for hydrocarbonand hence(E-1) is - 6.5. As S,,(E-1) is close to Sf in value, Method 2
85
will give slightly lower vent massflow ratewhereit is clear that basedon flame speeds,
Sf. The correctionfactorsfor expansibility and density as a function of vent areaand
Predshownthatthe correlationfitted well with the datain the regioncomparedto others.
Both explosion parameters are dependent on the design of the vessel. In general, the
values Of Predand dP/dt are changed because of alterations in the amount of heat loss
from the flame to the vessel walls. The geometrical term i. e. size/shape of the vessel is
generally written as the ratio of the vent area to some other representative area of the
system (the intemal surface area of the vented enclosure, A, or its volume to the 2/3
power). However, the vent coefficient, K method is the most often used for the
assessment of vented gas explosion and this parameter will be used in all graphs'
plotting. The size and shape of the enclosure influences the final pressure generated; for
instance higher pressures will be generated in long narrow enclosures than in more
nearly cubical ones but in practice, it is the characteristics of the explosion relief that
will determine whether or not a plant is damaged by an internal explosion (Cubbage and
Marshall, 1974). These two approaches in correlating the venting design do affect the
on the dependence of overpressure on the test vessel volume. If the spherical flame
86
The useof A, seemsto be appropriatedto link with the vessel'sshapeas if the vessel
shapeis substantiallydifferent from cubic or when UD >-2, the flame touchesthe
wall before there has been a significant pressurerise and sphericalflame propagation
cannotbe applied to model the explosion.In practise,maximum pressureusually but
not always occurs when the flame contactsa wall of the vessel and the flame area
begins to decrease(Nagy and Verakis, 1983). In non spherical vesselswith central
ignition, Ellis (Ellis, 1928)showedexperimentallythat the shapeof the flame tendsto
assumethe shapeof the vessel.This observationis agreedby Nagy et al (Nagy and
Verakis, 1983)with cylindrical vessel.In a large sphericalvessel,the buoyancyof the
hot gasesmay distort the flame development.If the ignition source is not centrally
located,the flame front will contactthe nearwall before combustionis completeas the
flame is allowed to expand in one direction resulting in an elongated flame with
increasinglylargersurfaceareaandhence,fasterexpansionthan centrally ignited flame.
Consequently,the maximum rate of pressureand the maximum pressurewill be less
than when ignition is central becauseof greaterheat losses.For central ignition, the
flame will expandsymmetricallyin all directionsin sphericalshape.After sometime,
the flame beginsto distort when first flames contactingthe vessel wall and stretching
towardsthe ventbeforechangingits shapeto a cylindrical flame.
From the data collected,it is observedthose cubic and sphericalvesselsare the most
prior to reaching the vessel wall can be expected to have a significant effect on the KC!
value in very small vessels. The reduced flame area after contacting with the walls
influences the mass burnt rate which finally decreasing the amount of mass flow rate
through a vent. Increasing L/D ratio, i. e. smaller D for constant L, less mass burnt rate is
obtained which results to decreasePred.Thus in large LID vessel, the spherical flame
portion of the flame propagation occupies a negligible part of the explosion event. This
is why the standard cubical relationship involved in Kr, or K, (dust deflagration
t
constant) parameters only apply to compact or near spherical vessels. For spherical
vessel, 14% of initial mass or at overpressureof 1.12 bars for maximum pressure inside
closed vessel is assumed to be 8 bars, has burnt when the flames first touch the vessel
wall compared to about 7.5 % initial mass ( overpressure= 600 mbars) burnt for cubic
vessel. It was found that spherical vessels give higher Predrather than cubic vessel due
to a lower quenching effect to the vessel wall and since the mass burnt rate is proportion
to rate of pressure rise (dP/dt), it gives rise to an accelerating flame front and an
increasedP,,d. This scenario is shown in Fig. 3.13.
V-1 rro
10 - CF-Wair V=22 m3
x,
A- ýe V-Mrrß
0 Vml m3
- Bartkrech
/-0.
ty Bracley &
c]
#0
Mtcheson
All .... Wkov
Mathod 1
1a8 Mathod 2
li
O.j ý,,,,, -x, -,-
0.1 1 10 100 1000
K
Figure 3.13 Influence of vessel's shapeon methane/air explosion o Sphere and o cubic
equationregardedto be applied for Prd < 200 mbars, it shown that the application of the
equation can be extended. This is shown by most of the data fall under the Swift's line
as previously mentioned in Section 3.4. In rectangular and cylinder vessel, the majority
of the flame propagation is non-spherical as the flame is elongated to follow the vessel
shapewhere - 2.5% of the mass has been burnt when spherical flames reaches the wall
and hence, there are large heat losses subjective to convection heat. Since quenching
rate is higher, the reduced flame area influences the mass burnt rate which finally
decreasesthe amount of mass flow rate through a vent (Andrews, 2004). Thus for most
practical situations, there will be large heat losses from the explosion and venting will
occur after the flame has touched the vessel wall as shown in Fig. 3.14,3.15 and 3.16
on non-spherical vessels data scattered outside the published correlation lines even at
low K.
From Fig. 3.16, it was found that cylinder vesseldata gavebit scatteredplots for K> I
from the proposedequationsand published correlations.In a cylindrical vessel, the
flame shapeis divided into two periods.In the first period immediately after ignition,
the flameswerehemisphericand then half-spheroidalwith sectionsparallel to the wall.
In a very short period of time, nearly 80 % of the flame areavanishesbecauseof the
10 V-4). 17 rml
0
0 V=O.18 n-LA
0 V60-5 n-L3
* V--35 nO
0 V=40 n LA
0 V--70 m3
* VL-81 ryLl
* V--2W n LI
0 W-203,8 rrB
- Ban knecht
....... SM ft
---- NbIkov
-- Wt hod I
0.1 6 MLAhod 2
Figure 3.14 The influence of rectangular vessels on P,, d for propane/air 0 rectangular.
10 Tavmgas/air
00,
'0'
13 13 V=0.23
V6-1.47
\162.8
dL
0 13 V= 1 M3
Barfimcht
0.1 00
Brachey& Mtd-eson
, '0 NUkav
----
/"--- Nbtx)d 1
6 M3tx)d 2
0.01
01 10 100
K
Figure 3.15 The influence of various geometries on Ped for town gas/air 0 rectangular, o
cubic, o sphere
90
1,11Acetorelä r
V-0.19M3
Z0
A
V-OMM3
J
V- ZO rn3
Bm*nf)cht
-i ,
0/ , Smft
----
al Bradleyand
----
Mtcheson
.... Nbikav
Nbthod 1
Nbüiod 2
1 10 100 1000
K
As been discussedabove on how the use of V25/A, (K) term is failed to give
satisfactoryresults in non-cubic and non-sphericalvessels,the term A, /A, will be
replaced the K term in order to investigatethe differencesin using both terms for
correlating the vessel's geometry in venting explosion. It should be noted that this
analysis is only involved the whole set experimentaldata of methane/airin various
geometries.
The data resultswill first be shown the P,,d vK for cubic, rectangularand cylinder
vesselsas illustratedin Fig. 3.17. Again, turbulent enhancementfactor, P is included
accordinglyfor all correlationlines. From Fig. 3.17,it canbe saidthat cubic vesselgave
satisfactoryresultsfor all correlationsbut not with rectangularvessel.The datascattered
aroundthe correlationlines, but mostly aroundMethod 2 and Bradley and Mitcheson's
correlation lines.This suggeststhat both Method 2 and Bradley and Mitchesoncan be
applied for safeventing guide in venting design purpose.However, the cylinder data
seemsto fit well below the correlationlines.
91
clvýw
10 10
CLNC R
,c
0 WOý2-1 A/
0 W. 4-1
0 W16 nil
A
W2.8-3
V -I m3 \611.2 -1
V =22 m3 \ý20A M3
t.
x V---*) m3 w2b. 6 -1 *"
*
Bank-chi
....... Swill
CIVar
1()
: cywlier
li
ä \ý-1 m3
NIKhod 2
0.1 -
When the same set of data was compared with P, d v AjAv plots, it apparently gave
better trend. The cubic data results seem not given good agreement when A,,/A, term is
applied for V= 30 m.3 All rectangular data fall below Bartknecht's correlation line,
suggesting that the use of A, term is more favourable in the case of non-cubic vessels as
been postulated above. It also showed the same trend when applied to cylindrical vessel.
From the analysis, it can be said that the K term is more favourable when it applied to
cubic and spherical vessels as it fits the spherical flame propagation scaling law. The
satisfactory results on correlating the influence of vessel's shape is given when AJA, is
applied for non-cubic and non-spherical vessels as been proved in Fig. 3.18. It can be
A, instead V 2/3term to get more precise results if the
recommended to use term of
vessel's shape influence have to be taken into account in venting design.
92
10.00 CFVair
Rellancjd- V- 02 n13
V-0.6 m3
0 V-2.8rn3
1. 0V- 112M3
00 V-20SM3
0* V-26.6m3
/,
Bw4macht
00
Smil
0
1110 Bradley and Mlcrmon
-------
00.... Wdkov
0
kftiod 1
.a08 000
0A M*od 2
(101
0.1 1.0 10.0 100.0 1000.0
NA,
The data shows that Predexponent should be 0.5 in the subsonic flow regime and unity
exponent in the sonic flow regime are supported by the data. Bartknecht's 0.58
exponent is a compromise between the two regimes which not a good correlating
exponent and does not fit the fluid mechanics of the vent flow. The main data scatter is
due to the additional influence of volume. It can be summarised that the proposed
equation (Method 2) for safe venting has more advantages compared to others. The
analysis is based on the mass and momentum balance, only employing an empirical
relationship to account for compressible flow and sonic flow. The use of S,, (E-1) and A,
terms for non-cubic and non-spherical vessels to correlate the important parameters
involved in venting gas explosions seems to agree well with the published experimental
data. The simplicity of this 'safe' venting correlation can be used for any required
The datafor eachvolumetestedover a rangeof K was plotted as PredvK and later Prýd
at fixed K of 2,4,10 and 30. This allowedthe influenceof volumeto be investigatedat
constantK. There is insufficient data to undertakethese plots for town's gas and
acetone.Figure3.20 showsthat for methanethereis a very significantvolumeinfluence
for K<4 with no consistenteffect at higher K. The higherK datais mainly in the sonic
flow regime.At low K the increasein Predwith volumeis probablydue to the flame self
accelerationeffect,which extensivelydiscussedin NFPA 68.
Roo
too
00
01 "SOW P401
pressureof I atm(Chippett,1984).
air leaks within structures-,despite this the published literature on this topic is quite
sparse.Current literature on stratified explosionsdeals predominantly with buoyant
gases(DeHaan,Crowhurst, Hoare, Bensilum and Shipp, 2001, Liebman, Corry and
Perlee,1970).
95
K=2
100 K=4
K=10
K=30
10
l
1 Ai:.....
.........
Methane-air
0.1 llý
0.1 1 10 100
V(ff
K--2
1001
K-4
K=l 0
K--X
e"'
rIr
d
0.1
10 100
V(M)
K=2
10
K=4
X--
, - x- A K=l0
ET x K=30
Hydrogen-air
0.1
10
v(M)
A significant flame self accelerationeffect for subsonicventing was shown for K<-5
and this effect is similar to vent induced turbulenceand could be accountedby the 0
term in the burning velocity equation.The turbulenceenhancementpredictedbasedon
Bartknecht'sequationandproposedMethod 2 were in a good agreementwith 0 derived
from tabulatedexperimentaldatabasedon Method 2 as seenin Table 3.5. It can be
said
that the 0 derivedwas perfectly reasonablevalue for P as usedby other experimenters
(Munday, 1963, Pasman,Groothuizen and Gooijer, 1974, Yao, 1974). The name
experiencesdecelerationeffect in larger volumeswhen pressureis high in larger K and
this effect hasneverbeenhighlightedpreviously. It is postulatedthat at high K and Pred
97
It can be suggestedthat the useof K term is more suitableto be applied on cubic and
sphericalvesselfor UD =2 but failed to give satisfactoryresultsfor non-cubicvessels.
The AdA, term is more favourableto correlatethe influenceof vessel'sgeometriesfor
non-cubicvessels.The data and figures shown in this work also illustrated that P,,d
exponentof 0.5 in the subsonicflow regime and unity in sonic flow regime are
supportedby the experimentaldata. The main data scatter is due to the additional
influence of volume. The Use Of Predexponent of 0.582 in Bartknecht's equation
comprornisesbetweensubsonicand sonic flow regime which gave most of the data
scatteredoutsidethe line.
4.0 Introduction
From previous chapter, it is shown that the published guidelines, particularly with
regardto the questionof appropriatescaling of the correlation parametersare having
their own limitations and validity to use in practise.It is crucial to note that there is
anotherimportant factors i. e. turbulent, self-acceleration,volume and geometryeffect
influencing the mechanism in venting explosion which subsequently giving the
publishedcorrelationa poor agreementwith the experimentalresults. For the effect of
burstingvent on maximum pressurein ventedexplosion,most of the work was limited
to stoichiometricfuel concentrations(Bartknecht,1993,Cooper,Fairweatherand Tite,
1986, Cubbageand Marshall, 1974, Cubbageand Simmonds, 1955, Cubbage and
99
The geometriesof Test vesselI and 2 usedin the presentstudy were cylindrical vessels
which give the volume of 0.2 and 0.0065m3respectively.The vesselwas closedat the
end andfitted at the other with a circular tube with diameterof 0.162 m, simulatingthe
vent (eitherinitially openor covered)to the 52 m3dumpvessel,which effectively gave
free vent dischargeconditions(its volume was 260 timesgreaterthan the ventedvessel)
but enablethe work to be carried out under laboratoryconditions.The length of Test
vesselI is Im while 0.315 m length is appliedto Test vessel2. For Test vessel1, all
experimentswere ignited either at end wall ignition or at the centre of the vessel
centrelinebut in a caseof Test vessel2, only end wall ignition wasconsidered.Pressure
developmentwas monitored at various locations along the length of the test vessels
(Referto Fig.2.2). The distance,x from the sparkdivided by the test vesseldiameter,D,
x/D of the pressuretransducerswas kept approximatelyconstantfor all test vessels.
For Test vessel1, the,pressurewas measuredat threepositions,Po,P, and P2on the test
vessel and at the dump vessel,pressureis measuredby P6. For Test vessel 2, two
piezoresistivepressuretransducerswere mountedalong the test vessel,namely POand
PI and P6 situatedat the dump vessel.For maximum reducedpressure,P,,,,,,this was
taken from P, pressuretransduceras it locatedat the centreof the vessel for both test
vessels.The flame front motion was determinedusingbarebeadthermocouplesarrayed
alongthe vesseland the tube centreline(symbolsas TI-T3 in Fig.2.2). The flame speed
100
to the 0.5 min diameterthermocouplebead that was used.This preventedits use for
arrival detection technique is that it can be used with hydrogen. The alternative
techniqueof ionisationprobe flame detectiondoesnot work with hydrogen names as
they haveno ionisation.As part of the experimentalprogramme,three,repeattestswere
performed at each condition and these demonstrated good consistency and
reproducibility,with peakpressuresvarying by less than±5% in magnitude.
TEST VE SSEL 1
Fuel/air (D K Vent type Ignition
Methane/air 0.84-1.43 16.4 Open End/Centre
Propane/air 0.8-1.5 16.4 Open End/Centre
Hydrogen/air 0.34-0.54 16.4- Open End/Centre
- Ethy- lene/air 0.6-1.6 16.4 Open End/Centre
TEST VESSEL 2
Fuel/air (D K Venttype Pv Ignition
(mbar)
Methane/air 0.84-1.43 1-0,2.1,3.3, Magazine 98 End
16.4 paper
Aluminium 178
foil
Melinex 209
membrane
I OOgA4 424
paper
Propane/air 0.8-1.5 16.4 Open 0 End
Hydrogen/air 0.34 0.54 16.4 Open 0 End
-
Ethylene/air 0.6 IL 16.4 Open 0 End
-
102
4.2 Generalexplosiondevelopment
The vertical dashedline is the time at which the flame exited the test vesselmarked as
Tm. From Fig. 4.1, the flame initially propagatedin a slow, laminar phasewith low
pressure and flame speed before the flame approachesthe vent which took
approximately110 ms and 72 ms to travel from the ignition point of end and central
respectively.The averageflame speedis 22.8 m/s for end ignition and 16.5 for central
ignition which was considerablygreaterthan the laminar sphericalflame speedof about
3.0 m/s. Basedupon previousresearchon ventedvessels(Chow, Cleaver,Fairweather
and Walker, 2000, Cooper, Fairweatherand Tite, 1986) and records of additional
thermocouplespositioned in the radial direction of the vessel, the flame initially
developedhernisphericallyfrom the point of ignition at the end wall. Then asthe vented
flow field was set up, the flame beganto elongatetowards the vent and during this
stage,unbumt gaseswere also being ventedout to dump vessel.For centrally ignited,
the flame is initially sphericalflame, progressivelystretchedon one side towards the
vent (flame acceleratetowardsthe vent) andalso at the oppositedirection of vent. From
this explanation,it suggeststhat end ignition had a much larger flame area than the
centralignition and this reflectswith end ignition gavehigher overpressurethat central
ignition. For Fig.4.2, it hasthe sametypical pressure-timeprofile for venting at initially
openvent but the flame accelerates very rapidly towardsthe vent in respectwith the one
illustrated in Fig.4.1. It can be anticipatedthat suction effect is influencing the flame
Figure 4.1 Typical pressure-time histories for centre and end ignition in 'rest vessel I
O= 1.06,CH44r
0.30 0.14
End igrition TOA
V 0.0065 0.12
0.25
0.10
0.20
0.08 E
0.15
0.06
W Flame position
CD 0.10
0.04
0.05 0.02
O. oo 0.00
0.02 0.04 0.06 0.08 0.1
Tirre(s)
Figure 4.2 Typical pressure-time histories and flame position In Test vessel 2.
As can be seen in Fig. 4.1 and 4.2, the peak I)I-CSSLII'e When tile fl, 1111C
OCCUITC(I ZIII-Cady
exited 1rom the vessel cither when ignited at the end WýIll01-Mltl-ýIl. It is considered that
it is the variation 0[ tile Mass bi.11-11
1-,
ItC ilild I'lanic speed of' the Flame approaching, tile
that has a Strong IIII'luence on the vent I'low and on tile SLIhSCLILICIII
VC111 C011A)LIStioll
hCh,IVIOLn,.
A major feature of' tile explosions Is that there arc suhstantial proportions of'
.
104
the original flammablemixture in the test vesselafter the flame has exited the vessel
andtrappedat the comerregion insidethe vessel.This is largerfor centralignition than
for endignition asit takeslongertime for the combustionto takeplacebefore reaching
the maximumpressureas suggestedfrom Fig. 4.L The directionof unburnedgas flow,
due to gas expansionbehind the flame front, was preferentiallyin the axial direction
towardsthe vent,wherethe unburnedgasesaredisplaced.The inducedflow through the
vent, aheadof the flame leadsto a significant increasein flame speedsand expansion
ratio in the main vessel.If the ignition is initiated at the end wall of the vessel, it
resultedin an elongationof the flame shapewith a correspondingincreaseof its area
andthus,increasingthe burning rate andflame speedeventually.In the caseof centrally
ignited, the flame will be in spherical vessel initially before progressively being
stretchedon onesidetowardsthe vent andthus,reducingthe flame area.This condition
has beenpostulatedby Ferraraet al (Ferrara,Benedetto,Salzanoand Russo,2006). In
caseof centralignition, there is an indication of higher quantity of residual unburned
mixture in the vesselwhereasalmostcompletecombustionoccurredin the caseof end
ignition that eventuallyleadsto higherpeakpressure.
ventingexplosions.Pm,,,,
is shownas a function of the equivalenceratio in Fig. 4.3 - 4.6
for methane,propane,ethyleneand hydrogenexplosionswith end and central ignition
respectively.
Figure 4.3 showedthat P,,, is higher at all equivalenceratios, (D for end ignition
correspondsto onefor central ignition. The sametrend showedthat end ignition results
higherP.,,, comparedto centralignition for propane/airandethylene/airexplosion to
up
(D = 1.13 and 1.0 respectively.However, the univocal trend is not consistent where
higher P seemsto favour the central ignition for rich concentrationsin
...ý, propane/air
andethylene/airi. e. (D> 1.3.This inconsistenttrendcan be explaineddue to the mixture
reactivity itself Cooper at al (Cooper,Fairweatherand Tite, 1986) reported that rich
propanemixtures have much higher susceptibility than methaneto develop surface
105
0.40 CH4/air
0.35 6 End ignition
Centreign.itio.n.,
0.30
0.25
.01
0.20
cz
E 0.15
CL
0.10
0.05
0.00
0.6 0.8 1.0 1.2 1.4 1.6
Equivalenceratio
1.8 (ý3Halajr
ig ticn
1.6
Certre I ti
1.4
1.2-,
0.8
0.6
0.4
0.2
0
-1
0.5 0.7 0.9 1.1 1.3 1.5 1.7
Eqjvalence ratio
In order to get insight about the reasonof such behaviour,it should be basedon the
3.5 C2H4/air
2.5
E"
2.0
1.5
E
CL
1.0
0.5
0.0 ý-
0.5 0.7 0.9 1.1 1.3 1.5 11.7
Equivalence ratio
180] C3H./ajr
160 -ý F-- - ---,
ej Endignibon 1
140 Centre grlfioný,
-I-
120-
100-
80-
60
40
20
0---
0.6 0.8 1 1.2 1.4 1.6
Eqd%derr,e ratio
3001 Q2Po/eir
250 A Endl
ceml
200
1501
100
50
0
0.5 0.7 0.9 1.1 1.3 1.5 1.7
EquKderre retio
propane/air(above)andethylene/air(below).
6.0 1 H2/air
.6 End Ignition
5.0 41
--. I...
4.0
-cm
Co
m 3.0
;:ý
cu
2,0
1.0
0.0 4--
0.30 0.35 0.40 0.45 0.50 0,55 0.60
Equivalence ratio
ignition positionsfor (D=0.34 and0.41 beforeit gavea rapid increasein Pn,,,,at centre
ignition i.e. Pmax
= 2.7 barg (D
at = 0.54 compared to end ignition at higher equivalence
ratio. Combustion times for central ignition are roughly half of the end ignition,
thereforeless is for
time available gases in the vesselto flow out before combustion is
It also can be concludedthat end ignition gave higher P,ý,,,,in lean to stoichiornetric
concentration in hydrocarbon/air mixtures and P.. in centrally ignited is highly
dominantcomparedto P.. at endignition in rich propane/airand ethylene/airmixtures
((D= 1.1to 1.6).
The faster flame speedwith end ignition can be explained by the fact that the burnt
gasesare only allowed to in
expand one direction which in turn will result in an
elongated flame with increasinglarger surface area (reaction front) and hence faster
expansionthan the centrally ignited flames.This accelerationof the flame towards the
vent is also associatedwith self-accelerationof the hydrocarbon-airmixtures inside the
vesselthroughthe formation of cellularflames.
109
CH4/ajr
25.017
Endigiltion
Cartrelgiticn
20.0
llao
1
19.0
s
5.0
0.0-
0.8 OL9 1 1.1 1.2 1.3 1.4 1.5
EquivaJenceratio
Figure 4.8 Flame speedsof methane/airat end and central ignition as a function of
equivalenceratio.
CHAlair
30.0
(U. 1.05
25.0
---- WPA 68 --I
10LO
5,0 -jj
0.0
0.0 OL2 0.6 0.8 1.0 1.2 1.4
DsLarm from spark
140.0-1 Vajr
i=0.542
120.0- NFPA 6B
-*-
imo ia End igrition
-,er- Calre Igritionl
80.0-
6(10-
40.0-
20.0-
0.0-90
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Distarce frcm spa* (M)
C, 1
3Pold
,&,; n, WPA 68
FromFig.4.9, all graphsshowa first flame speedmeasurement at 0.16m from the spark
of around3- 4 m/s and acceleratingrapidly at 0.5 m from the sparkto the flame speed
of 28 m/s for methane/air,42 m/s for propane/airand 124 m/s for hydrogen/air
respectively.This is about5 times higher than the sphericalflame speedbasedon data
excerptedfrom NFPA 68 for thosegases.For centralignition, Fig. 4.9 showsthat the
flame acceleratingis more significant after 0.17 m from the spark.It should be noted
that suctioneffect is 'alsoplayed a significant role in stretchingthe flame towards the
vent,resultinglargerflame areaand massburningrate,thushigh flame speedinside the
vesselin respectto its own laminar burning velocity. The net effect is flame self-
accelerationoccurred at half diameter of the test vessel for hydrocarbon/air and
hydrogen/air.It is also clear that for end ignition the distanceto the vent is greaterand
hencethe flame accelerationcontinuous over twice the distance as for the central
ignition. Bradley et al (Bradley, Cresswell and Puttock, 2001) have shown that the
However,the theory of end ignition gives higher flame speedis not supportedby Fig.
4.10 and 4.13 for propane/airand ethylene/airmixturesrespectively.At (D = 0.8 and
1.0,end ignition gavehigher flame speedsof II and 20 nVsbeforethe flame speedsof
27 m/s attainedfor centralignition comparedto 24 m/s at end ignition at (D= 1.12. As
beenstudiedpreviously,propaneexhibits spontaneous cell structureor flame cellularity
in rich mixtures (Tseng,Ismail and Faeth, 1993,Wingerdenand Zeeuwen, 1983a).In
the caseof centralignition wherethe flame travel is lesscomparedto end ignition, very
little unburnt gas mixture has been vented from the vessel; hence there is still
substantialamountof unburntgasestrappedinside the vesselasillustrated in Fig. 4.11.
112
30.0 CýFVadr
End igriti
25.0 Certre ig
2D.0
15.0
10.0-
5.0-
Figure 4.10 Flame speedsof propane/airat end and central ignition as a function of
equivalenceratio.
TimeatPrnax -7
1
400 12 RamaleftIhevessel r, Test w-ssel 1
C3FVad
Ramaat the comerregion
WO
x.. Time at Prnax
3M X--- Flameleftthevessel
250 11 Flameat comerregion
2DO -
150-,
4
1001
50
ý- --X-
- .... eK................. .x,
_Erdigrition
Certre Igrition
EqLjvalerK>--
rabo
Unburnt grises
lu 0uu0u0
Ignition 00 D0 00 P0 000c0 00
cp 800C
00 CE00 0 ýO "0 00
position
0' 0 0
a000001 8 cl n- 00
co 00 00C. OC8 CP
00 0 00 0808808
C. 0n `4 cp
Figure 4.12 Diagramon the flame movementat end ignition and the unburnt gasesleft
inside the vessel
4aO C,2HWajr
35.0 B-d lgrItion
30.0 Certre igpitionj
25.0
20.0
iso
10.0
&0
0.0
-1
0.4 0.6 0.8 1 1.2 1.4 1.6 1.8
EqLivcdence
rctio
90.0 -, FVajr
80.0
70.0
60.0
50.0
40.0
30.0
20.0
10.0
0.0 4-
0.30 0.35 0.40 0.45 0.50 0.55 0.60
EclAvaderce
ratio
0 -7
0.3D 0.35 0.40
...........
0.45 0.50 0.55 0.60
EqLivalance
rafio
within the vessel;as the burning velocity increases,the magnitudeof the maximum
pressurepeakincreases,whilst the time of occurrenceafter ignition decreases.It also
appearedthat maximumpeakpressureoccurredafter the flame exited the vesselwhich
suggestedthatthereis the variationof the massbum rate and flame speedof the flame
approachingthe vent which has a strong influence on the vent flow and on the
combustionbehaviour.
subsequent
2.0
TOUICH4
1.0
0.0
Figure4.16 Pressure-time
history for different fuels at (D= 1.0
plug with a 16A there is only a very low probability that a deflagraýtion
to detonation
will occur in a mannersimilar to that observedin elongatedtubes.It can be thoughtthat
the externalexplosionmight causethis phenomenato happen.However,if this occurred
then all the pressuresin vesseltest and the dump vesselwould increaseand eventually,
Fig. 4.17 showsthat this doesnot occur as there is no significant
pressuredifference
betweenpressureinsidethe vesselandpressureinside the dumpvessel.
118
6.0
5.0
'a 4.0
m
3.0
(n
2? Z.U
"
a.
1.0
0.0
0.005 0.01 0.015 0.02 0.025 0.03
-1.0 Time(s)
Figure 4.17 Pressuretime history for maximum pressureinside the vessel (PI) and
pressureinsidethe dumpvessel(P6)at stoichiometrichydrogen/air
Figure 4.18 compares the maximum pressure, P.. of all fuel/air mixtures used in the
tests as a function of equivalence ratio. It showed that the highest Pma,was recorded and
observedat (D= 1.0 for methanelairand propane/airand (1)= 1.18 for ethylene/airand
hydrogen/air.It can be said that the result hasa good agreementto other investigations
reported that at stoichiometric and slightly off-stoichiometric gas/air mixtures,
maximumpressureinside the vesselexhibits higher pressurecomparedto the near and
very rich concentrationas at theseconcentrations,laminar burning velocity is at the
highest values. It is interesting to note the highest P., for hydrogen/air explosion
recordedat 35 % concentrationby volume ((D= 1.18),not at 40 % concentration((D=
1.36) that beenreportedto have highest maximum overpressureand laminar burning
velocity (Andrewsand Bradley, 1972,Andrews and Bradley, 1973,Kumar, Dewit and
Greig, 1989).It shouldbe notedthat all maximumpressureobtainedfrom hydrogen/air
testswererecordedbasedon the peakpressuresnot the suddenpressurespike.
119
3-0 - CH4/edr
--X-
o C3HB/air
25-, C21-Wair
1ý
j --x- H2/är j
20
1.5
1.0
0.5
0.0
0.6 0.8 1.0 1.2 1.4 1.6 1.8 2-0 2-2
Eqdvalercerafio
1/2
1
As
cl,
05
P,,' E 1)
which C, ,.
Cd.e.20"
(4.2)
121
From the calculation,it is found that the flame areais about40 - 45 % of the vessel
surfaceareafor hydrocarbon/air when the pressuremarked its P,, but only 10 % of A.
a,,
in the case of hydrogen/air. This is the valid explanationfor the occurrenceof
detonationspikeobservedin the testvesselthat will be discussedlater.
Endigribon
60.0-, H,
1v0.0065 rrý 21air
50.0
4aD
30-01
Cl-Wair
10.0
0.0
0.0 0.5 1.0 1.5 2-0 25
EgAvadercerzdio
The increasein burning velocity ahead of the flame due to turbulence might be
attributedto the greatersurfaceareaof the outer boundaryof the flame. As turbulence
develops,the small scale turbulence begins to influence the mechanismsof flame
propagationand further enhancesthe burning rate (Abdel Gayed and Bradley, 1981).
Basedupon the previous investigationon ventedexplosion (Cooper,Fairweatherand
Tite, 1986),the flame initially developedhernisphericallyfrom the point of ignition at
the end wall. In the presenceof the obstacle(in this case,the vent), the flow sets up a
gradient filed leading to subsequentflame-area increase due to stretching, The
interactionof the flame front and inducedrapid turbulentgenerationassociatedwith the
instability mechanism.However,as beenobservedby McCann et al (McCann,Thomas
andEdwards,1985),cellular instabilities do not appearuntil the late stagesof the name
developmentin smaller vesseland thus, they have no significant effect on the initial
flame and pressuregrowth. The induced turbulent also triggers a rapid combustion of
122
From Fig. 4.20, the onsetof the denotationspike occurredat (D= 0.76 with Sf = 28.6
ao -i H2-air
7.0
ao
510
4.0
ao
20
1.0
0.0 i,
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2-0 Z2 2-4
EqLi\rafercerabo
4.0 -1 FVair
1 flý- 1-1+1-
3.51
1
3.01
2.5
2.0
1.51
i
1.0
0.5
0.0
0.04 14
ý
-0.5
Time(s)
0= 1. QC«2Wär
3.5 -
2.5-
,62.0-
:3
CL 1.0-
0.5-
0.0
0 0.02 0.04 0.06 O.OB 0.1 0.12 0.14
TIMM
comer region of the vesseland both top and bottom of the vessel.As shown in section
125
region of the Test vessel I at the spark end, where a flame arrival thermocouplewas
locatedas shownin Table 4.3. The resultsshowedthat the time of the flame arrival in
this comerregion was very closeto the time of the pressurespike's occurrencefor end
ignition but eventually not in the caseof centrally ignited. It is consideredthat the
comerregion is not a high turbulencezoneandhence,the auto ignition point or hot spot
is the bestdescribedfor the observation.Similar observationis reportedfor large scale
test using 35 % H2/air concentrationwherethe localisedexplosionoccurredinside the
vesselafter the outflow of flame in venting explosion (Dorofeev, Bezmelnitsin and
Sidorov, 1995). However, this phenomenais not observed to propane/air and
methane/airin both vessels.
4.4 Influence of vent coefficient, K, volume and burst vent pressure,P, on Pmax
A= 1
K= (4.3)
A, 1- BR
whereA, is the vent areaand A is the areaof the orifice plate. It is illustrated by Fig.
4.23 that the maximumoverpressurewill be reducedaccordinglywith the size of vent
areai.e. lower K will give lower Pmax in
which this case,Pn,,,,= 0.0178and 0.171 barg
for K =1 and 16.4respectivelywhich about 10timesdifferent in P,,,.. This is due to less
flow restrictionto spell out the unburnedmixtures from the vessel.Since all test has
doneon the initially openventing,there is a largedifferent on P"'a,betweenK=I and
K= 16.4.The maximumpressure,Pmaxof 0.0178achievedin the unrestrictedventing
case(K = 1) happenedat t= 57 ms and interestingly,secondpeak pressureoccurred
later at t= 140 ms. This secondpeak exhibited sustainedoscillatory pressuredue to
situation wherethe vent is relatively small and flame speedis higher in which dP/dt
remainspositiveafter the burnt gasesbeenexpelledfrom the vent(Harris, 1983).
0.25 1 CH4/air
o= 1.05
ü2D
K- 16.4
0. e
CLOC
-ao5 lak;*j
1111
out inside the vessel.The net effect is lesseramountof unbumt gasesis left inside the
vesselandeventuallywill give lower P.,,,,as shownin Fig. 4.24.
128
mixture.
H2/air
7.0 Tw
6.0
5.0
-ý3 4.0
ACo
3.0
2.0
1.0
0.0
0.01 0.02 0.03 0.04
Time(s)
other than K. However, there is not what the experiment show in the case of the same K
for different volumesas shownin Table4.4.
From the listed data in Table 4.4, it can be said that in larger volume, 0.2 M3 in the
presentstudy,it is obvious that self-accelerationis the important featurein increasing
the P,,. in Test vessel 1. Comparingthe ratio of Pmjtxj/Pmax2
in Table 4.4 with the
turbulent values calculatedfrom the original Bartknecht's work to the predicted Prd
from Method I and 2 given in Table 3.2, it can be consideredthat those values gave a
good agreement.This implies that the turbulent factor, 0 can be regardedas name self-
accelerationfactor occurred during venting. It can be postulated that the ratio of
PmaxIlPmaX2canbe consideredon how fast the flame acceleratesinside bigger vessel,and
self-accelerationis the main factor of the fast flame propagationandthe high pressure.
130
Table 4.4 Summaryof experimentalP.. for Test vessel I and 2 for K= 16.4. The
ignition position is end ignition.
betweenP,,,,,,,
and Sf,,,s in Test vessel I and 2 were about the samebut significantly
higherfor Sfý,,ratio for hydrogen/air.
g
and Edwards, 1985) that flame cellularity is appearedin the early stage of the explosion
in larger volume compared to the smaller volume and hence,influence the mass burning
rate and P.,,, inside the vessel. Further, the ratio Of SfavgJ/SfAvg2
is about 2 to 3 for
As discussedin Chapter3, the pressureat which the vent relief cover, Pvbeginsto fail
has a more significant effect on the magnitude of the P, and P2 (Figure 3.2). The
opening/breakingof the vent would delay or hinder the venting process,causing the
maximum burning rate increasingdue to the bulk flame areacompressedtowards the
vent and hence,the rates of pressurerise increasesas well as the pressureinside the
vesselin comparisonto the openvent mechanism.At low openingvent pressureP,, the
resistanceto flow expansionout of the chamberto the duct is less than at higher P, At
133
higher P, sincethe vent openedat a relatively late stage,whenthe total flame areahad
increasedsignificantlycomparedto the lower P, case,the rate of burnedgasproduction
exceededthe rate of unburnedgas venting which in turn causeda continuation in
pressurerise within the vessel(Chow,Cleaver,FairweatherandWalker, 2000).
0.9 -1 e.-Li
-
a8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
0 50 100 ý 150 2DO 250 30D 3w 400 450
P, (barg)
However, this condition is only satisfied the previous works for low K i. e. K =1 and 2.1
but not with larger K in this test configuration. It is apparent that a decrease in P,,,.,,
occurred at K= 16.4 from P, = 209 mbar to P, = 424 mbar, about 1.5 times different
and there was no significant increase in P,,,,,.for P, = 209 and 424 mbar respectively at
K=3.3. It should be noted that the concentration used (10 % methane/air) was the
slightly off stoichiometric which showed the higher P. in the tests conducted at
a.,
varying equivalence ratio. This similar phenomena has been observed in Bartknecht's
134
work (Bartknecht, 1993) and Harris and Briscoe (Harris and Briscoe, 1967). The
occurrenceof this peculiarbehaviouris dueto two opposingeffects.At higher P, even
if the flow accelerationis stronger,the explosionis vented at a later stage when the
flame is closer to the walls. Thus, there can only be small increase in flame area and
almost immediately followed by gas cooling and flame quenching at the nearby wall.
This results a small increase in burning rate and hence, the overpressure inside the
vessel (Ponizy and Leyer, 1999a, Ponizy and Leyer, 1999b). Similar result has shown in
Leeds work for larger cylinder. Heat losses after the flame has touched the wall make
the peak pressure decreases(Phylaktou and Andrews, 1991, Phylaktou, Andrews and
Herath, 1990). Experiments with turbulent gas mixtures at large vent showed that the
maximum pressure increased rapidly with high pressure of the bursting disc compared
to smaller vent (Harris and Briscoe, 1967). Turbulent could increase the burning
velocity and cause the increase in maximum overpressure. Hence, an explosion vented
in its early stageswill give a higher pressurethan one vented at a later stage.
Again, the sameK resultsthe sameP,,,,,for different volume as shownin Table 4.6 for
experimentalresultsfrom Test vesselI andTest vessel2. It is recommendedto usethe
As term to replaceK in order to provide preciseestimationon venting
area and thus,
Pma,,for given vesselvolumes.
136
Table 4.7 Percent error deviation for experimental and calculation result
Test vesselI
Gas/air Method I Method 2 BartknechtEq Swift Eq Bradley and Molkov
Mitcheson Eq Eq
C21-14/air
-96.9 -7.2 256.9 583.7 9.5 0.3
I I-
H2/air
- - -
I- I
Test vesseli -- -
K=1 -1
Table4.8 cont
vessel I in P,
comparisonwith ax in Test 2
vessel at the sameequivalenceratio. This is
justified by the reversedcalculationdoneusingMethod 2 equation.It shown that about
80-90 % of flame areahasbeenoccupyingthe vesseltotal surfacearea.It is confirmed
the observationreportedby McCann et al (McCann, Thomasand Edwards, 1985). In
their work, they said that the flame cellularity (self-acceleration)appearedin earlier
stage in larger volume and this give significant effect on the overpressureinside the
vessel.For smallervesseli. e. 0,0065 M3, experimentshave shownthat the presenceof
pressureoscillation and this coupling with the induced turbulence by the vent flow
increasedthe P... For centrally ignited explosion, the increaseintensity of the flame
cellularity during flame propagation produces accelerating flame front which later
interact with the vessel wall. Due to the rapid decelerationof the flame front as it
approachthe vesselwall, it results on the strong rare fraction waves which triggering
further combustionof a large amountof unburntgasesleft insidethe vesseland lead to
a significant P.,,, inside the vesselfor reactivegas/airmixturesasshown on ethylene/air
and hydrogen/airmixtures. The ratio Of Pmaxl/PmW
also illustratedthat there was about
139
2-7 times increasesin overpressurein Test vesselI comparedto Test vessel2 at the
sameequivalenceratio and K, suggestingthat flame is self-acceleratedin bigger vessel
andpropensityof cell cellularity is susceptiblein rich mixturesof hydrocarbon/air.
5.0 Introduction
This is a consequence
of the need to overcomethe inertia of the air within the duct,
beforean outwardflow of gasesfrom the vesselcan be establishedand pressurerelief
commence.In practicetherefore,a dischargeduct should be as short as possibleand
havea cross-sectionalareaat least as large as the areaof the relief vent. It should be
also preferablybe the sameshape,Since evenshort ducts can lead to doubling of the
overpressuregeneratedin the vesselbeing protected,it is suggestedthat whenevera
dischargeduct is usedin conjunctionwith the relief vent, the effective strengthof the
plant should be taken to be half of the value originally assumed.The strengthof the
dischargeduct itself shouldbe sufficient to withstandthe anticipatedexplosionpressure
andthepressureassociatedwith any explosionsthatoccurwithin the duct itself.
P,,
d with theduct= 0.779 ( Pred
1.013) '-'6'+1.013 Lt<3m (5.1)
withoutthe duct -
Pmd
with theduct -0,172(pred 1.013)*9" + 1.013 3m: 5Uý: 6m (5.2)
withouttheduct -
version NFPA 68 (NFPA68,2002), the only differencewas the value of the constant
used. The previous equations given for different duct length were provided by
Bartknecht(Banknecht,1993,Bartknecth,1981)asthe correlationsof his experimental
resultscarriedout in aI m3 explosionvesselwith a vent burstpressureof 150mbar and
a vent diameterof 0.2 rn or K= 33.3.The equationsare,
p-red 0.8614
with theduct -I 124 Pred withouttheduct Lt<3m (5.3)
P, 0.5165
d with theduct -2'48predwithout theduct 3m:5 Lt ý!6m (5.4)
barg,Eq. 5.1 predictsthat the additionof a vent pipe of the samesize asthe vent would
increasethe pressureto 1.27bar and to 1.04bar for Eq. 5.2. This is about0.56 bar for
Eq, 53 and 1.54given by Eq. 5.4. The presentwork investigatesaK of 16.4 and0.162
m diameter and 1.0 m long vent pipe (LJD of 6.2) for a 0.2 m3 cylindrical vessel.
However,the resultsthat will be shownlater give over conservativeagreementwith Eq.
5.2 and no agreementwith Eq. 5.3 but betteragreementwith Eq. 5.4. Therehasbeenno
validationof Eq. 5.1 and 5.2 for different K or vesselvolumeseven thoughthey are the
basisof the US andEuropeanStandardsfor the useof ductpipes.
The increasein P,,,,,,with the additionof a vent duct is so largethat vent ductscannotbe
used without increasingthe vent area and duct size to achieve a reduction in the
overpressure.However,there is insufficient design data for gasesto enablethis to be
done effectively and the physics of the processfor gas explosion venting is not well
understood.This contrastswith the situation for dust explosions,where a substantial
experimentaldata baseexists (Lunn, Crowhurst and Hey, 1988). Recently, Tamanini
and Fisher (Tamaniniand Fisher, 2003) derived a correlationto take into effect of the
duct when the duct explosionis occurring.All of thesecorrelationsexplicitly take into
accountthe dependence on the duct length,duct diameterandthe vesselvolume.
presenceof the duct to discharge the explosions products generally increases the
severity of the explosion comparedto the situation encounteredwith a vessel vented
directly to the atmosphere.It has beendemonstratedthat the deflagrationpressurecan
exceedby a factor of 10 or more than obtainedin experimentswith venting directly to
atmosphere(Bartknecth,1981,Molkov and Nekrasov, 1981).This is mainly due to the
complex interactionof the combustioninside the duct pipe and vessel. Despite of the
numerous experimentalstudies and their findings ( Bartknecht, 1993, DeGood and
Chartrathi, 1991, Ferrara, Willacy, Phylaktou, Andrews, Benedetto, Salzano, 2005,
Kasmani,Andrews,Phylaktou,Willacy, 2007, Kordylewski andWach, 1988,McCann,
Thomas, Edwards, 1985, Molkov and Nekrasov, 1981, Ponizy and Leyer, 1999a,
1999b,), the whole picture of the interactionbetweenthe gasdynamics inside the duct
144
and the vesselas well as the geometry and operating conditions on the maximum
pressureremainsunclear.
The additionalfrictional loss due to the presenceof the duct and the duct gas column
inertia are also consideredin reducing the venting flow rate, dependingon the duct
diameter and length (Ferrara, Benedetto, Salzano and Russo, 2006). Bartknecht
(Bartknecth.W, 1981) pointed out that the increaseof the explosion violence and the
peak pressurein the vesselare strongly affectedby the hydromechanicaldrag and gas
column inertia of the duct. The sameobservationhas been describedby Ponizy and
Leyer (Ponizy and Leyer, 1999a)that the frictional lossesand inertia are the key factors
in influencing the evolution of the explosion inside the vessel.This finding has been
confirmedby Ural (Ural, 1993)usingmathematicalmodel.
changesthe larninar burning velocity, S,,, of the flame and the mass burning rate of
sphericalflamesscalesas SU3 This directly influencesthe velocitiesin the ventpipe and
.
henceon the influenceof the vent pipe in the explosion.Thus a study of venting with
In the previous chapter, the details explanation for the physics and dynamics of the
simply vented gas explosion were given to get better understanding on the venting
mechanism. Figure 5.1 showed the pressure traces for simply vented and duct vented
explosion for methane/air at (D = 1.06 for end ignition. It illustrated that Pmaxfor vented
explosion with the duct fitted is about 4 times higher than for the corresponding Pmaxin
simply vented explosion. A similar trend is also showed for the maximum rate of
pressure rise, dP/dt, It is interesting to note that there are two peaks for dP/dt (ductvenwd)
traces. It shows that the maximum burning rate reaching the first peak after the flame
exited the duct pipe and it suggeststhat the combustion is still taking place after a rapid
decreasein dP/dt to give the second peak before dropping to barg s"' and finally
-80
deceasing to atmospheric condition inside the vessel. The maximum peak pressure for
both casesoccurred at the same time but at different magnitude. This can be proposed
that they have the same mechanism but the presence of the duct has triggering the
pressurerise inside the vessel.
146
1.4 O= 1.D6
End igrition 150
1.2
1.0
dP/cR SirToyv3"fted 100 m
08 _,
.
0.6 50
0.4
0
0.2
0.0
ý). 04 0.06 0.08 0.1 0.12 0..14 6 0.18 0.2 -50
_().
-0.4 Tin-u(s) -100
Figure 5.1 Pressure-time histories for simply vented and duct vented explosion for
17,
methane/airat 0) = 1.06 for end ignition (initially open venting). dP/dt traces were given
for both cases.
In order to get insight about the mechanism of vented duct eXPIOSiOll, tile flow and
cornhustion dynamic is presented in Fig. 5.2. Again, tile peak overpressure, P,,,,,, at
Pve. OCCUrsafter the flarne has exited the vent pipe, and it happened to all explosions
ýsej
with or without the vent bursting pressure, Pv observed from the experiments. The peak
pressure does not Occur WlIC11the flame is in tile duct as postulated by Lunn et al (Lunn,
Crowliurst and Hey, 1988). After the slow, larninar phase of flame propagation in tile
vessel, the flanic enters the duct (t = ti,,) Which Was followed by a Sudden increase in
pressure at the duct inlet caused by a very fast propagating flarne in the duct (Sjj,,,, = 350
in/s). This fast flame and associated sonic flow condition encounters the strong
tUrbuience field a few diameter of tile vena contracta which leads to a subsequent
violent combustion or known as 111.11-11-up.
The bL11-11-Up
phenomena or known as
secondary explosion lias been investigated
L,
by Bartknecht (BartknectliM, 1981), Poilizy
and Leyer (PonIzy and Leyer, 1999a, 1999h) and Molkov et al (Molkov, Dobashi,
SLIZUkiand Hinino, 1999). They identified the hurn-up phenomena and deduced Owt the
llow reversil resulting froin the action is tile main IIICCIIýInlsfn responsible for tile
This strong combustion inside the duct causes the massive reversal of the pressure
(negative pressure drop) towards the test vessel and this marked as AP2 This flow
ý.
pattern is associated with flow pressure (energy) losses which are characteristic of the
AP
(5.5)
PS9
2
study), it is of the order of 0.5 (Franzini and Finnemore, 1994). The maximum value of
this pressuredifference relating to the induced unburned gas flow was measured just
5.2. Substituting 150 mbar into Eq. 5.5 using density of p=1.2 kg/m 3 (ignoring any
small pressure rise and any compressibility effects), the calculated unburnt gas flow
is
velocity around 230 m/s. If this calculation is compared to the flame speed inside the
This means that there is a
duct of 350 m/s, the calculated Sgis 280m/s ( Sg= 0.8 Sfiame).
very large flow velocity into the duct pipe and thus will generate a very turbulent flow
field within the duct.
The ignition position clearly has not fundamentalinfluenceon the physicsof the duct
vented explosion,not does the static vent burst pressure,and although both do in
generalincreasethe peak overpressurewith end ignition having higher overpressures
than central and increasing the vent burst pressure usually increases the peak
overpressure.However,in somecasesthis doesnot occur.It is consideredthat it is the
variation of the massbum rate and flame speedof the flame approachingthe vent that
hasa stronginfluenceon the vent flow and on the subsequentcombustionbehaviour.A
major featureof the explosionsis that there are substantialproportionsof the original
flammable mixture in the test vesselafter the flame has exited the vent duct. This is
larger for centralignition than for endignition. It will be shownthat the initial vent flow
reachessonicconditionsandhencethe vent pipe is choked.Principally, the vent flow is
a linear function of the internal vesselpressureand from the relationship,the internal
vessel pressureincreasesuntil the mass of vented gasesreduced the vent flow to
subsonicandlowerpressureloss occur.
149
T. T.
2.2 A 2.2
1.4 1.4
A P sý6
1.0 1.0
0.2 0.2
0 02 0 04 O OB O OB ol 0 12 0 1r IO. 2
-0.2 . . . . . .
TmEXs)
AP3ý5
-0.6 AP2ý -0.6
Figure 5.2 Pressure records at selected positions along the test vessel for methane/air at
q) = 1.06 for end ignition. AP2-3IS the pressure difference at the vent entrance, AP3-5is
pressure difference inside the pipe and AP5-6is the duct exit pressure loss. Tin and T,,,,,
indicate the time flame enters and leaves the duct respectively
pressurewave which interacts with the flame front to distort it and henceincreaseits
surfaceareaand massburning rate. In this manner,the initial flow of gasesinto the duct
will be larger than for an initially open vent. Thesecombinedeffects ofturbulence and
pressurewaves createdby the vent bursting result in accelerationof the flame prior to
the vent duct and also increasethe flow velocity, turbulenceand flame speedin the vent
duct. Testswith ductedexplosion ventsgenerallydisplay Helmholtz oscillations, that is
the pocket of burned gas within the vessel undergoesbulk motion towards and away
from the vent opening, due to the massof the duct contentsand tile compressibility of
the gas in the primary enclosurethat acts as a spring (McCann, Thomas and Edwards,
1985).
In order to investigate the influence of the P, on four different vent covers which
P, ranged between 0 to 500 mbar were used. Comparison results between open venting
and closed venting was also studied to understand the physics of' vented explosion
mechanism on each cases. The variation ol'the maximurn over pressure, Pniax,with static
150
burstingpressure,P, is shownin Fig. 5.3 with centraland end ignition for (D=1.0 for
propane/airandmethane/airmixtures.
Figure 5.3 shows that all the present results were well below the correlation of
Bartknecht (Bartknecht, 1993) which has been adopted in NFPA 68 (NFPA68,2002)
and the draft European explosion venting guidance. The results at 100 mbar burst
pressure are at least 2 bars below those P,,,. predicted by the design correlation. Figure
5.3 also indicates that the present results for the influence of the vent static burst
s CF-WairEdend
14.0-
w CF-Wedr.
1 certre
6 C3HB/airat end
12-0-
?i MHB/dr. Edcerter
0-14, NFPA 68
10.0 C3-IB,WPA 68
8.0
4.0
zo
0.0
0 0.05 0.1 0.15 0.2 a25 0.3 0.35 0.4 0.45
PV(.
M)
much larger with a static vent burst pressurefor propaneand the effect is reversedfor
methaneonce a vent is in position. None of theseeffects are reflected in any vent
guidance(NFPA68,2002)or havebeenreportedby others.
The possibleexplanationis since the P, is high, so that the time the vent cover is
removed,the pressuredifferential acrossthe vent openingis larger and subsequently,
the pressureinsidethe vesselwill fall rapidly. In fact, the momentumof the gasesbeing
ventedmaybe sufficient to causethe pressureinsidethe vesselto fall temporarilybelow
atmosphericpressure(Harris, 1983).Molkov (Molkov, 1994)also observedthe same
trend in his work. At low P, values(P, < 0.2 MPa), the peakpressureincreaseswith
increasingP, but not at high P, (P, > 0.2 MPa). He explainedthat at high P, the
combustioninside the vesselis almostcompletedand as a consequence,
when the vent
opens,the turbulent 'micro explosion' in the duct does not have any effect on the
turbulenceinsidethe vesseldueto lack of 'Power'.
.........
..
5.0
0.0
0.00
0.05 0.10 0.15 0.20 0.25
R (barg)
0.30 0.35 0.40 0.45 0.50
Whilst at lowerP,, the mixture in the vesselis still not completelyburnedwhenthe vent
opensandthis eventuallyleadinga very strongturbulencein the vesseland combustion
intensity increaseddue to the back flow propagationwhen significant and fast energy
explosionsin theduct in the end casedueto higher terminal flame speedsin the vessel.
A simple duct flow pressureloss analysiswas used to explain the results, with the
critical event being the entry of the flame into the duct, which causeda dramatic
increasein the overpressuredue to the suddenincreasesin the unburnt gas velocity
aheadof the flame.
From Fig. 5.4 and 5.5, it is obvious that end ignition gavehigher P,,,,, in methane/air
and propane/airat slightly off-stoichiometric mixtures at the presenceof vent cover
compared to the central ignition. It can be elucidated that the different flame
propagationpatternscharacterizingend and centrally ignited n-dxturesare responsible
for different residualamountsof unburnt gasesleft inside the vessel at the time the
flame reachesthe vent and ignites fresh mixtures in the duct. As mentionedin previous
chapter,end ignition gives the maximumdistancebetweenthe flame and the vent. The
action of the vent is to distort the flame shapefrom hemisphericalas the flame develops
preferentially in the direction of the vent, where the unburntgasesare displaced.The
inducedflow throughthe vent duct, aheadof the flame leadsto a significant increasein
flame speedsand expansionratio in the main vessel, known as 'suction effect' in
previouschapter.
154
In the previous StUdý' (Bradley, Hicks, Haq, I, awcs, Sheppard and Woolley. 2003), the
inf'luence of' flarne stretch will also increase the turhulent velocity with a lower VaILIC01'
L.ewis, Le and Markstein no, Ma inside the vessel. Explosion I'lanies, particularly in rich
hydrocarbon IIlIXtUreS are conducive to tile development of' hydrodynamic and thernio
difl'us,ive effects. At sufficiently small Markstein nuniher, and particularly when they
are negative, the thernio diffusive el'l'cct,, (sinaller Lc,ývis no) are no longer stahilizing.
Thernio-diff-usive instabilities arising 1'rom the prel'erential dil'I'Lision of' reactants with
respect to thernial transport can lead to cellular flanic structures that augment I'lanic
propagation through an increase it) name surt'ace area (Bradley, Sheppard, Woolley,
Greenhalgh and Lockett, 2000). The flame sell' acceleration occurs afier a critical
laminar flame propagation distance and the available distance with end ignition is twice
I that with central ignition and hence self' acceleration is more likely. The net effect of'
the vent discharge and self acceleration are Ior the mass burning rate of' tile flame to
rapidly increase due to faster flarnes in end ignition, rather than due to the larger flatile
area of' the spherical flarne with central ignition. Further, f'roin Fig. 5.4 and 5.5, it is
worth noting that both figures Suggest that the enhanced cornbustion phase in the end
ignition is of very short duration while it stands for a quite longer tirne in centrally
ignited case. The latter case confirmed the hypothesis stated by Ferrara et al (Ferrara,
Benedetto, Salzano and Russo, 2006). They said that in case ot'central ignition, there is
an indication of' higher quantity of residual unburned mixture in the vessel whereas
almost complete combustion occurred in the case c4end ignition that eventually leads to
hig,her peak pressure,
1.0
0- 0.5
0.0 ý
0.04 0.06 0.08 0.1 0.12 0.14 0.
-0.5 T me(s)
However, this trend is not rnanifested in a very reactive mixtures i.e. ethylene/air and
hydrogen/air even at the lower equivalence ratio. It shows that central ignition gave
higher peak pressurecompared to end ignition. The same observation was illustrated at
n
high burning velocity mixtures for simply vented explosion previously. It should be
noted that in the present study, the vent coefficient I'm the rig configuration was 16.4
and due to the fast burning velocity ofhydrogen/air and ethylene/air, it creates higherjet
velocities at the vent opening and this sudden onset of the venting leading to the rapid
turbulent at the vena contraction region. In the case of sonic condition (P,,,,,,> 900
inbar), it causes the vent pipe to choke and theoretically, the vent flow is a linear
function of the internal vessel pressure. Since there is higher amount of unburnt gases
left in the vessel for the central ignition in respect of end ignition, it Increases the
cornbustion ratc dUe to the subsequentturhulisation (by the physical back-flow into tile
vessel) and by the interaction ol'the shock pressure waves with the flame 1'rame.Indeed,
the entrainment of' Freshgas pockets operated Lit high velocity hot gases is well known
to originate a violent ignition (Bradley, Emerson and GLI, 2003, Lee and Guirao, 1982)
and under certain conditions, leading to detonation its observed in Fig. 5.6 and 5.7.
150
8.0
6.0
4.0
2.0
0.0
0.04 0.05 0.06 0.07 0.08 0.09 0.1
-2.0
Tme(s)
central ignition showed much largocr influence of' P, on P,,,,,,, contradicting with the
the combustion time become shorter and less time is available for gases inside the
vessel to flow out before the combustion is completed, thereby reducing the
effectiveness of venting (Kumar, Dewit and Greig, 1989). Tile combustion times lor
central ignition are roughly hall' of' the end ignition and this explained the obtained
results.
8.0 C2Wai r
5.1.4 The acceleration of flame towards the vent duct and associated vent
velocities
respectively.For ethylene/airat (D = 0.8, the Lewis no and Markstein no are 1.02 and
4.34 respectively(Searbyand Quinard, 1990).However,there is no availabledata for
leanH2concentrationon Lewis no but -0.82 is Marksteinno at (D= 0.54.Markstein no,
Ma is a physicochernicalparameterthat expressesthe responseof flame to stretching
(Bradley,Lau and Lawes, 1992) while Lewis no is the ratio of conductiveto diffusive
fluxes where it describesthe relative rates of heat and mass transfer. Preferential
thermal diffusion is indicated by Le >I and Le <I is a preferential massdiffusion
(Clarke,2002).Lower value of Markstein no with Le <I indicatesthe propensityof the
flame becomeunstableandprone to cellular instability and self-fragmentation(Bradley,
Hicks,Haq, Lawes,Sheppardand Woolley, 2003).The net effect of the self-acceleration
and vent discharge'suction' effect is for the massburning rate of the flame to increase
due to fasterflames.This showsgreaterfor end ignition as the flame travel distanceto
158
the vent is twice that for centralignition. The increasein flame surfaceareadue to the
'suction' effectand longer flame travel distanceareto compensatefor the lower initial
massburning rate of end ignition with its 50 % reducedflame surfacearea and heat
lossesto the endflange.The net effect from this behaviouris higher overpressurewith
end ignition for methane/airand propane/airmixtures but not in the caseof reactive
mixturesi.e. ethyleneandhydrogenas shownin Table5.1.
In Table 5.1, simple momentum conservation and mass continuity were used to
calculate the downstream duct pipe velocity, Sg prior to the flame entry to the pipe.
Assuming ideal incompressible flow within a duct with no wall friction, mass continuity
gives,
where m is the massflow rate, A is the duct cross-sectionalarea and U is the flow
velocity. SubscriptsI and 2 refer to the duct measurementat two different point
positions.The gasvelocity is thengiven by,
(5.8)
pA
As the area of the circular pipe is a function of the diameter, we can reduce the
calculationfurtherby,
(A, )
so Sg,,
=
e. el A2
(5.9)
However,Eq 5.9 assumedthat the flame behavesasa piston, which is not the reality as
the flame is 'pear shapedand not flat' (Ellis, 1928).When comparedto the prediction
using the loss pressurecalculationusing Eq. 5.5 previously,the resultedgas velocities
from Eq. 5.9 are very high and greaterthan sonic (in somecasesfor propane/airand
Mach
ýRT
Co=ý
ME
Table 5.1 also indicatesthe little changein the peak final flame speedwith P, but it
should be rememberedthat the massburning rate will increasewith P, due to the flame
density term. For the highestvent burst pressure,this gives closeto a 30 % increasein
the massburningrate whenthe vent burstsaccordingto the Kci, t values.This doesnot
e
change the velocity of sound in the vent duct and hence the Mach number is not
affected.A methodof taking into accountthe bulk shapeof the flame is to look at the
rate of changeof P, with time just prior to the flame enteringthe duct. This is shownin
Table 5.1 with Kr,, This shows low values for central ignition at low P,, which
e,,t.
increasewith P, to be nearly 16timesthat for an openvent for methane/airand 10 times
for propane/air,For ethylene/airand hydrogen/air,it is about2 times and 1.4 times that
for an open vent respectively.With end ignition the effect is much lower apart from
propaneat the highest P, This may be due to the greaterinfluence of the increasing
160
The faster inass burn ratc approaching the vent as P, ilICI-CiSeSCMISeSSOI1iCFIOWIII the
vent in(] hence choked flow Vý'experiencLA This prevents there heing any oufflow 1'rom
the duct until the pressaire has HSCIIin the VCSSCItO drive the bLIl-Iltgases out. There then
l'ollows a period of' inixed hurnt gas and gas, venlino with micro explosions and
L- 1ý
detonations (propane/air real- ignition at P\ = 427 inhar and most cases in ethylene/air
and hydrogen/air) in the vent duct. This has heen detailed in section 4.2.
Figure 5.9 - 5.10 show the present flame speed meaSUITIllellISwere considerably higher
than the values expected on the hasis of' the development 01' Cellular flames, 1,01,
all
Studied fuel/mixtures in this present work. As been discussed in previous chapter to
tile Occurrenceof self-acceleration of the flame, tile sarne method was applied
COIII*il-l-n
based on the K(;,,,,,,adopted frorn NFPA 68. This enables a prediction from the NITA
68 data to be made of' the self-acceleration effect on flanne speed and this i,, shown in
Fig. 5.9 and 5.10.
25 CH4(NFPA)
C3H8(NFPA)
20
CH4, central, Pv--O
mb
end, Pv= 0 mb
15 -)K-CH4,
C31-18,end,Pv =0
mb
10 C31-18,central, Pv
0 mb
5
0
0.12 0.32 0.52 0.72 0.92 1.12
Distance from the spark (m)
FIgUre 5.9 Flanic spced a function of' tile I'lanic distance from (lie
spark I'm P, l'or
nictliane/air and propane/mr. (1)= 1.0 for hoth gas/mixtures.
At the K of 16.4,there was a high velocity createdin the unburnt gasflow through the
vent and this induces movement in the unburnt gas upstream of the vent, leading to
flame accelerationtowardsthe vent.The effect is to distort the flame into a pear shape,
with the apex of the pear passingthrough the vent first. This is
effect similar for
methaneand propane and the difference in the peak flame speedswere not large,
whereasthey should have been large if flame instabilities were the cause of the
acceleration.
speed is laminar for about 6 m/s at 0.36 m from the spark at end ignition for
C3118 Central 0.209 1.8 89.5 5.3 20 5.5 21.9 208 237 189
C3Hs End 0.209 2.8 201.2 5.5 35 12.4 49.6 472 431 345
CH4 Central 0.209 2.5 141.5 5.5 22 6.7 26.8 255 257 205
CH4 End 0.209 2.2 142.7 4.5 24 4.7 18.6 177 135 108
C3148 Central 0,424 2.9 117.0 5.5 30 5.2 20.8 198 216 173
4.0
*6.5
C3Hs End 0,424 296.5 5.5 55 20.2 81.0 771 666 533
CH4 Central 0.424 2,7 109.4 5.5 26 7.9 31.7 302 313 250
CH4 End 0.424 2.5 19011 5.5 34 5.5 21.8 208 290 232
164
62)2 Tfl/S
9.52
Sg
C2l14
Central 0 *7.7 469 2.7 17.1 3.4 13.6 130 181 187
CA
End 0 *5.0 323 7.9 16.8 3.4 13.5 128 272 331
H2 Central 0 *3.4 190 11.6 25,2 5.0 20.2 192 227 330
H2 End 0 *3.5 142 7.9 21.4 4.3 17.1 163 227 208
C2144 Central 0.098 *7.3 502 0.3 11.7 2.3 9.4 89 211 272
CA End 0.098 *6.1 364 8.2 17.2 3.4 13.7 131 156 296
H2 Central 0.098 *9.9 323 6.4 25.9 5.2 20.7 197 453 306
H2 End 0.098 *4.2 232 4.4 28.5 5.7 22.8 217 243 404
C2H4 Central 0.178 *9.7 580 1.8 15.2 3.0 12.1 116 234 385
CA End 0.178 *6.2 416 7.3 19.8 4.0 15.8 150 190 404
H2 Central 0.178 *8.8 354 12.4 27.4 5.5 21.9 209 509 410
H2 End 0.178 2.8 133 2.0 28.3 5.7 22.7 216 298 435
C2114 Central 0.209 *10.7 650 3.8 15.7 3.1 12.5 119 247 408
C2114 End 0.209 *6.9 480 7.4 22.8 4.6 18.3 174 235 459
H2 Central 0.209 *9.1 312 13.6 31.5 6.3 25.2 240 544 420
H2 End 0,209 3.2 181 10.5 30.5 -6.1 24.4 232 338 439
C2H4 Central 0.424 *7.4 461 5.1 10.0 2.0 8.0 76 252 467
CA 24,2 499
End 0.424 *6.9 519 7,7 4.8 19.4 184 243
H2 Central 0,424 *8.9 331 16.5 17.3 3.5 13.8 132 319 453
H2 End 0,424 2.0 104 5.7 37.2 7.4 29,8 284 341 499
50 C
--0- 2H4, end, Pv -0mb
45 --[4-C2H4, centre, Pv-Omb
40
35
30
25
20
a: 15
10
5
U
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Distancefrom the sr)ark Im)
Figure S.10 Flame speed as a function of distance from the spark for P, =0 for
ethylene/air ((D = 0.8) and hydrogen/air ((D = 0.54). Noted that only hydrogen line for
NFPA 68 was drawn as there is no available data for ethylene in the guidance.
To justify the influence of vent cover beenplaced during the presentwork, Fig. 5.11
and S.12 shownthe upstreamflame speedresultsfrom the experimentaltests.The mass
of the burnt massat the time of the ventburstingmay be calculatedfrom the assumption
that the pressurerise at the time of the vent burstingoccurswhenthe burnt massof gas
is P,/Piniti,,
1, wherePiiti,,,is the initial of
pressure I bar atmospheric If
pressure. the flame
pressure
(mbar)
Methane/air 0.16 m 0.20m 0.21 m 0.27 m
(Central
ignition)
Methane/air 0.21 m 0.25 m 0.27 m 0.34 m
(End ignition)
Propane/air 0.16 m 0.20 m 0.21 m 0.26 m
(Central
ignition)
Propane/air 0.20 m 0.24 m 0.26 m 0.33 m
(End ignition)
Hydrogen/air 0.15m 0.18 m 0.19 m 0.24 m
(Central
ignition)
Hydrogen/air 0.19 m 0.23 m 0.24 m 0.30 m
(End ignition)
Ethylene/air 0.14 m 0.17 m 0.18 m 0.23 m
(Central
ignition)
Ethylene/air 0.18 m 0.22 m 0.23 m 0.29
-m
(End ignition)
12 Tan 1.8
10 1.5
8 1.2
6 0.9
4 0.6 E
cl-
2 0.3
0 0
-2 -0.3
The reason for the greater flame movement than expected with end ignition was (ILICto
the greater expansion of the flame on the centrelinc compared with the edges, which
were subject to more cooling. This is a well known effect Coi large L/D vessel and for
very large L/D vessels results in an elongated U flame, with very fast centreline flarne
speeds(Phylaktou and Andrews, 1991). If' the distances in Table 5.2 are compared with
the flame speed results in Fig. 5.15 it can be seen that is soon as the vent burst there
was in increase in the flame speed. This occurs at a shomer flame radius for central
ignition, as predicted in Table 5.2. Also Fig. 5.9 and 5.10 shows that without a vent
cover, flarne acceleration above that (We to cellular flames, started at distances similar
to those in Table 5.2 for P, = 98mbar.
I(
25 1.5
2-0 1.2
0.9
1.0 0.6
0.5 0.3
0.0 0.0
-0.5 -0.3
0 0.05 0.1 0.15
Time(s)
Figure 5.12 Propane/airat (1)= 1.0 for end ignition. P, = 179 n1bar
7.5 tn 1.9
1.7
6.5
1.5-
5.5
1.31
4.5
3.5 0.9
2-5 0.7
0.5
1.5
0.3
0.5 0.1
-0.5 -0.1
0 0.04 0.08 0.12 0.16
Time(s)
3.5 2.1
tin totil P2
3.0 1.8
P,
2.5 1.5 E
P.)
-3
2.0 Pduct 1.2
1.5 P 5-6
.2
U.9- (-J)-,,,
0 CZ
1 Flame position P3-5 0.6 cc/)L
-0
0.5 0.3
U-
0.0 X- 0.0
1
-0.5 I! 1 -0.3
0.05 0.1 0.15
Time(s)
FigUre 5.14 Methane/air for (1)= 1.0 at central ignition. P, = 424 illbar
As mentioned earlier, the I'lame speed after the, vent burst increases rapidly for both
15
10
1,
5
---W- 5- -LL
00
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.6 1.0 1.2 1.4
Disfatý trcmftý spark (m) Di starce from ft sparx (m)
IT (r 0.0
x-5.0
10
0
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0
0.0
x
Figuire 5.15 Influence of Pv on flanlC SIVed LIPSti-CIlll01'thC VCIIJf0j' jjjetjjjjje/, jlj' ((I)
1.0) and propane/ali-((I)z--1.0) with end and central ignition.
170
60 1r-, Cortrelgrition
C2ýVaJr,
04-4, Erd Ignfion 7-z: Pv.
1!ýýPv
BBfft3w
16 :
so I- FV. 175ntw
ý Pv 178ntwl
FV . 209 rrbW 14 Pv: 209ntar
40 PV - 424 mber 424ftwi
12 _'',
'ýPv
lo
30
,
rL 20-, 6-
4,
10 A
2-
oo------ 0
O.D al a2 063 a4 0.5 0.6 0.7 0 0.05 al Ql 5 0.2 0.25 0.3 0.35 0.4
D"=fmmlhawark(rrO Cistarcefromthespark(ni
so I-po. ombff Wair.Erdloi6an
70 1-
Pv - 98 rvtw1 40 - 7-_pv: - WJair, Cove igritan
Orrim
ý-Pv- 178 =
60 -PV eafft)W
-N-208: 35
_P4, . ft. i-Pv . 178rrtw
50
("IA 68) 30 -Pv-2DDnter
ý-PV-424ntw!
,ii-W 25 11, Fe W 613
30 20
20
W 10
0 5
0 111 062 03 114 M5 Ole 0.7 0.8 0
DIOArtefromto spark(ni 0.2 03 0.4 0.5 O's
Dstw, wfmm*mopwk(rrj
The results of the maximum flame speedsin the vessel are shown as a function of
equivalenceratio in Fig. 5.17 - 5.20 for methane,propane,ethylene and hydrogen/air
with end andcentralignition. Thesefiguresreport the averageflame speedsmeasuredin
TV2andTV4in Fig,2.3).Theseflamespeeds
thesecondhalf of themainvessel(between
are considerablygreaterthan that for the laminar sphericalflame which 3m/s for
methane,3.6m/sfor propane,5.5m/sfor ethyleneand19.7m/sfor hydrogen/air.
entry into the vent duct, was 19 m/s for (D= 0.8 with no vent cover and increasedto 50
m./s for P, = 424 mbar. The laminarburning velocity at (D= 0.8 is approximately0.35
m/s and the adiabaticsphericalflame speedapproximately3 m/s. The measuredname
speedsare muchhigher thanthis andthe increaseis consideredto be due to flame self-
acceleration resulted from the formation of cellular flames
(Kasmani, Willacy, Phylaktou and Andrews, 2007). The maximum reactivity of
propane/airmixturesis at (D= 1.1,asthe peakflame temperatureoccurshere.
CH4/adr
Go
Centr.1 igNtion
50ýý 0 Pv=O
, -m- Pv = 98 nta
40- & Pv = 178 ntEu
x Pv = 209 rnbat
30 w Pv- 424 nba
1
201
gý
lz 10
0
0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4
EqLiWerr,
e atio
40-, CH4/ajr
36 Endlgribon
32
24, i
20-1
161 0 Pv- 0
12 " Pv- 9B n
1& Pv - 178
8
), Pv- Mg
I
)K Pv. 424
0 L-- -.- - -
0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4
Eclivalerre rclio
The resultsin Fig. 5.18 showthat the maximumflame speedwith no vent cover was 52
m/s, considerablyhigher than the 3.6 m/s for a sphericaladiabaticlarninar flame speed.
This increasedto 120 m/s with the 424 mbar vent burst pressure.However, for richer
173
mixtureswith (D= 1.35the flame speedincreasedto 234 m/s for the highestvent burst
pressureof 424 mbar. For no vent burstpressurethis mixture had a lower flame speed
than for (D = 1.1, as expectedon the basis of laminar flame speeds.For very rich
mixtures of (D = 1.65, which is closeto the rich flammability limit, the flame speeds
reducesharplyto valuescloseto thoseat (D= 0.8. Theseflamesspeedswere20 - 40 rn/s
andarevery fast for sucha nearlimit mixtures.
240 -j Q3Ft/air
ýýigribon
200,
-70-LPV-0-
-a Pv=98mbar
160 Pv-178mbar
x Pv=209mbar,
120 w Pv-424mbar
80
40-1
2W -IC, 3F6/ajr
Certreiarition
Pv- 0
160
Pv - 98 ntar
ý6 Pv- 178 rrbar!
120, ý
x Pv- 209 nbwý
w Pv - 424 atmý
80
40-,
m r: f-fEEE
0
0.6 0.8 1 1.2 1.4 1.6 1.8
Ec;ii\derr-a rdio
vent cover,the maximumflame speedwas 17 m/s and this increasedto 24 m/s with the
424 mbar vent burstpressureas picturedin Fig. 5.19.There is not that much significant
increasein flame speedin various P, for end and central ignition as the occurrenceof
flame cellularity is not great at this lean concentration stage but the turbulence
generatedby the backflow and flame front interactionwith pressurewave substantially
increasethe P.. as shownin Fig. 5.23 later.
175
30 Q2H14/ü r
Erd igribon
0 Pv -0 nter
1 la Pv - 9B ntar
20
Pv. 178 rrt)arý,
>
x Pv - 209 ntmr'
15 i
.9 NE Pv = 424 ntari
10
0,
0.55 0.6 0.65 0.7 0.75 0.8 0.85
EqLivAerr,e rdio
CýH4/edr
25-1
Certreigriton
20 s Pv= 0 ftar
Pvv
- 98 ntw
15-, 6 Pv= 178nteri
N= 209 ntarý
w Pv - 424 rrbTl
10-1
0, IrII11
45 hVair,
4o Erd igrffion
351
Mý E3 Pv- es n-tar ý
30-
16 Pv= 178 rrlmrýi
25 4m Pv= 209 rnbar
A 20 Pv= 424
15
10-
[Z e
5-1'
01A, I-
45
F61aj
r
40-1Certreigrit on
35 Pv= 0 mtxr
30 Pv - 98 nbar
Pv: 178 ntal
25 I&
Pv 209 nt w
20 Pv= 424
15
10-
i
5ýI-
The peakmeasuredflame speedsprior to the flame enteringthe duct pipe are shown in
Fig. 5.9,5,10,5.15-5.16,5.17-5.20as well as in Table 5.1. The associatedUnburntgas
velocity and the turbulent burning velocity upstreamof the vent, basedon this flame
speedand adiabaticburnt gas expansionare also shown in Table 5,1. The turbulence
burning velocitiesmeasuredwith central ignition give a ST/SLburning velocity ratio of
177
velocity in the ductjust prior to the flame entry to be predictedand theseare also shown
in Table 5.1. This calculationassumesthat the flame behavesas a piston, which is not
the reality as the flame is 'pear 'shaped and not flat. However, the resultant gas
velocitiesare very high and greaterthan sonic in a few cases(incompressibleflow was
assumedin the momentumequation).This indicatesthat nearsonic flow is likely in the
vent and the absolutepressureratio, P2/P6,acrossthe vent duct in Fig. S.II to 5.14
confirmedthis observation.It is importantto realisethat this duct flow is driven by the
fast flamesupstreamof the vent,shownin Fig. 5.9,5,10,5.15 -5.16 and 5.17-5.20.
The flame speedsin the vent duct were measuredand tabulatedin Table 5.1. These
were very fast flames even for an open vent. They were consistentlyhigher for end
ignition due to the fasterflames and the predictedfaster unbumt gas velocities in the
duct. The flame speedsare reasonablyclose to the predictedunburnt gas velocities
computed from the upstreamflame speeds.There are generally higher due to the
influenceof theturbulentburning velocity. However,it is clearthat the dominantcause
of the fast flames in the vent duct is the fast upstreamflames detailed above. The
unburnt gas velocity aheadof the flame in the duct has beenestimatedas 80 % of the
flame speedin Table5.1. For an adiabaticexplosionthe gasvelocity aheadof the flame
is (E-I)/E timesthe flame speedsandthis is 87 % of the flame speedfor stoichiometric
hydrocarbons.The lower value has been used to account for duct wall heat losses
(Phylaktou,FoleyandAndrews,1993).
In the caseof equivalenceeffect, the high flame speedupstreamof the vent induces
higher unburnt gas velocities ahead of the flame. These high vent pipe unburnt
velocitiesresultsin a very high pipe turbulencelevel and consequentlyvery high flame
speedin the duct pipe as shown of the worst case of propane/airin Fig. 5.21 and
hydrogen/airin Fig. 5.22. This showsa very large increasein flame speedinside the
vent duct with P, from 200 m/s to 700 m/s as the P, increasesfrom zero to 424 mbar in
the caseof propane/airat endignition for (D= 1.0.In the caseof hydrogen/air,the name
speedinside the vent duct was 300 m/s for P, =0 and increasessharply to 550 M/s for
Pv = 424 mbar for (D = 0.54. The causeof this is the large increase in flame speed
178
upstream of the vent in Fig. 5.20 from 23 to 40 m/s corresponding to the increase in Pv-
This high flame speed inside the duct pipe at lean concentration in hydrogen is
responsible to induce high turbulent flow leading to an intense mixing of cold and hot
wo , CArair
1 End ignition
7001
s Pv-o
Iz J
600 la Pv- 98 mbar
i tj Pv-. 178 barý
500
Pv=209mbari
m
4W - PV- 424 mbar
300
200 J'
100
0
0.6 0.8 1 1.2 1.4 1.6 1.8
Eclim3lerce rdo
Figure 5.21 Flame speedin the duct pipe for propane/airas a function of equivalence
ratio at endignition,
179
600 ] Wadr
! Certre lgribon
500-,
. A Flv ntw
=0
Pv- 9B ntar
4W - Pv- 178 rrbarý
x Pvý 209 rd",
300-, Pv- 424 n-bari
,
.9
200-1,
Ir 100 -,
0
0.25 0.3 0.35 0.4 0.45 o.5 0.55 0.6
EqLi\rcderr,
e rc-dio
Figure 5.22 Flamespeedin the duct pipe for hydrogen/airas a function of equivalence
ratio at endignition.
vent was alwaysthe lowest P.. for all equivalenceratios. Theselower overpressures
for 98 mbar and 178 mbar vent burst pressurescomparedwith open vents for the
propaneexplosionsare difficult to explain. The upstreamflame speedsresult in Table
5.1 doesshowinga higher flame speedfor the openvent for (D= 1.0 and at (D= 1.37 as
shown in Fig. 5.18. Propaneis propensityfor cellularity at rich concentration(Tseng,
Ismail and Faeth, 1993) and this would explain the higher P. in comparison of
"
methane/air.For propane,the onsetof cellular flamesoccursover a shorterdistanceand
the self accelerationof the flame is greater,which is a function of the distancefrom the
spark.The significant effect from this is to acceleratethe flame towardsthe vent more
quickly than for methane/airexplosions,as shownin Fig. 5.15. The effect is to stretch
the propaneflame so that it is more elongatedthan a methane/airflame. The action of
180
the staticburstpressureof the vent is, for low burstpressures,to delay the acceleration
of the flameto the vent. The flame doesnot know that the vesselis venteduntil the vent
burstsandwaspropagatingasin a closedvessel.For low ventburstpressurethis givesa
smallerflame areawhen the vent burstsandit is the effect of the vent burst pressureon
the flame bulk shapeand not on its centreline velocity that causesthe pressureto be
reduced. It is consideredthat this effect hasa greaterrelative influence on the propane
flamesdueto their greateraccelerationto the vent andthis is why propaneandmethane
flamesbehaveddifferently for low vent openingpressures.A higher initial flame speed
would indicatea higher massbum rate and hencea higher massflow rate in the vent
pipe and this would createmore turbulenceand hencehigher overpressuresdue to the
fasterflame speedsin the vent pipe.
Unfortunately,this was the only result that gavea higher flame speedprior to the vent
and hencethis mechanismcannotexplain all the results.Further,it is very significant
that only the propaneexplosionsdemonstrated
this effect. This is also indicatedthat the
quickeronsetof cellular flamesandthe greaterself accelerationof propaneflames may
be the causeof this difference.It is conjecturedthat P.,,,, was generatedby the rapid
turbulentcombustionof the unburntgasesleft upstreamof the vent in the main vessel.
The flame preferentially acceleratesin the direction of the vent once this is open,
leaving a large proportion of the unburnt mixture in the ventedexplosion vessel.The
gasesin the outerpart of the vesselare unbumtwhenthe flame entersthe vent. A high
flow back pressureis generatedby the high velocity unbumt g4sflow in the vent and
high flame speedsin the vent pipe, asshownin Fig. 5.21. This forcesthe
the subsequent
upstreamflow to reverse and ignite the outer unbumt gas flow in a fast burning
turbulent flame. This then acceleratesthe burnt gas flow out the vent pipe, which
increasingthe pressurerise.
181
CH,/air
CH4/air Co*e Igrition
Eridigriban 3.0 j-. *-Pv-ontEr
2151!
-6- Pv. 0 1
!
--o- FV. 98 fftw
- Pv. 96ntx
ao ý 2.5 -a- Pv .M rrbar,
Pv. 178 ntEr,
-G-
-. m- Pv - 209 rrberý
-, m- Pv- 209 ntrr
IT
20-,,
i --w-Pv»424 or,
x-
2.0-1
ts
or
a5i
;
ao 0.0
116 0.7 018 a9 1.1 1.2 1.3 1.4 0.6 (17 0.8 OL9 1.0 1.1 1.2 1.3 IA
Eco%derre refio Ectkderce rgdio
7.0 - C"air
End lgribon SO QkWir
4.5 Certrelgrit
-ä- Pv: 9Bn F'-Oýl 6
5.0- 4.0.
-a.- Pv 178 ý--A-PV-98
as
4.0 - -'»- - MB Pv- 171
*-PV. 4m ao- PV-2(X
Z5 - pv.
10-
1.5
to
US
0.0 0.0"
0.6 018 1 12 1.4 1.6 1.8 as 0.8 1.2 1.4 1.6 1.8
EcFJ%egerce
reMo EcpjKa6, ce rido
4.0
20
ao
(10 --I-- ----
OL85
a5o 0,55 OAD 0.65 0.70 OL75 a8o
0.50 0.55 (160 a65 VO CL75 (180 US
Ecdv,4erce rabo
EqLi\oWerre
mk
wat
45 kvalr 120 C"eigtifion
4.0 Erdignton
IýFV-Offbar Pv-Ontff
aI51
ýPV. 98rfbar Fv - OBrrtxr
3.0 1-*- Pv- 178Ma' - Pv. 178nbar,
Pv- 209ntivi 8.0
iý- Pw: 209ntwý
'Lý 2.5- FV- 424ft&
Pv 424nbý
201,
4.0
1.01
20
ao
a3D (135 0.40 a46 aso OM a60 0.0 i 10,
--
E4mmcerado (130 0.35 0.40 0.45 0.50 0.55 0.60
EcovJemsrado
Figure 5.27 showed the Predwith duct as a function of the reduced pressure obtained in
the same tests without the presence of duct. On the same figure, all experimental data
were compared to the correlation given in NFPA 68 for duct pipe length in between 0 to
3 m. As expected,the correlation given by NFPA 68 is under predicted the experimental
data results, suggestingthat there is more complex interaction between the vessels with
the presenceof duct with respect to the ductless vessel explosion i. e. simply vented. It is
also interesting to note that even though the same reduced pressure obtained in vented
explosion without the duct pipe attached, there was a variance Of Predin the presence of
duct pipe. This situation implied the correlation associated P,,d with the duct pipe and
without the duct pipe should not be taken as simple as given in NFPA 68 but other
parameters should be taken into account to accommodate the increase Of Predwhen the
duct pipe is attachedwith respect to simply vented explosion. As previously discussed,
ignition position, vessel volume, self-acceleration, induced turbulent, pressure loss as
well as secondaryexplosion should be included in the correlation in order to fit with the
results obtained experimentally. All published experimental data are listed in Appendix
A (Table A. 2).
184
10.0 -1
9*01 DeGoW and ChEtadit
8.0 t) Kbkovet d
c
7.0 ,& KtCarn et Ed
6.0- x Wharddr
x Prcpardar
5.0-,
o EthylenfMr
"": 4.0 0
+
3.0 NFPA68
2.0 - Badknoctt'se4jWicn,
1.0 ý
0.0 )r-I
0.0 20 4.0 6.0 8.0 10.0
ckict(bar)
PiedWthoLA
* MoCarnal al ( MEtheriEVair)
0.0
0 50 100 150 2DO
UD
L
The high overpressure due to the addition of a vent duct, to a 0.2 m3 closed vessel with
an LJD 2 K
of and of 16.4, occurs after the flame has exited the vent duct, but is not due
to an external explosion. The main effect of the ignition position is to give a greater
distance from the spark to the vent and greater flame acceleration. End ignition gave
higher Pmaxfor low burning velocities mixtures but central ignition is the worst case
scenario in Pnaxfor hydrogen/air and ethylene/air. It was found that substantial amount
of unburnt gases left inside the vessel after the vent burst is the leading factor in
increase of Prnaxfor high burning velocity mixtures for centrally ignited. The associate
gas velocities ahead of the flame create high unburnt gas flows conditions at entry to the
As P, increases, the distance of normal spherical flame propagation increases and there
is a further reduction in acceleration distance. This initially reduces the overpressure at
low P, The effect of the vent burst pressure is to increase the flow velocity in the duct
when the vent burst, as the flame has had more time to grow upstream of the vent and
this gives a higher vent duct flow velocity once the vent cover bursts. The effect of the
vent burst pressure is complex and non-linear and is not represented by the linear effect
in the correlation of Bartknecht (Bartknecht, 1993) used in NFPA 68 (NFPA68,2002).
186
6.0 Introduction
The influence of duct attachedto the main vesselon P,,,,,,has been investigatedin
earlier chapter.The investigation confirms the findings of the previous works and
highlights someconclusivefactors that lead to the severity of vented gas explosions
with the duct pipe attached. From the work, major influences on the behaviour of
duct
ventedexplosionswith attached pipe which leadsto a severePma,,
are;
validationof this procedureas this topic has been sparselyresearched.In the presence
of a vent duct, an increaseof venting area and duct diameterhas been found to not
alwaysresult in a decreasein the peak over pressure(Ponizy and Leyer, 1999a).This
returnsthe focus to the more realistic scenarioof whetherthe similar behaviourwould
188
be observedwhenthe bigger duct pipe is attachedin comparisonto the duct areais the
same as the vent.
The main reasonfor the increasein the overpressurewhenlong vent ductsare attached
to vents is due to the phasein the explosionwhen the flame is in the vent pipe with
unburnt gas mixture aheadof it. The expansionof the burnt gases in the vent pipe
the
greatly accelerates unburnt gas flow and this increases
the vent pipe friction, inlet
In this presentwork, only uncoveredvent casewas carried out with four different gases
(methane,propane,ethyleneandhydrogen)at different equivalenceratios. All mixtures
were ignited at the end wall and/orat the centreof the vesselcentrelineby an electrical
sparkwhich gives16 J energiesfor the gasexplosiontests.
Figure 6.1 shows, for the most reactive methane/air mixture, the pressure-time profile
inside the vessel at P, for the short 0.162m diameter pipe, which is effectively a free
discharge, and for the two Im long vent ducts of 0.162 and 0.315m diameter
respectively. The results show that the larger pipe diameter (duct area / vent area, AýA,
= 3.78) has little effect on the overpressure compared to that with the vent duct the same
diameter as that of the vent. This was not the expected result and did not agree with the
results of Nagy (2002, Nagy and Verakis) or Hey (Hey, 1991). However, this situation
is not been observed in hydrogen/air explosion as illustrated in Fig. 6.2 where vented
explosion with bigger pipe attached (Dp = 0.315m) shown the pressure peak is about 2
times higher than vented explosion with duct diameter of 0.162m. Interestingly, the
peak pressure for vented explosion with 0.162 m diameter pipe attached gave the same
peak pressure as the baseline case Le, simply vented. Interestingly, this trend is not
observed for end ignition where maximum pressure is lower for baseline case in respect
with duct vented explosion as shown in Fig. 6.3. Again, for centre ignition, the pressure
rise inside the vessel is directly influenced by the mixture reactivity left inside the
vessel after the flame exited the vessel and the presence of the duct gave little effect to
the final value of pressurein this case. This observation contradicts with the other works
190
FVadr O=0.54
8.01
I Certre igrition
7.0 1,
... Smplywnted
6.0
D_Jct"e nted DP_O. lr "2
5.0 Duct wnted, Dp 0.315 m
4.0
3.01'
ý.o
1.0
0.0
1) 0.02 0.04 0.06 0.08 0.1 0.12 0.14
-1.0
i Time (s)
-4u,
3.0 Simplyýented
2.0 ý
-Duct \.ented, Dp - 0.315m
1.5
La)
1.0
0.5
0.0
ý. 03 0.05 0.07 0.09 0.11
-0.
T ni--(s)
From Fig. 6.2 and 6.3, it can be deducedthat ignition position gives important inflUenCe
L-
on the final pressure inside the vessel. As discussed in previous chapters, the
combustion inside the vessel is almost completed when the flame exited into the duct
pipe, leaving only small amount 01'LinbUrntgasesinside the vessel which results M high
flame speedsand hence high unburnt burning velocity prior to duct entrance. Tile final
pressure inside the vessel is directly due to the high turbulent reversal flow from the
duct pipe to the vessel from a violent burn-up or secondary explosion and the flame
front interaction with the pressure/shock waves. In other words, there is a mutual
interaction between the vessel and duct during the explosion development marking to
the final pressure inside the vessel. For central ignition, the interaction between vessel
and the duct seemsnot playing an important role in order to determine the final pressure
hUt dUeto the substantial amount ofunbUrnt gasesleft inside the vessel.
TO FLINherInvestigate the overall trend of' vented gas explosion with bigger diameter
pipe attached, as a FunctiOll 01' e(lLIIVaIeIICCratio I*or all StUdied gases were
presented in Fig 6.4 - 6.7. For mellialle, at Q) = 0.84, the larger vent duct had only a
slightly larger overpressure than J'Orthe near 1'reevent condition and at (1)= 0.68, the
overpressures were the same. However, for tile 0.162 in Vent (ILICtthe OVerpreSS111-CS
were always much higher than for tile freC Vent it 111eLjLIiVIleIICCratios but slightly
lower in respect with larger vent p1peat (1)= 1.05 for end ignition. The saine trend is
192
shown in propane/airwhere the rich mixture concentrationgave high P., for larger
vent duct. However,in ethylene/airexplosion,it seemsfor the attached0.162 m vent
were higher than for the duct with bigger diameterattachedas
duct, the overpressures
shown in Fig. 6.6. From this data, it can be deducedthat it was only for the most
reactive mixtures that the larger vent duct did not solve the problem of the large
increasein the venting overpressurewhen the duct was the samediameteras the vent
andthereasonsfor this wereinvestigatedfurther.
CH4/adr
1.6 D-ctvertecL10 -
(1315rnendIgrition
D-xi vwtecL Dp
1.4 0.315MpentrEdotionon
Figure6.4Pmax
of methane/airasa function of equivalenceratio
Sn'PiyveMed,endlgrülon
a5 -9 C.3hWair
EI Srrptyverged,certre igrition
0.5
0.0
The overpressures
in Fig. 6.4 and 6.5 for the mostreactivemixture (P,,,,,,= 1.49bar for
(D = 1.08for methane/airat end ignition and P,,,,,,= 3.26 bar for (D = 1.375at centre
ignition for propane/air)cannot be predicted from recommendedcorrelations by
Bartknecht(Bartknecht,1993)as written in Eq. 5.3 and5.4. As mentionedin Chapter3
for simply ventedexplosion,Bartknecht(Bartknecht,1993)gavethe correlationfor 100
mbarstaticventburstpressureas;
The free vent overpressure is predicted by Eq. 6.1 to be 5.45 bar for methane/air
comparedwith the measuredvalue of 0.35 bar. It is 7.46 bar for propane/airfrom Eq.
6.1 and 1.62 bar experimentally.It is clear that Eq. 6.1 cannot be applied to smaller
vesselvolumes.Examinationof the venting dataof Bartknecht(2002) showsthat all of
his vesselvolumes had lower overpressuresthan for Eq. 6.1 and that this was the
correlationfor his datafor a 10 m3vessel.If the measuredfree vent overpressureof 0.35
bar for methaneis takenthen the effect of a vent duct < 3m long is predictedby Eq. 5.1
and 5.3 to be an increasein the overpressureto 1.23bar and only 0.5 bar respectively,
well below that measuredexperimentally.However,if the correlationin Eq. 5.4 for duct
length> 3m is usedthen the predictedincreasein the measuredoverpressureis 1.44bar
for methaneand3.18 bar for propane,which is closerto the measuredresultsin Fig. 6.4
and 6.5 for methane/airand propane/airrespectively.Yet, Eq. 5.2 given for duct length
in between3m and 6m seemsto give under predictedvalue of 1.03 bar for methane
and 1.44bar for propaneif comparedto Eq. 5.4. Summaryof the predictedvaluesbased
on Eq. 5.1-5.4andexperimentalresultsis shownin Table6.1.
194
o Smplywented,endigrýbon C,2Hýai r
9. 1
la 9mpdy\ented,contreigniton
8,0 -ý
7.0-1 6 Duct\ented, Dp-0.162m,
; end ignifion
EC) -ý
ý --o--Ductm. -nted, Op-o. 162rn,
5.0- centre ignißon
40 --X- DUC manted,Dp - 0.315 m,
-'
erd ignifion
NE Duct \ented, Dp=0.315m,
20 cantre ignibon
7
1.0-",
0.0
vent flow velocitieswould createless turbulencein the vent ductsand the presentlarge
back pressureswould be eliminated.The effect is essentiallythe sameas using leaner
nlixtures with the At
presentcase. lower K the velocities in the vent pipe would not be
ý
ao --4a-GnI3lyeritect cenre
lgrition
SO, Duct \Wecý Dp-0.162n
end igriticn
4.0-4 OW %Me Dp-0.162rr%
cergreigritlan
ao )E Dict writecý Dp. 0.315ril
end igrition
2-0 j1
ME D-Et NMe DP,431. IM
caire Igrition
1.0
0.0 ! fvw-'-
(130 0.35 0.40 0.45 a5o 0.55 0.60
Ec;Av.-derr,e redio
bum-up mechanism inside the duct can be ruled out to be the important role for the final
overpressure inside the vessel in the case of centre ignition. Conclusion can be drawn
based on this observation; the significant rapid rise of final pressure inside the vessel or
'detonation spike' depends on the flame instabilities (sudden venting can give rise to
flame instabilities and rarefaction wave (Dorofeev, Bezmelnitsin and Sidorov, 1995)),
pressure rise or 'spike detonation' observed in the hydrogen/air tests. The details of
deflagration-to-detonation mechanism are also discussed in Chapter 4.2 for simply
vented case. It then can be postulated that with central ignition, the bigger quantity of
fresh mixtures left for the residual combustion in the vessel is the main contribution of
the high magnitude of 'spike detonation' pressure in comparison with the pressure peak
attained at end ignition inside the vessel. However, as observed before, 'detonation
spiky pressure' is not exhibited in methane/air as methane does not auto-ignite easily
but it can undergo highly turbulence combustion.
Fi2/,-air
1.2
End igribon * Time of spike (simply%ented)
0.6 6 \A3sseloDmer(Dp=0,315m)
'o 0.4
X--
0.2
0.0
0.30 0.35 0.40 0.45 0.50 0.55 0.60
EqLi\ederceratio
The implication from this foregoing is that, if the flammable mixture is initially rich in
duct may be much severe than expected on consideration of the magnitude of simply
vented case, due to the high flame speeds and hence the high unbumt gas velocity ahead
of flame prior to duct entrance in the case of hydrocarbon/air. To get more insight to
this phenomena, the flow interaction between the main vessel and the duct were
CH4/bj r 4) = 1.08
2-1 Erd igribon k-tcx
for Dd = 0.315 m
1.7
0.9
0)
0.1
0.1
-0.3
-0.7 Time(s)
-1.1
The same trend observed where the maximum peak of Pma,occurred after the flame
exited the duct for both cases i. e. duct pipe diameter of 0.162 m and 0.315 m as shown
in previous chapters and there is no evidence for external explosion occurred in dump
vessel (indicated as P6 in Fig. 2.3 b). It is worth mentioning at this point that while
individual test was presented here for clarity, the results explained are indicative of the
trends displayed for this test vessel for all studied gas/air. At this point, methane/air at (D
= 1.08 would be the best representative to describe and explain the phenomena occurred
during the explosion development in comparison with the smaller duct pipe diameter
i. e. 0.162 m. The pressure difference between the explosion vessel and the vent pipe is
shown in Fig. 6.9 as a function of time, together with the vent duct explosion pressure.
198
This shows that when the flame was in the duct, there was a negative pressure
difference,which was higher for the large vent duct attached.The pressuredifference
betweenthe explosionvesseland the duct in the initial stageof the explosioncan be
to the
used compute meanvelocity of unburntgas into the duct,during the periodbefore
the flame enteredthe vent duct. The dynamicheadpressureloss for a pipe inlet for
incompressibleflow is 0.5. The pressuredifferencereachesa maximum of 0.2 bar,just
prior to the flame entry in the vent duct, as shownin Fig. 6.9. This correspondsto a
meanunburntgasvelocity in the vent of 258 m/s. However,the sharpedgeto the vent
will produce a vena contraction and the velocity at the vena contraction, using a
contraction'coefficient of 0.61, would be 423 m/s. This is closeto the speedof sound
(360 m/s) and the compressibleform of Bernoulli's equationshould be used.However,
this shows that very high unburnt gas velocities close to sonic conditions were
generatedat the vent and thesecreatehigh turbulenceconditionsin the vent pipe which
the shock waves will be generatedand this createsa high back pressureand the
observedreverseflow back into the explosionvessel(Kasmani,Andrews, Phylaktou,
Willacy, 2007, Ponizy and Leyer, 1999a,1999b).The details of the pressureloss and
unburntgasvelocity will be discussedlater in Section6.3.
The implication of this action is reverseflow causedby the larger duct diameterwould
result in largerscalemixing into the vesselas shownin Fig. 6.10, This phenomenahas
also been found by other workers (Kasmani, Andrews, Phylaktou, Willacy, 2007,
Ponizy and Leyer, 1999a,1999b).When the pressuredifferencebetweenthe vent duct
and the vesselwas high, the vent duct pressurebecamelow and this was the greatest
199
differencefor the large vent duct. For the small vent duct the staticpressurein the duct
to
wassimilar magnitude the difference
pressure between
the and
vessel the duct.
......... .
I ........ ..
...........
1.4
......... ...........
............ ",
bu'rned gases fresh mixture
..............
...........
:
--
..............
Figure 6.10 Flow regimesin vesselbefore and after explosion-likecombustionin duct,
effect of the flame front distortion-,
a) narrow ducts b) large duct- Reproducedfrom
(PonizyandLeyer, 1999b)
to very high flame speedsin the vent pipe. It is consideredthat this phenomenawas
occurringin thepresentwork with the large diameterventpipe and much smallervent.
However, using the bigger duct pipe diameter seemsto agreewith the theory stated
abovefor leanmixturesconcentrationwherethe overpressureinside the vesselis lower,
comparedto the pipe havingthe samediameterasthe vent asshownin ethylene/airat (D
200
= 0.8 for end ignition (Fig. 6.11). Apparently, the peak pressure for duct pipe of 0.162
m diameter reached its peak earlier after the flame exited the duct in comparison with
explosion with larger vent and simply vented. A possible explanation for this lies to the
Q2F4/,!r 2.7
9.0
End igrition
8.0 Sn-oy
2.4
\oened
2.1
7.0 - Duct wried, Dp = 0.1 E2 m
0.0 0
Figure 6.11 Pressure time histories and flame arrival for ethylene/air at (D = 0.8.
For leaner combustion at relatively small vent area, the outflow of unburnt gases into
the larger vent duct is less turbulent and hence, the combustion-] i ke-explosion event
inside the duct is less severe. Hence, the final explosion inside the vessel is reduced as
from the vessel but also to earlier arrival of the flame front at the duct entrance (refer to
Fig. 6.11 and 6.13). However, it is apparent that the combustion is still occurred inside
the 0.162 m duct pipe when the peak pressure inside the vessel marked its peak,
suggesting that there is still strong interaction between the vessel and the duct before the
flame exited the duct pipe and hence, producing a greatly increased pressure rise in
The study by lida et al (lida, Kawaguchi and Sato, 1985) mentioned that flame was
found to extinguish or hesitate in the channel before passing through in some cases,
depending on the equivalence ratio of the mixture, the channel with and the flame
inflow velocity. Other studies supported the above hypotheses by using relatively
201
C2H4/ar
1.2
e Lixt vened, Dp = 0.1 62m end
1.0 igriftion
9 ElictveriKDp=0.162rTi
cet14e igrition
0.8 D-ict vertad. Dp = 0.315 ni
x
Eyid igrýtion
0.6 w Duci veried, Dp = 0.315 m
ceritre igiitici
0.4
02
.
0.0
0.50 0.55 0.60 0.65 0.70 0.75 0.80 0.85
EquivaJencerclio
)( FIamearri\d, Dp=0.315m
0.2 1.0
-----------------------------------
Vessd region
0.1 0.5
-0.1 -0.5
Time(s)
FigUrc6.11 Pressuretime historiesand flanic arrival for propane/airat (1)= 0.8 for end
ignition.
203
On
--4- DuctVented,
mendignihon
Dp 0.315
CH4/aIr b)
Duelvented, I)p 0.316nk 600- CH41air
CentreIgnition -4- Dp - 0.315m, end
Ignition
-X- Ductvertled,Do. 0.162m -ý 500 Do - 0.315m, tentral
andignton I IRoo
milon,
-*-Duct vernedI)p- 0.162n a - 0.162 m, ond
CentreIllnition 400 - -a-
-0 Ignition
SwO vented, end Onown
-*- Dp - 0,162,conlrid
300 ignition
ý50" Vented,Caml,m)
200
100
0.4 0.5 0.6 0,7 0.8 0.9 1 1.1 1.2 0.4 0.5 0.6 0.7 0,11 0.9 1 111 1.2
Soulvalanceratio Equivalenceratio
and T3 as in Fig.2.3b) and (b) in the vent duct (betweenT4 and Tg) as a function of
equivalenceratio for methane/air.
The flame speedsin the main vesselapproachingthe vent were considerablyhigher than
for sphericallarninar flame speeds.This was due to two effects:firstly, self acceleration
of the flame fluough the cellular flame front mechanismand secondly,the suctioneffect
204
of the vent dischargeon the flame shapewhich would draw the flame expansion
preferentiallyin the direction of the vent (Kasmani,Andrews,Phylaktou and Willacy,
2007a, Kasmani,Andrews, Phylaktouand Willacy, 2007b). These effects were both
higher for end ignition as the distanceto the vent was double that for central ignition
and hence the two effects were enhanced.This was the main reason why the end
ignition gavethehighestoverpressures
comparedwith centralignition.
The flame speedsinside the duct were similar for both duct diameters,apart from the
equivalenceratio at (D= 1.08wherethe highestflame speedsof 490 m/s experiencedat
LJD = 3.18.The larger vent duct createda flow expansionfrom the vena contractionat
the inlet vent to the duct wall. This flow expansioncreatesa pressureloss that is the
sourceof the turbulencethat accelerates the flame which is largerwhen there is a larger
flow expansion.For the presentgeometry,the unbumt gas dynamic head pressureloss
in termsof the vent areadynamichead,was0.47 for the 0.162m vent pipe and 1.27 for
the 0.315 m diametervent pipe. This producedmore turbulenceand a greaterflame
accelerationof the flame inside the larger vent pipe, as shown in Fig. 6.14. Also the
lower meanvelocitiesin the larger pipe would enablea flame to propagatein regions
wheretherewaslocal turbulent quenchingin the smallervent pipe (Kasmani,Andrews,
Phylaktouand Willacy, 2007a,Kasmani,Andrews,Phylaktouand WillacY, 2007b) and
this would increasethe back pressure,as found experimentally.For leaner mixtures the
velocities were much lower and the turbulencegenerationwas significantly lower as
this is proportionalto the squareof velocity. Hence, the effect of the vent pipe was
muchlower for the slowerburning leanermixtures.
However, the trend for leaner concentrationresults in lower flame speedsin the duct
seemsnot applicablein the caseof hydrogen/airexplosion.As illustrated in Fig. 6.15,
the flame speedsinside the duct was higher for all equivalenceratio for vent duct the
samesize as the vent comparedto 0.315m diameterduct for end ignition and marking
its highestvalue at 484 m/s. This high flame speedsare due to the comparatively fast
flame from the vesseldue to large flame areaeffect. Studyby Aung et al (Aung, Hassan
and Faeth, 1997)showedthat hydrogenat lower equivalenceratio 0.3 < (D < 0.6 are in
unstablepreferential-diffusionconditions where developedchaotically irregular name
surfaces.
205
a) (b)
H21air
1 -4- " veriled.endgnam fiýair 600
100 Ductvented.
Dp- 0,162m,
end
Z
90
-"
Ignition
w4d'sawpow
500
-4-DUCIVOnted.00-0.142m,
'm I
voted, DP-0.162lkWd cenkoIgma"
-*-Dxl
70 4"Alo J? D-uclvented,on 0,31Sn end
400
Dict, ol-G. Dp- 0,162mcatwe ig*on
60 :p
Ductvented,
Dp 0.315m,
50 300
contisignAmm
40 1 ýDXIVWW np-0315M emit
200
a 30
20 100
r
w 10
0
0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.30 0.35 0.40 0.45 0.50 0.55 0.60
Equivalence
ralio EquWlertce
ratio
where SL_is the value of larninarburning velocity when the flame stretch,K=0, SL is
the Ian-tinarburning velocity, Ka is the Karlovitz number and Ma is the Markstein
numbergiven by equationsasbelow respectively,
Ka = K,61ý/ SL (6.3)
8D=
D. 'SL (6.4)
206
Ma = L1,5D (6.5)
K= flame stretch
Du = massdiffusivity
L= Marksteinlength
A similar effect would occur for smaller K or larger vent diameters.The vent flow
velocitieswould be reduced,lower turbulencelevels in the pipe would be generatedand
lower overpressure
increasesdue to the useof a vent would occur. It is significant that
much of the publishedwork on vent pipes, which all showsa large increasein the
overpressurecomparedwith a free vent discharge,was all carried out with relatively
large valuesof K. More works on vent pipes has beencarriedout for dust explosions,
but generallywith K> 10. It is consideredthat in view of the limited experimentaldata
on the impactof a ventduct on the overpressurefor gaseousexplosionand the potential
importanceof K and mixture reactivity (which determinesthe vent flow velocity) that
more works is requiredto understandthis type of venting phenomenaand to provide
morereliableventingdesignguidance.
As mentioned earlier in Section 6.1, the pressure loss prior to the duct pipe entrance
(indicated as AP2-3)and inside the pipe (AP3-6),will be used to calculate the unbumt gas
velocity prior to duct entrance as well as the gas velocity inside the duct to ascertain the
enhanced burning rate due to turbulisation and enhanced pressure drops due to the
geometry are playing the important role to the final pressure inside the vessel. This can
then leads to the comparison of the pressure loss for the vent pipe as being the same as
the vent orifice and larger vent pipe. This simple analysis can be evaluated using;
APi,,+ APd,,,,
+ AP.,,,=1 pSg'(Kin + 4f -ýL+ Ni + K,,,,,) (6.6)
2D
I ()1ý1'
tit ;I
ýI.,.
II
(III
. .1ýý
,I--I,
11-I'laMCdCVCIOjIIIICllt
Figurc 0.16 Sclifici-cii cinc photographS01'CCIIUI, CXIIIOSiOll
(IL11-ilig
of' H-,/air IIIIXILII'C,(1) = 0.26. Reproduced 1'rom Bradley and Harper (Bradley and
Harper, 1994).
209
Ki,, and K t are respectively the pressureloss coefficients for sudden flow area
..
(in
restriction/enlargement the caseof a flow to larger pipe from the vessel)and sudden
flow areaenlargement(pipe to dump vessel).f is the friction factor for he flow inside
the duct evaluatedfrom Darcy-Weisbachequation,Ni is the velocity heads lost in
fittings in this presentcasefor duct pipe of 0.162m,one gatevalve and two couplings
in the duct; valuesreportedin (Ferrara,Willacy, Phylaktou, Andrews, Benedettoand
Salzano,2005)wasused.Below is the list of parametervaluesused;
Kin = 0.5; K t=0.95; f=0.005; 1: N; = 0.25
...
1
APin =K psg (6.7)
2
O(geometry,
Re) (6.8)
where Re is the pipe Reynoldsnumber given by Re = pSRD /,a . For many practical
applications,the Re no is large enoughso for the flow inside the pipe where inertia
effect is more dominated rather than the viscous effects, is usually found that the
pressuredrops and head lossescorrelatedirectly with the dynamic pressure.Thus in
most cases,the loss coefficientsfor componentsare a function of geometryonly as Eq.
6.8 will be K= 0(geometry)(Munson, Young and Okiishi, 2006). Meanwhile, for
lossinsidethe pipe,
presSUTe
I
2 L+
PSv (4f Nj) (6.9)
2D
From Eq.6.6,6.7 and 6.9, unburnt gas velocity on each phasecan be evaluatedand
determinedaslisted in Table6.2.
210
From this analysis,it canbe confirined that the high unburritgasvelocity inside the duct
induceda very high turbulent level and thus a severesecondaryexplosion inside the
duct.The secondaryexplosionin the duct which in turnsaffectsthe residualcombustion
in the main vesseland showsmore violent for endignition as shownin Fig. 6.4 and 6.5.
However,the higher intensity of the secondaryexplosion does not affect literally in
practisethe final pressureinside the vesselin the caseof central ignition especiallyfor
high burning velocities mixtures i. e. hydrogen/airand ethylene/air.It is postulatedby
others (Ferrara, Willacy, Phylaktou, Andrews, Benedetto, Mkpadi, 2005, Ferrara,
Willacy, Phylaktou, Andrews, Benedetto,Salzano,Russo, 2008, Ponizy and Leyer,
212
6.4 Concludingremarks
Enlarged vent ducts, i.e. ducts with areas greater than the vent area gave vent
similar to those for free vents for lean mixtures,but for the maximum
overpressures
reactivity mixture ((D =1.08 for methane/airand (D = 1.375for propane/air)P ..a, was
similar to that with a vent ductthe samediameterasthe vent.
The causeof the large increasein overpressurefor both ducts in the case of rich
concentrationfor methaneand propane/airmixtures was due to the high induced
unburnt gas velocity into the vent and inside the vent duct itself. For the presentK=
16.4,this conditioncreatednearsonicflow conditionsat the vent venacontraction.The
arrival of the flame in the vent createdsonic flow in the vent duct and the high back
pressurecreateda reverse flow into the explosion vessel.This high turbulent flow
acceleratesthecombustionof remainingunburnt,mixture andthis further acceleratesthe
flow in the ductcreatingthe peakoverpressure.
impact of K, mixture reactivity and static burst pressureon vent design with vent
dischargeductsattachedif morereliabledesignguidanceis to be given.
CHAPTER 7
usedin Method 1. The net effect is as Sgis close to the flame speed,Sf in value, the
approachis only slightly lower ventmassflow ratethanthat basedon Sf.
A significant flame self accelerationeffect for subsonicventing was shown for K<-5
and this effect is similar to vent inducedturbulenceand could be accountedby the P
term in the burning velocity equation,The turbulenceenhancementpredictedbasedon
Bartknecht'sequationand proposedMethod 2 was in a good agreementwith P derived
from tabulatedexperimentaldata. It can be said that the 0 derived was perfectly
reasonablevaluefor P usedby other works (Munday, 1963,Pasman,Groothuizenand
Gooijer, 1974,Yao, 1974).However,the flame experiencesdecelerationeffect in larger
volumeswhenpressureis high in larger K and this effect has never beenhighlighted
previously. It is postulated that at high K and Predwith sonic venting during the
explosions, the self-accelerationis likely to have already occurred at the smaller
volumes.
It can be suggestedthat the use of K term is more suitableto be applied on cubic and
The high overpressuredue to the addition of a vent duct, in a 0.2 m3 cylindrical vessel
with an LID of 2 and K of 16.4,occursafter the flame has exited the vent duct, but is
not due to an external explosion. There is found that substantialamount of unburnt
gasesJeft insidethe vesselafter the vent burstsis the leadingfactor in increaseof P. ý,"
for high burning velocity mixtures at centrally ignited, The associategas velocities
aheadof the flame createhigh unburntgasflows conditionsat entry to the vent and this
give rise to high back pressures.
when the vent burst,as the flame has had more time to grow upstreamof the vent and
this givesa highervent duct flow velocity oncethe vent cover bursts,The effect of the
vent burst pressureis complex and non-linear and was not representedby the linear
effect in the correlationof Bartknecht(Bartknecht,1993)usedin NFPA 68 (NFPA68,
2002).
Enlarged vent ducts, i. e. ducts with areas greater than the vent area gave vent
similar to those for free vents for lean mixtures, but for the maximum
overpressures
reactivity mixture ((D =1.08 for methane/airand 0=1.375 for propane/air)P,,,,,was
),
similar to that with a vent duct the samediameterasthe vent.
The causeof the large increasein overpressurefor both ducts for rich concentrationin
low burningvelocity mixtures wasdue to the high inducedunburntgasvelocity into the
vent and inside the vent duct. For K= 16.4, this condition created near sonic flow
conditionsat the vent venacontraction.The arrival of the flame in the vent createdsonic
flow in the vent duct and the high back pressurecreated a reverse flow into the
218
It can be said that the ignition position play an important role in determiningthe final
pressureinside the vesselwhereend ignition gave highestPn,,,,with respectto central
ignition for methane/airandpropane/airin all cases.End ignition gives greaterdistance
towardsthe vent andhence,theflameacceleration overtwice thedistanceasfor
continuous
the centralignition,However,for high burning velocity mixturesi. e. hydrogen/air,the
trend is seemsnot to havea goodagreementwith centralignition gavehighestPmaxwith
a profoundly 'spike detonation' peak. The major finding is the occurrenceof auto-
ignition at the unburnedpocket of mixturestrappedinside the vesselwhich eventually
playing a significantrole in the hydrogen/airexplosionsbut not in methane/air.It then
can be postulatedthat at central ignition, the bigger quantity of fresh mixtures left
inside the vesselmarked as the main contribution of the high magnitudeof 'spike
detonation'pressurein comparisonwith the pressurepeakattainedat endignition inside
the vessel.
validate and quantify the A,/A, term to P,,,,, in venting design insteadof K for non-
cubic vessel.
The presentdesign correlationsfor explosion vents and with vent ducts, which are
basedon very limited experimentaldata, do not predict the presentresults and their
reliability for small vesselvolumes with high K is in doubt. In order to justify the
reliability and applicabilityof the currentdesigncorrelationon this subjectmatter,more
work on smaller vesselwith high K with different mixture reactivity are needed.The
useOf instead
Pdynamic of Pvis morefavourable in orderto correlatethe influenceof vent
cover to Pmax, dP/dt and flame speeds.It is recommendedto use commercial vent
coverfor future work in orderfor vent coverinertia effect to be taken into account.
more precise information regarding the mechanism developed during vented gas
explosion,Further work is requiredin the areaof the impact of K, mixture reactivity
and static burst pressureon vent design with vent dischargeducts attachedif more
reliabledesignguidanceis to be given.
modelof gasflow needto be introducedin orderto take into accountthe fluid dynamic
effect and the dynamicsof flame propagationand associatedpressureand temperature
during the ventedgasexplosion.Further,to mitigatethe severityof unbumt gasvelocity
at the duct, the useof flame arresterat the venacontractaregion would be the solution
to minimise the intensificationof secondaryexplosioninside the vesseland hence,less
physicalbackflow to the decreasing
vessel, the final inside the vessel.
overpressure
REFERENCES
221
RERERENCES
EuropeanStandard:GasexplosionventingguidanceEN 14994:2007.
Abdel Gayed,R. G. & Bradley, D. (1981) A two-eddy theory of premixed turbulent
flamepropagation.TransProceedingsof RoyalSocietyof London,A 301,1-25.
Andrews, G. E. (2004) Course notes on gas and dust explosionprotection design,
Universityof Leeds,United Kingdom.
Andrews, G. E. & Bradley, D. (1972) Determinationof burning velocity: A critical
review. Combustionand Flame, 18,133.
Andrews, G. E. & Bradley, D. (1973) Determinationof burning velocity by double
ignition in a closedvessel.Combustionand Flame,20,77.
Andrews, G. E., Bradley, D. & Lwakabama,S. B. (1975) Turbulenceand turbulent
flamepropagation-acritical appraisal.Combustionand Flame,24,285-304.
Anthony, E. J. (1978) The use of venting formula in the design and protection of
building and industrial plant from damageby gas or vapor explosions. Journal of
HazardousMaterials, 2,23-49.
Aung, K. T., Hassan,M. 1. & Faeth,G. M. (1997)Flame stretchinteractionsof laminar
prernixedhydrogen/airflames at normal temperatureand pressure,Combustionand
Flame, 109,1-24.
Bartknecht,W. (1985) Effectivenessof explosionventing as a protective measurefor
silos.PlantlOperationalprogress,4,4-12.
Bartknecht,W. (1993)Explosions-schultz, Berlin, Springer-Verlag.
Bartknecth, W. (1981) Explosions course prevention protection, Berlin,New York,
Springer-Verlag.
Bell, J. B., Cheng,R. K., Day, M, S. & Shepherd,I. G. (2007) Numerical simulation of
lewis numbereffectson leanpremixedturbulentflames.Proceedingsof the Combustion
Institute,31,1309-1317.
Bradley, D. (1999) Self accelerationof flames. 2nd Fire and F-TplosionsHazards
Seminar.
Bradley,D., Cresswell,T, M. & Puttock,J. S. (2001)Flame accelerationdue to flame-
inducedinstabilitiesin large scaleexplosions.Combustionand Flame, 124,551-559.
Bradley, D., Emerson,D. R. & Gu, X. J. (2003) Modes of reaction front propagation
from hot spots.Combustionand Flame, 133,63-74.
222
Ferrara, G., Benedetto, A. D., Salzano, E. & Russo, G. (2006) CFD analysis of gas
effects of changing the position of ignition source. Journal of Loss Prevention in the
Process Industries, 4,217-222.
225
Lunn, G., Crowhurst,D. & Hey, M. (1988) The effect of vent ducts on the reduced
explosionpressuresof venteddust explosions.Journal of Loss Prevention in Process
Industries,1,182-196,
Maisey, H. R. (1965) Gaseousand dust explosion venting- Part 1. Chemical and
ProcessEngineering,527-535.
Markstein,G. H. (1964)Nonsteadyflame propagation,New York, Macmillan.
McCann, D. P. J., Thomas,G. 0. & Edwards,D. H. (1985) Gasdynamicsof vented
explosionspart i: Experimentalstudies.Combustionand Flame,59,233-250.
Moen, 1.0., Bjerketvedt, D., Jenssen,A. & Thibault, P. A. (1985) Transition to
detonationin a large fuel-air cloud (brief communication).Combustionand Flame, 61,
285-291.
Moen, 1.0., Lee, J. H. S., Hjertager,B. H., Fuhre,K. & Eckhoff, R. K. (1982) Pressure
developmentdue to turbulentflame propagationin large-scalemethane-airexplosions.
Combustionand Flame,47,31-52.
Molkov, V., Dobashi,R., Suzuki, M. & Hirano, T. (2000) Venting of deflagrations:
Hydrocarbon-airand hydrogen-airsystems.Journal of Loss Prevention in the Process
Industries,13,397-409.
Molkov, V. V. (1994)Venting of deflagrations:Dynamic of the processin systemswith
materials,Al 16,1-10.
Molkov, V. V., Grigorash,A. V., Eber, R. M., Tarnanini,F. & Dobashi, R. (2004)
Vented gaseousdeflagrationswith inertial vent covers-State-of-the-artand progress.
ProcessSafetyProgress,23,29-36.
Molkov, V. V., Kokolchenko,A. & Alexandrov, S. (1997)Venting of deflagrationsin
buildings and equipment:Universalcorrelation.Fire SafetyScience-proceedings
of the
5th InternationalSymposium,1249-1260.
Molkov, V. V. & Nekrasov,V. P. (1981) Dynamicsof gas combustionin a constant
volume in the presenceof exhaust.Balashika. Translatedfrom Fizika Goreniya i
Vzryva,17,17-24.
Munday, G. (1963)The calculationof venting areasfor pressurerelief of explosionsin
vessel. Proceedingsof Symposiumon Chemical Process Hazards with Special
Referenceto Plant Design Series No-15. Manchester,The Institution of chemical
Engineers.
Munson,B. R., Young,D. F. & Okiishi, T. H. (2006)Fundamentalsoffluid mechanics:
Fifth edition,JohnWiley andSonsInc.
Nagy, J. & Verakis, H. C. (1983) Developmentand control of dust explosions,New
York, MarcelDekkerInc, New York andBasel.
NFPA68 (2002) Nfpa 68: Guide for venting of deflagrations:2002. National Fire
ProtectionAssociation.
Ng, H, D. & Lee,J. H. L. (2007)Commentson explosionproblemsfor hydrogensafety.
Jounal of Lossof Preventionin ProcessSafety.
Pasman,H. J., Groothuizen,T, M. & Gooijer, H. d. (1974) Design of pressurerelief
vents. Loss Preventionand Safety Promotion in the ProcessIndustries: Edited by
C.H.Buschmann,185-189.
Phylaktou, H. N. & Andrews, G. E. (1991) Gas explosionsin long closed vessels.
CombustionScienceand Technology,77,27-39.
Phylaktou,H. N., Andrews, G. E. & Iferath, P. (1990) Fast flame speedsand rate of
pressurerise in the initial period of gas explosionsin large l/d cylindrical enclosures.
Journal of Lossof Preventionin ProcessSafety,30355-364.
Phylaktou,H. N., Foley, M. & Andrews, G. E. (1993)Explosionsin a tube with a 90"
bend.Journal of Lossof Preventionin ProcessSafety,6,21-29.
228
Methane-air
10.4 0.2 1.93 Rec 1.5 Centre 0.03 0.02 0.03)15 RX I. /alosh. K 6. Vent covei is
alumill Itill)
10.4 0.2 1.93 Rec 1.5 Centre WO 0.00 0.0 172 K.( i. /alosh. KII. Vcnt cover is
alumillitini
10.4 0.2 1.9-1 Rec 1.5 Centre 0.0 0.21 0.0007S K.( J.Aalosh. K 28. Vent cover is
almilinimll
10.4 0.2 2.01 Rec I Centre 0.03 0.0 0.054 R.G. /alosh. K 5.8. Vent cover is
aluillillillill
10.4 0.2 2.01 Kec I Centre 0.01 0.06 0.031 R.Gzalosh. K 10. Vent cover is
alLill1illiLlIll
10.4 0.2 2.01 Rec I Centre 0.03 0.21 0.01(1 R.G./Aosh. K 20. Vent cover is
altilifillimil
10.4 0.2 2.01 Rec I Centre 0.03 0.69 0.0077 R.G./alosh. K-41. Vent co\,cr is
alunlinium
10.4 0.2 2.08 Rec I Centre 0.0.1 0.01 0.045 R.C. /alosh. K- 5.8. Will cover is
aluillinillill
10.4 0.2 2.08 Rec I Cclitrc 0.03 0.04 0.0". 3 R.O./alosh. K 8. Vent covvr is
10.4 0.2 2.08 Rec I Centre 0.03 0.14 0.015 luj. /alosh. K 18. Vent cover i's
allimilliulli
10.4 0.2 2.08 Rec I CCIIII-C 0.0.1 0-48 0.0001) KAi. /alosh. K 18. Vent covci Is
.1111111inium
activation I-)I.Csstll.c
NR 0.4 3.25 Rec I Centre 0.00.) 0.1 0.1, Clibbagc and Sililmonds 1: 101)
I I activation pressure
NR 0.4 3.25 Rec I Centre 0.004 0.13 0.3 Cubbage and Simmonds 1: Top
alurninium
10.4 0.6 4.38 Rec 1.5 Centre 0.03 0.05 0.049 R.G.Zalosli. K-8. Vent cover is
aluminium
I I
10.4 0.6 4.38 Rec 1.5 Centre 0.03 0.76 0.036 R.G.Zalosli. K=l 1. Vent cover is
aluminium
10.4 0.6 4.38 Rec 1.5 Centre 0.03 0.36 0.0098 R.G.Zalosh. K=40. Vent cover is
aluminium
9.5 1.0 5.5 Cylinder 1.554 Cent e 0.16 1.0 0.1 P.F.Thome et at. Polyethylene
l30um, K=4,0. lgnition at the
center of rear wall
9.5 1.0 5,5 Cylinder 1.554 Centre 0.32 2.0 0.05 S.Chippett
9.5 1.0 6.0 Cubic I Centre O'l 0.2 0.36 Bartknecht(Abb.2.223a)
9.5 1.0 6.0 Cubic I Centre 0.2 0.2 0.36 Bartknecht(Abb.2.223(a))
9.5 1.0 6.0 Cubic I Centre 0.2 0.4 0.36 Bartknecht(Abb.2,223b)
9.5 1.0 6.0 Cubic I Centre 0.5 0.64 0.36 Razus et at
9.5 1.0 6.0 Cubic I Centre 0.5 0.8 0,36 Bartknecht(Abb.2.223b)
9.5 1.0 6.0 Cubic I Centre 0.8 0.9 0.36 Bartknecht(Abb.2.223(a))
9.5 1.0 6.0 Cubic I Centre 1.0 1.2 0.36 Bartknecht(Abb.2.223b)
9.5 1.0 6.0 Cubic I Centre 1.0 1.3 0.36 Bartknecht(Abb.2.223(a))
9.5 Lo 6.0 Cubic I Centre 1.5 1.7 0.36 Bartknecht(Abb.2.223(a))
9.5 1.0 6.0 1 Cubic I Centre 2.0 2.4 0.36 Bartknecht(Abb.2.223(a))
9.5 1.0 6.0 Cubic I Centre 0.1 0.43 0.3 Razusetal
9.5 1.0 6.1 Cylinder 1,55 Centre 0.18 0.5 0.3 Paper as vent cover. One peak
9.5 1.0 6.1 Cylinder 1.55 Centre 0.18 1.2 0.2 Paper as vent cover. One peak
pressure.Pasmaii,I-I.J et at
9.5 1.0 6.1 Cylinder 1.55 Centre 0.20 0.1 0.2 Mylar as vent cover. Double peak
9.5 1.0 6.0 Cubic I Centre 0.5 1.2 0.16 same as above
9.5 1.0 6.0 Cubic I Centre 0.5 1.4 0,16 Razuset al
9.5 1.0 6.1 Cylinder 1.55 Centre 0.20 1.0 0.1 Mylar as vent cover. Double peak
I I I pressure.Pasman,113 et al
9.5 1.0 6.0 Cubic I Centre 0.5 1.95 0.1 Razuset al
9.5 1.0 6.1 Cylinder 1.55 Centre 0.20 2.0 0.05 Mylar as vent cover.One peak
pressure.Pasman,ti. 3 et al
9.5 1.0 6.0 Cubic I Centre 0.1 1.8 0.04 Bartknecht(Abb.2.223b)
9.5 1.0 6.0 Cubic I Centre 0.2 2.6 0.04 Bartkneclit(Abb. 2,223(a))
NR 2.8 12.46 Rec 1.8 Centre 0 0.014 1.14 Cubbage and Simmonds 1: Top
10.5 11.2 31.54 Rec 0.9 Centre 0.03 0.5 2.25 R.G.Zalosh.P2 is taken.
I Aluminium foil as a membrane
10.5 11.2 31.54 Rec 0.9 Centre 0.03 0.3 0.77 R.G.Zalosh.P2 is taken.
I Aluminium foil asa membrane
10 22.0 47.11 Cubic I Centre 0.012 0,016 7.85 P.F.Thorne et al. Polyethylene
130gm:K=I
10 22.0 47.11 Cubic I Centre 0.03 0.04 3,41 P.F.Thorne et al. Polyethylene
I 130gm;K=2.3
10 22.0 47.11 Cubic I Centre 0.05 0.02 3.41 P.F.Thorne et al. Polyethylene
130gm; K=2,3: Ignition at the
centreof front wall
10 22.0 47.11 Cubic I Centre 0.05 0.07 3.41 P.F.Thorne et al, Polyethylene
130i,im; K=2.3: Ignition at tile
centreof rear walI
10 22.0 47.11 Cubic I Centre 0.05 0.06 3.41 P.F.Thorne et al. Fibrcboard
I 12.5iim,K=2.3
10 22.0 47.11 Cubic I I Centre 0.05 0.16 3,41 FI.F.Thorne et al. Polyester50lim;
K-2.3
10 22.0 47.11 Cubic I Centre 0.04 0.07 1.96 P.F.Thorne et al. Polyethylene
130gm;K-4.0
10 22.0 47.11 Cubic I Centre 0.05 0.11 1.96 P.F.Thorne et al, Polyethylene
I I I I 1 130^ K-4.0* Ignition at tile
234
I I
centreof rearwall
10 22T 47.11 Cubic I Centre 0.09 0.6 0.98 P.F.Thorne et aL Fibreboard
12.5pm;K=8.0
10.0 26.6 56.1 Rec I Centre 0.02 0.07 0.50 Buckland(G.A.Lunn book)
10.0 26.6 56.1 Rec I Centre 0.06 0.11 0.50 Buckland(G.A, Lunn book)
10.0 26.6 56.1 1 Rec I Centre 0.06 0.05 0.50 Buckland(G.A, Lunn book)
10.0 26.C 56.1 Rec I Centre 0.07 0.06 0.50 Buckland(G.A. Lunnbook)
10.0 26.6 56.1 Rec I Centre 0.04 0.11 0.40 Buckland(G.A. Lunn book)
10.0 26.6 56.1 Rec I Centre 0.07 0.11 0.40 Buckland(G.A.Lunn book)
10.0 26.6 56.1 Rec I Centre 0.02 0.10 0.25 Buckland(G.A.Lunn book)
10.0 26.6 56.1 Rec 3 Centre 0.07 0.07 0.52 Buckland(G.A. Lunnbook)
10.0 26.6 56.1 Rec 4 Centre 0.08 0.13 0.63 Buckland(G.A.Lunn book)
10.0 26.6 56.1 Rec 2 Centre 0.09 0.07 0.40 Buckland(G.A.Lunn book)
10.0 26.6 156.1 Rec I Centre 0.09 0.22 0.20 Buckland(G.A.Lunn book)
10.0 26.6 56.1 Rec I Centre 0.12 0.22 0.20 Buckland(G.A.Lunn book)
10.0 26T 56,1 Rec - I Centre 0.007 0.08 0.13 Buckland(G.A.Lunn book)
10.0 26.6 56.1 Rec I Centre 0.04 0.06 0.13 Buckland(G.A,Lunn book)
9.5 30.0 57.93 Cubic I Centre 0.1 0.1 3.9 Bartknecht(Abb.
2.222).F= nx
I DN 1000
9.5 30.0 57.93 Cubic I Centre 0.1 0.2 3.4 Bartknecht(Abb.
2.222).F= nx
DN 1000
9.5 30.0 57.93 Cubic I Centre 0.1 0.4 2.5 Bartknecht(Abb.
2.222).F=n x
I I DN800
9.5 30.0 57.93 Cubic I Centre 0.1 0.5 2.0 2.222).F=n
Bartknecht(Abb. x
DN800
9.5 30.0 57.93 Cubic I Centre 0.1 1.0 L5 2.222).F=n
Bartknecht(Abb. x
I DN800
9.5 30.0 57.93 Cubic I Centre 0.1 0.9 1.4 Bartknecht(Abb.
2.222),F=n x
DN600
9.5 30.0 57.93 Cubic I Centre 0.1 1.4 1.0 Bartknecht(Abb.
2.222).F=n x
I DN800
9.5 30.0 57.93 Cubic I Centre 0.1 1.88 0.85 Bartknecht(Abb.
2.222).F=n x
DN600
9.5 30.0 57.93 Cubic I Centre 0.1 2.2 0.8 Bartkneclit(Abb.
2.222).F= 11 x
I DN 1000
9.5 30.0 57.93 Cubic I Centre 0.1 2.6 0.5 Bartkneclit(Abb.2.222).F=n x
DN800
9.5 30.0 57.93 Cubic I Centre 0.1 3.7 0.3 Bartknccht(Abb.
2.222).F-n x
I I I I I I I I IDN600
10.5 34.0 66.12 Rec 2.7 Centre 0.0 1 1.0 3 RA Ualosli. pl i.,, taken.
10.5 34.0 66.12 Rec 2.7 Centre 0-031 2.2 1 R.( ;,zalo."JI.P1 I,, taken.
-
Altumnium 1'()ilas a Ilicillbralle
10.5 34.0 66.12 Rec 2.7 Centre 0.011 1 25 KA I./alosh. P2 11.1 taken.
.3
Aluminlum 1'()il Illembralle
10.5 34.0 66.12 Rec 2.7 Centre 2.25 R.(,. ZaIosh.P2 is taken.
Aluminlum I'Oil Inell1hralic
(Tentre - 7 -
)-031 1.6 5 R.G.Zalosh. A 0.91-m dianicler
10.5 34.0 06.12 Rec 2T
and 2.7m long tank and an arraý ol
2.5cm pipc, suspended froll,
Propane/air
5.0 0.03 0.61 Vessel 1.54 Center 4.8 10.9 0.055 Cousin &, Cottons( 1051
5.0 0.03 0.61 Vessel 1.54 Center 29.0 32.1 0.055 l Cousin & Cottons( 1951
5.0 0.03 0.61 Vessel 1.54 Center 33.1 331.1 0.055 Cousin & Cottons( 1951).
5.0 0.03 0.61 Vessel 1.54 Center 16.4 20.7 0.055 Initial pressure is 4-5 psig(.'). 10
5.0 0.03 0.61 Vessel 1.54 Center 16.6 23.4 0.021 COLISin& Cottons( 1051
5.0 0.03 0.61 Vessel 1.54 Center 4.8 15.2 0.021 Cousin & Cottons( 1951
0.03 0.61 Vessel 1.54 Center 24.8 30.3 0.021 Cousin & Cottons( 1()51
-5.0
5.0 0.03 0.61 Vessel --1.54 -Uentci* 4.9 22.1 0.01 1 Cousin & Cottons( 11)51
ýý 54 Center 18.6 29 0.011 COLISin& Cottons( 1951
5.0 0.03 0.61 Vessel
.
5.0 0.03 0.61 Vessel 1.54 Center 20 33.1 0.011 cousin & cottons( 1951
5.0 0.085 1.134 Tank 2.3 Centre 1.5 2.8 0.0417 Cousin A, Cottons( 1951
5.0 0.085 1.134 Tank 2.3 Centre 9,0 8.0 0.0417 Cousin Colton.,"( 1951)
5.0 0.085 1.134 Tank 2.3 Centre 10.6 11.9 0.0417 Comill Cottons( 1951
5.0 0.085 1.1'14 Tank 2.3 Centre 0.0 0.1 0.0408 Cousin & Cottons( 1951)
5.0 0.085 1.134 Tank Centi 1.8 5.3 0.0 101 Cousin & Uottons(l 951
- -1 I .1
5.0 0.095 1.1'34 Tank 2.1 Centre 7.0 8.9 0.0 161 Cousin & Cottons( I )ý I
T-ank 2.3 Centre 12.7 13.3 0.0 161 CoUsin & Cottons( 195 1
236
5.0 0.085 1.1341 Tani_ _233 -Centre 0.0 1.1 0,0135 Cousin& Cottons(1951)
5.0 0.085 1.134 Tank 2.3 Centre 1.6 6.5 0.0080 Cousin & Cottons(I 95 1)
5.0 0.085 1.134 Tank 2.3 Centre 5.9 9.6 0.0080 Cousin & Cottons(I 95 1)
5.0 0.085 1.134 Tank 2.3 Centre 11.9 13.7 0.0080 Cousin & Cottons(I 95 1)
5.0 0.085 1.1341 Tank 2.3 Centre 0.0 2.0 0,0060 Cousin & Cottons(] 95 1)
5.0 0.085 1.1341 Tank 2.3 Centre 0.0 2.7 0.0036 Cousin & Cottons(l 95 1)
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.021 0.07 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.028 0.07 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 1.55 0.0525 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.048 0.0525 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 1.32 0.0438 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.0276 0.035 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.344 0.0314 R.G.Zalosh
4.60 0.17 1.9 Rec 1.5 Centre 0.03 0.028 0.0314 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.037 0.0314 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.28 0.0314 R.G.Zalosh
5.20 0.17 1.9 Rec 1.5 Centre 0.03 0.66 0.0314 R.G.Zalosh
1
5.20 0.17 1.9 Rec 1.5 Centre 0.03 1.61 0,0263 R.G.Zalosh
6.00 0.17 1.9 Rec 1.5 Centre 0.03 0.03 0.0175 R.G.Zalosh
4.80 0.65 2.9 Cylinder 1.198 Centre 0.86 2.75 0.0993 P.F.Thome et al
4.80 0.65 2.9 Cylinder 1.198 Centre 0.86 3.14 0.073 P.F.Tliorne et al
4.80 0.65 2.9 Cylinder 1.198 Centre 0.72 5.59 0.0324 P.F.Thorne et al
4.80 0.65 2.9 Cylinder 1.198 Centre 0.65 5.22 0.0182 P.F.Thome et al
4.80 0.65 2.9 Cylinder 1,198 Centre 0.65 5.31 0.0182 P.F.Thome et al
4.02 1.0 1.0 Cubic I Centre 0.1 0.28 0.38 Bartknecht data(Abb.2.219)
5.00 1.0 4.8 Sphere I Centre 0.1 0.3 0.31 C.Donat(l 973)
4.02 1.0 1.0 Cubic I Centre 0.1 0.55 0.23 Bartknecht data(Abb.2.219)
4.02 1.0 1.0 Cubic I Centre 0.1 0.92 0.18 Bartknecht data(Abb.2.219)
4.02 1.0 1.0 Cubic I Centre 0.2 2.0 0.1 Bartkneelit(Abb. 2.224(b))
4.02 1.0 1.0 Cubic I Centre 0.1 1.9 0.1 Bartknecht data(Abb.2.219)
4.02 1.0 1.0 Cubic I Centre 0.1 2.0 0.1 Bartkneclit(Abb. 2.224(b))
5.00 1.0 4.8 Sphere I Centre 0.2 2 0.095 Donat
5.00 1.0 4.8 Sphere I Centre ni C. Donat(l 973)
238
5.00 1.0 4.8 Sphere I Centre 0.52 2.52 1 0.065 C.Donat(I 973)
5.00 1.7 6.9 Sphere I Centre 0.52 1.02 0.308 1larris and Briscoc
5.00 1.7 6.9 Sphere I Centre 0.22 1.52 0.156 Harris and Briscoe
5.00 1.7 6.9 Sphere I Centre 0.22 0.52 0.139 Harris and Briscoe
5.00 1.7 6.9 Sphere I Centre 0.52 2.52 0.139 Harris and Briscoe
5.00 2.6 11.91 Rec 2.3 Centre 0.096 0.199 0.56 R.DeGood/K.Chatratlii: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Bottom 0.103 0.159 0.56 R.DeGood/K.Chatratlii: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.1 0.235 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Ree 2.3 Centre 0.1 0,314 0.56 R.DeGood/K. Chatrathi: Study of
factors affecting Pred.
I
5.00 2.6 11.91 Rec 2.3 Centre 0.086 0.185 0.56 R.DeGood7K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.086 0.311 0.56 R,DeGood/K. Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.086 0.385 0.56 R.DeGood/K. Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Bottom 0.086 1.007 0.56 R.DeGood/K.Chatrathi- Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.086 0.215 0.56 R.DeGood/K. Cliatrathi: Study of
Ifactors affecting Pred.
5.00 2.6 11.91 Ree 2.3 Centre 0.93 0.241 0.56 R.DeGood/K. Cliatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.9 0.172 0.56 R.DeGood/K. Chatrathi: Study of
I factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.9 0.16 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.9 0.144 0.56 1R.DeGood/K.Cliatratiii: Study of
-
239
5.00 2.6 11.91 Rec 2.3 Centre 0.103 0.262 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0,103 0.26 0.56 R.DeGood/K.Cjlatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.103 0.156 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Cen re 0.11 0.246 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
I
5.00 2.6 11.91 Rec 2.3 Centre 0.103 0.334 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0,103 0.325 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.103 0.29 0.56 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.103 0.441 0.56 R.DeGood/K.Chatrathi: Study of
I. factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.103 0.282 0.54 R.DeGood/K.Chatrathi: Study of
factors affecting Pred.
5.00 2.6 11.91 Rec 2.3 Centre 0.1 0.3 0.54 R.DeGood/K.Chatrathi: Study of
I I I factors affecting Pred.
4.45 30.4 61.4 Rec 0.363 Centre 0.4 0.7 0.58 P.F.Thome et at
4.3 35.0 67.4 Rec 1.14 Center 0.1 0.18 2.1 Solberg et al(1980)
4.3 35.0 67.4 Rec 1.14 Front 0.1 0.19 2.1 Solberg et a1(]980)
4.3 35.0 67.4 1 Rec 1.14 Rear 0.1 0.1 2.1 Solberg et al(1980)
4.3 35.0 67.4 Rec L 14 Front 0.1 0.3 1.6 Solberg et al(I 980)
4.3 35.0 67.4 Rec 1.14 Rear 0.1 0.18 1.6 Solberg et al(1980)
4.3 35.0 67.4 Rec 1.14 Center 0.1 1.05 1.1 ISolberg et al(1980)
4.3 35.0 67.4 Re 1,14 Front 0.1 0.71 1.1 Solberg et al(1980)
4.3 35.0 67.4 Rec 1.14 Rear 0.1 0.72 1.1 Solberg et al(I 980)
I vent
5.00 35.0 67.4 Rec 1.415 Centre 0.0 1.37 1 P.F,Thome et at
4.3 35.0 67.4 Rec 1.14 Center 0.1 1.55 0.53 Solberg et al(I 980)
4.3 35.0 67.4 Rec 1.14 Front 0.1 1.68 0.53 Solberg et a](1980)
6.0 40.0 73.7 Rec 'I Center 0.0137 0.0299 8.48 Bromma(Sweden)
6.0 40.0 73.7 Rec 'I Center 0,0196 0.0295 8.48 Bromma(Sweden)
5.0 60.0 74.2 Sphere I Centre 0.5 0.5 5.22 C.Donat(l 973)
5.0 60.0 74.2 Sphere I Centre 0.1 0.3 4.8 C.Donat(1973)
5.0 60.0 74.2 Sphere I Centre 0.2 0.5 3.76 C.Donat(I 973)
5.0 60.0 74.2 Sphere I Centre 0.1 0.5 3 IC.Donat(1973)
5.0 60.0 74.2 Sphere I Centre 0.5 1 2.22 C.Donat(I 973)
5.0 60.0 74.2 Sphere I Centre 0.1 1 1.2 C.Donat(I 973)
5.2 0.17 1.() Rec 0,3 Center 0.03 6 0.069 0.oý Z; iiosli
5.2 0.17 1.9 Rec 0.3, Center 0,0-36 0.0481 0A0 zalos11
1 1 1
5.2 0.17 1.9 Rec 0.3 Center 0.036 0.0552 0.10 Allosh
5.0 81.00 117.9 Rec 1 Center 0.014 0.049 18.72 flowards and Karabinis.
5.0 91.00 117.9 Rec 2 Cellter 0.0252 0.0462 18.72 Howards and Karahinis.
5.0 81.00 117.9 Rec 3 Cellter 0.0735 0.1008 18.72 1lowards and Karahinis.
5.0 81.00 117,9 Rec 4 Center 0.0351 0.0469 18.72 Howards and Karahinis.
5.0 91.00 117.9 Rec 5 Celltel, 0.0 16 0.035 19.72 1lowards and Karabinis,
_3
5.0 81.00 117.9 Rec 6 ('eiltet- 0.0136 0.043 18.72 1Imards and Karahinis-
5.0 81.00 117.9 Rec 7 Center 0.0189 0.0679 18.72 1lowards and Karabinis.
5.0 81.00 117.9 Rec 8 Center 0.014 0.021 18.72 1lowards and Karahinis.
(Vol ýýo) (III on) sliallu ratio (barg) (bill -g) (III)
10 1.47 7.76 Rec 4.48 Center 0.0005 0.01 0.84 Cubbage& Simmonds
15 1.47 7.76 Rec 4.48 Center 0.0005 0.029 0.84 Cubbage& Simmonds
20 1.47 7.76 Rec 4.48 Center 0.0005 0.069 0.84 Cubbage& Simmonds
25 1.47 7.76 Rec 4.48 Center 0.0005 0.096 0.84 Cubbage& Simmonds
30 1.47 7.76 Rec 4.48 Center 0.0005 0.041 0.84 Cubbage& Simmonds
35 1.47 7.76 Rec 4.48 Center 0.0005 0.022 0.84 Cubbage& Simmonds
25 1.47 7.76 Rec 4.48 Center 0.0008 0.069 0.84 Cubbage& Simmonds
25 1.47 7.76 Rec 4.48 Center 0.001 0.131 0.84 Cubbage& Simmonds
25 1.47 7.76 Rec 4.48 Center 0.003 0.248 0.84 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0012 0.076 1.24 Cubbage&, Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0013 0.083 1.24 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0026 0.141 1.24 Cubbage& Simmonds
25 2.8 11.92 , Rec 1.3 Center 0.0005 1 0.041 1.14 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0014 0.086 1.14 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.002 0.124 1.14 Cubbage& Simmonds
25 2.8 11.92 Rec 1,3 Center 0,0026 0.15 1.14 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0006 0.076 0.84 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.001 0,103 0.84 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.002 0.14 0.84 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0006 0.034 0.62 Cubbage& Simmonds
25 2.8 11.92 Rec 1,3 Center 0.0007 0.09 0.62 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0014 0.17 0.62 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0007 0.103 0.47 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.001 0.124 0.47 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.001 0.131 0.47 Cubbage& Simmonds
25 2.8 11.92 Rec 1.3 Center 0.0012 0.152 0.47 Cubbage&, Simmonds
I lydrogen/air
Vesse I -T-54 Centre 4.8 11-1 0.011 Cousin A, Cottons. Initial 1' 4S
40 0.01 0.004
psi (3.1 barg)
0.03 0.004 Vessel -ý-.54 -77ct-itrc 11.7 22.1 0.021 Cousill & Cottoll..". Initial 11 45
40
psh-,(.). I barvo
40 0.03 0.004 Vessel 1.54 Centre 27.6 0.021 COLISill & COU011S.llliliýd P -15
psig barg)
Vessel 1.54 C'ciltrc 4.8 11.7 O.Osý Cousin & Cottons. Initial 1) 4S
40 0.03 0.004
psig (3.1 barg)
40 0.03 0.004 Vessel -Fý4- -Z7- I1 10.5 0.055 Cousin & Cottons. hillial P 45
elltrC
psig (3.1 barg)
1 1
40 0.0.1 0.004 Vessel Lý4 Centre 22.8 2S.4 0.055 Cousin & Cottons. Initial P 4S
40 0.085 1.134 Tank 2-3) (7entre 0 3.1 0.14 Cousin & Cottons. Initial P0
pSig
40 0.085 1.134 Tank 2.3 CCIIII-C 5.() 11 0.008 Cousill Cottons 1111tialP I.S
.
1)Sig
0.085 1.134 Tank 2.3 Centre 11.9 12.o 0.008 Cousin A-, ('011011s. inilial PIS
140
O.W 1.114 '1ank 2.3 Centre 1.8 TO 0.0101 Cousin & Cottoll". Initial PIS
-10
psig
40 0.085 1.134 '1ank 2.3 Centre 7 10.3 0.0161 COLISill & COttOllS. 11116illP 1.5
245
I psig
40 0,085 1.134 Tank 2.3 Centre 11 11.9 0.0161 Cousin & Cottons. Initial P= 15
psig
40 0.085 1.134 Tank 3.3 Centre 1.5 4.6 0.0417 Cousin & Cottons. Initial P= 15
psig
40 0.085 1.134 Tank 4.3 Centre 6.5 8.2 0.0417 Cousin & Cottons. Initial P= 15
I psig
40 0.085 1.134 Tank 5.3 Centre 10.6 11.8 0.0417 Cousin & Cottons. Initial P= 15
Ipsig
40 0.22 0.5 Drum 1.41 Centre 0 4.5 0.09 Cousin & Cottons. Initial P= 0
psig
40 0.22 0.5 Drum 1.41 Centre 0 2.8 0.23 Cousin & Cottons. Initial P= 0
psig
40 0.22 0.5 Drum 1.41 Centre 0 1.5 0.43 Cousin & Cottons. Initial P= 0
psig
40 0.4 5.8 Pipe 22.1 Centre 0 4 0.05 Cousin & Cottons. Initial P= 0
psig
40 0.4 5.8 Pipe 22.1 Centre 0 2.9 0.09 Cousin & Cottons. Initial P= 0
psig
40 0.4 5.8 Pipe 22.1 Centre 0 2.6 0.2 Cousin & Cottons. Initial P= 0
I psig
40 0.4 5.8 Pipe 22.1 Centre 0 1.3 0.69 Cousin &, Cottons. Initial P= 0
1 1 psig
NR 0.95 4.673 Cylinder 2 Centre 0.075 1.3 0.2 1Razuset a]
NR 0.95 4.673 Cylinder 2 Centre 0.135 0.4 0.3 Razuset at
20 1 6 Cubic I Centre 0.1 0.5 0.34 Bartknecht(Abb.2.224(b))
20 1 6 Cubic I Centre 0.1 2 0.159 Bartknecht(Abb.2.224(b))
20 1 6 Cubic I Centre 0.2 0.5 0.36 Bartknecht(Abb.2.224(b))
20 1 6 Cubic I Centre 0.2 2 0.16 Bartknecht(Abb.2,224(b))
20 1 61 Cubic I Centre 0.5 2 0.18 Bartknecht(Abb.2.224(b))
10 6.37 16.6 Sphere I Centre 0.1 0.3 0.64 R.K. Kumar et.a).() 989)
10 6.37 16.6 Sphere I Centre 0.1 0.9 0.2 R.K. Kurnar et.al,(l 989)
10 6.37 16.6 Sphere I Centre 0.1. 1.7 0.071 R.K. Kumar et.al,(I 989)
20 6.37 16.6 Sphere I Centre 0.1 3.4 -U4 R.K. Kumar et.al.(1989)
20 6.37 16.6 Sphere I Centre 0.1 4.2 0.2 R.K. Kumar et,al,(I 989)
20 6.37 16.6 Sphere I Centre 0.1 4.8 0.071 R.K. Kumar et.al.(] 989)
12 6.37 16.6 Sphere II Centre 0.1 1 0,64 R.K. Kumar et.al.(I 9 89)
12 637 16.6 Sphere I Cent e 0.1 2 0.2- R.K. Kumar et.al.(] 989)
12 Sphere I Centre 0.1 2.5 0.071 R.K. Kurnar et.al.(I 9 89)
246
14 6.37- 16.6 Sphere I Centre 0.1 1,6 0.64 R.K. Kumaret. al.(1989)
14 6.37' 16.6 Sphere I Centre 0.1 4.9 0.2 R.K. Kurnar et.al.(l 989)
14 6.37 16.6 Sphere I Centre 0.1 5.1 0.071 R.K. Kumar et.al.(1989)
16 6.37 16.6 Sphere I Centre 0.1 2.2 0.64 R.K. Kumar et.al.(1989)
16 6.37 16.6 Sphere I Centre 0.1 3.4 0,2 R.K. Kumar et.al.(l 989)
16 6.37 16.6 Sphere I Centre 0.1 3.6 0.071 R.K. Kurnar et.al.(1989)
18 6.37 16.6 Sphere I Centre 0.1 3.9 0.2 R.K. Kumar et.al.(] 989)
18 6,371 16.6 1 Sphere I "ntre 0.1 5.1 0.071 R.K. Kumar et.al.(l 989)
Acetone/air
3.0 0.008 -0.19 Sphere 0 1.17 0.00013 Proceedingof I st & 2nd Sym. On Comb.Inst.
3.0 0.008 0.19 Sphere 0 0.48 0.00051 Proceedingof I st & 2nd Sym. On Comb.Inst.
3.0 0.008 0.19 Sphere 0 0.14 0.002 Proceedingof I st & 2nd Sym. On Comb. Inst.
3.0 0,008 0.19 Sphere 0 0.034 0.0082 Proceedingof I st & 2nd Sym. On Comb.Inst.
3.0 0.008 0.19 Sphere 0 0.034 0.018 Proceedingof I st & 2nd Sym. On Comb.Inst.
3.0 0.008 0.19 Sphere 0 0.034 0.031 Proceedingof I st & 2nd Sym. On Comb. Inst.
3.0 0.008 0.19 Sphere 0 0.034 0.0016 Proceedingof I st & 2nd Sym. On Comb, Inst.
3.0 0.008 0.19 Sphere 0 0.034 0.0036 Proceedingof I st & 2nd Sym. On Comb. Inst.
4.0 0.008 0.19 Sphere 0 2.76 0.00013 Proceedingof Ist & 2nd Sym. Oil Comb.Inst.
4.0 0.008 0.19 Sphere 0 1.38 0.00051 Proceedingof I st & 2nd Sym. On Comb.Inst.
4.0 0.008 0.19 Sphere 0 0.48 0.002 Proceedingof I st & 2nd Sym. On Comb.Inst.
4.0 1 0.19
0.008 Sphere 0 0.17 0.0082 Proceedingof I st & 2nd Sym. On Comb.Inst.
4.0 8 0.19
0.008 Sphere 0,034 0.018 Proceedingof I st & 2nd Syrn. On Comb.Inst.
4.0 8 0.19
0.008 Sphere 10 0.034 0.031 Proceedingof I st & 2nd Sym. On CombAnst.
4.0 0.008 0.19 Sphere 0 0.17 0.0016 Proceedingof I st & 2nd Sym. On Comb.I nst.
4.0 0.008 0.19 Sphere 0 0.034 0.0036 Proceedingof I st & 2nd Sym. On Comb.Inst.
5.0 0.008 0.19 Sphere- -0 4.48 0.00013 Proceedingof I st & 2nd Sym. On Comb.Inst.
5.0 8 0.19
0.008 Sphere 0 3.1 0.00051 Proceeding of I st & 2nd Sym. On Comb.Inst.
5.0 0.008 0.19 Sphere 0 0.83 0.002 Proceedingof Ist & 2nd Sym. On Cornb.lnst.
5.0 0.008 0.19 Sphere 0 0.28 0.0082 Proceeding of I st & 2nd Sym, On Comb.Inst.
5.0 0.008 0.19 Sphere 0 0.1 0.018 Proceeding of I st & 2nd Synn,On Comb.Inst.
5.0 0.008 0.19 Sphere 0 0.0344 0.031 Procecclingof I st & 2nd Syrn. On Comb-Inst.
5.0 0,008 0.19 Sphere 0 0.34 0.0016 Proceeding of I st & 2nd Sym. On Comb.Inst.
-
5.0 0.008 ý 0.19 1 Sphere 0 0.1 0.05367 Proceedingof I st & 2nd Sym. On Comb.Inst.
247
5.5 0.008 0.19 Sphere 0 3.38 0.00051 Proceedingof I st & 2nd Sym. On Comb.1nst.
5.5 0.008 0.19 Sphere 0 0.97 0.002 Proceedingof I st & 2nd Sym. On Comb.Inst.
5.5 0.008 0.19 Sphere 0 0.31 0.0082 Proceedingof I st 2nd Sym. On CornbAnst,
5.5 10.008 0.19 1 Sphere 0 0.14 0-018 Proceedingof I st 2nd Sym. On Comb.1nst.
.
5.5 0.008 0.19 Sphere 0 0.034 0.031 Proceedingof Ist & 2nd Sym. On CombAnst.
5.5 0.008 0.19 Sphere 0 0.41 0.0016 ProceedingofI st & 2nd Sym. On Comb.1nst.
5.5 0,008 0.19 Sphere 0 0.14 0.0036 Proceedingof I st &, 2nd Syrn. On Comb.] nst.
4.907 0.022 0.374 Cylinder 0.100 0.680 0.680 Molkov et al: Venting of Deflagration: IIC-air
I and H2-air system
4.907 2.000 7.677 Cylinder 0.150 2.270 2.270 Molkov et a]: Venting of Deflagration: HC-air
Tablc A. 2 List ol'publislied experimental data t'()]-vellwd gas explosion mth a pi-C.
"Clicc
Pillic.
()I'CILICI
A, P" Prod Ignition Prod Souces
Fuel L D Ud (M)
(M) (M) (M) (bar) (bar) without
vent (bar)
I'l-opillic ill] 0. () 10 0.0 1 00, (Io 1,01 lo I Ild \oI Polil"\ mid I -\cl
111casillud
4.0(ýi, 0.6 0.021 28.6 0.000140 0.00366 1.0 1 1-11d Not POlll/\ and LCN,
cl-
-1.18
measurccl
0.6 0.036 16.7 0.001017 0.003166 1.01 2.2 8 End Not Pollizý and Leyer
mcasurcd
1.1 0.016 68.8 0.000201 0.00366 1.01 2.81 Fnd Not Poili/ý and kcýer
111casill-cd
1.1 0.021 52.4 0.000146 0.00366 1.01 2.40 I.,nd Not Polli/) and Leyel.
measured
1.1 0.036 30.6 0.001017 0.00366 1.01 2,93 I-Ind Not Poliizý and Leyer
measured
2.6 0.016 162.5 0.000201 0.00366 1.01 2.93 Fild Not Polliz), and Leyer
measured
2.6 0.021 123.8 0.00031460.00366 1.01 2.56 F,nd Not Volliz) and Lever
measured
2.6 0.036 72.2 0.0010170.00366 1.01 2.4 Fnd Not Ponizýand Lcyer
111casill-cd
1.7 0.036 47.2 0.001017 0.00366 1.01 3.02 Central Not and Leyer
11011i/ý
nicasured
1.7 0.036 47.2 0.001017 0.00366 1.31 3.17 Centrill Not llollizý and Leyer
measured
1.7 0.036 47.2 0.001017 0.00366 1.92 3.67 Central Not Pollin and I, CNef,
illeasurcd
1.7 0.036 47.2 0,001017 0.00360 33 1 4.38 Cciltral Not Polii/ý and Leyel.
IIIC,tStll*CLI
1.7 0.036 47.2 0.001017 0.00360 1.01 1:11d Not 11onizyand I,cvef-
-1.77
1.7 0.036 47.2 0.001017 0.00306 1.33 2.99 1md Not Poni/y and I.c\,cr
111castil-cd
1.7 0.036 47.2 0.0010 17 0.00166 1.84 2.82 1 Ild Not Pollin and 1,c\ cl
249
Propane/air
5.0% 1 0.8446 1.2 0.56 2.6 1 11 1.2 Central 1.212 DeGood a
.
Chartrathi
2 0.8446 2.4 0.56 2.6 1 11 1.32 Central 1,212 DeGood a
.
Chartrath!
3 0.8446 3.6 0.56 2.6 1 11 1.4 Central 1.212 DeGood a
.
Chartrathi
3 0.8446 3.6 0.56 2.6 1 11 2.02 Bottorn 1.16 DeGood a
.
Chartrathi
Acetone/air
5.0% 25 0.5 50.0 0.19625 10 1.11 5.11 Central Not Molkov (1993)
measured
25 0.5 50.0 0.19625 10 1.06 3.81 Central Not MoIkov (1993)
measured
4 0.2 20.0 0.0314 2 1.16 5.31 Central Not Molkov (1993)
measured
10 0.2 50.0 0.0314 2 1.16 6.21 Central Not Molkov (1993)
measured
10 0.38 26.3 0.113354 2 IA 3.16 Central Not Molkov (1993)
measured
1.83 0.05 36.6 0.001963 0.027 1.21 6.01 Central 1.7 Molkov (1993)
2.35 0.05 47.0 0.001963 0.027 1.26 5.41 Central 1.7 Molkov (1993)
2.35 0,05 47.0 0.001963 0.027 1.26 4.51 Central Not Molkov (1993)
measured
2.35 0.05 47.0 0.001963 0.027 2.66 2.91 Central Not Molkov (1993)
measured
1.83 0.05 36.6 0.001963 0.027 2,43 5.41 Central Not Molkov (1993)
measured
Methane/air
9.5% 0.11 0.1129 1.0 0.01 0.00564 1.01 1.11 Not 1.092 McCann(1985)
available
0.3 0.1129 2.7 0.01 0,00564 1.01 1.2 Not 1.092 McCann(1985)
available
0.52 0.1129 4.6 0.01 0.00564 1.01 1.18 Not 1.092 McCann(1985)
available