Vortex Flow Research Articles

Effect of oyster shell filling in artificial reefs on flow field environment and assessing the potential of carbon fixation

Artificial reefs are basic fishery facilities for marine habitat conservation and construction. In order to enhance fishery resources and restore the ecological environment in Laizhou Bay, where is an important fishing ground in China, an assembly-type oyster reef was designed based on the biological community, water depth and sea current in Laizhou Bay. This paper studied the effect of filling oyster shell in the assembly-type oyster reef on the flow field distribution by computational fluid dynamics (CFD), optimized the structure of reef, and presented the construction deployment of reef habitat and assessed potential of carbon fixation of oyster reefs. Indexes of upwelling, slow flow and vortex were chosen to describe the flow field effect of oyster reefs. The distribution characteristics of flow field under three types of filling oyster shell were analyzed: filling no shell (UAR), filling shells with an 83.6 void ratio (OAR), and filling shells with a 0 void ratio (FAR). The results showed that the upwelling, vortex and slow flow efficiency indicators of the OAR had obvious advantages compared with the other two filling methods, and the efficiency indicators of OAR with the spacing between basic elements of 280 m and deployment angle of 0° were higher than the others. Finally, according to the study result of filling method, spacing and deployment of oyster reef, the assessment showed that reefs could fix 2178.9 t carbon by themselves on the basis of national marine ranching demonstration area regulations in China.

Gaussian smoothed particle hydrodynamics: A high-order meshfree particle method

Smoothed particle hydrodynamics (SPH) has attracted significant attention in recent decades, and exhibits special advantages in modeling complex flows with multiphysics processes and complex phenomena. Its accuracy depends heavily on the distribution of particles, and will generally be lower if the particles are distributed non-uniformly. A high-order SPH scheme is proposed in the present work for simulating both compressible and incompressible flows. It uses a Gaussian quadrature rule to perform the particle approximation of SPH by introducing Gaussian nodes. Unfortunately, the Gaussian nodes hardly overlap with SPH particles due to the Lagrangian feature, and thus we use a high-order interpolation method to obtain the corresponding physical quantities at the Gaussian nodes. The accuracy and robustness of the proposed Gaussian SPH are demonstrated by several numerical tests, including the Sod problem, Poiseuille flow, Couette flow, cavity flow, Taylor–Green vortex and dam break flow, and a convergence analysis is also conducted to evaluate the effects of particle resolution and distribution for reconstructing a given function. The simulation results for each test case are in good agreements with the available analytical, experimental or numerical results, showing that the proposed Gaussian SPH method is accurate and reliable but expensive for simulating compressible and incompressible flow problems.

Numerical investigation of energy dissipation and vortex characteristics on the S-shaped region of a reversible pump-turbine

In order to meet the regulation needs of power system, the pump-turbine is prone to operate in the S-shaped region during frequent start-up, leading to unit vibration and grid connection failure. In this paper, the S characteristic experiment was carried out, and the entropy production, unstable flow and vortex structure under turbine, runaway, turbine breaking, reverse pump conditions with optimum guide vane opening are detailed analyzed. The relationship between the energy dissipation and flow field characteristics is clarified. The results show that the runner entropy production is the largest in all components, and the wall entropy production is greater than the fluctuating entropy production under turbine condition. With the decrease in discharge, the fluctuating entropy production and wall entropy production first increase and then decrease, reaching the maximum values in the runaway condition. The circumferential movement trend of water enhances, and forms circumferential flow in the vaneless region. The separation vortex formed on the blade suction surface extends from the blade inlet to the middle of the flow channel, inducing a high entropy production. Under the turbine braking condition, the tangential velocity increases rapidly in the vaneless region. The circumferential flow develops into water ring, blocking the water flows into the runner. Under the reverse pump condition, the water flows in the opposite direction, with negative tangential and streamwise velocities. The absolute flow angle becomes obtuse, and a large-scale reverse vortex is formed. The runner inlet is blocked by the water ring, and the flow channel is almost completely filled with the low-speed region. The proportion of the guide vanes entropy production reaches maximum. Unstable flow leads to the increase in vorticity. Separation vortex, bypassed vortex, and reverse flow vortex leads to an increase in VST, forming a high VST region at corresponding positions. Due to the rotational motion accompanying the development of separated vortex, it leads to the increase in CFT. In addition, the circumferential flow and water ring have higher tangential velocities, which can lead to high CFT regions.

