首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 359 毫秒
1.
Using solid obstacles to accelerate the deflagration to detonation transition (DDT) process induces additional thrust loss, and fluidic obstacles can alleviate this problem to a certain extent. A detailed simulation is conducted to investigate the effects of multiple groups of fluidic obstacles on the flame acceleration and DDT process under different initial velocities and gas types. The results show that, initially, the propagation of reflected shock wave formed by jet impingement is opposite to the flame acceleration direction, thus increasing the initial jet velocity will hinder the flame acceleration. Later, the vortex structure and enhanced turbulence can promote flame acceleration. As the flame accelerates, the virtual blockage ratio of the fluidic obstacles decreases, and increasing initial jet velocity or using reactive jet gases both affect the virtual blockage ratio. Further, increasing initial jet velocity or using reactive jet gases can shorten the detonation initiation time and distance. Compared with solid obstacles, it is concluded that fluidic obstacles can achieve faster detonation initiation with a smaller blockage ratio. Overall, the detonation phenomena in this study are all triggered by hot spots formed by the interaction between reflected waves and distorted flame, but the formation of reflected waves varies.  相似文献   

2.
Numerical simulations were performed to study flame acceleration and deflagration-to-detonation transition (DDT) in hydrogen-air mixture in a channel with a two-dimensional array of cylindrical obstacles. A high-order numerical algorithm with adaptive mesh refinement was applied to solve reactive Navier-Stokes equations. The effect of obstacle layout was examined by considering three layouts at a fixed blockage ratio of 0.5: staggered, inline-concentrated, and inline-scattered. Three blockage ratios, 0.33, 0.5, and 0.67, were used for the case of staggered obstacles to explore the influence of blockage ratio. The results show that both obstacle layout and blockage ratio have significant effects on flame acceleration and DDT occurrence, although the basic mechanism of detonation initiation is consistent for all the cases involving shock focusing. In the staggered case, the head-on collisions of flame and pressure waves with obstacles greatly promote the growth of flame surface area and thus lead to the fastest flame acceleration and shortest detonation onset time. In the inline-concentrated case, flame propagates slower than that in the staggered case due to smaller flame surface area. However, compared to the zigzag path in the staggered case, the straight passages parallel to flame propagation direction in the inline-concentrated case are more conducive to producing strong shock focusing and thus result in the shortest DDT distance. In the inline-scattered case, the straight passage along the centerline of channel facilitates the early acceleration of flame, but it has the slowest flame propagation in lateral directions and thereby the longest DDT time and distance. For the staggered obstacles at different blockage ratios, flame acceleration rate increases with increasing blockage ratio. The occurrence of DDT is hindered by the most congested obstacles, because the shock focusing is insufficiently strong to initiate detonation after passing through the excessively narrow gaps.  相似文献   

3.
《Combustion and Flame》2007,148(1-2):4-47
This paper summarizes a 10-year theoretical and numerical effort to understand the deflagration-to-detonation transition (DDT). To simulate DDT from first principles, it is necessary to resolve the relevant scales ranging from the size of the system to the flame thickness, a range that can cover up to 12 orders of magnitude in real systems. This computational challenge resulted in the development of numerical algorithms for solving coupled partial and ordinary differential equations and a new method for adaptive mesh refinement to deal with multiscale phenomena. Insight into how, when, and where DDT occurs was obtained by analyzing a series of multidimensional numerical simulations of laboratory experiments designed to create a turbulent flame through a series of shock–flame interactions. The simulations showed that these interactions are important for creating the conditions in which DDT can occur. Flames enhance the strength of shocks passing through a turbulent flame brush and generate new shocks. In turn, shock interactions with flames create and drive the turbulence in flames. The turbulent flame itself does not undergo a transition, but it creates conditions in nearby unreacted material that lead to ignition centers, or “hot spots,” which can then produce a detonation through the Zeldovich gradient mechanism involving gradients of reactivity. Obstacles and boundary layers, through their interactions with shocks and flames, help to create environments in which hot spots can develop. Other scenarios producing reactivity gradients that can lead to detonations include flame–flame interactions, turbulent mixing of hot products with reactant gases, and direct shock ignition. Major unresolved questions concern the properties of nonequilibrium, shock-driven turbulence, stochastic properties of ignition events, and the possibility of unconfined DDT.  相似文献   

