Numerical Simulation of Interaction between Large-Scale Congestion and Vent during the Natural Gas Explosion in a Kitchen

'e influence of large-scale congestion on a confined natural gas explosion in a typical Chinese kitchen was studied using the computational fluid dynamics technology. It was found that opening the explosion venting surface promotes the development of turbulence, flame propagation velocity, and multipeak overpressure in the explosion flow field. Large-scale congestion can significantly strengthen the influence of the explosion venting surface on the flow field; the congestion and the explosion venting surface have a synergistic effect on the explosion flow field. At the moment of gas explosion, the flow fields in each area of the kitchen exhibit different distribution characteristics. A flow field near small-scale congestion is more likely to produce greater turbulence, combustion rate, and flame speed. 'e obstruction effect of large-scale congestion perpendicular to the flame propagation direction is dominant. 'e indoor flame propagation speed and overpressure development speed increase and the peak combustion rate and indoor peak overpressure decrease with an increase in obstacle blockage. Increases in the large-scale volume congestion rate and volume blockage in the kitchen induce changes in the indoor flame propagation mode and increase the external explosion overpressure. 'is paper investigated the correlation behavior between large-scale congestion and vent surface in a typical Chinese civil kitchen during natural gas explosion process and provided important support for understanding the mechanism of congestion on gas explosion process and the distribution of explosion hazards in a kitchen.


