Comparison of explosion parameters for Methane–Air mixture in different cylindrical flameproof enclosures

Flameproof enclosures having internal electrical components are generally used in classified hazardous areas such as underground coalmines, refineries and places where explosive gas atmosphere may be formed. Flameproof enclosure can withstand the pressure developed during an internal explosion of an explosive mixture due to electrical arc, spark or hot surface of internal electrical components. The internal electrical component of a flameproof enclosure can form ignition source and also work as an obstacle in the explosion wave propagation. The ignition source position and obstacle in a flameproof enclosure have significant effect on explosion pressure development and rate of explosion pressure rise. To study this effect three cylindrical flameproof enclosures with different diameters and heights are chosen to perform the experiment. The explosive mixture used for the experiment is stoichiometric composition of methane in air at normal atmospheric pressure and temperature.It is observed that the development of maximum explosion pressure (P_max) and maximum rate of explosion pressure rise (dp/dt)_ex in a cylindrical flameproof enclosure are influenced by the position of ignition source, presence of internal metal or non-metal obstacles (component). The severity index, K_G is also calculated for the cylindrical enclosures and found that it is influenced by position of ignition source as well as blockage ratios (BR) of the obstacles in the enclosures.     

Effects of significant damage of flame gap surfaces in flameproof electrical apparatus on flame gap efficiency

Electrical apparatus for use in the presence of explosive gas atmospheres has to be specially designed to prevent the apparatus from igniting the gas. Flameproof design is one of several options, and one requirement is then that any holes and slits in the enclosure wall be designed to prevent a possible gas explosion inside the enclosure from being transmitted to an explosive gas cloud outside it. Current standards (IEC) require that joint surfaces have a surface roughness of <6.3 μm. Any damaged joint surface has be restored to this quality. The present investigation has demonstrated that flame gap surfaces in flameproof electrical apparatuses can suffer considerable mechanical and corrosive damage before the flame gaps no longer function satisfactorily. In some cases very significant mechanical surface damage in fact improves the gap performance. This indicates that current high costs of repairing and replacing flameproof electrical apparatus in process plants offshore and onshore can be reduced considerably without any increase of the explosion risk.     

Dynamic response and effect of apertures on explosion parameters of flameproof apparatus

Explosion parameters for closed flameproof apparatus are changed when apertures like gap (e.g. push button) and porous structures (breathing element) are introduced on the cover or wall of the flameproof enclosures. Similarly, an interconnecting tube between two enclosures, results in significant change in explosion parameters. It is observed that the maximum explosion pressure, maximum rate of pressure rise and severity index are higher for enclosures with apertures on cover or body than that of enclosures without apertures. In case of two interconnected identical enclosures, the explosion parameters are increased in the secondary enclosure and higher than that of primary enclosure and also of isolated enclosure.     

Effects of obstacles and deposited coal dust on characteristics of premixed methane–air explosions in a long closed pipe

An experimental system including pressure transducer, electric spark ignition device, data acquisition and control unit was set up to investigate methane–air explosions in a horizontal pipe closed at both ends with or without the presence of obstacles and deposited coal dust. The experimental results show that explosion characteristics depended on the methane content, on the layout of obstacles, and on the deposited coal dust. Pressure fluctuation with a frequency of 150 Hz appeared in its crest when the methane content was close to the stoichiometric ratio (9.5% methane percentage by volume). The pressure rise rate increased locally when a single obstacle was mounted in the pipe, but it had little effect on the pressure peak. Repeated obstacles mounted in the pipe caused the pressure to rise sharply, and the mean maximum explosion overpressure increased with the increase of the obstacle’s number. The amplitude of pressure fluctuation was reduced when deposited coal dust was paved in the bottom of the pipe. However, when repeated obstacles were arranged inside, the maximum overpressures were higher with coal dust deposited than pure gas explosions.     

Modified porosity approach and laminar flamelet modelling for advanced simulation of accidental explosions

The production of oil and gas is an inherently hazardous task. Therefore it is crucial to provide reliable estimates of the risks involved. The major contributors to the risk level of an offshore installation, for example, arise from accidents involving explosion and fire. Computational Fluid Dynamics (CFD) can be a powerful tool to help with the calculation of accidental explosion scenarios. In this context, the present work suggests a novel implementation of a model based on a modified Porosity Distributed Resistance (MPDR) approach within an unstructured 3D Navier–Stokes solver. The model operates by representing parts of the filtered geometry from the original model through porosity values attributed to an unstructured tetrahedral mesh. Extra resistance terms are added in the momentum equation as well as extra sources of turbulence. Two extra sources of turbulence are modelled. The first of these is due to the shear layers of the non-resolved obstacles, whilst the second is due to the presence of wakes behind the non-resolved obstacles. Results for 2D and 3D test cases are compared against both experimental data and simulations with fully-resolved geometry and good agreement is observed.     

