Discrete element and macro-particle model coupled simulation method for liquid metal particle filter clogging

Abstract

The application discloses a kind of liquid metal particle filter clogging discrete element and macroscopic particle model coupling simulation method, establishes the geometry model of flow calculation domain inside filter and divides grid;Calculate fluid velocity field, adopt discrete element method to track the motion of particle in fluid and the interaction between it and fluid, calculate the collision, adhesion and deposition behavior between particle and particle, particle and wall surface;Adopt macroscopic particle model to characterize the feedback effect of particle volume on local flow field, update velocity field;The method can realize the simultaneous simulation of important variables such as particle behavior and flow field change during filtration process, realize the dynamic characterization of particle migration, capture, clogging and filter cake formation process, and can provide calculation basis for filter structure design and operation optimization.

Description

Technical Field

[0001] This invention belongs to the field of nuclear reactor thermal-hydraulic calculation technology, specifically relating to a coupled simulation method of discrete element and macroscopic particle model for filtration blockage of liquid metal particles. Background Technology

[0002] During the operation of a liquid metal cooling circuit, solid particulate impurities are generated due to oxidation reactions and corrosion of structural materials. These solid particulate impurities cannot be stably dissolved in liquid metal and can only flow with the coolant in particulate form, migrating, colliding, depositing, and accumulating in localized areas such as filters. The continuous accumulation of particles leads to the contraction of filter channels, increased flow resistance, and increased pressure drop. In severe cases, it can cause filter blockage, affecting the safe and stable operation of the liquid metal circuit.

[0003] Particle clogging in filters is not a simple adhesion process, but rather involves complex behaviors such as particle migration, inter-particle collisions, particle-wall collisions, adhesion, bridging, and filter cake formation. In the initial stages of filtration, larger particles are preferentially intercepted by the filter medium, while smaller particles can still penetrate the filter pores. As particles accumulate, the filtration process gradually transitions from deep filtration to surface filtration, with a stable filter cake layer gradually forming on the filter surface. This leads to further contraction of local flow channels, increasing pressure drop, and ultimately filter failure. Numerical simulations of particle filtration and clogging processes in filters must consider not only the interaction between particles and fluid, but also inter-particle collisions, particle-wall contact adhesion, the evolution of particle packing structures, and the feedback effects of clogging on the local flow field and pressure drop distribution. Summary of the Invention

[0004] To overcome the problems existing in the prior art, the present invention aims to provide a coupled simulation method of discrete element method and macroscopic particle model for clogging of liquid metal particle filters. This method can combine the discrete element method with the macroscopic particle model to perform coupled simulation calculations of the migration, collision, stress, deposition and clogging processes of particles in liquid metal within the filter, and realize bidirectional transmission and dynamic updating between particle information and flow field information, thereby achieving synchronous characterization of particle motion behavior and local flow field changes, and providing a calculation basis for the structural design and operation analysis of filters in liquid metal.

[0005] To achieve the above objectives, the present invention adopts the following technical solution: This invention provides a coupled simulation method for liquid metal particle filter clogging using discrete element method (DEM) and macroscopic particle model. The DEM is used to calculate the particle position, velocity, forces acting on the particles, and collision processes between particles. The macroscopic particle model is used to characterize the particle volume occupancy effect and the feedback effect of particles on the local flow field. By establishing a geometric model of the internal flow computational domain of the filter, performing mesh generation, initializing the flow field, and conducting flow field calculations, particle and flow field information are transferred between the DEM and the macroscopic particle model. The fluid velocity at the particle contact mesh is calculated, and a momentum source term is added at the contact mesh to update the flow field, thereby achieving a coupled simulation of particle motion behavior and local flow field changes.

