Two-phase Flow Field Research Articles

Most Suitable Blade Inclination Angle for Multiphase Flow in Soybean Milk Machine

The blade inclination angle of soybean milk machine is a key geometric parameter for efficient crushing. For the purpose of obtaining optimal design, the gas-liquid two-phase flow field inside a soybean milk machine is simulated. The gas holdup from simulation is in agreement with the experiment. The simulation result shows that the lower blade A has a great influence on the internal flow field of soybean milk machine, while the upper blade B has a small influence on the flow field. As the angle l αA l increases, the peak value of radial velocity decreases and moves to the interior of the cavity, so does the total pressure. When αA changes from -24° to -26°, the velocity vector at the bottom of the cavity changes from the connected state to the separated state, and the pressure difference between the up and the bottom surface of blade A becomes large. When αA = 24°, the flow field has the strongest turbulent kinetic energy and dissipation. When αB =28°, the pressure difference reaches the maximum. In summary, the best inclination angles are αAopt =-24° ∼ -26° and αBopt =28°, respectively.

Research on multiphysics coupling fields in electrochemical trepanning of lateral flow

Electrochemical trepanning (ECTr) is an effective method for machining the ruled surface parts. Generally, the forward flow mode is used in ECTr, under which the streamlines at the outlet are divergent, resulting in the obvious flow patterns at the outlet and the instability of the machining process. In ECTr of a diffuser with a special structure, the lateral flow mode is adopted to improve the uniformity of the flow field, thereby improving the surface quality at the hub. ECTr is a complicated multiphysics coupling process. To investigate the distributions of electric field, two-phase flow field and thermal field in ECTr with lateral flow, a multiphysics coupling field model was established. In this model, a coupling relationship was formed between the various physical fields through the change of the electrolyte conductivity. And a coupled multiphysics simulation was performed to show how the gas bubbles volume fraction, electrolyte temperature, electrolyte conductivity and current density change along the flow path. Compared with the inlet of the electrolyte, the gas bubbles volume fraction and the temperature at the outlet increased by 38.8% and 6.3 K, respectively. Under the combined influence, the conductivity decreased by 7.227 S/m at the outlet, resulting in a decrease of 57.81 A/cm2 in the current density. The simulation results were then verified against those from a corresponding experiment involving lateral flow ECTr. Along the flow path, the thickness of the machined blade gradually increased, varying from 2.09 mm to 2.76 mm. The surface quality gradually deteriorated along the flow path and the surface roughness varied from Ra 0.72 μm to Ra 1.05 μm. Combining the simulation and the experiment, the correctness and the effectiveness of the multiphysics coupling model and simulation were confirmed. The results can be applied to other ECM processes.

Multi-Physics Coupling Modeling and Experimental Investigation of Vibration-Assisted Blisk Channel ECM.

Due to its advantages of good surface quality and not being affected by material hardness, electrochemical machining (ECM) is suitable for the machining of blisk, which is known for its hard-to-machine materials and complex shapes. However, because of the unstable processing and low machining quality, conventional linear feeding blisk ECM has difficulty in obtaining a complex structure. To settle this problem, the vibration-assisted ECM method is introduced to machine blisk channels in this paper. To analyze the influence of vibration on the process of ECM, a two-phase flow field model is established based on the RANS k-ε turbulence model, which is suitable for narrow flow field and high flow velocity. The model is coupled with the electric field, the flow field, and the temperature field to form a multi-physics field coupling model. In addition, dynamic simulation is carried out on account of the multi-physics field coupling model and comparative experiments are conducted using the self-developed ECM machine tool. While a shortcut appeared in the contrast experiment, machining with vibration-assisted channel ECM achieved fine machining stability and surface quality. The workpiece obtained by vibration-assisted channel ECM has three narrow and straight channels, with a width of less than 3 mm, an aspect ratio of more than 8, and an average surface roughness Ra in the hub of 0.327 μm. Compared with experimental data, the maximum relative errors of simulation are only 1.05% in channel width and 8.11% in machining current, which indicates that the multi-physics field coupling model is close to machining reality.

Flow visualisation in real-size optical injectors of conventional, additised, and renewable gasoline blends