Introduction
As high-quality green energy, natural gas is widely used in industrial production of raw materials and for cooking and heating. With the rapid development of the global economy, the demand for natural gas has increased significantly. However, with the aging of pipelines and the irregular use of gas facilities, gas explosion accidents occur more frequently, with significant casualties and economic losses. According to China News, there were 548 gas explosion accidents in 2020, resulting in 84 deaths and 670 injuries, of which 218 were indoor gas explosion accidents, mostly in kitchens [1]. Vented explosions induced by kitchen doors and windows, coupled with the turbulence of indoor cupboards, cabinets, stoves, and other items causing congestion, promote the development of indoor explosion flow fields and explosion hazards, increasing the complexity of the vented explosion process.
As a general representation of indoor congestion, volume blockage is widely used to study the effects of congestion on explosive flow fields. To theoretically explore the effect of congestion on the acceleration mechanism of gas explosion flames and the influence of explosion hazards, researchers have studied the influence of regularly arranged slender congestion objects with different obstacle blockages on vented explosions. Peter et al. [2] used numerical methods to study the effect of a fixed arrangement of vertical cylinders (with an outer diameter of 2 in, which is 0.0508 m) on the vented explosion process of propane gas and found that the indoor peak overpressure increased with an increase in obstacle volume blockage. Kuznetsov et al. [3] investigated the influence of vertical wooden beams with lengths of 1 m and cross-sections of 0.04 m × 0.04 m, 0.04 m × 0.05 m, and 0.05 m × 0.05 m in a regular arrangement on a hydrogen vented explosion in an explosion chamber with dimensions of 1 m × 0.96 m × 0.98 m.
e results showed that greater obstacle blockage produced a greater flame combustion speed and overpressure. Using numerical methods, Zhang et al. [4] vertically arranged slender columnar obstacles of different sizes in a room with dimensions of 4.6 m × 4.6 m × 3 m to study their influence on the vented explosion hazard characteristics of LPG-air mixtures. It was found that the explosion temperature of the LPG-air mixture changed little with obstacle blockage, but the explosion pressure was sufficient to damage brick and mortar walls with a thickness of 370 mm. Lian et al. [5] used numerical methods to study damage to the building structure caused by the explosion of methane, hydrogen, and liquefied petroleum gas (LPG) in the presence of small square columns evenly arranged in a device with dimensions of 4.6 m × 4.6 m × 3 m. It was found that as the obstacle blockage increased from 0.24 to 0.36, the maximum flight distance of explosive fragments generated by the explosion vent surface increased from 6.18 m to 18.78 m, with an expansion factor of 3.04. Congestion is an important factor affecting indoor explosion hazards. Tomlin et al. [6] studied the effect of cylinders with a length of 4.5 m and a diameter of 0.18 mm on the explosion of natural gas in a 9 m × 4.5 m × 4.5 m explosion chamber. It was found that the overpressure is caused by flame acceleration resulting from indoor congestion and the explosion venting surface. Under normal circumstances, when the size of the explosion vent is fixed, the indoor peak overpressure increases with increasing obstacle blockage. Some research has studied the influence of obstacle shape on gas explosions. Masri et al. [7][8][9] studied an explosion chamber with dimensions of 0.195 m × 0.195 m × 0.545 m. e effect of circular, triangular, and rectangular congestion cross-sections on a propane/ air explosion was studied to determine the influence of different obstacle blockages and congestion shapes on the flame shape and speed. For rectangular cross-sections, the volume blockage had a significant effect on flame acceleration; the effect was relatively small for triangular and circular cross-sections.
To determine the explosion hazard characteristics in actual accident scenarios, the influence of gas cylinder congestion on the vented explosion process was studied using actual conditions at a combustible gas cylinder storage site. Schiavetti et al. [10,11] analyzed the influence of congestion with different cylinder sizes on the hydrogen explosion process with explosive devices of 1.15 m 3 and 25 m 3 . It was found that congestion increases the maximum peak overpressure in the chamber and a larger volume blockage produces a greater flame area and combustion rate. Skjold et al. [12] conducted hydrogen explosion experiments in ISO 20GP containers using gas cylinders for congestion to study the effect of obstacle blockage on a hydrogen-vented explosion; the results showed that the explosion pressure increased with an increase in gas concentration and congestion. However, when the gas concentration increased from 12% to 15%, the maximum explosion pressure decreased with an increase in obstacle blockage. Chen et al. [13] conducted a hydrogen-air mixture deflagration experiment in a 27 m 3 square explosion chamber and studied the effects of gas concentration, explosion venting area, and gas cylinder congestion on vented explosions. It was found that both the peak overpressure and flame velocity increased with an increase in the volume blockage. Li et al. [14] used a 1 m 3 rectangular container to study the effect of gas cylinder congestion parallel to the direction of the explosion vent in a hydrogen-air vented explosion. It was found that when the volume blockage increased from 0 to 0.44, the indoor and outdoor peak overpressures gradually increased and the peak overpressure and peak overpressure arrival time gradually decreased. ese studies provide important support and a theoretical basis for understanding the influence mechanism of congestion on the vented gas explosion process and the distribution of explosion hazards. However, the existing research has usually considered only small-scale congestion in regular shapes and uniform arrangements, or the structure of the congestion was simplified, making it difficult to understand the actual accident process to accurately reveal the nature of the explosion hazard. e influence of kitchen congestion characteristics and layout on the explosion flow field is often ignored, especially for gas explosion accidents in residential kitchens. e existing research results are somewhat inadequate for guiding safety design for gas explosion prevention and control in residential kitchens. At the same time, gas explosion may also cause damage to the building structure. erefore, some scholars [15] have used robust design criteria to conduct related studies on structural damage in order to obtain a more robust seismic design and reduce explosion hazards.
In terms of technical standards, NFPA68 [16] in the United States provides a simplified model for the maximum vented explosion overpressure, which can be used in large spaces with a length-to-diameter ratio not greater than 5. However, the surface areas of indoor equipment and facilities are omitted when calculating the surface area of the surrounding body, and the enhancement effect of the builtin congestion on the flame is ignored. e European Union EN14994 [17] provides a calculation model for the maximum vented explosion overpressure suitable for explosion chambers with a length-to-diameter ratio not greater than 2 and a volume not greater than 1000 m 3 . e turbulence enhancement effect in the model is based mainly on the complexity factor assigned to the explosion chamber. e model focuses mainly on small-scale continuous obstacles in industrial plants; the four assignments of complexity factors in the model are not applicable to residential buildings. Although the "Code for Fire Protection Design of Buildings" in China [18] includes vented explosion requirements for industrial plants and warehouses, it does not include vented explosion safety design methods and technical requirements for residential buildings with typical large-scale congestion.
Based on the study of the explosive flow field in a residential kitchen by Yang et al. [19], this study uses computational fluid dynamics technology to further investigate 2 Advances in Civil Engineering the correlation between large-scale congestion and vent surfaces in typical Chinese kitchens during a natural gas explosion and determine the influence of volume blockage on indoor flow field parameters.