[0006] A coupled discrete element method and macroscopic particle model simulation method for clogging of liquid metal particle filters includes the following steps: Step 1: Establish a geometric model of the computational domain for the liquid metal fluid flow inside the filter, and then perform mesh generation based on the geometric model; Step 2: Set the initial conditions of the flow field, initialize the flow field, and obtain the initial velocity field and temperature field of the liquid metal fluid; Step 3: Based on the liquid metal fluid velocity field obtained in Step 2, calculate the drag force on the particles. If collisions occur between particles or between particles and the wall, calculate the forces acting on the particles during the collision, and calculate the particle velocity in the current time step based on the forces acting on the particles and the particle velocity in the previous time step. By using the discrete element method to calculate the drag force, collision contact force, and damping force on each particle, the migration, collision, stress, deposition, and accumulation behavior of impurity particles during liquid metal filtration can be described simultaneously. This overcomes the problem that traditional continuous medium particle models cannot characterize discrete collisions and local accumulation evolution between particles, and is especially suitable for filtration clogging processes with high particle volume fraction and significant influence of particle size on the local flow field. The specific steps are as follows: Step 3-1: Based on the liquid metal fluid velocity field obtained in Step 2, the drag force on the particles from the fluid is calculated as follows: (1) in: —The particles are subjected to the drag force of the fluid, N; — Fluid viscosity, Pa·s; —Fluid density, kg / m³ 3 ; —Particle size, in meters; —Drag coefficient; —Fluid Reynolds number; —Fluid velocity, m / s; —Particle velocity, m / s; Step 3-2: When a particle collides with another particle or the particle wall, calculate the normal contact force during the collision process. Normal phase damping force Tangential contact force , and tangential damping force The calculation of the normal contact force is as follows: (2) (3) (4) in: — Normal contact force of the particle, N; —Equivalent Young's modulus, GPa; —Normal overlap, m; , —Young's modulus of each component, GPa; , —Poisson's ratio of each component, GPa; —Equivalent radius, m; , —Radii of each component, in meters; The normal phase damping force is calculated as follows: (5) (6) in: —Normal damping force of the particle; —Coefficient of normality; —The normal component of the relative velocity, in m / s; —Normal stiffness, N / m; —Equivalent mass, kg; The tangential contact force is calculated as follows: (7) (8) in: — Tangential contact force of the particles, N; —Normal stiffness, N / m; —Equivalent shear modulus, GPa; — Tangential overlap, m; The tangential damping force is calculated as follows: (9) in: — Tangential damping force of the particle, N; — Tangential coefficient of restitution; — Tangential component of relative velocity, m / s; Based on the forces acting on the particle and the particle velocity in the previous time step, the formula for calculating the particle velocity in the current time step is as follows: (10) in: —The particle velocity at time step n+1, in m / s; —The particle velocity at time step n, in m / s; —Time step, s; —Particle mass, kg; —Acceleration due to gravity, m / s² 2 ; Step 4: Based on the particle position and velocity obtained in Step 3, the velocity of the fluid elements in contact with the particles is corrected using a macroscopic particle model. For each particle and its contact fluid element, the mixed velocity of the particle's current velocity and the fluid element's velocity in the previous time step is obtained by weighting according to the volume fraction. Based on the obtained mixed velocity, the velocity of the fluid elements in contact with the particle is extrapolated to the center of the fluid element, and the velocity of these fluid elements is matched with the particle velocity using a momentum source term. By introducing a macroscopic particle model and using a volume fraction weighting method to handle the momentum exchange between particles and fluid elements, the effects of particle accumulation, pore occupancy, and local blockage on the flow field can be equivalently fed back to the fluid computational unit without the need for complex body-fitted mesh reconstruction of each particle surface. This method retains the detailed description capability of the discrete element method for particle motion and collision behavior, while improving the feasibility and efficiency of fluid-structure interaction calculations under high particle volume fraction conditions. It overcomes the technical obstacles of inconsistent particle scale and fluid mesh scale and excessive computational cost when the number of particles is large.