Effect of rusting and mechanical damage of gap surfaces on efficiency of flame gaps in flameproof electrical apparatus

Potentially incendiary electrical apparatus for use in the presence of explosive gas atmospheres have to be specially designed to prevent the apparatus from igniting the gas. Flameproof design is one of several options, and one requirement is then that any holes and gaps in the enclosure wall be designed to prevent a gas explosion inside the enclosure from being transmitted to an explosive gas cloud outside it. Current standards (IEC) require that flame gap surfaces have a surface roughness of <6.3 μm. Any damaged flame gap surface has to be restored to this quality. The present investigation has demonstrated that flame gap surfaces in flameproof electrical apparatuses can suffer considerable corrosive and mechanical damage without any reduction of gap performance. In some cases very significant mechanical surface damage in fact improves gap performance. Possible physical reasons for this are discussed. These findings indicate that current high costs of repairing and replacing flameproof electrical apparatus in process plants offshore and onshore can be significantly reduced without any increase of explosion risks.     

纵向风和阻塞比对隧道双火源温度分布的影响*

为研究城市公路隧道内纵向通风和障碍物对双火羽流行为的影响,采用数值模拟方法分析双火羽流纵向烟气温度变化规律。研究结果表明:随着风速的增加,顶棚下方最高温度不断降低,烟气逆流现象逐渐减弱至消失;随着阻塞比的增加,下游火源一直向下游倾斜,而上游火源逐渐由向下游倾斜转变为向上游倾斜;基于流体力学理论,引入阻塞比修正无障碍物时的弗劳德数,进而建立适用于隧道内有障碍物的双火羽流顶棚最高温升分段预测模型,研究结果可为隧道火灾防治提供基础数据和理论参考。     

Flame acceleration in tube explosions with up to three flat-bar obstacles with variable obstacle separation distance

The effect of obstacle separation distance on the severity of gas explosions has received little methodical study. It was the aim of this work to investigate the influence of obstacle spacing of up to three flat-bar obstacles. The tests were performed using methane-air (10% by vol.), in an elongated vented cylindrical vessel 162 mm internal diameter with an overall length-to-diameter, L/D, of 27.7. The obstacles had either 2 or 4 flat-bars and presenting 20% blockage ratio to the flow path. The different number of flat-bars for the same blockage achieved a change of the obstacle scale which was also part of this investigation. The first two obstacles were kept at the established optimum spacing and only the spacing between the second and third obstacles was varied. The profiles of maximum flame speed and overpressure with separation distance were shown to agree with the cold flow turbulence profile determined in cold flows by other researchers. However, the present results showed that the maximum effect in explosions is experienced at 80 to 100 obstacle scales about 4 times further downstream than the position of maximum turbulence determined in the cold flow studies. Similar trends were observed for the flames speeds. In both cases the optimum spacing between the second and third obstacles corresponded to the same optimum spacing found for the first two obstacles demonstrating that the optimum separation distance does not change with number of obstacles. In planning the layout of new installations, the worst case separation distance needs to be avoided but incorporated when assessing the risk to existing set-ups. The results clearly demonstrate that high congestion in a given layout does not necessarily imply higher explosion severity as traditionally assumed. Less congested but optimally separated obstructions can lead to higher overpressures.     

An integrated model for vent area design of hydrocarbon-air mixture explosion inside cubic enclosures with obstacles

This study aims to develop an integrated model - NFPA-68-BRANN model, which can be used to calculate the vent areas of cubic enclosures with obstacles. Seven experiments regarding vented explosion inside the obstructed enclosure are reviewed and applied to check the accuracy of two existing standards, i.e. the NFPA-68 2018 and the BS EN 14994:2007. Accordingly, the parameters to describe the flame development in the NFPA-68 2018 are amended by adopting the Bauwens model. Bayesian Regularization Artificial Neuron Network (BRANN) model presenting the non-linear relationship between the turbulent flame enhancement factor X and its affecting factors is subsequently developed. Eventually, the NFPA-68-BRANN model is generated by incorporating the BRANN model into the modified NFAP-68 2018. The accuracy of the NFPA-68-BRANN model is validated by using a series of the New Baker Test data.     

