A numerical simulation method for gas-solid two-phase flow in granular beds based on the Euler-Euler model

Abstract

This invention discloses a numerical simulation method for gas-solid two-phase flow in a granular bed based on the Euler-Euler model. To address the problem that existing processes for recycling and reusing spherical metal powders easily damage the vacuum system of vacuum induction melting furnaces or result in low production efficiency, this invention proposes a numerical simulation method for gas-solid two-phase flow in a granular bed based on the Euler-Euler model. This invention fully considers the influence of vacuuming speed on particle trajectory, and determines the maximum flight speed of particles at different vacuuming speeds through curve fitting. Furthermore, it studies the particle trajectory at different vacuuming speeds within the induction melting furnace, achieving improved production efficiency without damaging the vacuum system. The simulation results provide guidance for the vacuuming process before powder recycling and remelting.

Description

Technical Field

[0001] This invention relates to a numerical simulation method for gas-solid two-phase flow in a granular bed. Background Technology

[0002] As the most widely used raw material in powder metallurgy, spherical metal powder is commonly produced through processes such as vacuum induction melting with inert gas atomization, plasma rotating electrode atomization, and plasma wire atomization. Due to the limitations of powder preparation processes, the particle size range of powder particles exhibits a normal distribution (0.1 μm to 200 μm). Different powder metallurgy processes require different particle size ranges; powders that are too small or too large cannot be used in powder metallurgy, resulting in approximately 20% to 40% waste. Therefore, the recycling of spherical metal powders with particle sizes that do not meet actual production requirements is 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 a master alloy, followed by inert gas atomization and powder spraying of the melt, is a novel powder recycling production method in China. To ensure that the melt is not oxidized before the atomization process, the mixture must be induction heated in a vacuum environment.

[0004] Previous studies investigating the motion of gas-solid two-phase flow in granular beds have not considered the influence of the coupling effect between complex cavities and gas flow on particle flow. In vacuuming processes, excessively high vacuuming speeds lead to increased gas flow velocity within the induction melting furnace. Fine powder particles accumulated on the surface of the granular bed in the melting crucible are dragged out of the crucible by the gas drag force, and some powder particles enter the vacuum system with the gas, thus disrupting the vacuum system. Therefore, excessively high vacuuming speeds can damage the vacuum system, while excessively low vacuuming speeds reduce production efficiency. No numerical simulation methods for gas-solid two-phase flow in granular beds during the vacuuming process before powder recycling and remelting have been reported in the current technology. Summary of the Invention

[0005] To address the problems of potential damage to the vacuum system or low production efficiency in existing spherical metal powder recycling processes, this invention proposes a numerical simulation method for gas-solid two-phase flow in a particle bed based on the Euler-Euler model. This invention fully considers the influence of vacuuming speed on particle trajectory and determines the maximum flight speed of particles at different vacuuming speeds through curve fitting. Furthermore, it studies the particle trajectory at different vacuuming speeds within the induction melting furnace, thereby improving production efficiency without damaging the vacuum system. The simulation results provide guidance for the vacuuming process before powder recycling and remelting.

[0006] The numerical simulation method for gas-solid two-phase flow in granular beds based on the Euler-Euler model of this invention is carried out 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 region of the particle bed in the induction melting furnace; establish the 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 external dimensions of the induction melting furnace and the size of the vacuum port of the induction melting furnace.

[0008] Step 2: Use mesh generation software to mesh the computational domain of the numerical simulation of gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace, and refine the mesh in the region where the gas and the particles in the particle bed are adjacent.

[0009] Step 3: Use 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 setup includes: using pressure-based transient solution; defining gravity and its direction; performing numerical simulation calculations using an Euler-Euler two-phase flow model combined with an SST turbulence model; defining the phases, material properties, and interphase interactions in each computational domain; initializing the physical conditions of the computational domain by setting the vacuum port pumping velocity boundary conditions according to the vacuuming process; using the SIMPLE algorithm for pressure-velocity coupling; and setting the convergence residual, computation time step, and total computation time.