Numerical Models.
e computational fluid dynamics (CFD) technique is used to solve the transient flow field in confined space. e model mainly uses the finite volume method to solve a series of equations including the massconservation equation, momentum-conservation equation, and energy-conservation equation. With the help of the Cartesian tensor notation, it can be expressed as follows [20]. e mass-conservation equation is zρ zt e momentum-conservation equation is e energy-conservation equation is where x is the space coordinate, t is the time coordinate, ρ is the density, u is the velocity, p is the static pressure, and i and j are the coordinate directions. E is the specific internal energy; its expression is where C V is the constant-volume specific heat, T is the temperature, m fu is the mass fraction of the fuel, and H c is the heat of combustion. e turbulent-diffusion coefficient is Γ * � μ t /(σ) * , where (σ) * is the default value of the turbulent Prandtl constant. τ ij is the viscous stress tensor, and its expression is where the turbulence-viscosity coefficient μ t � C μ ρk 2 /ε, where k and ε are the turbulent kinetic energy and its dissipation rate, δ ij is the Kronecker delta, and the model constant C μ � 0.09 m 2 /s. e turbulence in the vented gas explosion process is described by the standard k-ε turbulence model. It consists of two conservation equations for k and ε, respectively: where C 1 and C 2 are constants (C 1 and C 2 � 1.44 and 1.79). A one-step reaction model is used to describe the combustion process of the fuel/air mixture. is is mathematically formulated as a conservation equation for the fuel mass fraction: e volume combustion rate R fu � C t ρ (s t 2 /Γ fu ) R min , where R min is the minimum of the mass fractions for fuel, oxygen, and the reaction product, and C t is the combustionmodel constant (C t � 40). e turbulence combustion velocity is expressed as where u t is the turbulence intensity, L t is the characteristic length of the turbulence, S l is the specific laminar combustion velocity, and v is the fluid dynamic viscosity. e effective laminar burning velocity S b can be expressed as where S l is the specific laminar combustion velocity, R f is the spherical flame radius, F s is the laminar flame acceleration coefficient, and the value is 0.15.

Verification of Grid Independence and Experimental
Validation of the Model. e numerical model uses an eightnode structural grid to discretize the space and subgrid technology to describe the effect of small-scale entities on the flow field [21]. e physical model and the scenario simulated in this study is similar to that in [22]. To verify the correctness of the model, explosion overpressures calculated in [22] were compared with the experimental results of Bauwens et al. [23]. e variation trend of the two was similar, which proved the applicability of the numerical method [24]. Meanwhile, the effect of the grid size has been verified also. e grid size (0.1 m × 0.1 m × 0.1 m) used in this study is same as those in [22].

Research Plan
Based on the "Design Code for Residential Buildings" [25] and the "National Building Standard Design Atlas--Residential Kitchen" (14J913-2) [26] in China, this study selected an apartment kitchen with a type II structural layout as the research object, as shown in Figure 1. e model kitchen size was set as 3.4 m (length) × 2.1 m (width) × 2.4 m (height). e layout of the objects in the room and the main simulation parameters are consistent with the description in [22], including the following: the window in Figure 1 was regarded as a square vent with dimensions of 0.8 m × 0.8 m; the lower edge was 1 m from the floor, and the opening pressure was 20 kPa. To accurately describe the impact of an external explosion on the indoor explosion flow field, the calculation domain was extended to five times the length of the room along the direction of the vent. e extended part Advances in Civil Engineering of the calculation domain was set as a free-flow boundary.
e calculated domain size was 20.4 m × 2.1 m × 2.4 m, and the total number of grids generated by a grid size of 0.1 m is 112,750. e kitchen was filled with premixed methane gas/ air with a volume concentration of 10% in a static state. e initial pressure and temperature of the environment in the calculation domain were set to 1.01325 × 102 kPa and 300 K, respectively. A spherical ignition source with a radius of 0.015 m was placed at the geometric center position, 0.05 m away from the back wall of the kitchen. A physical model of a type II kitchen is shown in Figure 2. To accurately describe the distribution and sudden changes in the explosion flow field characteristics, several measuring points were set up. Table 1 provides a detailed description of the measuring points.
Based on the type II layout of the kitchen, the core congestion characteristics in the explosion flow field induced by the explosion venting surface were investigated. e volume blockage was changed by changing the size of the cabinet and wall cabinet. e simulation conditions are presented in Table 2.

