Numerical study of dust lifting in a channel with vertical obstacles

In the paper, several results of numerical computation of multiphase flows in a channel with complex geometry are considered. The objective of the research was to study the dust lifting process from a layer behind a shock wave in a rectangular channel with vertical obstacles in the upper part of the tube. It is to be shown that that kind and also any sort of geometry may crucially change the whole phenomena of dust enhancement and of combustion. This is very important for safety in, for example, coal mines where channels are usually of more sophisticated structure than is usually assumed by most researchers.     

Experimental investigation on spontaneous ignition caused by pressurized hydrogen suddenly release into an S-shaped tube

An experimental investigation on the effects of continuous semicircular curved structure on spontaneous ignition during pressurized hydrogen suddenly release was conducted. An S-shaped tube with 700 mm in length and 10 mm in diameter was used in our experiments, and a straight tube with the same configuration was adopted for comparison. The results show that the continuously generated rarefaction waves and reflected shock waves make the pressure curves in the S-shaped tube more complicated. Meanwhile, the mean velocity and intensity of the leading shock wave undergo considerable attenuation when it propagates in the S-shaped structure. By comparing with the straight tube, the minimum critical pressure condition for spontaneous ignition in the S-shaped tube is slightly difficult to reach, but the difference is not huge. Nevertheless, the S-shaped structure can effectively promote hydrogen-air mixing and make combustion more intense. A secondary overpressure peak detected by the pressure transducer near the nozzle occurs in the spontaneous ignition cases and no such pressure increase is caught in the non-ignition cases. The transition from spontaneous combustion flame to a jet flame at the nozzle and the complete out-tube jet flame development process are captured and discussed.     

气固两相流下球阀磨损特性研究

旋塞球阀是钻柱内防喷系统中的关键设备,在气固两相流下球阀易受磨损而失效,并造成严重的井喷事故。为此,将计算流体动力学理论与冲蚀磨损理论相结合,运用FLUENT软件对球阀壁面在气固两相流下的磨损分布情况进行研究,并进一步分析了球阀结构参数对于球阀壁面磨损的影响规律。结果表明:当气固两相流流经球阀时,固体颗粒会与气流分离,并在壁面上产生三处磨损集中区;随着球阀开度的减小,球阀壁面磨损量会急剧增大,且阀球内通道壁面上的磨损集中区由块状逐步转化为带状,而球阀出口处的磨损集中区则会逐渐向下移动;球阀流道直径的减小也会使得壁面磨损量增加,但磨损集中区的分布基本不变。研究结果可为进一步优化球阀流道结构以减轻其壁面磨损提供理论依据。     

油水两相流弯管处安全分析

针对油水两相流经过弯管时的流向改变会导致流体速度和压力发生突变,造成发生静电事故和腐蚀事故风险上升 的不利影响,提出了RSM模型和Mixture模型相结合的安全分析方法。该方法对不同入口速度和含水率的油水两相流进行 数值模拟,并用Origin软件拟合了进口最大允许流速与管径及含水率的经验关系。结果表明,在含水率和入口速度一定 时,随着管径的增加,弯管处的最大速度呈现逐渐减小的趋势:当管径和入口速度一定时,随着含水率的增加,弯管处 的最大速度也逐渐减小。最大压力出现在弯管外拱壁处,最小压力出现在弯管内侧拱壁处。在实际生产中,增加弯管下 游直管段内侧壁的壁厚,可有效防止空化腐蚀所造成的危害;通过含水率来确定安全流速,可有效降低静电事故的风险 。     

泄爆容器中粉尘爆炸的试验研究与数值模拟

为了解泄爆容器中粉尘爆炸的发展过程,采用试验和数值模拟相结合的方法对玉米淀粉在圆柱形容器内的泄爆过程进行研究。数值模型采用欧拉–拉格朗日方法模拟粉尘爆炸的两相流问题,通过求解非稳态的湍流两相反应流守恒方程对试验进行二维仿真。试验和模拟结果表明,点火位置对爆炸发展过程有明显影响,点火位置离泄爆口越远,容器中的最大泄爆压力Pred,max越高。在粉尘爆炸的安全防护设计中,应把点火位置作为重要影响因素之一加以考虑。     