4.
Flame acceleration and deflagration-to-detonation transition (DDT) in a channel containing an array of staggered cylindrical obstacles and a stoichiometric hydrogen-air mixture were studied by solving the fully-compressible reactive Navier-Stokes equations using a high-order numerical algorithm and adaptive mesh refinement. Four different longitudinal spacings (ls) of the neighboring obstacle rows (i.e., ls = 15.28, 19.1, 25.4, and 38.2 mm, corresponding to 1.2, 1.5, 2 and 3 times of obstacle diameter, respectively) were used to examine the effect of obstacle spacing on flame acceleration and DDT. The results show that the main mechanisms of flame acceleration and transition to detonation in all the cases studied are consistent. While the flame acceleration is caused by the growth of flame surface area in the initial stage, it is governed by shock-flame interactions in the later stage when shock waves are generated. The focusing of strong shocks at flame front is responsible for the initiation of detonation. It was found that the flame propagation speed and the DDT run-up distance and time are highly dependent on ls. Specifically, the flame acceleration declines as ls increases, since a larger ls leads to less disturbance of flow by obstacles per unit channel length. For detonation initiation, both the run-up distance and time increase with the increase of ls. It is interesting to note that the DDT distance and time increase significantly as ls increases from 19.1 mm to 25.4 mm. This is related to the slowdown of the increase rate of energy release over a period before DDT occurs under large ls condition, because every time the flame passes over an obstacle row the shock-flame interaction is delayed and numerous isolated pockets of unburned gas material are formed.  相似文献   

5.
Combustion of hydrogen can take place in different modes such as laminar flames, slow and fast deflagrations and detonations. As these modes have widely varying propagation mechanisms, modeling the transition from one to the other presents a challenging task. This involves implementation of different sub-models and methods for turbulence-chemistry interaction, flame acceleration and shock propagation. In the present work, a unified numerical framework based on OpenFOAM has been evolved to simulate such phenomena with a specific emphasis on the Deflagration to Detonation Transition (DDT) in hydrogen-air mixtures. The approach is primarily based on the transport equation for the reaction progress variable. Different sub-models have been implemented to capture turbulence chemistry interaction and heat release due to autoignition. The choice of sub-models has been decided based on its applicability to lean hydrogen mixtures at high pressures and is relevant in the context of the present study. Simulations have been carried out in a two dimensional rectangular channel based on the GraVent experimental facility. Numerical results obtained from the simulations have been validated with the experimental data. Specific focus has been placed on identifying the flame propagation mechanisms in smooth and obstructed channels with stratified initial distribution. In a smooth channel with stratified distribution, it is observed that the flame surface area increases along the propagation direction, thereby enhancing the energy release rate and is identified to be the key parameter leading to strong flame acceleration. When obstacles are introduced, the increase in burning rate due to turbulence induced by the obstacles is partly negated by the hindrance to the unburned gases feeding the flame. The net effect of these competing factors leads to higher flame acceleration and propagation mechanism is identified to be in the fast deflagration regime. Further analysis shows that several pressure pulses and shock complexes are formed in the obstacle section. The ensuing decoupled shock-flame interaction augments the flame speed until the flame coalesces with a strong shock ahead of it and propagates as a single unit. At this point, a sharp increase in propagation speed is observed thus completing the DDT process. Subsequent propagation takes place at a uniform speed into the unburned mixture.  相似文献   

6.
The accidental release of hydrogen into enclosures can result in a flammable mixture with concentration gradients and possible deflagration-to-detonation transition (DDT). This numerical study aims to investigate the effect of obstacle spacing and mixture concentration on the DDT in a homogeneous and inhomogeneous hydrogen-air mixture. The paper focuses on the mechanisms behind the DDT in two mixtures with an average hydrogen concentration of 15% and 30%. Unlike the near-stoichiometric mixture, in the lean mixture, DDT only occurs in the inhomogeneous mixture. Depending on obstacle spacing, three different regimes of DDT were observed in the near-stoichiometric inhomogeneous mixture: i) Detonation was ignited when a strong Mach stem formed and propagated between the obstacles; ii) two explosion centers appeared when incident shock and Mach stem reflected from upper and lower obstacles, respectively; iii) Mach stem did not form but DDT occurred behind the flame front at the top of the obstacle.  相似文献   