Solution adaptive CFD simulation of premixed flame propagation over various solid obstructions

This paper presents the results of a number of calculations carried out in order to simulate combustion past obstacles of different shape and blockage ratio. The obstacle shapes considered are circles, squares, triangles and flat plates. Two-dimensional simulations are carried out with the McNEWT code. The code solves the reacting flow field with a laminar flamelet model on an unstructured mesh. Adaptive mesh refinement is applied so that the flame front is accompanied by mesh refinement throughout the calculation domain. A transition from laminar to turbulent combustion induced by passage past the obstacle is seen in the simulations. Evidence for the transition is found in the change in flame shape, flame speed and pressure. The simulations are compared with experimental data and there is good agreement between experiment and simulation.     

The acceleration of flames in tube explosions with two obstacles as a function of the obstacle separation distance

The separation distance (or pitch) between two successive obstacles or rows of obstacles is an important parameter in the acceleration of flame propagation and increase in explosion severity. Whilst this is generally recognised, it has received little specific attention by investigators. In this work a vented cylindrical vessel 162 mm in diameter 4.5 m long was used to study the effect of separation distance of two low blockage (30%) obstacles. The set up was demonstrated to produce overpressure through the fast flame speeds generated (i.e. in a similar mechanism to vapour cloud explosions). A worst case separation distance was found to be 1.75 m which produced close to 3 bar overpressure and a flame speed of about 500 m/s. These values were of the order of twice the overpressure and flame speed with a double obstacle separated 2.75 m (83 characteristic obstacle length scales) apart. The profile of effects with separation distance was shown to agree with the cold flow turbulence profile determined in cold flows by other researchers. However, the present results showed that the maximum effect in explosions is experienced further downstream than the position of maximum turbulence determined in the cold flow studies. It is suggested that this may be due to the convection of the turbulence profile by the propagating flame. The present results would suggest that in many previous studies of repeated obstacles the separation distance investigated might not have included the worst case set up, and therefore existing explosion protection guidelines may not be derived from worst case scenarios.     

Thermal evaluation of junction and connection boxes in explosion protection

To reliably avoid potential ignition sources and thus ignition of the potentially explosive atmosphere in junction and connection boxes of type of protection Increased Safety ‘e', the self-heating shall not exceed a specified level depending on the temperature class. For the conformity assessment of such products, complex thermal tests are necessary due to the great variety of mounting types and arrangements of terminal blocks in the enclosure, depending on the enclosure size. These tests are very time-consuming for the manufacturer of junction and connection boxes or for the testing laboratory and are therefore associated with considerable costs. To reduce this effort and to ensure a uniform assessment in the conformity evaluation by a certification body, it is therefore essential to create fitting charts depending on the enclosure size.This work introduces a calculation tool by means of systematic investigations on different enclosure sizes and with different assemblies, which enables the calculation of fitting charts with justifiable effort. For this purpose, the maximum temperature in the enclosure was determined as a function of the current and the assembly at 20 ± 5 °C ambient temperature. From this, electrical and enclosure-specific constants such as the maximum permissible current per conductor and a conductor specific factor were determined and combined with an exponential dependence of the power dissipation. It is shown that this relationship is valid for an overtemperature of 40 K for compliance with temperature class T6 up to an ambient temperature of +40 °C. Finally, to verify the reliability of the calculation tool, the results are compared with the enclosure-specific rated value of the maximum power dissipation according to IEC 60079-7, 5.7.     

Crowd dynamics discrete element multi-circle model

A Discrete Element Method (DEM) technique for modelling crowd dynamics is developed. Each person is represented by three overlapping circles with a position, orientation and velocity in 2D. Contact forces between elements are modelled as well as psychological forces and motive forces and moments. Motion is then modelled in a Newtonian manner with a numerical integration time-stepping scheme. The model is tested on a single enclosure entry scenario where some model parameters are scaled. Motion is generally realistic, although areas for improvement are identified. Flow and contact forces are monitored under different scenarios. The model is then used on a multi-enclosure entry scenario. It is shown that stewarding significantly reduces contact forces and hence risk of injury. The potential for further applications is demonstrated on hypothetical scenarios on the London Underground.It is concluded that the model is reasonably realistic for dense crowd flow scenarios, but more complex situations like a supermarket would require aspects of Artificial Intelligence in the model, which is a feasible development. The three-circle representation gives a reasonable model of the 2D geometry.     

The effect of barriers on reducing thermal heat fluxes from a hydrocarbon pool fire