[0007] Step 5: Introduce the momentum source term calculated in Step 4 into the fluid control equation, re-solve the fluid flow process in the filtration area, and obtain the new liquid metal fluid velocity field distribution results, thereby realizing the dynamic characterization of the influence of particle accumulation and blockage behavior on flow field changes.

[0008] Step 6: In the next time step, the liquid metal fluid velocity field calculated in Step 5 is transferred back to the discrete element method for subsequent calculations of particle forces, collisions, and motion states. After each round of particle-fluid phase coupling calculation is completed, the current simulation time is determined. If the preset termination time has not been reached, the particle motion, mixing velocity calculation, and flow field update process continue. If the termination time is reached or exceeded, the calculation ends and the particle distribution, flow field distribution, and related evolution results during the particle filtration clogging process are output. By cyclically executing particle motion calculation, fluid element velocity correction, and flow field re-solution within each time step, the dynamic simulation of the entire filtration clogging process can be achieved, continuously characterizing the coupling evolution relationship between particle migration, collision, deposition, accumulation, clogging, and flow field redistribution.

[0009] Preferably, step 1 is as follows: Step 1-1: Use geometric modeling software to create a geometric model of the computational domain for the liquid metal fluid flow inside the filter; Step 1-2: Based on the geometric model of the liquid metal fluid flow calculation domain inside the filter obtained in Step 1-1, the filter is structured and meshed using mesh generation software, and the wall mesh is locally refined.

[0010] Preferably, step 2 is as follows: Step 2-1: Based on the calculated operating conditions, set the initial temperature and flow velocity at the inlet of the liquid metal fluid domain inside the filter, and set the pressure at the outlet of the fluid domain and the fluid domain wall as a symmetrical boundary; Step 2-2: Initialize the flow field based on the mass, momentum, and energy equations of the liquid metal fluid to obtain the initial velocity and temperature fields of the liquid metal fluid. The mass, momentum, and energy equations are shown below: (11) (12) (13) in: —Fluid density, kg / s; —Fluid velocity, m / s; —Fluid pressure, m / s; —Viscous stress tensor, Pa; —All volumetric forces and drag sources received by the fluid, N / m 3 ; — Enthalpy of fluid, J / kg; Thermal conductivity, W·m -1 ·K -1 ; —Fluid temperature, K; Liquid metal velocity field results are derived from the velocity of the liquid metal. u Composition, the temperature field results of liquid metal are determined by the temperature of liquid metal. T composition.

[0011] Preferably, step 4 is as follows: The mixing velocity of the particle's current velocity with the velocity of the fluid element in the previous time step is calculated as follows: (14) in: —The mixing velocity at time step n+1, in m / s; —The velocity of the fluid element at time step n, in m / s; —The particle velocity at time step n, in m / s; —Volume fraction, dimensionless; Volume fraction The calculation method is as follows: (15) in: —The overlapping volume of the particle and fluid unit, m 3 ; —Fluid unit volume, m 3 ; The formula for calculating the momentum source term is as follows: (16) in: —The momentum source term corresponding to the (n+1)th time step, N / m 3 ; —Fluid density, kg / s; —The mixing velocity at time step n+1, in m / s; —The velocity of the fluid element at time step n, in m / s; —The velocity of the fluid element at time step n, in m / s; —Time step, s; Compared with the prior art, the present invention has the following beneficial effects: 1) This invention achieves bidirectional transmission and real-time updating of particle information and flow field information through the cyclic coupling calculation of discrete element method and macroscopic particle model, so that particle movement, deposition and accumulation and the resulting local flow field changes can be solved synchronously under the same calculation framework, thereby enabling a more realistic simulation of the interaction between particle filtration clogging process and fluid flow. 2) This invention can simultaneously characterize the migration, collision, stress, deposition, accumulation and clogging behavior of particles during liquid metal filtration. It is especially suitable for complex particle filtration problems with high particle volume fraction, significant influence of particle size on flow field and continuous change of local pore structure during filtration, thus improving the completeness and applicability of numerical simulation of filtration clogging process. 3) This invention can be used for numerical analysis of the filtration behavior, clogging mechanism and filter structure design of impurity particles in liquid metal reactors, providing more accurate calculation support for the optimization of filter pore size, structure, operating flow rate and maintenance cycle, thereby improving the reliability of filter design and operation optimization. Attached Figure Description

