Granular bed gas-solid two-phase flow numerical simulation method based on Euler-Euler model

Abstract

The invention discloses a particle bed gas-solid two-phase flow numerical simulation method based on an Euler-Euler model, and relates to a particle bed gas-solid two-phase flow numerical simulation method. The particle bed gas-solid two-phase flow numerical simulation method based on the Euler-Euler model aims to solve the problem that a vacuum system of a vacuum induction melting furnace is prone to being damaged or production efficiency is low in the existing spherical metal powder recycling and reusing process, and the particle bed gas-solid two-phase flow numerical simulation method based on the Euler-Euler model is provided. According to the method, the influence of the vacuumizing speed on the particle motion trail is fully considered, and meanwhile, the maximum flight speed of the particles at different vacuumizing speeds is determined through curve fitting; and the motion trails of the particles at different vacuumizing speeds in the induction melting furnace are researched, so that the production efficiency is improved on the premise of not damaging a vacuum system, and a simulation result has guiding significance on a vacuumizing process before powder recovery and remelting.

Description

Technical Field

[0001] The invention relates to a numerical simulation method for gas-solid two-phase flow in a particle bed. Background Art

[0002] As the most widely used raw material in powder metallurgy, the common production processes of spherical metal powders are vacuum induction melting inert gas atomization, plasma rotating electrode atomization, plasma wire atomization, etc. Restricted by the principle of powder preparation process, the particle size range of powder particles is normally distributed (particle size range is 0.1μm~200μm). Different powder metallurgy processes require different particle size ranges. Powders with too small and too large particle sizes cannot be used in powder metallurgy processes. About 20% to 40% of powders are wasted. Therefore, the recycling and reuse of spherical metal powders whose particle sizes do not meet actual production needs has become an urgent problem to be solved. The recycling and reuse of spherical metal powders can effectively reduce production costs.

[0003] The process of directly induction heating a mixture of spherical metal powder and master alloy, and then atomizing and spraying the melt with inert gas is a new type of powder recycling production method in China. In order to ensure that the melt is not oxidized before the atomization process, the mixed material must be inductively heated in a vacuum environment.

[0004] Previous studies did not consider the effect of the coupling between the complex cavity and the gas flow on the particle flow when exploring the movement law of the gas-solid two-phase flow in the particle bed. In the vacuum process, excessive vacuum speed will increase the gas flow rate in the induction melting furnace. The fine powder particles accumulated on the surface of the particle bed in the melting crucible will be dragged out of the crucible under the action of the gas drag force. Some powder particles will enter the vacuum system with the gas, thereby destroying the vacuum system. Therefore, too fast a vacuum speed will destroy the vacuum system, while too slow a vacuum speed will reduce production efficiency. There is no report in the prior art on the method of numerical simulation of the gas-solid two-phase flow of the particle bed during the vacuum process before the powder recovery remelting process. Summary of the invention

[0005] In order to solve the problem that the existing recycling and reuse process of spherical metal powders easily causes damage to the vacuum system of the vacuum induction melting furnace or low production efficiency, the present invention proposes a numerical simulation method of gas-solid two-phase flow in a particle bed based on the Euler-Euler model. The present invention fully considers the influence of the vacuum pumping speed on the particle motion trajectory, and determines the maximum flight speed of the particles at different vacuum pumping speeds through curve fitting; and studies the particle motion trajectory at different vacuum pumping speeds in the induction melting furnace, so as to improve the production efficiency without damaging the vacuum system. The simulation results have guiding significance for the vacuum pumping process before powder recycling and remelting.

[0006] The present invention is based on the Euler-Euler model of the particle bed gas-solid two-phase flow numerical simulation method according to the following steps:

[0007] Step 1: Use computer-aided design software to establish a two-dimensional model of the gas-solid two-phase flow area of ​​the particle bed in the induction melting furnace; establish a numerical simulation calculation domain of the gas-solid two-phase flow of the particle bed in the vacuum process of the induction melting furnace according to the location of each component in the induction melting furnace, the outer dimensions of the induction melting furnace and the vacuum port dimensions of the induction melting furnace;