Detonation flows in aluminium particle gas suspensions,inhomogeneous in concentrations

A physical and mathematical model of the reduced kinetics is presented describing heterogeneous detonation in suspensions non-uniform in particle concentration. The model is based on the heterogeneous media approaches, semi-empirical laws of ignition and combustion, and data on the dependence of the detonation velocity on particle concentration. Formation of suboxides and incomplete combustion of aluminum are taken into account integrally. The dependence of the heat release of chemical reactions and the fraction of unburnt particles on the initial composition is determined from the solution of the stationary problem of the structure of the detonation wave. In the calculations of unsteady detonation flows, it is supposed to solve an additional equation for the spatial distribution of initial concentrations. The problems of initiation and development of cellular detonation in flat channels in suspensions of micron-sized aluminum particles are studied. Dependences of the cell size on particle concentration in uniform suspensions are determined. The flow patterns of cellular structures, the forms of the leading front, and the propagation velocities in channels with longitudinal or transversal gradients of particle concentration are analyzed.     

Benchmarking of two-phase flow through safety relief valves and pipes

This paper compares calculated results for two-phase flows through safety relief valves and pipes using the TPHEM, CCFLOW and RRERSP computer programs. These studies were conducted to locate errors in the programs as well as to further our understanding of how each program worked. For most low-to-moderate viscosity flow examples, and for the frozen flow examples, the program results agreed despite differences in the calculation methods. Current thinking is that the TPHEM computer program gets better results for high viscosity flows through safety relief valves (nozzles) and for pipe down flows. This is because this program successively iterates nozzle flows to achieve results consistent with the choke pressure, temperature and quality and successively iterates pipe flows to achieve results consistent with the pipe length. Safety relief valve flows from the RRERSP computer program are reduced by the inlet pipe non-recoverable pressure loss. This effect is significant for several cases.     

Numerical modelling of dust-layered detonation structure in a narrow tube

The problems of lifting and dispersing of a dust layer behind the propagating shock wave as well as ignition, combustion of coal particles and dust-layered detonation formation in a tube are numerically investigated. The layered detonation is formed at large distance from the place of the primary shock wave initiation (～100 diameters of the tube). The strong oblique transverse shocks caused by combustion zone were discovered. The acceleration of leading shock wave and dust-layered detonation formation are connected with increasing and intensification of combustion zone which strongly depends on arising system of the oblique waves due to the development of the dust layer instabilities and vice versa. In the applied model, the moving medium is treated as a two-phase, two-velocity and two-temperature continuum with mechanical and thermal interphase interaction. The numerical procedure is based on the finite-volume approach and is implemented for parallel computing. The results obtained are of interest for applications in predictive modelling of accidents in industrial systems with reactive dust.     

Experimental investigation on shock waves generated by pressurized gas release through a tube

An experimental investigation on the flow structures and the strength of shock waves generated by high-pressure gas release through a tube into air was conducted. The results demonstrated that a leading shock wave was generated in front of the compressed gas jet and the shock wave speed increased firstly, then decreased and finally kept constant with an increase of the propagation distance in the tube. The experimentally measured Mach numbers of shock waves were close to those calculated from the theory of ideal shock tube flow. After spouting out of the tube, the normal shock quickly developed into a hemispherical shape. The Mach disk was observed in the under-expanded jet. For high-pressure combustible gas release, the concept of theoretical critical pressure of ignition was introduced and several theoretical critical pressures of common gaseous fuels were obtained.     

Modeling the initial flame acceleration in an obstructed channel using large eddy simulation