Research on renewable and alternative fuels is crucial for improving the energy and environmental efficiency of modern gasoline internal combustion engines. To highlight the influence of fuel rheological and thermodynamic properties on phase change and atomisation processes, three types of gasoline blends were tested. More specifically, the campaign comprised a reference gasoline, an ethanol/gasoline blend (10% v/v) representative of renewable fuels, and an additised gasoline sample treated with viscoelasticity-inducing agents. High-speed imaging of the transient two-phase flow field arising in the internal geometry and the near-nozzle spray region of gasoline injectors was performed employing Diffuse Backlight Illumination. The metallic body of a commercial injector was modified to fit transparent tips realising two nozzle layouts, namely a two-hole real size model resembling the Engine Combustion Network spray G injector and an enraged replica with an offset hole. Experiments were conducted at realistic operating conditions comprising an injection pressure of 100 bar and ambient pressures in the range of 0.1–6.0 bar to cover the entire range of chamber pressures prevailing in Gasoline Direct Injection engines. The action of viscoelastic additives was verified to have a suppressive effect on in-nozzle cavitation (6% reduction in cavitation extent) , while also enhancing spray atomisation at flash-boing conditions, in a manner resembling the more volatile gasoline/ethanol blends. Finally, persisting liquid ligaments were found to form after the end of injection for the additised sample, owing to the surfactant nature of the additives.

Hydrodynamic Analysis of Self-Propulsion Performance of Wave-Driven Catamaran

The wave-driven catamaran is a small surface vehicle driven by ocean waves. It consists of a hull and hydrofoils, and has a multi-body dynamic structure. The process of moving from static state to autonomous navigation driven by ocean waves is called “self-propulsion”, and reflects the ability of the wave-driven catamaran to absorb oceanic wave energy. Considering the importance of the design of the wave-driven catamaran, its self-propulsion performance should be comprehensively analysed. However, the wave-driven catamaran’s multi-body dynamic structure, unpredictable dynamic and kinematic responses driven by waves make it difficult to analyse its self-propulsion performance. In this paper, firstly, a multi-body dynamic model is established for wave-driven catamaran. Secondly, a two-phase numerical flow field containing water and air is established. Thirdly, a numerical simulation method for the self-propulsion process of the wave-driven catamaran is proposed by combining the multi-body dynamic model with a numerical flow field. Through numerical simulation, the hydrodynamic response, including the thrust of the hydrofoils, the resistance of the hull and the sailing velocity of the wave-driven catamaran are identified and comprehensively analysed. Lastly, the accuracy of the numerical simulation results is verified through a self-propulsion test in a towing tank. In contrast with previous research, this method combines multi-body dynamics with computational fluid dynamics (CFD) to avoid errors caused by artificially setting the motion mode of the catamaran, and calculates the real velocity of the catamaran.

Experimental study on the flow and heat transfer behavior of polymer solutions in the closed liquid ring vacuum pump system

The closed liquid ring vacuum pump system, which mainly includes the liquid ring vacuum pump and the heat exchanger, is widely used for gas pumping in various industries. In this paper, the energy saving effect in the liquid ring vacuum pump and the flow and heat transfer behavior in the shell-and-tube heat exchanger are investigated for aqueous solutions of rigid xanthan gum (XG) and flexible Polyacrylamide (PAM), and the effects of the concentration, Reynolds number, polymer conformation and flow geometry on the flow behavior of polymer solutions are studied. Results show that, despite the significant differences in molecular conformation and rheological properties between XG and PAM, the similar energy saving rate, about 1.5 ∼ 17.5%, is obtained in the liquid ring vacuum pump. However, in the shell-and-tube heat exchanger, XG solutions exhibit the more obvious pressure drop enhancement and heat transfer reduction than that of PAM solutions for the same concentration. Furthermore, due to the polymer shear degradation, the shaft power in the liquid ring vacuum pump decreases first and then stabilizes, conversely, the flow resistance decreases and the flow rate increases in the heat exchanger. The great difference in the flow condition between the intense rotating gas–liquid two-phase flow field and the gentle liquid pipe flow causes the different flow behavior for XG and PAM solutions.

Three-dimensional simulation and experimental investigation of electrolyte jet machining with the inclined nozzle