The radiant heat flux from a pool fire is frequently calculated using the solid flame model, where the flame envelope is approximated as a stationary cylinder whose surface emits thermal radiation at a constant rate. Radiant heat flux calculations using the solid flame model assume the target to be at a given elevation, typically at ground level, and to have an unobstructed view of the fire. The presence of obstacles (e.g., walls, buildings, etc.) or terrain features that would create shaded areas and provide shielding of a target from the fire is typically neglected in these calculations: this is a conservative approach, but it is not accurate. This paper presents a methodology to utilize the solid flame model to calculate the heat flux to a target while taking into account the presence of an obstruction between the target and the fire. The shielded solid flame method can quantitatively account for the presence of obstacles as a passive mitigation measure and allows project developers or designers to optimize their facility layout to meet safety requirements. The methodology presented in this paper uses the same correlations found in currently used solid flame models (e.g., LNGFIRE3), therefore, it remains consistent with current regulatory requirements for LNG facilities in the U.S.     

Development of a Computationally Efficient Full Human Body Finite Element Model

Introduction: A simplified and computationally efficient human body finite element model is presented. The model complements the Global Human Body Models Consortium (GHBMC) detailed 50th percentile occupant (M50-O) by providing kinematic and kinetic data with a significantly reduced run time using the same body habitus.

Methods: The simplified occupant model (M50-OS) was developed using the same source geometry as the M50-O. Though some meshed components were preserved, the total element count was reduced by remeshing, homogenizing, or in some cases omitting structures that are explicitly contained in the M50-O. Bones are included as rigid bodies, with the exception of the ribs, which are deformable but were remeshed to a coarser element density than the M50-O. Material models for all deformable components were drawn from the biomechanics literature. Kinematic joints were implemented at major articulations (shoulder, elbow, wrist, hip, knee, and ankle) with moment vs. angle relationships from the literature included for the knee and ankle. The brain of the detailed model was inserted within the skull of the simplified model, and kinematics and strain patterns are compared.

Results: The M50-OS model has 11 contacts and 354,000 elements; in contrast, the M50-O model has 447 contacts and 2.2 million elements. The model can be repositioned without requiring simulation. Thirteen validation and robustness simulations were completed. This included denuded rib compression at 7 discrete sites, 5 rigid body impacts, and one sled simulation. Denuded tests showed a good match to the experimental data of force vs. deflection slopes. The frontal rigid chest impact simulation produced a peak force and deflection within the corridor of 4.63 kN and 31.2%, respectively. Similar results vs. experimental data (peak forces of 5.19 and 8.71 kN) were found for an abdominal bar impact and lateral sled test, respectively. A lateral plate impact at 12 m/s exhibited a peak of roughly 20 kN (due to stiff foam used around the shoulder) but a more biofidelic response immediately afterward, plateauing at 9 kN at 12 ms. Results from a frontal sled simulation showed that reaction forces and kinematic trends matched experimental results well. The robustness test demonstrated that peak femur loads were nearly identical to the M50-O model. Use of the detailed model brain within the simplified model demonstrated a paradigm for using the M50-OS to leverage aspects of the M50-O. Strain patterns for the 2 models showed consistent patterns but greater strains in the detailed model, with deviations thought to be the result of slightly different kinematics between models. The M50-OS with the deformable skull and brain exhibited a run time 4.75 faster than the M50-O on the same hardware.

Conclusions: The simplified GHBMC model is intended to complement rather than replace the detailed M50-O model. It exhibited, on average, a 35-fold reduction in run time for a set of rigid impacts. The model can be used in a modular fashion with the M50-O and more broadly can be used as a platform for parametric studies or studies focused on specific body regions.     

Numerical study of the effect of obstacles on the spontaneous ignition of high-pressure hydrogen

A numerical simulation of the spontaneous ignition of high-pressure hydrogen in a duct with two obstacles on the walls is conducted to explore the spontaneous ignition mechanisms. Two-dimensional rectangular ducts are adopted, and the Navier–Stokes equations with a detailed chemical kinetic mechanism are solved by using direct numerical simulations. In this study, we focus on the effects of the initial pressure of hydrogen and the position of the obstacles on the ignition mechanisms. Our results demonstrate that the presence of obstacles significantly changes the spontaneous ignition mechanisms producing three distinct ignition mechanisms. In addition, the position of the obstacles drastically changes the interaction of shock waves with the contact surface, and spontaneous ignition may take place at a relatively low pressure in some obstacle positions, which is attributed to the propagation direction and interaction timing of two reflected shock waves.     

Consequence modeling of a real rupture of toluene storage tank