[0008] Step 2: meshing the calculation domain of the gas-solid two-phase flow numerical simulation of the granular bed in the vacuum process of the induction melting furnace using meshing software, and encrypting the mesh in the area where the gas and the granular bed particles are adjacent;

[0009] Step 3, using a computational fluid dynamics solver to solve the numerical simulation of the gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace;

[0010] The solution settings include: using pressure-based transient solution; defining gravity and direction; using the Euler-Euler two-phase flow model combined with the SST turbulence model for numerical simulation calculations; defining the phases, material properties, and phase interactions of each calculation area; setting the vacuum port exhaust velocity boundary conditions according to the vacuum pumping process to initialize the physical conditions of the calculation domain; using the SIMPLE algorithm for pressure-velocity coupling; setting the convergence residual, calculation time step, and total calculation time;

[0011] In step 3, the number of Euler phases of the Euler-Euler two-phase flow model is 2, the primary phase is set to air, and the secondary phase is set to granular material;

[0012] In step 3, the material properties include particle diameter, particle viscosity parameter, and packing density; the particle diameter is 10 μm and 100 μm; the particle viscosity adopts the Lun-et-al model; the packing density is 0.55;

[0013] The pressure-velocity coupling in step 3 uses the SIMPLE algorithm and sets the convergence residual to 1e -5 , the calculation time step is 1e -4 s, calculate the total time 10s;

[0014] The initialization setting of the physical conditions in the computational domain in step 3 includes: setting the boundary condition to a fixed wall with a no-slip wall, setting the inlet to a velocity inlet, and setting the volume fraction value of air to 1;

[0015] Step 4: Use computational fluid dynamics post-processing software to output numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed; the numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed include a volume flow velocity cloud map, a particle motion velocity cloud map, and a particle maximum flight velocity fitting curve.

[0016] The beneficial effects of the present invention are:

[0017] In the prior art, the experimental method is usually used to study the gas-solid two-phase flow of the granular bed in the vacuum process in the induction melting furnace, which is easy to damage the vacuum system and has a high economic cost. The present invention uses the Euler-Euler two-phase flow model in the tracking of the interface of the granular bed gas-solid two-phase flow, which has the advantages of low computational cost and high computational efficiency compared with the Euler-Lagrange method. Combined with the SST turbulence model, more accurate transient gas flow characteristics can be obtained, including the position and structure of gas vortices, the instantaneous velocity distribution of the gas flow field, etc., while overcoming the problem of huge computing power requirements brought about by direct numerical simulation due to the need to solve all turbulent information. The numerical simulation method of gas-solid two-phase flow in a granular bed based on the Euler-Euler model provided by the present invention can accurately simulate the movement of particles in the vacuum process, clarify the influence of factors such as vacuum speed and particle diameter on the particle movement trajectory, and the simulation process is close to the actual situation. Compared with the research method of experimental in-situ measurement, it has the advantages of not destroying the vacuum system, low economic cost and high research efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] Figure 1 This is a grid division diagram of the calculation domain for numerical simulation of gas-solid two-phase flow in a granular bed in a vacuum process in an induction melting furnace in Example 1;

[0019] Figure 2 This is the gas flow field velocity cloud diagram when the boundary condition of the vacuum port exhaust velocity is 0.2 m / s in Example 1;

[0020] Figure 3 This is the gas flow field velocity cloud diagram when the boundary condition of the vacuum port exhaust velocity is 2 m / s in Example 1;

[0021] Figure 4 This is the particle flow field velocity cloud diagram when the boundary condition of the vacuum port exhaust velocity is 0.2 m / s in Example 1;

[0022] Figure 5 This is the particle flow field velocity cloud diagram when the boundary condition of the vacuum port exhaust velocity is 2 m / s in Example 1;

[0023] Figure 6 This is a fitting curve diagram of the maximum flight speed of particles in Example 1;

[0024] Figure 7 This is a velocity cloud diagram of the particle flow field in the vacuuming process when the particle diameter is 100 μm in Example 1;

[0025] Figure 8 This is a particle flow field velocity cloud diagram of the vacuuming process in Example 1 with a particle diameter of 10 μm. DETAILED DESCRIPTION

[0026] The technical solution of the present invention is not limited to the specific implementation modes listed below, but also includes any reasonable combination of the specific implementation modes.