Correlation Mechanism of Vented Explosion Surface and
Large-scale Congestion in the Explosion Flow Field 4.1.1. Turbulence Characteristics. Figure 3 shows the turbulent kinetic energy-time curves at the central location in an airtight room, an empty constrained explosion room, and a typical congested room. e turbulent kinetic energy in the airtight room is maintained at a low level. e turbulent kinetic energy in the constrained room first increases and then decreases with time, but the turbulent kinetic energy in the empty room has a lower growth rate, reaching its peak at approximately 0.6 s. Before the explosion venting surface opens, the turbulent kinetic energy in the congested room is almost unaffected by congestion. e explosion increases rapidly after opening, reaching a peak value at approximately 0.33 s. e indoor peak turbulent kinetic energy with congestion increases by two orders of magnitude compared with the empty room. It is observed that turbulence is weakest in a closed room. e opening of the explosion vent promotes the gas in the room to accelerate outward axially, producing a relatively small turbulence. Although the largescale congestion has little effect on the turbulence, the opening of the explosion vent causes the turbulence of the flow field to increase; the explosion vent and congestion have a synergistic effect on the indoor turbulence. Figure 4 shows the peak turbulent kinetic energy-distance curves for the three working conditions. e turbulence energy in the closed room is weak. In the empty vent room, opening the explosion vent induces a weak turbulent flow, which is more obvious near the explosion vent. In the congested room, the coupling of congestion with different scale characteristics and the explosion vent induces fluctuations in the indoor turbulence. Small-scale obstacles such as range hoods promote turbulent flow; large-scale congestion objects such as refrigerators hinder gas flow and inhibit turbulent energy.

Combustion Rate and Flame.
e turbulence characteristics of the flow field affect the gas combustion rate and flame propagation speed. e variation trend of the indoor peak combustion rate along the axial direction of the room is shown in Figure 5. Figure 6 shows the flame cloud map on the z � 1.5 m slice for the three working conditions. When there is no congestion, the indoor peak combustion rate curve is almost the same; there is a significant difference near the explosion vent, similar to the turbulent kinetic energy change trend. With congestion, the indoor combustion rate increases earlier and almost linearly. It is also observed in Figure 6 that before the explosion venting surface was opened, the flame structure was hemispherical in the early stage of explosion, and the flame propagation speed was relatively low. With congestion, although the flame passes through congestion objects such as cupboards and cabinets, its forward shape does not seem to change significantly. A possible reason is that large-scale congestion objects are arranged on the edge of the chamber wall, and their arrangement direction is the same as the flame propagation direction; there is not much turbulence on the surface of the congestion. e shape of the flame changed significantly after it passed through the first sudden change in congestion size. In the middle of the flame, an elliptical area appeared where a cloud of combustible gas burned. With increasing time, the flame on the side of the room appeared in a depression area near the water heater with rapid gas burning, caused by the sudden change in the congestion cross-section and rapid turbulent combustion induced by the small-scale water heater. e variation trend of the flame propagation velocity along the axial direction of the room is shown in Figure 7.
e variation trend is similar to that of the peak turbulent kinetic energy, as shown in Figure 4, indicating a positive feedback mechanism between the forward at a relatively stable low speed. e action of the vent surface has little influence on the flame because the opening pressure of the vent surface is 20 kPa. With congestion, the change in flame propagation velocity depends mainly on the effect of congestion on the flow field flame propagation and turbulence in the flow field. When there is no congestion, the flame propagates [27]. e flame exhibits two stages of propagation with congestion. In the first stage, when the flame interacts with the congestion, the increased flame surface area promotes the flame propagation speed. In the second stage, the flame front and the vortex behind the congestion interact, inducing flame deceleration. is may be due to the unburned gas flowing in front of the flame, causing the convective flame tip to expand and decelerate [28,29].