Development of the full Lagrangian approach for modeling dilute dispersed media flows (a review)

Continuum models of media with zero pressure are widely used in various branches of physics and mechanics, including studies of a dilute dispersed phase in multiphase flows. In zero-pressure media, the particle trajectories may intersect, “folds” and “puckers” of the phase volume may arise, and “caustics” (the envelopes of particle trajectories) may appear, near which the density of the medium sharply increases. In recent decades, the phenomena of clustering and aerodynamic focusing of inertial admixture in gas and liquid flows have attracted increasing attention of researchers. This is due to the importance of taking into account the inhomogeneities in the impurity concentration when describing the transport of aerosol pollutants in the environment, the mechanisms of droplet growth in rain clouds, scattering of radiation by dispersed inclusions, initiation of detonation in two-phase mixtures, as well as when solving problems of two-phase aerodynamics, interpretation of measurements obtained by LDV or PIV methods, and in many other applications. These problems gave an impetus to a significant increase in the number of publications devoted to the processes of accumulation and clustering of inertial particles in gas and liquid flows. Within the framework of classical two-fluid models and standard Eulerian approaches assuming single-valuedness of continuum parameters of the media, it turns out impossible to describe zones of multi-valued velocity fields and density singularities in flows with crossing particle trajectories. One of the alternatives is the full Lagrangian approach proposed by the author earlier. In recent years, this approach has been further developed in combination with averaged Eulerian and Lagrangian (vortex-blob method) methods for describing the dynamics of the carrier phase. Such combined approaches made it possible to study the structure of local zones of accumulation of inertial particles in vortex, transient, and turbulent flows. This article describes the basic ideas of the full Lagrangian approach, provides examples of the most significant results which illustrate the unique capabilities of the method, and gives an overview of the main directions of further development of the method as applied to transient, vortex, and turbulent flows of “gas-particle” media. Some of the ideas discussed and the results presented below are of a more general interest, since they are also applicable to other models of zero-pressure media.

An experimental study of vortex evolution around an offshore stationary dual-chamber oscillating water column converter

The Oscillating Water Column (OWC) is a promising Wave Energy Converter (WEC) due to its practicality and efficiency. Conventional single-chamber OWCs exhibit constraints in extracting energy over diverse wave periods, while dual-chamber OWCs demonstrate enhanced adaptability, leading to improved energy capture rates. This study investigates an offshore stationary dual-chamber OWC with an underlying transverse channel designed to facilitate wave ingress and augment energy extraction. Experiments across various wave periods reveal that dual-chamber OWCs with transverse channels achieve superior wave energy extraction compared to single-chamber counterparts. Time-resolved particle image velocimetry (TR-PIV) is employed for detailed flow field measurement to elucidate the complex hydrodynamic processes. Analysis of the water column oscillation and flow dynamics reveals an intensified interaction between hydrodynamic processes in the dual chambers, especially with elongated transverse channels. The Ω method is used to examine vortex evolution, indicating heightened complexity in vortex formation within dual-chamber OWCs compared to single-chamber OWCs. Notably, large-scale vortices emerging in the forward chamber impede oscillation and diminish energy extraction efficiency. Optimizing the design of transverse channel entries, frontal walls, and internal chamber corners is crucial to mitigate flow separation and vortex generation, thereby enhancing overall energy extraction performance.

Leidenfrost Effect-Induced Chaotic Vortex Flow for Efficient Mixing of Highly Viscous Droplets.

Efficiently mixing highly viscous liquids in microfluidic systems is appealing for green chemistry such as chemical synthesis and catalysis, but it is a long-standing challenge owing to the unfavorable diffusion kinetics. In this work, a new strategy is explored for mixing viscous droplets by harnessing a peculiar Leidenfrost state, where the substrate temperature is above the boiling point of the liquid without apparent liquid evaporation. Compared to the control experiment where the droplet stays at a similar temperature but in the contact boiling regime, the mixing time can be reduced significantly. Moreover, it is demonstrated that the liquid mixing originates from the chaotic convection flow in the Leidenfrost droplet, characterized by the internal vortex motion evidenced by the microscale visualization. A correlation between mixing time and droplet volume is also proposed, showing a good agreement with experimental results. It is further shown that Leidenfrost droplets can be used to synthesize nanoparticles of the desired morphology, and it is anticipated that this simple and scalable fabrication approach will find applications in the biological, pharmaceutical, and chemical industries.