[0011] The Euler-Euler two-phase flow model described in step 3 has 2 Euler phases, with the main phase set as air and the secondary phase set as particulate material.

[0012] In step 3, the material properties include particle diameter, particle viscosity parameters, and bulk density; the particle diameters are 10 μm and 100 μm; the particle viscosity adopts the Lun-et-Al model; and the bulk density is 0.55.

[0013] Step 3 describes pressure-velocity coupling using the SIMPLE algorithm, with the convergence residual set to 1e. -5 The calculation time step is 1e. -4 s, Calculate the total duration as 10s;

[0014] Step 3 involves initializing the physical conditions of the computational domain, including setting the boundary conditions to a fixed wall with no slippage, setting the inlet to a velocity inlet, and setting the volume fraction of air to 1.

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

[0016] The beneficial effects of this invention are:

[0017] In existing technologies, experimental methods are typically used to study the gas-solid two-phase flow in a particle bed during vacuuming processes in induction melting furnaces. This approach is prone to damaging the vacuum system and incurs high economic costs. This invention employs an Euler-Euler two-phase flow model for interface tracking in the particle bed gas-solid two-phase flow, which offers advantages in computational cost and efficiency compared to the Euler-Lagrange method. Combined with the SST turbulence model, more accurate transient gas flow characteristics can be obtained, including the location and structure of gas vortices and the instantaneous velocity distribution of the gas flow field. Simultaneously, it overcomes the significant computational demands of direct numerical simulation, which requires solving for all turbulence information. The particle bed gas-solid two-phase flow numerical simulation method based on the Euler-Euler model provided by this invention can accurately simulate particle motion during the vacuuming process, elucidating the influence of factors such as vacuuming velocity and particle diameter on particle trajectories. The simulation process closely approximates actual conditions. Compared to experimental in-situ measurements, this method offers advantages such as not damaging the vacuum system, low economic cost, and high research efficiency. Attached Figure Description

[0018] Figure 1 This is the mesh partitioning diagram of the computational domain for the numerical simulation of the gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace in Example 1.

[0019] Figure 2 This is a velocity contour map of the gas flow field when the boundary condition for the vacuum port pumping speed in Example 1 is 0.2 m / s;

[0020] Figure 3 This is a velocity contour map of the gas flow field when the boundary condition for the vacuum port pumping speed is 2 m / s in Example 1.

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

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

[0023] Figure 6 This is a curve showing the maximum flight speed of the particles in Example 1.

[0024] Figure 7 This is a particle flow field velocity cloud map of the vacuum process for particles with a diameter of 100 μm in Example 1;

[0025] Figure 8 This is a velocity cloud map of the particle flow field during the vacuuming process of particles with a diameter of 10 μm in Example 1. Detailed Implementation

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

[0027] Specific Implementation Method 1: This implementation method, based on the Euler-Euler model, uses a numerical simulation method for gas-solid two-phase flow in granular beds, and proceeds 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 region of the particle bed in the induction melting furnace; establish the 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 external dimensions of the induction melting furnace and the size of the vacuum port of the induction melting furnace.

[0029] Step 2: Use mesh generation software to mesh the computational domain of the numerical simulation of gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace, and refine the mesh in the region where the gas and the particles in the particle bed are adjacent.

[0030] Step 3: Use 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 setup includes: using pressure-based transient solution; defining gravity and its direction; performing numerical simulation calculations using an Euler-Euler two-phase flow model combined with an SST turbulence model; defining the phases, material properties, and interphase interactions in each computational domain; initializing the physical conditions of the computational domain by setting the vacuum port pumping velocity boundary conditions according to the vacuuming process; using the SIMPLE algorithm for pressure-velocity coupling; and setting the convergence residual, computation time step, and total computation time.

[0032] The Euler-Euler two-phase flow model described in step 3 has 2 Euler phases, with the main phase set as air and the secondary phase set as particulate material.

[0033] In step 3, the material properties include particle diameter, particle viscosity parameters, and bulk density; the particle diameters are 10 μm and 100 μm; the particle viscosity adopts the Lun-et-Al model; and the bulk density is 0.55.