Electrolyte jet machining (EJM) is an innovative form of electrochemical machining (ECM) with excellent machining flexibility and accuracy. Recently, the nozzle inclination is introduced into EJM to further expand its machinability. In this work, to clearly reveal the mechanism behind this EJM process, a three-dimensional (3-D) multiphysics model of the stationary EJM process accounting for two-phase flow field and electric field was developed for the first time and experiments were also performed for validation. With the inclined nozzle, the distribution of flow velocity and thickness of electrolyte film around the jet exhibits uneven, thus resulting in non-circular hydraulic jump and non-axisymmetric anodic current density distribution, respectively. Besides, the effects of nozzle inclination angle and electrolyte ejecting flow velocity were investigated. Larger inclination angle leads to increased asymmetry of flow velocity and current density distribution, which results in more seriously non-circular hydraulic jump and asymmetric machined dimple shape. The electrolyte ejecting flow velocity affects the jet shape and current density distribution that eventually influences the machined results with the inclined nozzle, but this effect does not exist with the conventional vertical nozzle. Moreover, this work was extended to the translating EJM process to study the effect of nozzle translating direction with the inclined nozzle. It is demonstrated that through selecting translating directions, grooves with different shapes and surface finish can be machined even using a certain nozzle inclination. The mechanism for these phenomena was clarified by simulation, paving the way for design of the EJM process with the inclined nozzle.

Pressure calculation for flows with moving surface boundaries from particle tracking velocimetry (PTV)

The examples of flow conditions, where an object of a fixed or deformable body moves in a fluid, or the interface between the flow phases instantaneously changes its topology, are numerous in industry and natural sciences. The advent of particle image velocimetry (PIV) [1] and particle tracking velocimetry (PTV) [2] enabled the measurement of the instantaneous velocity fields in these types of complicated flow fields. As a next step, several methodologies have been developed in the past decade to calculate the pressure fields from PIV or PTV data [3,4]. These methods were developed based on the assumption of a stationary flow domain, with surface boundaries that are fixed and independent of time. This makes the current pressure calculation methods inapplicable to a flow domain with deformable moving surface boundaries. Also, for most of the two-phase flows, the capillary forces are significant and the pressure drop over the two-phase interface must be considered. Therefore, the current pressure calculators require an improvement in the formulation of the algorithms to account for the deformable volume conditions and the effect of the surface tension force. For the calculation of pressure from sparse PTV velocity data, firstly, a tessellation method is required to interconnect the irregularly spaced vectors in the flow field using a highquality mesh grid. The mesh must be dynamic and adjust itself to the moving boundaries. This tessellation method has already been developed by the current authors [5]. As the next step, equations of motion for a deformable C.V. need to be coupled with the tessellation method to calculate the instantaneous pressures in a two-phase flow field, with a moving interface, which will be the ultimate goal of the current study.

Comprehensive particle image velocimetry measurement and numerical model validations on the gas–liquid flow field in a lab-scale cyclonic flotation column

The investigation of flow field can provide strong theoretical support for improving the flotation performance of fine particles in a cyclonic flotation column. In this study, particle image velocimetry (PIV) is combined with endoscopic measurement and phase discrimination technique to measure the cross section and axial section of gas–liquid two-phase flow field in a lab-scale cyclonic flotation column. The axial, radial, and tangential velocity of both gas and liquid phases are obtained. Based on the PIV experimental data, computational fluid dynamics (CFD)-based numerical models are validated. The results show the good prediction ability of Eulerian–Eulerian multiphase model, Reynolds stress model (RSM), and Tomiyama drag model. Then the simulated results with validated models are displayed to indicate the flow characteristics of the flotation column. The gas concentration is observed to be higher in the center and low on the sides. The gas phase always moves upward, while the liquid phase moves upward at the axis, downward between the axis and the wall, and upward again near the wall. Overall, this study provides an innovative approach for PIV measurement of complex multiphase flow field, and a useful reference for appropriate numerical models selection.

Experimental investigation of solute transport in variably saturated porous media using x-ray computed tomography

Solute transport through variably saturated porous media is ubiquitous in multiple subsurface flows, piquing the geoscience community's interest. This study adopts a novel experimental approach using microfocus x-ray computed tomography for real-time imaging of a three-dimensional NaI tracer plume in a partially saturated packing column. A stabilized two-phase flow field is achievable through continuous co-injection of two-phase fluids: NaCl solvent and pump oil. Thus, the critical role of the NaCl saturation Sw and Péclet number on dispersion can be fully studied by controlling the NaCl fractional flow rate and the total flow rate from the Buckley–Leverett theory. Furthermore, we study solute transport behavior based on statistical moments, the dispersion coefficient, the dilution index, and the mean scalar dissipation rate. Experimental results indicate that the solute transport is Fickian for high Sw ≥ 0.34. In contrast, anomalous transport behavior is found for Sw < 0.34, where the concentration distribution is initially left-tailed and leptokurtic before reaching a well-dispersed regime. The dispersion coefficient is 2–10 times larger for partially saturated cases compared with the fully saturated case and shows a non-monotonical dependency on Sw. Finally, the analysis of the dilution index indicates that the overall mixing strength increases when Sw decreases, whereas the mean scalar dissipation rate reveals that the time scaling of transverse mixing is the largest at an intermediate Sw. The results can be used to elucidate the solute transport behavior in a two-phase system.