[0012] Figure 1 This is a schematic diagram of the forces acting on a particle when it collides with a wall.

[0013] Figure 2 This is a schematic diagram of the particle-contact mesh coupling relationship in a macroscopic particle model.

[0014] Figure 3 This is a schematic diagram of coupled calculation.

[0015] Figure 4 This is a flowchart of the present invention. Detailed Implementation

[0016] The following combination Figure 4 The flowchart shown illustrates the invention in further detail, using numerical simulation of the filtration and clogging process of particles in a filter screen as an example.

[0017] This invention provides a coupled simulation method for discrete element method and macroscopic particle model of liquid metal particle filter clogging, comprising the following steps: Step 1: To perform numerical simulation of the filtration and clogging process of particles in the filter screen, it is necessary to establish a geometric model of the computational domain for the liquid metal fluid flow inside the filter screen, and then perform mesh generation based on the geometric model. The specific steps are as follows: Step 1-1: Using geometric modeling software, a geometric model of the computational domain for the liquid metal fluid flow inside the filter was established. A square flow channel model with four filter holes was created, with a side length of 0.4 mm and a length of 2 mm. The filter hole thickness is 100 μm and the pore diameter is 100 μm.

[0018] Step 1-2: Based on the geometric model of the liquid metal fluid flow calculation domain inside the filter obtained in Step 1-1, the mesh generation software is used to perform structured mesh generation, and the wall mesh is refined. The number of boundary layers is set to 3.