[0034] Step 3 describes pressure-velocity coupling using the SIMPLE algorithm, with the convergence residual set to 1e. -5 The calculation time step is 1e. -4 s, Calculate the total duration as 10s;

[0035] Step 3 involves initializing the physical conditions of the computational domain, including setting the boundary conditions to a fixed wall with no slippage, setting the inlet to a velocity inlet, and setting the volume fraction of air to 1.

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

[0037] This embodiment has the following beneficial effects:

[0038] In existing technologies, experimental methods are typically used to study the gas-solid two-phase flow in a particle bed during the vacuuming process in an induction melting furnace. This approach is prone to damaging the vacuum system and is costly. This embodiment employs an Euler-Euler two-phase flow model for interface tracking in the particle bed gas-solid two-phase flow, which offers advantages in computational cost and efficiency compared to the Euler-Lagrange method. Combined with the SST turbulence model, more accurate transient gas flow characteristics can be obtained, including the location and structure of gas vortices and the instantaneous velocity distribution of the gas flow field. This also overcomes the significant computational demands of direct numerical simulation, which requires solving for all turbulence information. The particle bed gas-solid two-phase flow numerical simulation method based on the Euler-Euler model provided in this embodiment can accurately simulate particle motion during the vacuuming process, elucidating the influence of factors such as vacuuming velocity and particle diameter on particle trajectory. The simulation process closely approximates actual conditions. Compared to experimental in-situ measurements, this method offers advantages such as not damaging the vacuum system, low cost, and high research efficiency.

[0039] Specific Implementation Method Two: This implementation method differs from Specific Implementation Method One in that the components inside the induction melting furnace described in step 1 include the furnace body, vacuum port, crucible, and granular raw materials.

[0040] Specific Implementation Method 3: This implementation method differs from Specific Implementation Method 1 or 2 in that the computer-aided design software mentioned in step 1 is Ansys SpaceClaim.

[0041] Specific Implementation Method Four: This implementation method differs from Specific Implementation Methods One to Three in that: the external dimensions of the induction melting furnace described in step 1 are 2m and 1.8m, and the diameter of the vacuum port of the induction melting furnace is 0.8m.

[0042] Specific Implementation Method 5: This implementation method differs from one of the specific implementation methods one to four in that the mesh generation software mentioned in step 2 is Ansys meshing.

[0043] Specific Implementation Method Six: This implementation method differs from one of the specific implementation methods one to five in that the computational fluid dynamics solver mentioned in step 3 is Ansys Fluent.

[0044] Specific Implementation Method Seven: This implementation method differs from Specific Implementation Methods One to Six in that the computational fluid dynamics solver described in step 3 is based on the Euler-Euler model of the particle bed gas-solid two-phase flow mathematical model established by the turbulence equation and the double Euler control equation.

[0045] Specific Implementation Method Eight: This implementation method differs from Specific Implementation Methods One through Seven in that the turbulence equation is:

[0046]

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

[0048] Specific Implementation Method Nine: This implementation method differs from Specific Implementation Methods One through Eight in that the double Euler governing equations include a continuity equation and a momentum equation.

[0049] The continuity equation is:

[0050]

[0051] In the formula: α q ρ is the phase volume fraction; q The q-phase density vector; The q-phase velocity vector; The mass transfer mass from phase p to phase q is the mass transfer mass between phase p and phase q. t represents the interphase mass transfer mass from phase q to phase p; t represents time.

[0052] The momentum equation is:

[0053]

[0054] In the formula: α q ρ is the phase volume fraction; q The q-phase density vector; Pq is the velocity vector of phase q; Pq is the pressure term. For phase q pressure; Let q be the corresponding force strain tensor; For mass force; It is an interphase force; Let p be the interphase mass transfer rate from phase p to phase q; Let be the interphase mass transfer rate from phase q to phase p; This is the additional force for phase q; The mass transfer mass from phase p to phase q is the mass transfer mass between phase p and phase q. Let be the mass transfer mass between the q phase and the p phase.