[0027] Specific implementation method 1: The numerical simulation method of gas-solid two-phase flow in a particle bed based on the Euler-Euler model in this implementation method is carried out according to the following steps:

[0028] Step 1: Use computer-aided design software to establish a two-dimensional model of the gas-solid two-phase flow area of ​​the particle bed in the induction melting furnace; establish a numerical simulation calculation domain of the gas-solid two-phase flow of the particle bed in the vacuum process of the induction melting furnace according to the location of each component in the induction melting furnace, the outer dimensions of the induction melting furnace and the vacuum port dimensions of the induction melting furnace;

[0029] Step 2: meshing the calculation domain of the gas-solid two-phase flow numerical simulation of the granular bed in the vacuum process of the induction melting furnace using meshing software, and encrypting the mesh in the area where the gas and the granular bed particles are adjacent;

[0030] Step 3, using a computational fluid dynamics solver to solve the numerical simulation of the gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace;

[0031] The solution settings include: using pressure-based transient solution; defining gravity and direction; using the Euler-Euler two-phase flow model combined with the SST turbulence model for numerical simulation calculations; defining the phases, material properties, and phase interactions of each calculation area; setting the vacuum port exhaust velocity boundary conditions according to the vacuum pumping process to initialize the physical conditions of the calculation domain; using the SIMPLE algorithm for pressure-velocity coupling; setting the convergence residual, calculation time step, and total calculation time;

[0032] In step 3, the number of Euler phases of the Euler-Euler two-phase flow model is 2, the primary phase is set to air, and the secondary phase is set to granular material;

[0033] In step 3, the material properties include particle diameter, particle viscosity parameter, and packing density; the particle diameter is 10 μm and 100 μm; the particle viscosity adopts the Lun-et-al model; the packing density is 0.55;

[0034] The pressure-velocity coupling in step 3 uses the SIMPLE algorithm and sets the convergence residual to 1e -5 , the calculation time step is 1e -4 s, calculate the total time 10s;

[0035] The initialization setting of the physical conditions in the computational domain in step 3 includes: setting the boundary condition to a fixed wall with a no-slip wall, setting the inlet to a velocity inlet, and setting the volume fraction value of air to 1;

[0036] Step 4: Use computational fluid dynamics post-processing software to output numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed; the numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed include a volume flow velocity cloud map, a particle motion velocity cloud map, and a particle maximum flight velocity fitting curve.

[0037] This implementation has the following beneficial effects:

[0038] In the prior art, the experimental method is usually used to study the gas-solid two-phase flow of the particle bed in the vacuum process in the induction melting furnace, which is easy to damage the vacuum system and has a high economic cost. This embodiment uses the Euler-Euler two-phase flow model in the tracking of the gas-solid two-phase flow interface in the particle bed, which has the advantages of low computational cost and high computational efficiency compared to the Euler-Lagrangian method. Combined with the SST turbulence model, more accurate transient gas flow characteristics can be obtained, including the position and structure of gas vortices, the instantaneous velocity distribution of the gas flow field, etc., while overcoming the problem of huge computing power requirements brought about by direct numerical simulation due to the need to solve all turbulent information. The numerical simulation method of gas-solid two-phase flow in a particle bed based on the Euler-Euler model provided in this embodiment can accurately simulate the movement of particles in the vacuum process, clarify the influence of factors such as vacuum speed and particle diameter on the particle movement trajectory, and the simulation process is close to the actual situation. Compared with the research method of experimental in-situ measurement, it has the advantages of not destroying the vacuum system, low economic cost and high research efficiency.

[0039] Specific implementation method 2: This implementation method is different from specific implementation method 1 in that: the components in the induction melting furnace described in step 1 include a furnace body, a vacuum port, a crucible and granular raw materials.

[0040] Specific implementation method three: This implementation method is different from specific implementation methods one or two in that: the computer-aided design software described in step 1 is Ansys SpaceClaim.

[0041] Specific implementation method 4: This implementation method is different from any one of specific implementation methods 1 to 3 in that: the outer dimensions of the induction melting furnace in step 1 are 2m and 1.8m, and the vacuum port diameter of the induction melting furnace is 0.8m.