Local flow dynamics in the motion of slug bubbles in a flowing mini square channel

For a slug flow flowing in a square capillary, the surrounding liquid flows through a thin film between the slug bubble and the wall and the corners of the channel and generates an asymmetric two-phase flow field with significant spatial variations. The current literature mainly discusses the bulk features of such a flow field, while local phenomena in these flows are not well known yet. In this study, the local dynamics of the liquid film and corner flows in a flowing square channel with a width of 3 mm were investigated by experiments and numerical simulations. Instantaneous velocity fields in the liquid flow were measured at the channel's central plane using particle tracking velocimetry (PTV). The obtained data were utilized to calculate the axial variations of the film and corner velocities. Three-dimensional numerical simulations with adaptive mesh refinement of the interface were performed for a series of flow conditions. Comparing the film velocity obtained from PTV and simulations demonstrates an overall good agreement between the results, with discrepancies close to the solid walls, where the PTV data is sporadic. Film and corner flows show rapid velocity changes adjacent to the minimum film area. Both experiments and simulations elucidate the backflow phenomena in the analyzed flow fields, where the liquid locally flows in the opposite streamwise direction. The backflow region shrinks as the liquid flow rate and speed of the bubble increase. Film flow shows significant pressure drops as it passes through the minimum film area. On the contrary, the corner flow does not demonstrate any rapid pressure fluctuations. The vortex structures’ identification revealed dominant vortex loops inside the gas and small-scale vortex filaments along the channel walls, generated from the liquid film rupture.

Interaction and influence of a flow field and particleboard particles in an airflow forming machine with a coupled Euler-DPM model.

Particleboards are widely used in the artificial board market, which can be constructed from a variety of raw materials and require small amounts of energy to be produced. In the particleboard production process, forming machines play an important role as the key equipment for achieving continuous production. In recent years, airflow forming machines have received increasing attention in particleboard production lines because of their strong separation ability and low price. However, the internal flow field is complex and difficult to control, which affects the surface quality and strength of the particleboard. The most pressing technical difficulty is controlling the flow field characteristics of the airflow paver. At present, the research on this subject is conducted primarily through repeated experiments, which entail long research periods and high processing costs. To reduce human and financial costs, in this study, Computational Fluid Dynamics (CFD) is employed to investigate the flow field and the gas-solid two-phase flow field coupled with particle movement of an airflow forming machine. The accuracy of the calculation model is verified by comparing characteristic point velocities obtained from experimental analysis and a simulation. The simulation results show that in practical production, the frequency of a negative pressure fan should be greater than 27 Hz. It is necessary to set the shoulder properly, and the slab smoothness can be improved by moving the shoulder back on the premise of meeting the strength requirements of the box. The distance between the shoulders of the box body should be less than 2570 mm, and particles with uniform diameter should be added to the paving box to reduce the turbulence effect, improve the quality of particle forming and provide actual particleboard production with a solid theoretical foundation.

Parametric study of gas−liquid two-phase flow field in horizontal stirred tank

Based on Fluent software, the gas−liquid two-phase flow in the horizontal stirred tank was simulated with SST k−ω turbulence model, Eulerian−Eulerian two-fluid model, and multi-reference flame method. The mixing process in the tank was calculated by tracer method. The results show that increasing the rotating speed or gas flow is conducive to a more uniform distribution of the gas phase and accelerates the mixing of the liquid phase. When the rotating speed exceeds 93 r/min, the relative power demand remains basically constant. The change in the inclination angle of the upper impeller has minimal effect on the gas phase distribution. When the inclination angle is 50°, the relative power demand reaches the maximum. An appropriate increase in the impeller distance from the bottom improves the gas holdup and gas phase distribution but increases the liquid phase mixing time.

The Influence of Different Blade Outlet Angle and Wedge Groove on the Performance of the Double Suction Pump