[0055] Specific Implementation Method 10: This implementation method differs from Specific Implementation Methods 1 to 9 in that the computational fluid dynamics post-processing software is CFD-POST.

[0056] Example 1:

[0057] This embodiment uses the Euler-Euler model to perform numerical simulation of gas-solid two-phase flow in granular beds, following these steps:

[0058] Step 1: Use computer-aided design software to establish a two-dimensional model of the gas-solid two-phase flow region of the particle bed in the induction melting furnace; establish the 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 external dimensions of the induction melting furnace and the size of the vacuum port of the induction melting furnace.

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

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

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

[0062] Step 2: Use mesh generation software to mesh the computational domain of the numerical simulation of gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace, and refine the mesh in the region where the gas and the particles in the particle bed are adjacent.

[0063] The mesh generation software is Ansys Meshing;

[0064] Step 3: Use 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 setup includes: using pressure-based transient solution; defining gravity and its direction; performing numerical simulation calculations using an Euler-Euler two-phase flow model combined with an SST turbulence model; defining the phases, material properties, and interphase interactions in each computational domain; initializing the physical conditions of the computational domain by setting the vacuum port pumping velocity boundary conditions according to the vacuuming process; using the SIMPLE algorithm for pressure-velocity coupling; and setting the convergence residual, computation time step, and total computation time.

[0066] The computational fluid dynamics solver mentioned 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 main phase is set as air, and the secondary phase is set as particulate material.

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

[0069] Step 3 describes pressure-velocity coupling using the SIMPLE algorithm, with the convergence residual set to 1e. -5 The calculation time step is 1e. -4 s, Calculate the total duration as 10s;

[0070] Step 3 involves initializing the physical conditions of the computational domain, including setting the boundary conditions to a fixed, non-slip wall, setting the inlet to a velocity inlet, and setting the air volume fraction to 1. The initialization conditions are set according to the actual working conditions.

[0071] The computational fluid dynamics solver described in step 3 is based on the Euler-Euler model of the granular bed gas-solid two-phase flow mathematical model established by the turbulence equation and the double Euler control equation.

[0072] The gas flow velocity is relatively high and the flow domain structure is complex during the vacuuming process. Therefore, the turbulence equation is used to solve the gas flow behavior. The turbulence equation is as follows:

[0073]

[0074]

[0075] In the formula: α g For gas phase fraction, ρ g U is the gas phase density; k is the turbulent kinetic energy; ω is the specific dissipation rate; U g Γ represents the gas phase velocity. k Γ represents the effective diffusion rate of turbulent kinetic energy. ω G represents the effective diffusion rate relative to the dissipation rate. kW is the term representing the generation of turbulent kinetic energy. k S is the dissipation term of turbulent kinetic energy. k G is the source term for turbulent kinetic energy. ω W is the generating term for the specific dissipation rate. ω S is the dissipation term for the specific dissipation rate; ω The source term for the specific dissipation rate;

[0076] The bi-Euler governing equations include a continuity equation and a momentum equation:

[0077] The continuity equation is:

[0078]

[0079] In the formula: α q ρ is the phase volume fraction; q The q-phase density vector; The q-phase velocity vector; The mass transfer mass from phase p to phase q is the mass transfer mass between phase p and phase q. t represents the interphase mass transfer mass from phase q to phase p; t represents time.

[0080] The momentum equation is:

[0081]

[0082] In the formula: α q ρ is the phase volume fraction; q The q-phase density vector; Pq is the velocity vector of phase q; Pq is the pressure term. For phase q pressure; Let q be the corresponding force strain tensor; For mass force; It is an interphase force; Let p be the interphase mass transfer rate from phase p to phase q; Let be the interphase mass transfer rate from phase q to phase p; This is the additional force for phase q; The mass transfer mass from phase p to phase q is the mass transfer mass between phase p and phase q. The mass transfer mass from phase q to phase p is the mass transfer mass between phases q and p.

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