7.
Computational Fluid Dynamics solvers are developed for explosion modelling and hazards analysis in Hydrogen air mixtures. The work is presented in two parts. These include firstly a numerical approach to simulate flame acceleration and deflagration to detonation transition (DDT) in hydrogen–air mixture and the second part presents comparisons between two approaches to detonation modelling. The detonation models are coded and the predictions in identical scenarios are compared. The DDT model which is presented here solves fully compressible, multidimensional, transient, reactive Navier–Stokes equations with a chemical reaction mechanism for different stages of flame propagation and acceleration from a laminar flame to a highly turbulent flame and subsequent transition from deflagration to detonation. The model has been used to simulate flame acceleration (FA) and DDT in a 2-D symmetric rectangular channel with 0.04 m height and 1 m length which is filled with obstacles. Comparison has been made between the predictions using a 21-step detailed chemistry as well as a single step reaction mechanism. The effect of initial temperature on the run-up distances to DDT has also been investigated.  相似文献   

8.
A numerical approach has been developed to simulate flame acceleration and deflagration to detonation transition in hydrogen-air mixture. Fully compressible, multidimensional, transient, reactive Navier–Stokes equations are solved with a chemical reaction mechanism which is tuned to simulate different stages of flame propagation and acceleration from a laminar flame to a turbulent flame and subsequent transition from deflagration to detonation. Since the numerical approach must simulate both deflagrations and detonations correctly, it is initially tested to verify the accuracy of the predicted flame temperature and velocity as well as detonation pressure, velocity and cell size. The model is then used to simulate flame acceleration (FA) and transition from deflagration to detonation (DDT) in a 2-D rectangular channel with 0.08 m height and 2 m length which is filled with obstacles to reproduce the experimental results of Teodorczyk et al.The simulations are carried out using two different initial ignition strengths to investigate the effects and the results are evaluated against the observations and measurements of Teodorczyk et al.  相似文献   

9.
The present study aims to test the capability of our newly developed density-based solver, ExplosionFoam, for flame acceleration (FA) and deflagration-to-detonation transition (DDT) in mixtures with concentration gradients which is of important safety concern. The solver is based on the open source computational fluid dynamics (CFD) platform OpenFOAM® and uses the hydrogen-air single-step chemistry and the corresponding transport coefficients developed by the authors. Numerical simulations have been conducted for the experimental set up of Ettner et al. [7], which involves flame acceleration and DDT in both homogeneous hydrogen-air mixture as well as an inhomogeneous mixture with concentration gradients in an obstucted channel. The predictions demonstrate good quantitative agreement with the experimental measurements in flame tip position, speed and pressure profiles. Qualitatively, the numerical simulations have reproduced well the flame acceleration and DDT phenomena observed in the experiment. The results have revealed that in the computed cases, DDT is induced by the interaction of the precursor inert shock wave with the wall close to high hydrogen concentration rather than with the obstacle. Some vortex pairs appear ahead of the flame due to the interaction between the obstacles and the gas flow caused by combustion-induced expansion, but they soon disappear after the flame passes through them. Hydrogen cannot be completely consumed especially in the fuel rich region. This is of additional safety concern as the unburned hydrogen can be potentially re-ignited once more fresh air is available in an accidental scenario, resulting in subsequent explosions.  相似文献   

10.
Deflagration-to-detonation transition of a flame ring circularly expanding in a 260 μm gap filled with stoichiometric ethylene/oxygen mixtures initially at atmospheric pressure and temperature has been experimentally visualized. The results show that DDT can occur under the influence of wall confinement even for an expanding flame. DDT could be observed at a distance as short as ~70 mm from the ignition spot, which corresponds to ~130 μs after the ignition spark voltage breakdown. Velocity overshoot of reaction front velocity exceeding Chapman–Jouguet velocity was characterized. Cell structures were observed on the reaction fronts after DDT occurred. The visualizations also showed that smooth circular flame developed right after ignition quickly evolved into wrinkled flame as the flame ring propagated outwards. Flame propagating velocity was accelerated from ~600 m/s to ~1000 m/s during the wrinkled flame stage. A series of local explosion on the flame ring was observed during the DDT process, and resulted in an abrupt surge on reaction front propagation velocity.  相似文献   