This study aims to significantly improve the performance of the centrifugal pump by designing a wedge grooves on the outside of the impeller cover, and investigating the method for determining the blade outlet angle that plays important roles in the energy conversion. The solidliquid two-phase flow fields with the designed pump and an original pump were numerically simulated in the ANSYS-FLUENT. Compared with the original pump, reducing the blade outlet angle can effectively improve the sand-water wear on the blade surface and opening wedge grooves can slow down impeller ring wear. The head of the improved model pump decreases by about 1.5% and the efficiency increases by about 1.5%.

Inferring incompressible two-phase flow fields from the interface motion using physics-informed neural networks

In this work, physics-informed neural networks are applied to incompressible two-phase flow problems. We investigate the forward problem, where the governing equations are solved from initial and boundary conditions, as well as the inverse problem, where continuous velocity and pressure fields are inferred from scattered-time data on the interface position. We employ a volume of fluid approach, i.e. the auxiliary variable here is the volume fraction of the fluids within each phase. For the forward problem, we solve the two-phase Couette and Poiseuille flow. For the inverse problem, three classical test cases for two-phase modeling are investigated: (i) drop in a shear flow, (ii) oscillating drop and (iii) rising bubble. Data of the interface position over time is generated by numerical simulation. An effective way to distribute spatial training points to fit the interface, i.e. the volume fraction field, and the residual points is proposed. Furthermore, we show that appropriate weighting of losses associated with the residual of the partial differential equations is crucial for successful training. The benefit of using adaptive activation functions is evaluated for both the forward and inverse problem.

On the three-dimensional features of a confined slug bubble in a flowing square capillary

The motion of long bubbles in tubular capillaries has typically been described by bulk characteristics. However, the dynamics of slug bubbles in square capillaries are more complex due to a corner flow and a thin film flow. The physics can be correctly explained by elucidating local 3D features of the two-phase flow field. To this aim, an experimental study based on particle tracking velocimetry (PTV) and a numerical simulation based on the volume-of-fluid method were conducted to investigate the dynamics of slug bubbles rising in a flowing square capillary with a cross-sectional area of 3 × 3 mm2. To precisely analyze the phases' interaction, interfacial flow data were mapped onto a radial-tangential coordinate system on central and diagonal planes. The simulated interface topology and velocity fields show a good agreement with the experimental PTV data on the central plane, with an absolute error of less than 1.2% for terminal bubble speed. Tangential speeds show their maxima occurring in the channel corners, where pressure is maximum. The thin liquid film flow that occurs where the bubble approaches the wall applies noticeable shear stress on the channel walls, where high and low-pressure regions are generated. Structures of vortices inside the bubble were identified using isosurfaces of the Q-criterion, and their cores were detected based on the parallel vector method. Results reveal a dominant vortex ring adjacent to the liquid film flow and two oblique vortex tubes close to the bubble's nose.

Coaxial helical gas assisted laser water jet machining of SiC/SiC ceramic matrix composites

• A novel coaxial helical gas atmosphere is firstly introduced to extend the stable length of water jet as well as expel the accumulated water layer during the machining of LWJ. • The influences of gas component and pressure on the stable length of water jet and accumulated water layer status on substrate surface are revealed. • Scribing and through cutting of CMCs substrates with thickness of 3 mm are realized by CGALWJ machining without the drawing of SiC fibers, formation of recast layer and delamination. SiC/SiC ceramic matrix composites (CMCs) offer an excellent combination of properties at high temperature such as high specific strength, chemical inertness and irradiation tolerance. Those superior properties make CMCs beneficial for use in high-temperature structural applications that are exposed to extreme environments such as aerospace and nuclear energy. However, well machining qualities can hardly be achieved by conventional machining techniques owing to these properties. Laser water jet (LWJ) machining is a promising solution, which is capable of ablating materials with less/no heat defects, well machining precision and consistency. Nevertheless, the machining capacity of LWJ is still limited by the stability of water jet to a great extent. A water layer may form on substrate surface during the impingement of LWJ, which also sets up obstacle for sufficient ablation. Therefore, a novel coaxial helical gas atmosphere is introduced to promote the machining capacity of LWJ in this paper. A theoretical model is established to describe the gas-water two-phase flow field during the ejection and impingement of coaxial gas assisted LWJ (CGALWJ). The influences of gas component and pressure on the stable length of water jet and surface water layer status are analyzed based on numerical simulations and experiments. Scribing experiments are further carried out on CMCs substrates with thickness of 3 mm. Groove with maximum depth-to-width ratio of 13.6 as well as through cutting are realized without the drawing of SiC fibers, formation of recast layer and delamination. The theoretical and experimental results provide solid foundation for the high-quality machining of ceramic matrix composites and other hard-to-process materials.