[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 velocity of the gas flow field inside the furnace increases, and a stable gas vortex is formed between the upper wall of the vacuum port and the outer wall of the crucible. However, this has little impact on the structure and position of the gas flow field vortex. With increasing vacuuming speed, the gas velocity inside the furnace increases, the drag force of the gas on the particles strengthens, and the particle velocity increases. When the vacuuming speed exceeds 2 m / 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 granular bed based on the Euler-Euler model, characterized in that: The numerical simulation method for gas-solid two-phase flow in granular beds based on the Euler-Euler model is carried out according to the following steps: Step 1: Use computer-aided design software to establish a two-dimensional model of the gas-solid two-phase flow region of the particle bed in the induction melting furnace; establish the 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 external dimensions of the induction melting furnace and the size of the vacuum port of the induction melting furnace. Step 2: Use mesh generation software to mesh the computational domain of the numerical simulation of gas-solid two-phase flow in the particle bed of the vacuum process in the induction melting furnace, and refine the mesh in the region where the gas and the particles in the particle bed are adjacent. Step 3: Use 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 setup includes: using pressure-based transient solution; defining gravity and its direction; performing numerical simulation calculations using an Euler-Euler two-phase flow model combined with an SST turbulence model; defining the phases, material properties, and interphase interactions in each computational domain; initializing the physical conditions of the computational domain by setting the vacuum port pumping velocity boundary conditions according to the vacuuming process; using the SIMPLE algorithm for pressure-velocity coupling; and setting the convergence residual, computation time step, and total computation time. The Euler-Euler two-phase flow model described in step 3 has 2 Euler phases, with the main phase set as air and the secondary phase as particulate material. The material properties described in step 3 include particle diameter, particle viscosity parameters, and bulk density; the particle diameter is 10 μm and 100 μm; the particle viscosity adopts the Lun-et-Al model; and the bulk density is 0.

55. Step 3 describes pressure-velocity coupling using the SIMPLE algorithm, with the convergence residual set to 1e. -5 The calculation time step is 1e -4 s, Calculate the total duration as 10s; Step 3 involves initializing the physical conditions of the computational domain, including setting the boundary conditions to a fixed wall with no slippage, setting the inlet to a velocity inlet, and setting the volume fraction of air to 1. Step 4: Output the numerical simulation results of the gas-solid two-phase flow motion characteristics of the particle bed using computational fluid dynamics post-processing software; the numerical simulation results of the gas-solid two-phase flow motion characteristics of the particle bed include the volume flow velocity cloud map, the particle motion velocity cloud map, and the fitting curve of the maximum particle flight speed.

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

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

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

5. The numerical simulation method for gas-solid two-phase flow in granular beds based on the Euler-Euler model according to claim 1, characterized in that: The mesh generation software used in step 2 is Ansys Meshing.

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

7. The numerical simulation method for gas-solid two-phase flow in granular beds based on the Euler-Euler model according to claim 1, characterized in that: The computational fluid dynamics solver described in step 3 is based on the Euler-Euler model of the granular bed gas-solid two-phase flow mathematical model established by the turbulence equation and the double Euler control equation.

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

9. The numerical simulation method for gas-solid two-phase flow in granular beds based on the Euler-Euler model according to claim 1, characterized in that: The bi-Euler governing equations include a continuity equation and a momentum equation: The continuity equation is: In the formula: α q ρ is the phase volume fraction; q The q-phase density vector; The q-phase velocity vector; The mass transfer mass from phase p to phase q is the mass transfer mass between phase p and phase q. t represents the interphase mass transfer mass from phase q to phase p; t represents time. The momentum equation is: In the formula: α q ρ is the phase volume fraction; q The q-phase density vector; Pq is the velocity vector of phase q; Pq is the pressure term. For phase q pressure; Let q be the corresponding force strain tensor; For mass force; It is an interphase force; Let p be the interphase mass transfer rate from phase p to phase q; Let be the interphase mass transfer rate from phase q to phase p; This is the additional force for phase q; The mass transfer mass from phase p to phase q is the mass transfer mass between phase p and phase q. Let be the mass transfer mass between the q phase and the p phase.

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

Images