4.1.3.
Overpressure. e overpressure-time curves of the central measurement points in different rooms are shown in Figure 8. e pressure development trend in the rooms was almost the same before the explosion venting surface was opened. An increase in pressure is caused by the expansion of combustion products [13]. Unburned gas produces turbulence around congestion [28], and the combustion rate is faster.
us, the pressure of a congested room increases        Advances in Civil Engineering faster, and the opening pressure of the explosion vent surface is reached earlier. After the vented explosion surface is opened, a large amount of gas is discharged from the vented explosion port, forming an overpressure peak P 1 . Because the closed room is in an ideal adiabatic environment, the overpressure increases monotonically with time, and there is no energy loss. e overpressure peak P 2 may be the result of the combined effect of the peak indoor maximum flame area and the external secondary explosion. Pang et al. [30] also had the same explanation about the relevant content. ey thought that when the combustion rate is relatively larger and the peak combustion rate changes rapidly, an external secondary explosion may have occurred. e secondary explosion wave can propagate back into the room and block the vent at the same time; therefore, a stronger external explosion peak appears in the room. Because the turbulence induced by congestion increases the flame area and increases the combustion rate in the room, peak P 2 is higher and its arrival time is earlier. e pressure increase in the empty room was relatively slow, the peak value was lower, and the peak value was reached later, indicating that congestion has a significant effect on the development of explosive overpressure, especially for the second peak. Figure 9 shows the peak overpressure-distance curves on the central axis in different rooms. e peak overpressure at each measuring point in the closed room was almost unchanged and much greater than in the other two conditions. e overpressure at each position in the congested room was significantly greater than in the empty room. In the vented explosion room, the peak overpressure near the explosion vent gradually decreased; the overpressure near the explosion vent in the congested room decreased faster than in the other rooms.

Turbulence Characteristics.
To study the dual effects of the distribution characteristics of indoor congestion on the explosion flow field, the typical congestion room (VB � 33.72%) was divided into five regions to analyze the explosion process, as shown in Figure 10. e turbulent kinetic energy of the intermediate measuring points on the central axis in different areas of the room changes with time, as shown in Figure 11. e indoor turbulent kinetic energy first increases and then decreases, but the peak value and arrival time in different regions are different due to the difference in congestion distribution. e turbulence kinetic energy in region I is the greatest, followed by region IV. e turbulence kinetic energy in other regions is lower, differing by a factor of more than 150 in peak turbulent kinetic energy. Region IV reaches the peak first and region V reaches the peak last. Figure 12 shows the peak turbulent kinetic energydistance curves along the axial direction of the room in different regions. Figure 13 shows the velocity vector distribution on slices at different heights. It is observed that region I is directly opposite the explosion vent and is most affected by it. Opening the explosion vent accelerates the discharge of indoor gas; the gas in this region has a higher flow velocity, causing the flow field turbulence to increase  further. e turbulence generated by the small-scale congestion (cooker hood) in region IV also increases the peak turbulent flow energy in region I. Regions II and III are greatly affected by the roof and the ground. In Figures 13(a) and 13(c), it is observed that the gas flow at the front and rear walls is obstructed; reflux occurs, hindering further development of turbulence. e cooker hood in region IV enhances the turbulent flow; the vortex generated by the gas Advances in Civil Engineering flow near the cooker hood is also observed in Figure 13(b). e turbulent flow in region V is hindered mainly by largescale congestion (refrigerator), and the turbulent flow is significantly weakened. e effect of congestion on turbulence is observed in two forms: promotion and hindrance. Small-scale congestion can generate stronger turbulent eddies and promote an increase in turbulent flow energy.
Large-scale congestion perpendicular to the direction of flame propagation affects turbulence; the obstructive effect is dominant. e arrangement of cupboards, cabinets, and other congestion along the flame propagation direction has little effect on indoor turbulence, which is consistent with the conclusions of Bjerketvedt et al. [31]. e difference in turbulence in different regions is also due to the influence of   ere is no congestion in region I, although the turbulence is the greatest, indicating that the influence of congestion on the flame is much stronger than on the volume of unburned gas [32]. Figure 14 shows the variation trend of the indoor peak combustion rate along the axial direction of the room. It is observed in the figure that the peak combustion rate in each region gradually increases with increasing distance, but the development trend is 14 Advances in Civil Engineering increase. Region V is mainly obstructed by congestion, manifested by the relatively slow increase of the combustion rate. e flame propagation speed along the length of the room is shown in Figure 15. e flame propagation speed first increases and then decreases with the distance. e increase in region II is the most obvious; the increase in region V is the slowest; region I has a secondary increase near the explosion vent. From analysis of the flame cloud maps on different height slices of the room, it is observed in Figure 16(a) that the flame reaches the vicinity of the cooker hood in approximately 290 ms at a vertical height of 2.1 m and is affected by small-scale congestion to produce a burning red ellipsoid area. It continues to increase, passing through multiple congestion objects; the flame velocity in region II has a large peak, but the flame propagation velocity near the explosion vent is blocked by the wall and rapidly decreases. Figure 16(c) shows the flame distribution at a vertical height of 0.4 m. e congestion in this region is mainly a large cabinet arranged along the flame propagation direction, which has a small shearing effect on the flame; the flame grows slowly in region V. e vertical height in Figure 16(b) is most affected by the explosion vent and is 1.4 m directly in front of its center. e turbulence near the vent causes a secondary increase in the flame velocity in region I. Figure 17 shows the overpressure-time curve of the indoor center measuring point. e trends of the overpressure time curve in each region are similar, but the P 1 and P 2 peaks on each curve and their arrival times