[0019] Step 2: Set the initial conditions for the computational domain of the liquid metal fluid flow inside the filter, specifically including parameters such as inlet velocity, inlet temperature, outlet pressure, and wall boundary conditions. Based on the physical properties of the liquid metal, initialize the fluid state within the computational domain, thereby calculating the velocity and temperature field distributions of the flow computational domain inside the filter at the initial moment. This provides the initial conditions for subsequent iterative flow field calculations and the coupled simulation of particle motion, stress, and clogging processes. The specific steps are as follows: Step 2-1: Based on the calculated operating conditions, set the inlet of the liquid metal fluid domain inside the filter as a velocity inlet and set the initial inlet temperature and inlet velocity; set the fluid domain outlet as a pressure boundary; and set the fluid domain wall as a symmetrical boundary. Step 2-2: Initialize the flow field based on the mass, momentum, and energy equations of the liquid metal fluid to obtain the initial velocity and temperature fields of the liquid metal fluid. The mass, momentum, and energy equations are shown below: (1) (2) (3) in: —Fluid density, kg / s; —Fluid velocity, m / s; —Fluid pressure, m / s; —Viscous stress tensor, Pa; —All volumetric forces and drag sources received by the fluid, N / m 3 ; — Enthalpy of fluid, J / kg; Thermal conductivity, W·m -1 ·K -1 ; — Fluid temperature, K; Step 3: Based on the liquid metal fluid velocity field obtained in Step 2, calculate the drag force exerted on the particles by the fluid. If collisions occur between particles or between particles and the wall, calculate the forces acting on the particles during the collision process, such as... Figure 1 As shown, Figure 1 This is a schematic diagram of the collision process between particles and the wall, where and These are the velocities of the fluid and the particles, respectively. The particles are subjected to the drag force of the fluid; , , and This involves considering the normal contact force, normal damping force, tangential contact force, and tangential damping force experienced by the particle during its collision with the wall. The particle velocity at the current time step is then calculated based on the forces acting on the particle and its velocity from the previous time step. The specific steps are as follows: Step 3-1: Based on the liquid metal fluid velocity field obtained in Step 2, calculate the drag force on the particles caused by the fluid using the non-spherical drag force model proposed by Haider and Levenspiel. : (4) in: —The particles are subjected to the drag force of the fluid, N; — Fluid viscosity, Pa·s; —Fluid density, kg / m³ 3 ; —Particle size, in meters; —Drag coefficient; —Fluid Reynolds number; —Fluid velocity, m / s; —Particle velocity, m / s; Step 3-2: When particles collide with each other or with the particle wall, the Hertz-Mindlin contact model is used to calculate the normal contact force during the collision process. Normal phase damping force Tangential contact force , and tangential damping force The calculation of the normal contact force is as follows: (5) (6) (7) in: — Normal contact force of the particle, N; —Equivalent Young's modulus, GPa; —Normal overlap, m; , —Young's modulus of each component, GPa; , —Poisson's ratio of each component, GPa; —Equivalent radius, m; , —Radii of each component, in meters; The normal phase damping force is calculated as follows: (8) (9) in: —Normal damping force of the particle; —Coefficient of normality; —The normal component of the relative velocity, in m / s; —Normal stiffness, N / m; —Equivalent mass, kg; The tangential contact force is calculated as follows: (10) (11) in: — Tangential contact force of the particles, N; —Normal stiffness, N / m; —Equivalent shear modulus, GPa; — Tangential overlap, m; The tangential damping force is calculated as follows: (12) in: — Tangential damping force of the particle, N; — Tangential coefficient of restitution; — Tangential component of relative velocity, m / s; Based on the forces acting on the particle and the particle velocity in the previous time step, the formula for calculating the particle velocity in the current time step is as follows: (10) in: —The particle velocity at time step n+1, in m / s; —The particle velocity at time step n, in m / s; —Time step, s; —Particle mass, kg; —Acceleration due to gravity, m / s² 2 ; Step 4: Based on the particle position and velocity obtained in Step 3, the velocity of the fluid elements in contact with the particles is corrected using a macroscopic particle model. For each particle and the fluid element in contact with it, the mixed velocity of the particle's current velocity and the velocity of the fluid element in the previous time step is obtained by weighting according to the volume fraction, such as... Figure 2 As shown. Figure 2In the diagram, white grids represent fluid grids, light gray areas represent particles, and dark gray grids represent areas where particles contact or overlap with the fluid grids. Within the contact grids, the mixing velocity is calculated as follows: (14) in: —The mixing velocity at time step n+1, in m / s; —The velocity of the fluid element at time step n, in m / s; —The particle velocity at time step n, in m / s; —Volume fraction, dimensionless; Volume fraction The calculation method is as follows: (15) in: —The overlapping volume of the particle and fluid unit, m 3 ; —Fluid unit volume, m 3 ; Based on the obtained mixing velocity, the velocities of the fluid units in contact with the particle are extrapolated to the center of the fluid unit, and the momentum source term is used to uniformly match the velocities of these fluid units with the particle velocity. The formula for calculating the momentum source term is as follows: (16) in: —The momentum source term corresponding to the (n+1)th time step, N / m 3 ; —Fluid density, kg / s; —The mixing velocity at time step n+1, in m / s; —The velocity of the fluid element at time step n, in m / s; —The velocity of the fluid element at time step n, in m / s; —Time step, s; Step 5: Introduce the momentum source term calculated in Step 4 into the fluid control equation, resolve the fluid flow process in the filtration area, and obtain the updated velocity field distribution of the liquid metal fluid, thereby realizing the dynamic characterization of the influence of particle accumulation and blockage behavior on the local flow field.