11.
Shock-initiated ignition is studied numerically for single step Arrhenius kinetics. This is relevant to hydrogen safety, because hydrogen detonates easily, and detonation in shocked mixture may occur in deflagration to detonation transition scenarios, due to shock reflections on obstacles subsequent to flame acceleration. Simulation of ignition behind a shock moving into combustible mixture is difficult because of the singular nature of the initial conditions. The solution method includes two components. First, space as an independent variable is replaced by the ratio space over time, and second, initial conditions at a small non-zero time are used, obtained in closed form from short time asymptotics. This way, the initial singularity is effectively removed and the early process is well resolved. This method was used to study how the leading shock strength, and the heat release, affects the shock-initiated ignition process. The Essentially Non-Oscillatory algorithm used captures ignition, hot spot growth, birth of a secondary shock and transition into a detonation. Results show that for weaker shock cases as well as for higher heat release the evolution is more rapid, that a secondary shock forms closer to the contact surface and quickly becomes a detonation wave.  相似文献   

12.
Ensuring hydrogen safety has become of great significance nowadays, whose leakage can possible result in the deflagration to detonation (DDT). The numerical study aims to explore the effect of solid obstacle distribution on the DDT in a homogeneous hydrogen-air mixture. Results show that the detonation initiation process can be classified into two types: i) local spherical detonation caused by the coupling of flame surface and high-pressure region; ii) detonation triggered by the interaction between the upper wall and multiple compression waves in front of the flame. Also, this study finds that the flame acceleration (FA) experiences two periodic “acceleration-deceleration” processes before the detonation initiation, and the initiation distance and time are the shortest when the obstacles are symmetrically distributed. Further, the higher the unilateral blockage ratio, the more unfavorable DDT occurs. The present results highlight the effect of different solid obstacle distribution patterns on the FA and DDT process.  相似文献   

13.
Flame acceleration was investigated in an obstructed, square-cross-section channel. Flame acceleration was promoted by an array of top and bottom surface mounted obstacles that were distributed along the entire channel length at an equal spacing corresponding to one channel height. This work is based on a previous investigation of the effects of blockage ratio on the early stage of flame acceleration. This study is focused on the later stage of flame acceleration when compression waves, and eventually a shock wave, form ahead of the flame. The objective of the study is to investigate the effect of obstacle blockage on the rate of flame acceleration and on the final quasi-steady flame-tip velocity. Schlieren photography was used to track the development of the shock-flame complex. It was determined that the interaction between the flame front and the reflected shock waves produced from contact of the lead shock wave with the channel top, channel bottom, and obstacle surfaces govern the late stage of flame acceleration process. The shock-flame interactions produce oscillations in the flame-tip velocity similar to that observed in the early stage of flame acceleration, but only much larger in magnitude. Eventually the flame achieves a globally quasi-steady velocity. For the lowest blockage obstacles, the velocity approaches the speed of sound of the combustion products. The final quasi-steady flame velocity was lower in tests with the higher obstacle blockage. In the quasi-steady propagation regime with the lowest blockage obstacles, burning pockets of gas extended only a few obstacles back from the flame-tip, whereas burning pockets were observed further back in tests with the higher obstacle blockage.  相似文献   