16
Advances in Civil Engineering are slightly different, indicating that the congestion has an influence on the time and space distribution of the pressure field. Combining the opening pressure arrival time and the overpressure double-peak interval time in the five regions, as shown in Figure 18, it is found that region I, which is most affected by the turbulence of congestion in all directions, increases to the opening pressure first. However, after the explosion vent is opened, the peak overpressure P 2 in this region is the smallest. In region IV, with the most small-scale congestion, peak overpressure P 2 is the greatest; the interval between the double peaks of overpressure in this region is the shortest. Figure 19 shows the peak overpressure-distance curve on the central axis of the explosion vent. Similar to the development of the turbulent kinetic energy and combustion rate, region I is most affected by the explosion vent. Near the explosion vent, the peak overpressure decreases rapidly; region IV is most affected by small-scale congestion, and the peak overpressure increases slightly near the cooker hood.

Peak Combustion Rate and Flame. Figures 20 and 21
show the curves of the peak combustion rate-distance and flame speed-distance on the central axis of the room with different volume blockages. It is observed in the figures that the peak combustion rate on the indoor axis shows an increasing trend with distance. e maximum peak combustion rate appears near the explosion vent and decreases with an increase in the volume blockage. e variation trend of the indoor flame propagation speed is similar to that of the turbulent kinetic energy; accelerated flame propagation occurs at the junction of the sudden change in the congestion cross-section and the small congestion and near the explosion vent. From analysis of the flame cloud map at a height of 1.4 m with different volume blockages in Figure 22, it is observed that because the flame development space decreases with an increase in volume blockage, the flame propagates to the two side walls faster in the first half and is constrained by the congestion. Forward propagation accelerates, thereby accelerating the effect of the cross-sectional mutation cabinet, small-scale cooker hood, and explosion vent on the propagation path, ultimately leading to faster flame propagation. In the second half, there may be more combustible gas in the room with a small volume blockage. An increase in the gas cloud released to the outside leads to a more violent external explosion. e explosion congestion causes more gas to burn in the room; the flame propagation speed fluctuates near the explosion vent, as shown in Figures 21 and 22. Figure 23 shows the indoor and outdoor overpressure and combustion rate-time curves with different volume blockages. Figure 24 shows the maximum peak overpressure-distance curves on the central axis with five volume blockages. Figure 25 shows the curves of the opening pressure arrival time and double peak interval of the overpressure for different volume blockages. It is observed in the figure that with different volume blockages, the indoor overpressure curves all present a bimodal structure. Peak P 1 is essentially the same; its value is approximately equal to the opening pressure of 20 kPa, but its arrival time is earlier with an increase in the volume blockage. Peak P 2 is the maximum indoor peak. Generally, the distribution trend of the maximum indoor peak overpressure is essentially the same with different volume blockages, decreasing with an increase in volume blockage. e time between the overpressure double peaks decreases first and then increases with an increase in the volume blockage. It is observed in Figure 23 that external explosions occur with different volume blockages. With an increase in volume blockage, the peak outdoor overpressure is larger, and the peak indoor overpressure is smaller.   Although the increase in volume blockage reduces the combustible gas inside and outside, it increases the indoor flame propagation speed, which leads to a stronger external explosion. Peak P 2 did not increase with increased intensity of the external explosion. It may be that the maximum peak overpressure in the room results from the combined effects of the external explosion and the maximum flame area [33]. e flame development area in a room with a low volume blockage is larger, thus the peak indoor combustion rate and maximum peak overpressure are greater, as shown in Figures 20 and 24.