The propagation and acceleration of a flame surface past obstructions in a closed square channel was investigated using large eddy simulation. The dynamic Smagorinsky–Lilly subgrid model and the Boger flame surface density combustion model were used. The geometry is essentially two-dimensional with fence-type obstacles distributed on the top and bottom surfaces, equally spaced along the channel length at the channel height. Flame propagation, however, is three dimensional as ignition occurs at a point at the center of the channel cross-section. The effect of obstacle blockage ratio on the development of the flame structure was investigated by varying the obstacle height. Three-dimensional cases were simulated from the initiation of a combustion kernel through spark ignition to the acceleration of the flame front at speeds up to 80 m/s. The transition from laminar flame propagation to turbulent flame propagation within the “thin reaction zone” regime was observed in the simulations. By analyzing the development of the three dimensional flame surface and unburned gas flow field, the formation of several flame structures observed experimentally are explained. Global quantities such as the total flame area and centerline flame velocity were ascertained and compared to the experimental data. High amplitude oscillations in the centerline flame velocity were found to occur from a combination of the unburned gas flow field and fluctuations in the volumetric burning rate.     

A numerical approach to investigate the maximum permissible nozzle diameter in explosion by hot turbulent jets

The ignition of a combustible environment by hot jets is a safety concern in many industries. In explosion protection concepts, for a protection of the type “flameproof enclosures” a maximum permissible gap is of major importance. In this work a numerical framework is described to investigate the ignition processes by a hot turbulent jet which flows out from such gaps. A Probability Density Function (PDF) method in conjunction with a reaction-diffusion manifold (REDIM) technique is used to model the turbulent reactive flow. In this paper the ignition of a stoichiometric mixture of hydrogen/air gas by a hot exhaust turbulent jet is examined. The impact of the nozzle diameter on the ignition delay time is investigated, too. The method is used to explore the maximum nozzle diameter for specific boundary conditions for which there is no ignition.     

高速气固两相流作用下的煤化工管材冲蚀磨损行为研究