[0042] Specific implementation method five: This implementation method is different from any one of specific implementation methods one to four in that: the meshing software described in step 2 is Ansys meshing.

[0043] Specific implementation method 6: This implementation method is different from any one of specific implementation methods 1 to 5 in that: the computational fluid dynamics solver in step 3 is Ansys fluent.

[0044] Specific implementation method seven: This implementation method is different from any one of specific implementation methods one to six in that: the computational fluid dynamics solver described in step 3 establishes a mathematical model of gas-solid two-phase flow in a particle bed of the Euler-Euler model based on the turbulence equations and the dual Euler control equations.

[0045] Specific implementation eight: This implementation differs from any one of specific implementations one to seven in that: the turbulence equation is:

[0046]

[0047] Where: α g is the gas phase fraction, ρ g is the gas phase density; k is the turbulent kinetic energy; ω is the specific dissipation rate; U g is the gas phase velocity; Γ k is the effective diffusion rate of turbulent kinetic energy; Γ ω is the effective diffusion rate relative to the dissipation rate; G k is the generation term of turbulent kinetic energy; W k is the dissipation term of turbulent kinetic energy; S k is the source term of turbulent kinetic energy; G ω is the generation term of specific dissipation rate; W ω is the dissipation term of the specific dissipation rate; S ω is the source term of the specific dissipation rate.

[0048] Specific implementation method 9: This implementation method is different from any one of the specific implementation methods 1 to 8 in that: the dual Euler control equation includes a continuity equation and a momentum equation:

[0049] The continuity equation is:

[0050]

[0051] Where: α q is the phase volume fraction; ρ q is the q-phase density vector; is the q-phase velocity vector; is the interphase mass transfer quality from phase p to phase q; is the mass transfer from phase q to phase p; t is the time;

[0052] The momentum equation is:

[0053]

[0054] Where: α q is the phase volume fraction; ρ q is the q-phase density vector; is the q-phase velocity vector; Pq is the pressure term; is the q-phase pressure; is the force-strain tensor corresponding to q; is the mass force; is the interphase force; is the mass transfer rate from phase p to phase q; is the interphase mass transfer rate from phase q to phase p; is the additional force of phase q; is the interphase mass transfer quality from phase p to phase q; is the interphase mass transfer mass from phase q to phase p.

[0055] Specific implementation method ten: This implementation method is different from any one of specific implementation methods one to nine in that: the computational fluid dynamics post-processing software is CFD-POST.

[0056] Embodiment 1:

[0057] The numerical simulation method of gas-solid two-phase flow in a particle bed based on the Euler-Euler model in this embodiment is carried out according to the following steps:

[0058] Step 1: Use computer-aided design software to establish a two-dimensional model of the gas-solid two-phase flow area of ​​the particle bed in the induction melting furnace; establish a numerical simulation calculation domain of the gas-solid two-phase flow of the particle bed in the vacuum process of the induction melting furnace according to the location of each component in the induction melting furnace, the outer dimensions of the induction melting furnace and the vacuum port dimensions of the induction melting furnace;

[0059] The components in the induction melting furnace include a furnace body, a vacuum port, a crucible and granular raw materials;

[0060] The computer-aided design software is Ansys SpaceClaim;

[0061] The outer dimensions of the induction melting furnace are 2m and 1.8m, and the vacuum port diameter of the induction melting furnace is 0.8m;

[0062] Step 2: meshing the calculation domain of the gas-solid two-phase flow numerical simulation of the granular bed in the vacuum process of the induction melting furnace using meshing software, and encrypting the mesh in the area where the gas and the granular bed particles are adjacent;

[0063] The meshing software is Ansys meshing;

[0064] Step 3, using a computational fluid dynamics solver to solve the numerical simulation of the gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace;

[0065] The solution settings include: using pressure-based transient solution; defining gravity and direction; using the Euler-Euler two-phase flow model combined with the SST turbulence model for numerical simulation calculations; defining the phases, material properties, and phase interactions of each calculation area; setting the vacuum port exhaust velocity boundary conditions according to the vacuum pumping process to initialize the physical conditions of the calculation domain; using the SIMPLE algorithm for pressure-velocity coupling; setting the convergence residual, calculation time step, and total calculation time;