14.
Experiments on shock waves propagation, spontaneous ignition, and flame development during high-pressure hydrogen release through tubes with symmetrical obstacles (O1-1) and asymmetrical obstacles (O1-2) are conducted. The obstacle's side is triangular with a length of 4 mm, a height of 3.6 mm, and its width is 15 mm. In the experiments, a reflected shock wave generates and propagates both upstream and downstream when the leading shock wave encounters the obstacle. At the same burst pressure, the reflected shock wave intensity in tube O1-1 is significantly greater than that in tube O1-2. Moreover, the presence of obstacles in the tube can induce spontaneous ignition. The minimum burst pressures for spontaneous ignition for tubes O1-1 and O1-2 are 2.84 MPa and 3.28 MPa respectively, lower than that for the smooth tube. Furthermore, both the initial ignition position and ignition time are greatly advanced in obstruction tubes, mainly affected by obstacle positions and burst pressures. Finally, the flame separation process near the obstacle is observed. After passing the obstacle, the flames grow rapidly in radial and axial directions on the tube sidewalls. And at the same burst pressure, the flame convergence time in tube O1-2 is usually longer than that in tube O1-1.  相似文献   

15.
The effect of blockage ratio on the early phase of the flame acceleration process was investigated in an obstructed square cross-section channel. Flame acceleration was promoted by an array of top-and bottom-surface mounted obstacles that were distributed along the entire channel length at an equal spacing corresponding to one channel height. It was determined that flame acceleration is more pronounced for higher blockage obstacles during the initial stage of flame acceleration up to a flame velocity below the speed of sound of the reactants. The progression of the flame shape and flame area was determined by constructing a series of three-dimensional flame surface models using synchronized orthogonal schlieren images. A novel schlieren based photographic technique was used to visualize the unburned gas flow field ahead of the flame front. A small amount of helium gas is injected into the channel before ignition, and the evolution of the helium diluted unburned gas pocket is tracked simultaneously with the flame front. Using this technique the formation of a vortex downstream of each obstacle was observed. The size of the vortex increases with time until it reaches the channel wall and completely spans the distance between adjacent obstacles. A shear layer develops separating the core flow from the recirculation zone between the obstacles. The evolution of oscillations in centerline flame velocity is discussed in the context of the development of these flow structures in the unburned gas.  相似文献   

16.
By utilizing a newly designed constant volume combustion bomb (CVCB), turbulent flame combustion phenomena are investigated using hydrogen–air mixture under the initial pressures of 1 bar, 2 bar and 3 bar, including flame acceleration, turbulent flame propagation and flame–shock interaction with pressure oscillations. The results show that the process of flame acceleration through perforated plate can be characterized by three stages: laminar flame, jet flame and turbulent flame. Fast turbulent flame can generate a visible shock wave ahead of the flame front, which is reflected from the end wall of combustion chamber. Subsequently, the velocity of reflected shock wave declines gradually since it is affected by the compression wave formed by flame acceleration. In return, the propagation velocity of turbulent flame front is also influenced. The intense interaction between flame front and reflected shock can be captured by high-speed schlieren photography clearly under different initial pressures. The results show that the propagation velocity of turbulent flame rises with the increase of initial pressure, while the forward shock velocities show no apparent difference. On the other hand, the reflected shock wave decays faster under higher initial pressure conditions due to the faster flame propagation. Moreover, the influence of initial pressure on pressure oscillations is also analyzed comprehensively according to the experimental results.  相似文献   

17.
A new level-set algorithm is developed that distinguishes the interior from the exterior of complex geometry by parameterizing the correlation between the stereolithography (STL) information and grid points in order to embody three-dimensional boundary interfaces in the Cartesian coordinate system. The method is applied to the flame and shock wave interactions in hydrocarbon fuel blends including the transition to a detonation in various three-dimensional geometric complexities for predicting the maximum chamber pressure. The locations of the strong shock-flame interactions and the onset of the hot spot initiation in the regions of the strongest interaction are investigated. Moreover, the study is intended for providing safety assessments of the pressure confinements that include bends, steps, and narrow channels connecting cavities by utilizing the present algorithm that can identify the triggering effects of the transition to detonation.  相似文献   