[0020] Step 6: In the calculation of the next time step, the flow field information obtained in Step 5 is transferred back to the discrete element method for subsequent calculations of particle forces, collisions, and motion states. For example... Figure 3 As shown, within each macroscopic time step, the current fluid velocity is first calculated. The fluid velocity is then passed to the discrete element method, and the fluid velocity is subsequently used to calculate the fluid force acting on the particles. And by combining the contact forces between particles and between particles and the wall, the particle velocity is updated. And particle position. After completing the discrete element method calculation, the updated particle position, velocity, and particle feedback effect on the fluid are transferred to the macroscopic particle model for calculating and updating the flow field for the next time step. Figure 3 It is evident that the discrete element method and the macroscopic particle model have advantages in... t n 、t n+1 、t n+2 The coupling is performed alternately: the Discrete Element Method (DEM) is mainly responsible for updating particle forces, collisions, and motion, while the macroscopic particle model is mainly responsible for updating the flow field and calculating the particle feedback effect on the fluid. The DEM and the macroscopic particle model exchange particle and fluid information through a data transfer script, thereby achieving bidirectional coupling between the particle and fluid phases.

[0021] The coupled calculation process in steps 3, 4, and 5 will continue iteratively. In each iteration, the particle system and the fluid domain will be updated alternately, and the interaction between particles and fluid will be continuously calculated and transmitted. After each round of coupled calculation of the particle phase and fluid phase is completed, the current simulation time is judged. If the preset termination time has not been reached, the particle motion calculation, mixing velocity calculation, and flow field update process will continue. If the preset termination time is reached or exceeded, the calculation will end, and the particle distribution, flow field distribution, and evolution results during the particle filtration clogging process will be output.

Claims

1. A coupled simulation method of discrete element method and macroscopic particle model for filtration clogging of liquid metal particles, characterized in that: Includes the following steps: Step 1: Establish a geometric model of the computational domain for the liquid metal fluid flow inside the filter, and then perform mesh generation based on the geometric model; Step 2: Set the initial conditions of the flow field, initialize the flow field, and obtain the initial velocity field and temperature field of the liquid metal fluid; Step 3: Based on the liquid metal fluid velocity field obtained in Step 2, calculate the drag force on the particles. If collisions occur between particles or between particles and the wall, calculate the force on the particles during the collision. Then, calculate the particle velocity in the current time step based on the force on the particles and the particle velocity in the previous time step. The specific steps are as follows: Step 3-1: Based on the liquid metal fluid velocity field obtained in Step 2, the drag force on the particles from the fluid is calculated as follows: (1) in: —The particles are subjected to the drag force of the fluid, N; — Fluid viscosity, Pa·s; —Fluid density, kg / m³ 3 ; —Particle size, in meters; —Drag coefficient; —Fluid Reynolds number; —Fluid velocity, m / s; —Particle velocity, m / s; Step 3-2: When particles collide with each other or with the wall, calculate the normal contact force during the collision process. Normal phase damping force Tangential contact force , and tangential damping force The calculation of the normal contact force is as follows: (2) (3) (4) in: — Normal contact force of the particle, N; —Equivalent Young's modulus, GPa; —Normal overlap, m; , —Young's modulus of each component, GPa; , —Poisson's ratio of each component, GPa; —Equivalent radius, m; , —Radii of each component, in meters; The normal phase damping force is calculated as follows: (5) (6) in: —Normal damping force of the particle; —Coefficient of normality; —The normal component of the relative velocity, in m / s; —Normal stiffness, N / m; —Equivalent mass, kg; The tangential contact force is calculated as follows: (7) (8) in: — Tangential contact force of the particles, N; —Normal stiffness, N / m; —Equivalent shear modulus, GPa; — Tangential overlap, m; The tangential damping force is calculated as follows: (9) in: — Tangential damping force of the particle, N; — Tangential coefficient of restitution; — Tangential component of relative velocity, m / s; Based on the forces acting on the particle and the particle velocity in the previous time step, the formula for calculating the particle velocity in the current time step is as follows: (10) in: —The particle velocity at time step n+1, in m / s; —The particle velocity at time step n, in m / s; —Time step, s; —Particle mass, kg; —Acceleration due to gravity, m / s² 2 ; Step 4: Based on the particle position and velocity obtained in Step 3, the velocities of the fluid units in contact with the particles are corrected using a macroscopic particle model. For each particle and the fluid unit in contact with it, the mixed velocity of the particle's current velocity and the velocity of the fluid unit in the previous time step is obtained by weighting according to the volume fraction. Based on the obtained mixed velocity, the velocities of the fluid units in contact with the particle are extrapolated to the center of the fluid unit, and the momentum source term is used to match the velocities of these fluid units with the particle velocity. Step 5: Introduce the momentum source term calculated in Step 4 into the fluid control equation, re-solve the fluid flow process in the filtration area, and obtain the new liquid metal fluid velocity field distribution results, thereby realizing the dynamic characterization of the influence of particle accumulation and blockage behavior on flow field changes. Step 6: In the calculation of the next time step, the velocity field of the liquid metal fluid obtained in Step 5 is transferred back to the discrete element method for subsequent calculation of particle force, collision and motion state; after each round of particle phase and fluid phase coupling calculation is completed, the current simulation time is determined. If the preset termination time has not been reached, the particle motion, mixing velocity calculation and flow field update process continue to be executed; when the termination time is reached or exceeded, the calculation ends and the particle distribution, flow field distribution and related evolution results of the particle filtration clogging process are output.