[0066] The computational fluid dynamics solver in step 3 is Ansys fluent;

[0067] Step 3: The number of Euler phases in the Euler-Euler two-phase flow model is 2, the primary phase is set to air, and the secondary phase is set to granular material;

[0068] In step 3, the material properties include particle diameter, particle viscosity parameter, and packing density; the particle diameter is 10 μm and 100 μm; the particle viscosity adopts the Lun-et-al model; the packing density is 0.55;

[0069] The pressure-velocity coupling in step 3 uses the SIMPLE algorithm and sets the convergence residual to 1e -5 , the calculation time step is 1e -4 s, calculate the total time 10s;

[0070] The initialization settings of the physical conditions in the computational domain described in step 3 include: setting the boundary conditions as a fixed wall using a no-slip wall, setting the inlet as a velocity inlet, and setting the volume fraction value of the air to 1; the initialization conditions are set according to the actual working conditions.

[0071] In step 3, the computational fluid dynamics solver establishes a particle bed gas-solid two-phase flow mathematical model of the Euler-Euler model based on the turbulence equations and the double Euler control equations;

[0072] The gas flow rate during the vacuuming process is fast and the flow domain structure is complex, so the turbulence equation is used to solve the gas flow behavior. The turbulence equation is:

[0073]

[0074]

[0075] Where: α g is the gas phase fraction, ρ g is the gas phase density; k is the turbulent kinetic energy; ω is the specific dissipation rate; U g is the gas phase velocity; Γ k is the effective diffusion rate of turbulent kinetic energy; Γ ω is the effective diffusion rate relative to the dissipation rate; G kis the generation term of turbulent kinetic energy; W k is the dissipation term of turbulent kinetic energy; S k is the source term of turbulent kinetic energy; G ω is the generation term of specific dissipation rate; W ω is the dissipation term of the specific dissipation rate; S ω is the source term of specific dissipation rate;

[0076] The dual Euler control equations include the continuity equation and the momentum equation:

[0077] The continuity equation is:

[0078]

[0079] Where: α q is the phase volume fraction; ρ q is the q-phase density vector; is the q-phase velocity vector; is the interphase mass transfer quality from phase p to phase q; is the mass transfer from phase q to phase p; t is the time;

[0080] The momentum equation is:

[0081]

[0082] Where: α q is the phase volume fraction; ρ q is the q-phase density vector; is the q-phase velocity vector; Pq is the pressure term; is the q-phase pressure; is the force-strain tensor corresponding to q; is the mass force; is the interphase force; is the mass transfer rate from phase p to phase q; is the interphase mass transfer rate from phase q to phase p; is the additional force of phase q; is the interphase mass transfer quality from phase p to phase q; is the interphase mass transfer quality from phase q to phase p;

[0083] Step 4: Use computational fluid dynamics post-processing software to output numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed; the numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed include a volume flow velocity cloud map, a particle motion velocity cloud map, and a particle maximum flight velocity fitting curve.

[0084] The computational fluid dynamics post-processing software is CFD-POST.

[0085] Depend on Figure 2-6It can be seen that as the vacuuming speed increases, the instantaneous speed of the gas flow field in the furnace increases, and the gas forms a stable gas vortex on the upper wall of the vacuum port and the outer wall of the crucible, but the structure and position of the gas flow field vortex are less affected. As the vacuuming speed increases, the gas flow rate in the furnace increases, the drag force of the gas on the particles increases, and the movement speed of the particles increases. When the vacuuming rate exceeds 2m / s, the particles will be dragged out of the crucible by the gas. Figure 7-8 It can be seen that the larger the particle size, the weaker the dragging effect of the gas on the particles, and the less likely the particles are to be dragged out of the crucible.

Claims