Experimental Study on the Performance of the Gas Centrifugal Seal with Continuous Injection and Discharge of Sealing Fluid

The centrifugal seal has the advantages of simple structure, long service life and the cost is far lower than the dry gas seal. But under high-speed gas conditions, the temperature of the sealing fluid rises too high and low sealing pressure. This paper studies the gas centrifugal seal with continuous injection and discharge of sealing fluid, establishes a numerical model, simulates the gas-liquid two-phase flow field and temperature field in the sealing chamber. Revealing the effects of the speed and the sealing fluid flux on the sealing performance by comparing the changes in the sealing fluid pressure, temperature, and sealing capacity under different operating conditions. And the reliability of its application to gas mechanical seals is verified through experiments. The results show that at high speeds, the pumping compression effect of the sealing fluid is enhanced, and the sealing capability is improved. Within a certain range, increasing the flux of the sealing fluid can effectively reduce the temperature and prevent the sealing fluid from vaporizing, the continuous injection and discharge structure has a significant effect.

Numerical and Experimental Analysis of Transient 2D Flow around Ship Spoiler.(Dept.M)

An experimental analysis and numerical computation were conducted for flow around ship spoiler mounted at the bottom of the ship and subjected to air injection The air was injected beside the spoiler and at a distance equals the spoiler length to stabilize cavities behind spoilers. This system was proposed to control the roll, pitch motion and speed of the ship, with injecting exhaust gas behind spoilers. The forces and moments acting on the ship's bottom depend on spoilers' inclination, extensions and positions relative to the ship's center of mass. To improve the understanding of the parameters affecting the flow field aroundsuch spoiler, experimental images of bubble formation are triggered and compared with the computed flow field at different spoiler inclination angle, and air injection position. Two dimensional Navier-Stokes code is used to model the two-phase flow field around a ship spoiler with the free surface simulation in Piecewise Linear Interface Construction (PLIC) method. The governing equations are discretized on a structured grid using an upwind difference scheme. For different conditions, the bubbles shape, the two-dimensional flow field around the spoiler body and the pressure variation on the wake of the spoiler body are computed, Flow field and bubble formation around the spoiler are recorded experimentally at different condition and time sequence with scientific video camera The Froude number was about 1.6 based on spoiler length. The computational and experimental results show that there is a good matching of bubble formation in the spoiler with 90', but in the condition of 45° the experimental images show that bubble filled the space between the bubble and the body of the ship. This is may be due to the two-dimension computational model consideration. So, a water layer was isolated between the bubble and the ship body. In the case of spoiler with 45° inclination angle, the flow is circulated inside the bubble. The bubble shape is longer which gives more negative pressure on larger area. The moment around the 450 spoiler fixationpoint is seven times of that around the 90° spoiler. In double injection condition beside and at a distance equal the horizontal spoiler length downstream the spoiler, the bubbles are more stable and the moment is higher by 7.9 and 17.8 times after 1.2 and 2 seconds, respectively, than that around the 90° spoiler. Thus, changing inclination spoiler angle and position relative to injection angle will give different bubble shapes which give different forces to control the roll, pitch motion and speed of the ship.

Research Progress on Numerical Simulation of Two-phase Flow in the Gas-solid Fluidized Bed

It is difficult to accurately measure the parameters of solid particles in the experiment of the gas-solid fluidized bed. The numerical simulation plays an important role to accurately describe flow characteristics in the fluidized bed. Combined with the research work of the research group, this paper analyzes the application of numerical simulation of fluidized bed from the aspects of gas-solid coupling algorithm, drag model, flow characteristics, and reaction characteristics based on the previous studies. The specificity improvement of the gas-solid coupling algorithm and the regional application of the drag model is the trend of the recent development of numerical simulation. Previous studies mainly focus on the gas-solid two-phase flow field characteristics in the traditional fluidized bed, but few on the complex flow characteristics such as gas-solid reverse flow and the coupling with reaction characteristics. It is of great significance for designing a novel fluidized bed reactor to realize gas-solid continuous reaction to establish and improve the numerical simulation method of gas-solid non-catalytic reaction.