Applications of partial differential equations to transient heat and mass transfer in MHD flow over a porous medium

The flow problem discussed in this research is an investigation of compressible viscous fluid in an unstable 2D MHD laminar driven vortex and stable boundary layer flow, passing a heated stretched sheet at varying speeds via a porous medium exposed to a magnetic field with viscous dissipation, thermal diffusion, thermal radiation, chemical reaction and Dufour effects. With the aid of similarity transformations, PDEs that come out in governing equations of this formulations are converted into non-linear sets of ODEs. Well-defined bvp4c technique used to numerically solve these converted equations. With use of various graphs, the consequences of different factors on the distributions of velocity, temperature and concentration were examined numerically and graphically. Physical parameters for this issue are also discussed.

Turbulence approaches for numerical predictions of vehicle-like afterbody vortex flows

The numerical prediction of aerodynamic characteristics for vehicles is crucial to both industry and academia, with various numerical approaches playing a critical role in accurately resolving flow fields. This study aims to evaluate the effectiveness of three typical numerical approaches, including RANS, IDDES, and LES in predicting the afterbody vortex flows of a generic model, specifically a slanted-base cylinder. This study involved analyzing aerodynamic coefficients, time-averaged surface flow, time-averaged surrounding flow and transient flow, revealing the capabilities of each approach. RANS offers acceptable accuracy in predicting time-averaged aerodynamic coefficients and surface flow patterns, though it falls short in capturing time-varying physical quantities. LES, despite its higher computational cost, provides a more accurate prediction for both time-averaged and transient flow behaviors, particularly in capturing flow instabilities and multi-scale fluctuations. IDDES can be prioritized when a rough understanding of transient characteristics is sufficient. This study highlights the unique strengths and limitations of three typical numerical approaches in predicting vehicle-like afterbody vortex flows, guiding the selection of appropriate methods based on specific research needs.

Steady state regimes in lid driven cavity flows of magnetic fluids in the presence of an external magnetic field

In this work, we study the flow of magnetic fluids in a square top lid driven cavity in the presence of an applied magnetic field. The magnetic field is generated by a steady electric current flowing through a wire placed underneath the cavity. The magnetization of the magnetic fluid is described by an evolution equation that combines magnetization relaxation, advection and interaction with the flow via the vorticity. The governing equations are solved via a vorticity-streamfunction formulation implemented in a finite difference scheme. The main governing physical parameters of the examined cavity flow are the Reynolds number, ranging from 50 to 500, the magnetic pressure coefficient, which measures the relative importance between magnetic and inertial forces, ranging from 0 to 1000, and Péclet number which is O(1) and measures the ratio between the magnetization relaxation time and the flow time scale. We have examined the different mechanisms which cause changes in the local fluid magnetization, i.e., flow vorticity and advection. In particular, we examine how vorticity of the flow intensifies deviations of the magnetization of the fluid with respect to the local equilibrium of magnetization. In addition, we identify a complex topological structure of the flow inside the cavity characterized by the magnetic effects that dominate the generation of secondary structures of vorticity induced by the effect of field gradient, orientation and a non-uniformity in the fluid magnetization inside the cavity.

Numerical study of the impacts of stochastic forcing on the vortex in fluid flow

This paper focuses on a numerical study about the stochastic Navier–Stokes equations. Unlike previous studies, this paper focuses on studying these equations from the perspective of vortices. The vorticity–stream function method was proposed to deal with incompressible fluid flow. And a Crank–Nicolson Fourier pseudo-spectral method was put forward to solve the formulation of stream function equation. In addition, we have conducted some numerical experiments to observe the effects of random forcing on vortices in the fluid flow by utilizing stochastic solution.

The effect of Reynolds number on the separated flow over a low-aspect-ratio wing