1. A numerical simulation method for gas-solid two-phase flow in a particle bed based on the Euler-Euler model, characterized in that: The numerical simulation method of gas-solid two-phase flow in a particle bed based on the Euler-Euler model is carried out in the following steps: Step 1: Use computer-aided design software to establish a two-dimensional model of the gas-solid two-phase flow area of ​​the particle bed in the induction melting furnace; establish a numerical simulation calculation domain of the gas-solid two-phase flow of the particle bed in the vacuum process of the induction melting furnace according to the location of each component in the induction melting furnace, the outer dimensions of the induction melting furnace and the vacuum port dimensions of the induction melting furnace; Step 2: meshing the calculation domain of the gas-solid two-phase flow numerical simulation of the granular bed in the vacuum process of the induction melting furnace using meshing software, and encrypting the mesh in the area where the gas and the granular bed particles are adjacent; Step 3, using a computational fluid dynamics solver to solve the numerical simulation of the gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace; The solution settings include: using pressure-based transient solution; defining gravity and direction; using the Euler-Euler two-phase flow model combined with the SST turbulence model for numerical simulation calculations; defining the phases, material properties, and phase interactions of each calculation area; setting the vacuum port exhaust velocity boundary conditions according to the vacuum pumping process to initialize the physical conditions of the calculation domain; using the SIMPLE algorithm for pressure-velocity coupling; setting the convergence residual, calculation time step, and total calculation time; In step 3, the number of Euler phases of the Euler-Euler two-phase flow model is 2, the primary phase is set to air, and the secondary phase is set to granular material; The material properties in step 3 include particle diameter, particle viscosity parameter, and packing density; the particle diameter is 10 μm and 100 μm; the particle viscosity adopts the Lun-et-al model; the packing density is 0.55; The pressure-velocity coupling in step 3 uses the SIMPLE algorithm and sets the convergence residual to 1e -5 , the calculation time step is 1e -4 s, calculate the total time 10s; The initialization setting of the physical conditions in the computational domain in step 3 includes: setting the boundary condition to a fixed wall with a no-slip wall, setting the inlet to a velocity inlet, and setting the volume fraction value of air to 1; Step 4: Use computational fluid dynamics post-processing software to output numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed; the numerical simulation results of the motion characteristics of the gas-solid two-phase flow in the particle bed include a volume flow velocity cloud map, a particle motion velocity cloud map, and a particle maximum flight velocity fitting curve.

2. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The components in the induction melting furnace in step 1 include a furnace body, a vacuum port, a crucible and granular raw materials.

3. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The computer-aided design software described in step 1 is Ansys SpaceClaim.

4. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The outer dimensions of the induction melting furnace in step 1 are 2m and 1.8m, and the vacuum port diameter of the induction melting furnace is 0.8m.

5. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The meshing software described in step 2 is Ansys meshing.

6. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The computational fluid dynamics solver in step 3 is Ansys Fluent.

7. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The computational fluid dynamics solver described in step 3 establishes a mathematical model of gas-solid two-phase flow in a particle bed based on the Euler-Euler model based on the turbulence equations and the dual Euler control equations.

8. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The turbulence equation is: Where: α g is the gas phase fraction, ρ g is the gas phase density; k is the turbulent kinetic energy; ω is the specific dissipation rate; U g is the gas phase velocity; Γ k is the effective diffusion rate of turbulent kinetic energy; Γ ω is the effective diffusion rate relative to the dissipation rate; G k is the generation term of turbulent kinetic energy; W k is the dissipation term of turbulent kinetic energy; S k is the source term of turbulent kinetic energy; G ω is the generation term of specific dissipation rate; W ω is the dissipation term of the specific dissipation rate; S ω is the source term of the specific dissipation rate.

9. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The dual Euler control equations include the continuity equation and the momentum equation: The continuity equation is: Where: α q is the phase volume fraction; ρ q is the q-phase density vector; is the q-phase velocity vector; is the interphase mass transfer quality from phase p to phase q; is the mass transfer from phase q to phase p; t is the time; The momentum equation is: Where: α q is the phase volume fraction; ρ q is the q-phase density vector; is the q-phase velocity vector; Pq is the pressure term; is the q-phase pressure; is the force-strain tensor corresponding to q; is the mass force; is the interphase force; is the mass transfer rate from phase p to phase q; is the interphase mass transfer rate from phase q to phase p; is the additional force of phase q; is the interphase mass transfer quality from phase p to phase q; is the interphase mass transfer mass from phase q to phase p.

10. The method for numerical simulation of gas-solid two-phase flow in a particle bed based on the Euler-Euler model according to claim 1, characterized in that: The computational fluid dynamics post-processing software is CFD-POST.

Images