Conclusion
(1) e congestion and explosion vent have a synergistic effect on the explosion flow field. Opening explosion venting surfaces such as building doors and windows increases the indoor turbulent kinetic energy, combustion rate, and flame propagation speed to some extent, but the impact is small. A confinement explosion can significantly reduce the indoor peak overpressure. Compared with a vented explosion in an empty room, congestion objects of different scales such as furniture and home appliances significantly increase the turbulence, combustion rate, flame speed, and peak overpressure in a vented explosion room. All types of congestion in the kitchen and the explosion venting surface coordinate to increase the severity of natural gas explosion accidents. (2) e indoor gas explosion flow field distribution is asymmetric under the action of congestion objects of different sizes and layouts. Small-scale congestion such as the lampblack machine in the kitchen promotes the formation of turbulence in the adjacent area and increases the flame propagation speed and combustion rate. Large-scale congestion such as a refrigerator arranged in the vertical flame propagation direction hinders the development of turbulence and flame propagation and reduces the combustion rate. Large-scale cabinets and cupboards arranged along the flame propagation direction have little influence on the indoor turbulence. e overpressure distribution in indoor regions is relatively balanced, mostly between 28 and 29 kPa. e explosion vent and its wall have a greater impact on the overpressure distribution.
(3) An increase in volume blockage from increased large-scale indoor congestion has a significant impact on the vented explosion process. Increasing the size of the congestion along the flame propagation direction promotes a change in the flame propagation mode and a sudden change in the congestion cross-section, leading to an increase in the flame propagation speed and acceleration of the increase in internal pressure. An external explosion is greater in a room with a high congestion rate. e indoor peak combustion rate and peak overpressure decrease as the volume blockage increases because there are more combustible gases in a room with low congestion and a larger flame area during complete combustion.  Temperature (K) m fu : e fuel mass fraction H c : Heat of combustion (kJ/mol) C μ : Constant of k-ε model (m 2 /s) R fu : Volume combustion rate (kg/m 3 /s) R min : Minimum mass fraction C t : Combustion-model constant u: Velocity (m/s) u t : Turbulence intensity (m/s) L t : Turbulent length scale (m) S l : Laminar combustion velocity (m/s) R f : Spherical flame radius F s : Laminar flame acceleration coefficient Γ * : Turbulent-diffusion coefficient (m 2 /s) τ ij : Viscous stress tensor μ t : Turbulence-viscosity coefficient ε: Dissipation rate of turbulent kinetic energy (m 2 /s 3 ) δ ij : e Kronecker delta (σ) * : Default value of the turbulent Prandtl constant v: Fluid dynamic viscosity (m 2 /s) ρ: Gas density (kg/m 3 ).
Data Availability e data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.