At high incidence, low-aspect-ratio wings present a unique set of aerodynamic characteristics, including flow separation, vortex shedding and unsteady force production. Furthermore, low-aspect-ratio wings exhibit a highly impactful tip vortex, which introduces strong spanwise gradients into an already complex flow. In this work, we explore the interaction between leading-edge flow separation and a strong, persistent tip vortex over a Reynolds number range of $600 \leq Re \leq 10{\,}000$ . In performing this study, we aim to bridge the insight gained from existing low-Reynolds-number studies of separated flow on finite wings ( $Re \approx 10^2$ ) and turbulent flows at higher Reynolds numbers ( $Re \approx 10^4$ ). Our study suggests two primary effects of the Reynolds number. First, we observe a break from periodicity, along with a dramatic increase in the intensity and concentration of small-scale eddies, as we shift from $Re = 600$ to $Re = 2500$ . Second, we observe that many of our flow diagnostics, including the time-averaged aerodynamic force, exhibit reduced sensitivity to Reynolds number beyond $Re = 2500$ , an observation attributed to the stabilising impact of the wing tip vortex. This latter point illustrates the manner by which the tip vortex drives flow over low-aspect-ratio wings, and provides insight into how our existing understanding of this flow field may be adjusted for higher-Reynolds-number applications.

Changes of intracardiac flow dynamics measured by HyperDoppler in patients with aortic stenosis.

Assessment of intracardiac flow dynamics has recently acquired significance due to the development of new measurement methods based on echocardiography. Recent studies have demonstrated that cardiac abnormalities are associated with changes in intracardiac vortical flows. Yet, no previous study assessed the impact of aortic stenosis (AS) on intracardiac vortices. This study aims to explore the clinical potential of additional information provided by quantifying intracardiac flow dynamics in patients with AS. One hundred and twenty patients with severe AS, sixty patients with concentric ventricular remodelling (VR), and hundred controls (CTRL) were prospectively included and underwent non-invasive evaluation of intracardiac flow dynamics. In addition to standard echocardiography, fluid dynamics were assessed by means of HyperDoppler. Vortex depth (P < 0.001), vortex length (P = 0.003), vortex intensity (P < 0.001), and vortex area (P = 0.049) were significantly increased in AS compared with CTRL. In addition, mean energy dissipation was significantly higher in AS compared with CTRL (P < 0.001) and VR (P = 0.002). At receiver operating characteristic analysis, vortex depth showed the best discrimination capacity for AS (P < 0.001). Changes in fluid dynamics-based HyperDoppler indices can be reliably assessed in patients with AS. Significant changes in vortex depth and intensity can selectively differentiate AS from both concentric remodelling and healthy CTRLs, suggesting that the assessment of intracardiac flow dynamics may provide complementary information to standard echocardiography to better characterize patients' subsets.

Parametric study of a vortex-enhanced supersonic inductive plasma torch

The feed gas injection configuration in radio-frequency (RF) inductively coupled plasma (ICP) torches plays a critical role in discharge stability, gas heating, and device thermal management: particularly if a supersonic nozzle is used to subsequently accelerate the hot gas. A novel injection configuration is the bidirectional vortex, which segments the internal ICP flow field into two counter-propagating vortices that can significantly enhance gas heating and reduce heat losses. The diameter of the interface between the vortices (known as the mantle) is expected to be an important dimensional parameter affecting torch operation, especially relative to the nozzle size. In this work, we investigate the effect of nozzle throat diameter on the behaviour and performance of a vortex-enhanced supersonic ICP torch. The system is operated at RF powers and argon mass flow rates between 200–1000 W and 0–400 mg s−1 respectively, and different nozzle diameters ranging from 1.5 to 4 mm are explored. Because of the high-temperature environment, and to prevent disruption of the vortex flow fields, non-invasive diagnostics are used to measure the gas temperature and plasma density, and to infer the torch thermal efficiency and achievable gas specific enthalpy change. The maximum temperature is between 8500–9500 K with the 1.5 mm nozzle giving the highest temperature for a given power and mass flow rate, while plasma densities vary between 1020–1021 m−3 depending on the operating conditions. The thermal efficiency increases from 29% for the 1.5 mm nozzle to just above 70% for the 4 mm nozzle with a similar maximum specific enthalpy of around 1.5 MJ kg−1. These results demonstrate the important coupling between torch properties, and how system optimization can lead to tailored performance of potential interest to several ground and space-based applications.

Properties and Model of Pore-Scale Methane Displacing Water in Hydrate-Bearing Sediments