2. The method for coupled simulation of discrete element and macroscopic particle model of liquid metal particle filter clogging according to claim 1, characterized in that: Step 1 is as follows: Step 1-1: Use geometric modeling software to create a geometric model of the computational domain for the liquid metal fluid flow inside the filter; Step 1-2: Based on the geometric model of the liquid metal fluid flow calculation domain inside the filter obtained in Step 1-1, the filter is structured and meshed using mesh generation software, and the wall mesh is locally refined.

3. The method for coupled simulation of discrete element and macroscopic particle model of liquid metal particle filter clogging according to claim 1, characterized in that: Step 2 is as follows: Step 2-1: Based on the calculated operating conditions, set the initial temperature and flow velocity at the inlet of the liquid metal fluid domain inside the filter, and set the pressure at the outlet of the fluid domain and the fluid domain wall as a symmetrical boundary; Step 2-2: Initialize the flow field according to the mass, momentum, and energy equations of the liquid metal fluid to obtain the initial velocity and temperature fields of the liquid metal fluid; the mass, momentum, and energy equations are shown below: (11) (12) (13) in: —Fluid density, kg / s; —Fluid velocity, m / s; —Fluid pressure, m / s; —Viscous stress tensor, Pa; —All volumetric forces and drag sources received by the fluid, N / m 3 ; — Enthalpy of fluid, J / kg; Thermal conductivity, W·m -1 ·K -1 ; —Fluid temperature, K; Liquid metal velocity field results are derived from the velocity of the liquid metal. u Composition, the temperature field results of liquid metal are determined by the temperature of liquid metal. T composition.

4. The method for coupled simulation of discrete element and macroscopic particle model of liquid metal particle filter clogging according to claim 1, characterized in that: Step 4 is as follows: The mixing velocity of the particle's current velocity with the velocity of the fluid element in the previous time step is calculated as follows: (14) in: —The mixing velocity at time step n+1, in m / s; —The velocity of the fluid element at time step n, in m / s; —The particle velocity at time step n, in m / s; —Volume fraction, dimensionless; Volume fraction The calculation method is as follows: (15) in: —The overlapping volume of the particle and fluid unit, m 3 ; —Fluid unit volume, m 3 ; The formula for calculating the momentum source term is as follows: (16) in: —The momentum source term corresponding to the (n+1)th time step, N / m 3 ; —Fluid density, kg / s; —The mixing velocity at time step n+1, in m / s; —The velocity of the fluid element at time step n, in m / s; —Time step, s.

Images