18.
The whole evolution of flame propagation in a confined combustion chamber was firstly experimentally observed in a newly designed experimental apparatus equipped with a perforated plate. The effect of the flame-flow/acoustic/shock wave interaction on the flame propagation was studied. The experiment was conducted with a stoichiometric hydrogen-air mixture. According to the flame morphology and the flame tip velocity, the whole evolution of flame propagation in the experimental apparatus was classified into the following three stages: laminar flame, jet flame and turbulent flame. In the present work, different flame propagation modes were obtained in different conditions. Depending on the initial pressure, three different flame propagation modes were observed. At an initial pressure of 1 bar, the flame propagation after perforated plate was mainly controlled by the interactions of the flame and combustion-generated flow ahead of the flame front. As initial pressures went up to 3 bar and 5 bar, shock waves were clearly observed ahead of the flame, which played a significant role on the flame propagation. The flame decelerated sharply and even propagated backwards, induced by the flame-shock wave interactions. Depending on the intensity of the shock wave, the backward-propagation velocity was higher at 5 bar with a stronger shock wave. In addition, the pressure oscillation at different initial pressures was discussed.  相似文献   

19.
To investigate the effects of the geometry of downstream pipes on the shock ignition and the formation of the shock waves during high-pressure hydrogen sudden expansion, a series of bench-mark experiments were designed and high-pressure hydrogen were released into five types of pipes with different angles (60, 90, 120, 150 and 180°). It was found that the geometry of downstream pipes had a significant influence on the shock ignition of hydrogen. The incident shock wave would be reflected at the corner of the pipes with angles of 60, 90, 120 and 150°. The intensity of the reflected shock wave is higher if the angle is smaller. In addition, the average velocity of the leading incident shock wave would decrease when it passed the corner of the pipe. Using a pipe with smaller angle significantly increases the likelihood of shock ignition and lowers the minimal required burst pressure for shock ignition. The overpressure of the incident shock waves inside the exhaust chamber (for the cases with the angles of 60, 90, 120 and 150°) decreases sharply. There are three flame propagation behaviors inside the exhaust chamber: flame quenching, flame separation and no flame separation. The results of this study have implications concerning designs for storage safety of hydrogen energy and may help get better understanding of shock ignition mechanism of high pressure hydrogen and effect of pipeline geometry on ignition.  相似文献   

20.
The effect of nonspherical (i.e. cylindrical) bomb geometry on the evolution of outwardly propagating flames and the determination of laminar flame speeds using the conventional constant-pressure technique is investigated experimentally and theoretically. The cylindrical chamber boundary modifies the propagation rate through the interaction of the wall with the flow induced by thermal expansion across the flame (even with constant pressure), which leads to significant distortion of the flame surface for large flame radii. These departures from the unconfined case, especially the resulting nonzero burned gas velocities, can lead to significant errors in flame speeds calculated using the conventional assumptions, especially for large flame radii. For example, at a flame radius of 0.5 times the wall radius, the flame speed calculated neglecting confinement effects can be low by ∼15% (even with constant pressure).A methodology to estimate the effect of nonzero burned gas velocities on the measured flame speed in cylindrical chambers is presented. Modeling and experiments indicate that the effect of confinement can be neglected for flame radii less than 0.3 times the wall radius while still achieving acceptable accuracy (within 3%). The methodology is applied to correct the flame speed for nonzero burned gas speeds, in order to extend the range of flame radii useful for flame speed measurements. Under the proposed scaling, the burned gas speed can be well approximated as a function of only flame radius for a given chamber geometry - i.e. the correction function need only be determined once for an apparatus and then it can be used for any mixture. Results indicate that the flow correction can be used to extract flame speeds for flame radii up to 0.5 times the wall radius with somewhat larger, yet still acceptable uncertainties for the cases studied. Flow-corrected burning velocities are measured for hydrogen and syngas mixtures at atmospheric and elevated pressures. Flow-corrected flame speeds in the small cylindrical chamber used here agree well with previously reported flame speeds from large spherical chambers. Previous papers presenting burning velocities from cylindrical chambers report performing data analysis on flame radii less than 0.5 or 0.6 times the wall radius, where the flame speed calculated neglecting confinement effects may be low by ∼15 or 20%, respectively. For cylindrical chambers, data analysis should be restricted to flame radii less than 0.3 times the wall radius or a flow correction should be employed to account for the burned gas motions.With regard to the design of future vessels, larger vessels that minimize the flow aberrations for the same flame radius are preferred. Larger vessels maximize the relatively unaffected region of data allowing for a more straightforward approach to interpret the experimental data.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号