The flow characteristics of methane and water in sedimentary layers are important factors that affect the beneficial exploitation of marine hydrates. To study the influencing factors of methane drive-off water processes in porous media, we constructed nonhomogeneous geometric models using MATLAB 2020a random distribution functions. We developed a mathematical model of gas–water two-phase flow based on the Navier–Stokes equation. The gas-driven water processes in porous media were described using the level-set method and solved through the finite element method. We investigated the effects of the nonhomogeneous structure of pore media, wettability, and repulsion rate on gas-driven water channeling. The nonhomogeneity of the pore medium is the most critical factor influencing the flow. The size of the throat within the hydrophilic environment determines the level of difficulty of gas-driven water flow. In regions with a high concentration of narrow passages, the formation of extensive air-locked areas is more likely, leading to a decrease in the efficiency of the flow channel. In the gas–water drive process, water saturation changes over time according to a negative exponential function relationship. The more hydrophilic the pore medium, the more difficult the gas-phase drive becomes, and this correlation is particularly noticeable at higher drive rates. The significant pressure differentials caused by the high drive-off velocities lead to quicker methane breakthroughs. Instantaneous flow rates at narrow throats can be up to two orders of magnitude higher than average. Additionally, there is a susceptibility to vortex flow in the area where the throat connects to the orifice. The results of this study can enhance our understanding of gas–water two-phase flow in porous media and help commercialize the exploitation of clean energy in the deep ocean.

Impingement of vortex dipole on heated boundaries and related thermal plume dynamics

The profound influence of an externally induced vortex dipole on thermal plume dynamics is numerically studied for varying Rayleigh numbers (Ra) employing the Bhatnagar–Gross–Krook collision model-based lattice Boltzmann method with a double distribution function approach. This study is extended to vortex dipole impingement with different types of heated bottom boundaries of two-dimensional domain, such as flat, “V-shaped,” and “inverted-V-shaped.” The vortex dipole impingement with the heated boundaries generates secondary vortices, which in turn produce vortex-driven thermal plumes, thereby advancing plume generation. The subsequent merging of the plumes enhances heat transport and leads to a continuous plume ascent. The presence of convex corners facilitates flow separation and also gives rise to the formation of secondary vortex dipoles, thereby significantly impacting the continuous generation of jet-like plumes when compared to concave configurations. The lack of an external vortex in pure buoyancy-driven flows produces less pronounced jet-like plumes and a relatively low Nusselt number. The boundary types and Ra significantly influence the vorticity production, resulting in higher enstrophy and palinstrophy for convex boundaries compared to flat and concave ones. A lower Prandtl number increases secondary vortices and corner rolls, leading to larger velocity gradients, higher thermal diffusivity, resulting in increased kinetic energy and thermal dissipation rates. The increased cell height enhances heat transfer at the top boundary due to improved heat convection from the slanted boundary and influence of early dipole impingement. Furthermore, kinetic energy dissipates in the dipole-driven flows and increases in the buoyancy-dominated flows.

Topology optimization with vortex-intensity-based objective function enables novel film cooling hole with superior performance

Film cooling is widely adopted in the thermal protection of modern aerospace engines. The optimization of the film hole is an important issue in thermal designs. In this work, a density-based topology optimization is conducted to obtain a high-performance film cooling hole geometry. Unlike the traditional wall-temperature-based optimization objectives, a novel objective function based on the Q-criteria is proposed by taking the vorticity dynamics and flow mixing into account. Finally, a unique high-performance film hole is obtained, which shows superior performance compared to the standard cylindrical hole, particularly at high blowing ratios. The film cooling effectiveness can be elevated by 51 folds in the case of a high blowing ratio (M = 1.5). Analyses of the vortices and velocity fields suggest that the presence of the additional secondary vortex pair enhances the attachment of cooling air to the wall.

Scaling Up Perovskite Solar Cell Fabrication: Antisolvent‐Controlled Crystallization of Printed Perovskite Semiconductor

Scaling up perovskite solar cells stands as one of the frontiers in advancing this rapidly growing technology. Yet, controlling perovskite thin‐film crystallization during and post‐printing differs significantly from lab‐scale processes that have yielded record device efficiencies. This study investigates antisolvent treatment for slot‐die‐coated perovskite solar cells using in situ optical spectroscopy and comparing among multiple antisolvents. The antisolvent bath used in slot‐die coating affects the perovskite crystallization and film quality differently when comparing to the established spin‐coating antisolvent treatment process. A novel dynamic antisolvent method, employing either vortex or laminar flow, is developed. It outperforms steady‐bath techniques in generating high‐quality, haze‐free films. Optimization studies identify critical treatment times. Implementing this novel antisolvent treatment leads to a peak average power conversion efficiency of 15.62% and the highest device efficiency of 18.57%, an excellent performance for slot‐die‐coated MAPbI3 devices printed and tested under ambient conditions. The method is validated for an alternative perovskite composition, FA0.9Cs0.1PbI3, and printing technique, blade coating. This research highlights the importance of in situ analysis for enhancing perovskite film quality and introduces scalable approaches for controlling large‐area film crystallization kinetics, driven by the demand for efficient and scalable manufacturing processes in the field of perovskite solar cells.