Chemicals used in paint industry are generally comprised of toxic and inflammable substances; they may cause serious problems for humans and environment if safety conditions are neglected. In 2012, rupture of a toluene tank in a paint factory in Iran claimed lives of two workers near the tank. The present study is formulated in order to model the incident. The results showed the area prone to flash fire occurrence and extension. Furthermore, radiation was measured in terms of distance to the tank. A bund wall is also assumed around the vessel in order to determine its influence on evaporation in flash fire area, radiation, and maximum concentration of vapor cloud. The results provided important hints on minimum and required distance of electrical equipment or electrical enclosure (PLC) installation from vessels as well as suitable distance between the vessels. Finally, preventive recommendations were proposed to reduce the risk of potential accidents.     

Ammonia large scale atmospheric dispersion experiments in industrial configurations

A programme of large-scale experiments for atmospheric dispersion was carried out by INERIS over a period extending from December 1996 to April 1997. The objectives of the test campaign were to measure anhydrous ammonia concentrations in a range of few meters to 2 km from the release, in order to generate data to be used to improve 2-phase discharge and dispersion modelling.

The discharges were released from a 6-tonne storage tank of pressurised liquid ammonia and through a discharge device with an outlet diameter of 2 in. Fifteen trials were carried out with various release configurations corresponding to industrial situations (impinging jets on the ground and on a wall at various distances, release through a flange without seal…). The quantity of ammonia discharged from the liquid phase varied according to the tests, from 1.4 to 3.5 tons for durations between 7 and 14 min and, therefore, at flow rates between 2 and 4.5 kg/s. Approximately 200 sensors were settled downwind to measure ammonia concentrations and temperature in the plume. These tests showed that for discharges with identical flow rates the distances corresponding to the same concentration vary a lot according to the configurations. These distances tend to be reduced by the presence of obstacles or retention dikes that collected liquid ammonia. In the paper, the main experimental results are presented. In order to enable the comparisons with numerical predictions, more detailed information are given in [Bouet R. (1999). Ammoniac—Essais de dispersion atmosphérique à grande échelle. INERIS rapport, ref INERIS-DRA-RBo-1999-20410 (available at http://www.ineris.fr/recherches/recherches.htm).     

LES – DFSD modelling of vented hydrogen explosions in a small-scale combustion chamber

Accidental explosions are a plausible danger to the chemical process industries. In the event of a gas explosion, any obstacles placed within the path of the flame generate turbulence, which accelerates the transient flame and raises explosion overpressure, posing a safety hazard. This paper presents numerical studies using an in-house computational fluid dynamics (CFD) model for lean premixed hydrogen/air flame propagations with an equivalence ratio of 0.7. A laboratory-scale combustion chamber is used with repeated solid obstacles. The transient compressible large eddy simulation (LES) modelling technique combined with a dynamic flame surface density (DFSD) combustion model is used to carry out the numerical simulations in three-dimensional space. The study presented uses eight different baffle configurations with two solid obstructions, which have area blockage ratios of 0.24 and 0.5. The flame speed, maximum rate of pressure-rise as well as peak overpressure magnitude and timing are presented and discussed. Numerical results are validated against available published experimental data. It is concluded that, increasing the solid obstacle area blockage ratio and the number of consecutive baffles results in a raised maximum rate of pressure rise, higher peak explosion overpressure and faster flame propagation. Future model development would require more experimental data, probably in a more congested configuration.     

障碍物对油气-空气混合气体泄压爆炸火焰传播特性影响

为了研究障碍物对油气泄压爆炸火焰传播特性的影响规律,进行了不同数量障碍物工况下的对比实验,并利用纹影仪和高速摄影仪记录了火焰传播过程,针对障碍物对火焰形态、火焰锋面位置及火焰传播速度的影响规律进行了研究,结果表明:圆柱体障碍物会导致油气泄压爆炸火焰形态产生褶皱和弯曲变形,诱导层流火焰向湍流火焰转变,加速火焰的传播,对油气泄压爆炸火焰的初始传播形态有显著影响;随着障碍物数量的增多,火焰锋面传播距离点火端的最大距离增大,但到达最远距离的时间减少;障碍物能够增强火焰的传播速度,尤其对障碍物下游火焰影响最为显著,随着障碍物数量的增多,火焰传播的最大速度也随之增大,但达到最大火焰传播速度的时间却随之减少;障碍物的存在增大了油气泄压爆炸过程外部爆炸压力,并且随着障碍物数量的增多,外部爆炸压力峰值增长幅度增大。