为解决煤化工业中节流阀突扩口高速气固两相流对管壁材质的冲蚀磨损问题,利用基于激波管原理驱动的气固两相流冲蚀实验装置,试验研究冲击角度、温度对煤化工管材(10#、AISI304)的冲蚀磨损规律。研究结果表明:10#、AISI304管材的冲蚀率将随着冲击角度的增加而先增大后减小;室温下,10#、AISI 304钢的最大冲蚀率均出现在15°~30°区间;随着温度的升高,10#的最大冲蚀率出现在30°~45°区间,AISI304最大冲蚀率出现在30°。10#在30°,45°冲击角度下冲蚀磨损率会随温度上升显著上升,在15°冲击角度下冲蚀磨损率反而会随温度上升而下降。AISI 304在15°,30°,45°冲击角度下,冲蚀磨损率均会随温度上升而上升;在特定条件下,10#管材的冲蚀性能将优于AISI304。     

Study of the explosion parameters of vapor–liquid diethyl ether/air mixtures

The knowledge of the vapor–liquid two-phase diethyl ether (DEE)/air mixtures (mist) on the explosion parameters was an important basis of accident prevention. Two sets of vapor–liquid two-phase DEE/air mixtures of various concentrations were obtained with Sauter mean diameters of 12.89 and 22.90 μm. Experiments were conducted on vapor–liquid two-phase DEE/air mixtures of various concentrations at an ignition energy of 40.32 J and at an initial room temperature and pressure of 21 °C and 0.10 MPa, respectively. The effects of the concentration and particle size of DEE on the explosion pressure, the explosion temperature, and the lower and upper flammability limits were analyzed. Finally, a series of experiments was conducted on vapor–liquid two-phase DEE/air mixtures of various concentrations at various ignition energies. The minimum ignition energies were determined, and the results were discussed. The results were also compared against our previous work on the explosion characteristics of vapor–liquid two-phase n-hexane/air mixtures.     

Numerical simulation of flame acceleration and deflagration-to-detonation transition of ethylene in channels

Ethylene (C₂H₄) is a hydrocarbon fuel and widely used in chemical industry, however, ethylene is highly flammable and therefore presents a serious fire and explosion hazard. This work is initiated by addressing the hazard assessment of ethylene mixtures in different scale channels (d = 5 mm, 10 mm and 20 mm) from the aspect of flame acceleration (FA) and deflagration-to-detonation transition (DDT) by using large eddy simulation (LES) method coupled with the artificially thickened flame (ATF) approach. The fifth order local characteristics based weighted essentially non-oscillatory (WENO) conservative finite difference scheme is employed to solve the governing equations. The numerical results confirm that flame velocity increase rapidly at the beginning stage in three channels, and the flame acceleration rate is slower in the subsequent stage, afterwards, the flame velocity has an abrupt increase, and the onset of detonation occurs. Due to the fact that wall effect is significant in the narrow channel (e.g.,5 mm), especially in the ignition stage of the flame, flames have different shapes in wider channels (10 mm and 20 mm) and narrow channel (5 mm). Both the pressure and temperature profiles confirm DDT run-up distances are 0.251 m, 0.203 m and 0.161 m in 20 mm, 10 mm and 5 mm channels, respectively, which indicates that a shorter run-up distance is required in narrower channel. The cellular detonation structures for the ethylene-air mixture in different channels indicate that multi-headed detonation structures can be found in 20 mm channel, as the channel width decreases to 10 mm, detonation has a single-headed spinning structure, as the width is further reduced to 5 mm, only large longitudinal oscillation of the pressure can be observed.     

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.     

弯道明渠内危化品泄漏扩散的数值模拟

自由水面流动和危化品水体迁移扩散过程属有界性问题,利用两相流理论建立三维瞬态VOF动力学模型,在其基础上提出危化品水体迁移扩散的浓度压缩性微分方程,有限容积法离散方程,高分辨率Gamma离散格式适用于解决有界性问题,用Gamma格式离散VOF和危化品迁移扩散方程的对流项。对Chang实验的弯曲水槽的流动及危化品扩散进行数值模拟,模拟结果与实验结果较吻合。笔者建立的数学模型能准确地模拟水流三维特性突出的弯道水槽的流场和危化品水体迁移扩散过程。     

A simple,explicit formula for the critical pressure of homogenous two-phase nozzle flows

The critical pressure ratio of the homogeneous two-phase nozzle flow model known as the Omega method is expressed in function of the Omega Parameter as the exact numerical solution of a transcendental equation. A well fitting, easy to use, explicit approximation for flashing and non-flashing flows is presented here. The validation against the exact numerical solution proves that this new formula is better fitting than the other ones in the technical literature for both single and two-component flows.     

密闭管道瓦斯爆炸冲击波冲量及压力上升速率的实验研究

针对矿井瓦斯爆炸破坏模式主要在压力破坏和冲量破坏的研究,实验分析瓦斯在密闭管道发生爆炸时瓦斯浓度对冲击波冲量及压力上升速率的影响,利用管道中距离点火源不同位置的压力传感器测试了不同浓度瓦斯的爆炸压力,对冲击波冲量及压力上升速率进行分析,为防爆抑爆提供依据。研究结果显示:在管道中距离点火源的不同位置上,当浓度为9.5%时,瓦斯爆炸冲击波冲量及压力上升速率最大;由于超压衰减和传播距离的增加,在距离点火源4m和8m时压力冲量较大;在瓦斯浓度较低的范围内瓦斯爆炸时,其压力上升速率增长较快,而随着浓度的增加在较宽的浓度范围内,能较稳定地维持在高位值。     

燃烧管中烟煤粉粉尘云传播参数的试验研究

利用自行设计的长29.6 m,内径199 mm,配有特殊扬尘装置的大犁卧式燃烧爆炸管道试验系统,对弱点火条件下烟煤粉与空气两相悬浮流中的爆炸过程进行了试验研究,用压电传感器测量了管内各测点的压力信号,观测到快速爆燃的状态稳定,分析了爆燃波稳定传播机理.结果表明:在煤粉浓度为300g/m~3及弱点火条件下,悬浮烟煤粉粉尘云中爆燃波能够稳定传播,且稳态传播距离持续20 m以上,峰值超压和波速平均值分别约为70 kPa和430m/s.