Effects of submergence ratio on flow around a circular pier in combined wave–current flows using particle image velocimetry

The present study illustrates an experimental investigation of flow hydrodynamics in the vicinity of a submerged circular pier across various submergence levels under only current and wave–current combined flow conditions. The instantaneous velocity data are collected using particle image velocimetry for three distinct frequencies of waves to determine the influence of wave superimposition on the current-induced turbulence parameters. The distribution of phase averaged turbulence quantities, such as mean velocities, Reynolds shear stress, turbulent kinetic energy, flow patterns, and vorticity analysis by Q-criterion, are presented. The results provide insight into the impacts of wave frequency and submergence ratio on the formation of horseshoe vortices, trailing vortex, and reverse flow zones. From the results it is observed that a decrease in the submergence level of the structure causes the formation of strong horseshoe vortices and reverse flow zones in the absence of waves. Also, an increase in wave frequency intensifies the turbulence kinetic energy at the upstream of the pier and eddy generation behind the pier. The present findings highlight the effects of pier submergence and wave characteristics, such as frequency, wave height, and wave period, on flow patterns and turbulent flow characteristics and aid in the design of marine structures. Furthermore, experimental data serve as a valuable resource for validating theoretical or mathematical models related to combined wave–current environments.

Flow structure around a microswimmer at fluid–fluid interface

Living organisms such as bacteria and algae often form biofilms at air–liquid and/or liquid–liquid interfaces. Therefore, it is important to understand the hydrodynamic interaction between the fluid–fluid interface and microorganisms. In the present study, a squirmer model is employed to study the flow field structures around a symmetrically trapped spherical microswimmer at an interface separating two fluids with viscosity contrast. In these simulations, Reynolds ( Re ) and Capillary ( Ca ) numbers are very small, and hence the contribution of inertia and interface deformations are neglected. The squirmer model is validated against the analytical solution obtained by considering a source dipole. It is observed that the flow structure and vorticity distribution around the microswimmers are strongly influenced by the squirmer parameter (β) and viscosity contrast (λ). Furthermore, the interplay between force-dipole and source-dipole along with viscosity contrast leads to a range of flow structures such as symmetric four-lobe and asymmetric quadrupolar flow fields. This flow structure asymmetry around the swimmer is quantified with different steady state orientations (φ). The effect of flow profile on the passive tracer particle transport is explored in terms of trajectories. Specifically, larger values of φ result in more profound curly trajectories.

Numerical study on flame acceleration and deflagration-to-detonation transition: Spatial distribution of solid obstacles

This study conducts a detailed numerical investigation on the spatial distribution of solid obstacles using the large eddy simulation method. It is discovered that although flame acceleration induced by solid obstacles is dominated by factors such as flow field disturbances, vortices and recirculation zones, turbulence, flame surface areas and combustion heat release rates, etc., the characteristics of the leading shock wave are key to detonation initiation. Specifically, the intensity of the leading shock wave, its formation time, and its distance from the flame front significantly affect detonation initiation. Depending on the state of the shock wave, the detonation initiation process may occur through various mechanisms such as shock reflection, shock focusing. Overall, the types of detonation initiation in this study all belong to the shock detonation transition. However, the detonation initiation process can be further classified into two categories: (I) Detonation induced by shock wave reflection; (II) detonation triggered by shock wave focusing. Despite certain disparities in the detonation initiation process, all detonation initiation processes conform to the gradient theory, and the flame evolution processes in all cases consistently follow three stages: the laminar slow-ignition stage; the turbulent deflagration stage; the detonation initiation stage. Furthermore, the study further discerns that, compared to positioning obstacles on the wall, placing obstacles inside the combustion chamber can further augment the detonation-assisting effect. However, excessively sparse or dense spatial distributions of solid obstacles fail to yield the optimal detonation effect. An optimal distribution exists, which triggers the fastest detonation initiation.