A simulation method and system for dynamic distribution of inclusions in continuous casting billets

Abstract

The present invention relates to the technical field of continuous casting production of steel, and in particular to a simulation method, system, and computer program product for the dynamic distribution of inclusions in a continuous casting billet, comprising: constructing a three-dimensional model of the continuous casting billet, importing FLUENT software, and simulating the movement of fluid in the continuous casting billet; simulating the growth process of the solidified billet shell until complete solidification; injecting and tracking inclusions in a discrete phase model; obtaining the position of the inclusions in the continuous casting billet at each time step, and judging whether each inclusion meets the criteria of being captured by the solidified billet shell; if so, clearing the speed of the inclusion and assigning a new speed to the inclusion based on its position in the continuous casting billet; and deriving inclusion information at the outlet of the continuous casting billet. The simulation results of the dynamic distribution of inclusions obtained by the present invention are more in line with actual production, can accurately reflect the distribution of inclusions in the continuous casting billet, and effectively guide the production of the casting billet.

Description

Technical Field

[0001] The present invention relates to the technical field of steel continuous casting production, and in particular to a simulation method, system and computer program product for dynamic distribution of inclusions in a continuous casting billet. Background Art

[0002] With the rapid development of society and economy, the demand for high-quality steel is increasing. Currently, most steel products are produced by continuous casting, and the level of continuous casting technology is directly related to the quality and output of steel products.

[0003] Inclusions in continuous casting billets refer to various non-metallic impurity particles present in the molten steel during the continuous casting process. They are a key factor affecting billet quality. Inclusions have become a primary concern for continuous casting billet quality and are currently a focus of attention within the metallurgical community.

[0004] In actual production, due to the complexity of the steelmaking process, it is difficult to maintain constant continuous casting parameters. Fluctuations in the continuous casting parameters within a certain range can significantly affect the subsequent distribution of inclusions. Especially during the unsteady-state casting process, inclusions in the billet undergo a repetitive cycle of decreasing, increasing, and then decreasing again. Therefore, numerically simulating the distribution of inclusions within the continuous casting billet and accurately identifying and shortening steel plates with high inclusion concentrations is crucial for improving metal yield and the quality of cast steel products, thereby promoting the manufacture of high-end steel products.

[0005] In the existing FLUENT-based simulation method for inclusion distribution in continuous casting billets, inclusions remain stationary after being captured. This does not conform to the actual production situation in which inclusions continue to move with the billet after leaving the crystallizer. It cannot reflect the dynamic distribution law of inclusions and cannot accurately reflect the actual situation of inclusion distribution in the continuous casting billet at different casting stages. Summary of the Invention

[0006] Therefore, the technical problem to be solved by the present invention is to overcome the problem in the prior art that the actual situation of inclusion distribution cannot be accurately reflected, which is not conducive to improving the quality of cast steel products.

[0007] To solve the above technical problems, the present invention provides a simulation method for the dynamic distribution of inclusions in a continuous casting billet, comprising:

[0008] Construct a three-dimensional model of the continuous casting billet and import it into FLUENT software to simulate the fluid movement inside the continuous casting billet;

[0009] The energy equation and solidification-melting model are enabled. Under the condition that the flow velocity of the entire flow basin is set as the pulling velocity, the growth process of the solidified shell is simulated until it is completely solidified.

[0010] After the shell is completely solidified, the energy equation and solidification-melting model are turned off, the flow equation and VOF model are turned on, and inclusions are injected and tracked in the discrete phase model.

[0011] At each time step, the position of the inclusions in the continuous casting billet is obtained to determine whether each inclusion meets the criteria for being captured by the solidified billet shell. If so, the inclusion's velocity is removed and a new velocity is assigned to the inclusion based on its position in the continuous casting billet, so that the inclusion's velocity in the horizontal X direction is 0 and the total velocity is equal to the casting speed. If not, no action is taken.

[0012] Export inclusion information at the exit of the continuous casting slab calculation domain.

[0013] Preferably, in FLUENT software, a large eddy simulation model is used to simulate the fluid movement in the continuous casting billet.

[0014] Preferably, the standard for complete solidification of the shell is that the solidification front of the shell does not move within a preset time.

[0015] Preferably, the inclusions injected into the discrete phase model are multiple types of inclusions with the same number, and each type of inclusions has different mass and size.

[0016] Preferably, the criteria for inclusions to be captured by the solidified shell are: the inclusions move to the shell solidification interface, and the solid phase fraction of the shell solidification interface is greater than or equal to 0.67.

[0017] Preferably, assigning a new velocity to the inclusion according to its position in the continuous casting strand comprises:

[0018] According to the position of the inclusions in the continuous casting billet, it is judged whether the inclusions are located in the vertical section, curved section or horizontal section of the continuous casting billet;

[0019] If the inclusion is located in the vertical section of the continuous casting billet, the speed of the inclusion in the horizontal X and Y directions is set to 0, and the speed in the vertical Z direction is set to the pulling speed;

[0020] If the inclusion is located in the bending section of the continuous casting billet, the horizontal X-direction speed of the inclusion is set to 0, and the horizontal Y-direction speed is set to , the vertical Z direction speed is set to ,in is the pulling speed, angle , and are the Z and Y coordinates of the inclusion, and are the Z and Y coordinates of the center of the curved segment respectively;

[0021] If the inclusion is located in the horizontal section of the continuous casting billet, the speed of the inclusion in the horizontal X direction and the vertical Z direction is set to 0, and the speed in the horizontal Y direction is set to the pulling speed.

[0022] Preferably, the size of a single grid set for simulating the fluid movement in the continuous casting billet is smaller than the moving distance of the inclusion in one time step.

[0023] Preferably, inclusion information is derived at the outlet of the continuous casting strand, wherein the inclusion information includes the size, mass, velocity and component velocities in various directions of the inclusions.

[0024] The present invention also provides a simulation system for the dynamic distribution of inclusions in a continuous casting billet, comprising:

[0025] Model building module, used to build a three-dimensional model of the continuous casting billet, import it into FLUENT software, and simulate the fluid movement inside the continuous casting billet;

[0026] The shell solidification calculation module is used to start the energy equation and solidification melting model. Under the condition that the flow velocity of the entire flow basin is set as the pulling velocity, the solidification shell growth process is simulated until it is completely solidified.

[0027] The inclusion injection module is used to inject and track inclusions in the discrete phase model after the shell is completely solidified, close the energy equation and solidification-melting model, open the flow equation and VOF model;

[0028] The inclusion velocity setting module is used to obtain the position of inclusions in the continuous casting billet at each time step and determine whether each inclusion meets the criteria for being captured by the solidified billet shell. If so, the inclusion velocity is cleared and a new velocity is assigned to the inclusion based on its position in the continuous casting billet, so that the inclusion velocity in the horizontal X direction is 0 and the total velocity is equal to the casting speed. If not, no processing is performed.

[0029] The inclusion information export module is used to export inclusion information at the exit of the continuous casting slab calculation domain.

[0030] The present invention also provides a computer program product, comprising a computer program, which, when executed by a processor, implements the steps of the simulation method for dynamic distribution of inclusions in a continuous casting billet.

[0031] The above technical solution of the present invention has the following beneficial effects compared with the prior art:

[0032] The present invention provides a simulation method for the dynamic distribution of inclusions within a continuous casting billet. Inclusions are injected after the shell within the continuous casting billet is completely solidified, and their movement distribution is tracked. When the inclusions are captured by the solidified shell, different speeds are assigned to the inclusions according to their actual movement within the continuous casting billet and their different positions within the billet, thereby resolving the problem of inclusions remaining stationary in the bending section after being captured by the shell in existing methods. This method provides a simulation result of the dynamic distribution of inclusions that moves with the continuous casting billet. The simulation result of the dynamic distribution of inclusions obtained by the present invention is more consistent with actual production, can accurately reflect the distribution of inclusions within the continuous casting billet, effectively guide the production of the continuous casting billet, and improve the quality of cast steel products. Furthermore, the present invention consumes a small amount of computing resources to achieve the desired effect of inclusion movement, and has greater scalability. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to make the content of the present invention more clearly understood, the present invention is further described in detail below based on specific embodiments of the present invention in conjunction with the accompanying drawings, wherein:

[0034] Figure 1 This is a flow chart of a simulation method for dynamic distribution of inclusions in a continuous casting billet according to the present invention;

[0035] Figure 2 is a schematic diagram of the complete solidification of the shell, where Figure 2 (a) is a schematic diagram of the solidification of the shell with a calculation time of 600s. Figure 2 (b) is a schematic diagram of the solidification of the shell with a calculation time of 650s;

[0036] Figure 3 It is a schematic diagram for calculating the inclusion velocity based on the position of the inclusion in the continuous casting billet;

[0037] Figure 4 is a schematic diagram of the movement of inclusions in the vertical section of the continuous casting billet, where Figure 4 (a) is a schematic diagram of the initial capture position of inclusions in the vertical section of the continuous casting billet. Figure 4 (b) is a schematic diagram of the calculation results of the position reached by the inclusion in the vertical section of the continuous casting billet after t time;

[0038] Figure 5 is a schematic diagram of the movement of inclusions in the bending section of the continuous casting billet, where Figure 5 (a) is a schematic diagram of the initial capture position of inclusions in the bending section of the continuous casting billet. Figure 5 (b) is a schematic diagram of the calculation results of the position reached by the inclusion in the bending section of the continuous casting billet after t time;

[0039] Figure 6 is a schematic diagram of the movement of inclusions in the horizontal section of the continuous casting billet, where Figure 6(a) is a schematic diagram of the initial capture position of inclusions in the horizontal section of the continuous casting billet. Figure 6 (b) is a schematic diagram of the calculation results of the position reached by the inclusions after t time when they move downward in the horizontal section of the continuous casting billet. DETAILED DESCRIPTION

[0040] The present invention will be further described below with reference to the accompanying drawings and specific embodiments so that those skilled in the art can better understand the present invention and implement it. However, the embodiments are not intended to limit the present invention.

[0041] Reference Figure 1 As shown, the present invention provides a simulation method for the dynamic distribution of inclusions in a continuous casting billet, comprising:

[0042] S1: Construct a three-dimensional model of the continuous casting billet, import it into FLUENT software, and simulate the fluid movement inside the continuous casting billet.

[0043] In this embodiment, Solidworks was used to construct a three-dimensional model of the continuous casting billet. In FLUENT software, the most advanced large eddy simulation model was used to simulate the fluid movement in the continuous casting billet.

[0044] S2: Turn on the energy equation and solidification-melting model, and simulate the solidification shell growth process until it is completely solidified, setting the flow velocity of the entire flow basin as the pulling velocity.

[0045] The growth process of the solidified shell is as follows: after the molten steel is injected, the ingot solidifies from the outermost layer to the inner part, and the solidified part forms the shell. The shell is completely solidified to form the ingot, usually in the horizontal section of the continuous casting ingot.

[0046] In this embodiment, the criterion for determining whether the solidified shell is completely solidified is that the solidification front of the shell does not change with time, that is, the solidification front does not move within a preset time, which is considered to be completely solidified. In this embodiment, the preset time is 50 seconds.

[0047] The solidification front of the shell refers to the front end of the unsolidified molten steel. Figure 2 As shown, Figure 2 (a) is a schematic diagram of the solidification of the shell with a calculation time of 600s. Figure 2 (b) is a schematic diagram of the solidification of the billet shell with a calculation time of 650 s. The solidification of the billet is radial, and the point where the front disappears is the complete solidification section. The solidification front does not move, that is, the billet shell is completely solidified into the billet at this point.

[0048] S3: After the shell is completely solidified, the energy equation and solidification-melting model are turned off, the flow equation and VOF (Volume of Fluid) model are turned on, and inclusions are injected and tracked in the discrete phase model (DPM).

[0049] The Discrete Phase Model (DPM) is an advanced model in Fluent software for simulating the motion of discrete phase particles in fluids. It allows users to analyze the behavior of particles, droplets, or bubbles in fluid flow, detailing the interactions between particles and fluids, as well as particle-particle collisions, mergers, and breakup processes. DPM has a wide range of applications in engineering, including but not limited to industrial combustion, aerosol spraying, particulate material processing, and the petrochemical industry.

[0050] In this embodiment, the inclusions injected into the discrete phase model are multiple types of inclusions of equal quantity, each type of inclusion has different mass and size, so as to simulate the inclusions in the actual production process. The initial position of the inclusion injection is the inlet of the molten steel.

[0051] S4: Obtain the position of inclusions in the continuous casting billet at each time step and determine whether each inclusion meets the criteria for being captured by the solidified billet shell; if so, remove the velocity of the inclusion and assign a new velocity to the inclusion according to its position in the continuous casting billet, so that the velocity of the inclusion in the horizontal X direction is 0 and the total velocity is equal to the casting speed; if not, no processing is performed.

[0052] The horizontal X direction is the width direction of the continuous casting billet, and the horizontal Y direction is the thickness direction of the continuous casting billet.

[0053] Preferably, the criteria for inclusion capture by the solidifying shell are: the inclusion moves to the shell solidification interface, and the solid fraction of the shell solidification interface is greater than or equal to 0.67. Once captured by the shell, the inclusion will no longer move horizontally in the X direction, but will only move horizontally in the Y and vertical Z directions with the cast strand until it leaves the continuous casting strand.

[0054] According to the size of the continuous casting billet, it is divided into three ranges: vertical section, curved section and horizontal section. Taking the size of the continuous casting billet used in this embodiment as an example, Figure 3 As shown, the vertical segment is the portion where the Z coordinate is less than 3.305, the curved segment is the portion where the Z coordinate is greater than 3.305 and the Y coordinate is less than 4.914, and the horizontal segment is the portion where the Y coordinate is greater than 4.914.

[0055] Specifically, a new velocity is assigned to the inclusion according to its position in the continuous casting billet, including:

[0056] According to the position of the inclusions in the continuous casting billet, it is judged whether the inclusions are located in the vertical section, curved section or horizontal section of the continuous casting billet;

[0057] If the inclusion is located in the vertical section of the continuous casting billet, the speed of the inclusion in the horizontal X and Y directions is set to 0, and the speed in the vertical Z direction is set to the pulling speed; Figure 4 is a schematic diagram of the movement of inclusions in the vertical section of the continuous casting billet, where Figure 4(a) is a schematic diagram of the initial capture position of inclusions in the vertical section of the continuous casting billet. Figure 4 (b) is a schematic diagram of the calculation results of the position reached by the inclusion in the vertical section of the continuous casting billet after t time;

[0058] If the inclusion is located in the bending section of the continuous casting billet, the horizontal X-direction speed of the inclusion is set to 0, and the horizontal Y-direction speed is set to , the vertical Z direction speed is set to ,in is the pulling speed, angle , and are the Z and Y coordinates of the inclusion, and are the Z and Y coordinates of the center of the curved segment respectively; Figure 5 is a schematic diagram of the movement of inclusions in the bending section of the continuous casting billet, where Figure 5 (a) is a schematic diagram of the initial capture position of inclusions in the bending section of the continuous casting billet. Figure 5 (b) is a schematic diagram of the calculation results of the position reached by the inclusion in the bending section of the continuous casting billet after t time;

[0059] If the inclusion is located in the horizontal section of the continuous casting billet, the speed of the inclusion in the horizontal X direction and the vertical Z direction is set to 0, and the speed in the horizontal Y direction is set to the pulling speed; Figure 6 is a schematic diagram of the movement of inclusions in the horizontal section of the continuous casting billet, where Figure 6 (a) is a schematic diagram of the initial capture position of inclusions in the horizontal section of the continuous casting billet. Figure 6 (b) is a schematic diagram of the calculation results of the position reached by the inclusions after t time when they move downward in the horizontal section of the continuous casting billet.

[0060] The basic unit used by FLUENT software to simulate fluid motion is the grid. The numerical calculations for each step are the result of individual changes in grid parameters. If the distance an inclusion moves within a time step is smaller than one grid, the calculation of the inclusion's movement distance will be inaccurate, and it may also cause inclusions to move suddenly after accumulating multiple time steps. Therefore, in this embodiment, the single grid size used to simulate fluid motion within the continuous casting billet is smaller than the distance an inclusion moves within a time step to ensure the continuity and regularity of the inclusion's movement.

[0061] Numerical calculations were performed using FLUENT software, with data saved at regular intervals. The distribution of inclusions was observed through cloud plots until inclusions were observed to move normally in vertical, curved, and horizontal sections, with the calculated error between the combined velocity and the casting speed no greater than 20%. The criteria for determining whether inclusion distribution was uniform and regular were that the inclusions moved in a directional manner on the inner or outer arc side, without significant deviation. The wide surface of the continuous casting billet closest to the center of the curve was defined as the inner arc side, while the wide surface farther from the center of the curve was defined as the outer arc side.

[0062] Based on the user-defined functions in the FLUENT software, secondary development is carried out to determine whether inclusions are captured by the billet shell and to assign new speeds to inclusions. By utilizing the loading characteristics of the user-defined functions in the FLUENT software, the function of the inclusion moving with the billet is realized through the combination of several macro functions. On the basis of ensuring the simulation accuracy, the calculation of the present invention has higher stability, faster convergence speed, and is more in line with production practice. In addition, the technical route and concept of the present invention are not limited to the simulation of inclusion movement in continuous casting billets, but can also be extended to simulation application scenarios such as the removal of inclusions inside continuous casting billets.

[0063] S5: Export inclusion information through REPORT at the exit of the continuous casting slab calculation domain, including: inclusion size, mass, velocity and component velocity in each direction.

[0064] In order to verify the effect of the simulation method for dynamic distribution of inclusions in a continuous casting billet according to the present invention, two groups of experiments were conducted in this embodiment, wherein the first group of experiments included step S4 and the second group of experiments did not include step S4.

[0065] In the first set of experiments, the large eddy model was used in the FLUENT software to simulate the fluid flow. The size of the continuous casting billet was 170×1570mm, including the left nozzle, the right nozzle and the bottom nozzle. The pulling speed was set to 1.6m / min, the time step was 0.05s, and there was no blowing condition. The energy equation and the solidification and melting model were turned on to perform solidification calculations. When the calculation reached 600s, it was found that the solidification front hardly fluctuated with time, and it was considered to be completely solidified. Therefore, the energy equation and the solidification and melting model were turned off, the flow equation and the VOF model were turned on, the flow field calculation was performed, and inclusions were injected. In this set of experiments, the sizes of the injected inclusions were 0.0005mm, 0.005mm, and 0.01mm, respectively, and the injection time was 1s. The total mass of each type of inclusion injected was 1×10 -7 After the inclusions were injected, the dynamic distribution of the inclusions was observed after the calculation time reached 1000 s. It was observed that the inclusions were captured and remained stationary, without moving with the shell.

[0066] The second set of experiments used the large eddy model in FLUENT software to simulate fluid flow. The size of the continuous casting billet was 170×1570mm, including the left nozzle, the right nozzle and the bottom nozzle. The pulling speed was set to 1.6m / min, the time step was 0.05s, and there was no blowing condition. The energy equation and the solidification and melting model were turned on to perform solidification calculations. When the calculation reached 600s, it was found that the solidification front hardly fluctuated with time, and it was considered to be completely solidified. Therefore, the energy equation and the solidification and melting model were turned off, the flow equation and the VOF model were turned on, the flow field calculation was performed, and inclusions were injected. In this set of experiments, the sizes of the injected inclusions were 0.0005mm, 0.005mm, and 0.01mm, respectively, and the injection time was 1s. The total mass of each type of inclusion injected was 1×10 -7 kg. The position of inclusions within the continuous casting billet is determined at each time step to determine whether each inclusion meets the adsorption criteria for capture by the billet shell. If so, the velocity of the inclusion is removed and a new velocity is assigned to the inclusion based on its position within the billet. A monitoring surface is established at the bottom outlet to record information about inclusions that leave the billet with the billet and export it via REPORT. When the calculation time reaches 1000s, the dynamic distribution of the inclusions is observed to determine whether it has a stable dynamic distribution. It is found that the inclusions move with the billet shell and have a stable dynamic distribution.

[0067] In summary, the simulation method for the dynamic distribution of inclusions in a continuous casting billet described in the present invention injects inclusions after the shell in the continuous casting billet is completely solidified, and tracks its movement distribution; when the inclusions are captured by the solidified shell, different speeds are assigned to the inclusions according to the actual movement of the continuous casting billet and the different positions of the inclusions in the continuous casting billet, so as to solve the problem of the existing method that the inclusions remain stationary in the bending section after being captured by the shell, and obtain the dynamic distribution simulation results of the inclusions moving with the continuous casting billet. The simulation results of the dynamic distribution of inclusions obtained by the present invention are more in line with production practice, can accurately reflect the distribution of inclusions in the continuous casting billet, effectively guide the production of the continuous casting billet, and improve the quality of cast steel products; and the present invention consumes a small amount of computing resources to achieve the expected effect of inclusion movement, and has stronger generalizability.

[0068] Based on the above-mentioned simulation method for dynamic distribution of inclusions in a continuous casting billet, this embodiment further provides a simulation system for dynamic distribution of inclusions in a continuous casting billet, comprising:

[0069] Model building module, used to build a three-dimensional model of the continuous casting billet, import it into FLUENT software, and simulate the fluid movement inside the continuous casting billet;

[0070] The shell solidification calculation module is used to start the energy equation and solidification melting model. Under the condition that the flow velocity of the entire flow basin is set as the pulling velocity, the solidification shell growth process is simulated until it is completely solidified.

[0071] The inclusion injection module is used to inject and track inclusions in the discrete phase model after the shell is completely solidified, close the energy equation and solidification-melting model, open the flow equation and VOF model;

[0072] The inclusion velocity setting module is used to obtain the position of inclusions in the continuous casting billet at each time step and determine whether each inclusion meets the criteria for being captured by the solidified billet shell. If so, the inclusion velocity is cleared and a new velocity is assigned to the inclusion based on its position in the continuous casting billet, so that the inclusion velocity in the horizontal X direction is 0 and the total velocity is equal to the casting speed. If not, no processing is performed.

[0073] The inclusion information export module is used to export inclusion information at the exit of the continuous casting slab calculation domain.

[0074] This embodiment further provides a computer program product, including a computer program, which, when executed by a processor, implements the steps of the above-mentioned simulation method for dynamic distribution of inclusions in a continuous casting billet.

[0075] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware. Furthermore, the present application may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0076] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0077] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1The function specified in one or more boxes.

[0078] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0079] Obviously, the above embodiments are merely examples for clarity of explanation and are not intended to limit the implementation methods. Those skilled in the art will appreciate that other variations or modifications can be made based on the above description. It is not necessary and impossible to enumerate all implementation methods here. Obvious variations or modifications arising therefrom remain within the scope of protection of the present invention.

Claims

1. A simulation method for dynamic distribution of inclusions in continuous casting billets, characterized in that: include: Construct a three-dimensional model of the continuous casting billet and import it into FLUENT software to simulate the fluid movement inside the continuous casting billet; The energy equation and solidification-melting model are enabled. Under the condition that the flow velocity of the entire flow basin is set as the pulling velocity, the growth process of the solidified shell is simulated until it is completely solidified. After the shell is completely solidified, the energy equation and solidification-melting model are turned off, the flow equation and VOF model are turned on, and inclusions are injected and tracked in the discrete phase model. Obtain the position of inclusions in the continuous casting billet at each time step and determine whether each inclusion meets the criteria of being captured by the solidified billet shell; If so, the speed of the inclusion is removed and a new speed is assigned to the inclusion according to its position in the continuous casting billet, so that the speed of the inclusion in the horizontal X direction is 0 and the total speed is equal to the pulling speed, including: According to the position of the inclusions in the continuous casting billet, it is judged whether the inclusions are located in the vertical section, curved section or horizontal section of the continuous casting billet; If the inclusion is located in the vertical section of the continuous casting billet, the speed of the inclusion in the horizontal X and Y directions is set to 0, and the speed in the vertical Z direction is set to the pulling speed; If the inclusion is located in the bending section of the continuous casting billet, the horizontal X-direction speed of the inclusion is set to 0, and the horizontal Y-direction speed is set to , the vertical Z direction speed is set to ,in is the pulling speed, the angle between the direction of inclusion movement and the horizontal direction , and are the Z and Y coordinates of the inclusion, and are the Z and Y coordinates of the center of the curved segment respectively; If the inclusion is located in the horizontal section of the continuous casting billet, the speed of the inclusion in the horizontal X direction and the vertical Z direction is set to 0, and the speed in the horizontal Y direction is set to the pulling speed; If not, no action will be taken; Export inclusion information at the exit of the continuous casting slab calculation domain.

2. The method for simulating the dynamic distribution of inclusions in a continuous casting billet according to claim 1, characterized in that: In FLUENT software, the large eddy simulation model is used to simulate the fluid movement in the continuous casting billet.

3. The method for simulating the dynamic distribution of inclusions in a continuous casting billet according to claim 1, characterized in that: The standard for complete solidification of the shell is that the solidification front of the shell does not move within a preset time.

4. The method for simulating the dynamic distribution of inclusions in a continuous casting billet according to claim 1, characterized in that: In the discrete phase model, the injected inclusions are multiple types of inclusions with the same number, and each type of inclusion has different mass and size.

5. The method for simulating the dynamic distribution of inclusions in a continuous casting billet according to claim 1, characterized in that: The criteria for inclusions to be captured by the solidifying shell are: the inclusions move to the shell solidification interface, and the solid phase fraction of the shell solidification interface is greater than or equal to 0.

67.

6. The method for simulating the dynamic distribution of inclusions in a continuous casting billet according to claim 1, characterized in that: The single grid size set to simulate the fluid movement in the continuous casting billet is smaller than the moving distance of the inclusion in one time step.

7. The method for simulating the dynamic distribution of inclusions in a continuous casting billet according to claim 1, characterized in that: Inclusion information is derived at the outlet of the continuous casting strand, and the inclusion information includes the size, mass, velocity and component velocities in various directions of the inclusions.

8. A simulation system for the dynamic distribution of inclusions in continuous casting billets, characterized in that: include: Model building module, used to build a three-dimensional model of the continuous casting billet, import it into FLUENT software, and simulate the fluid movement inside the continuous casting billet; The shell solidification calculation module is used to start the energy equation and solidification melting model. Under the condition that the flow velocity of the entire flow basin is set as the pulling velocity, the solidification shell growth process is simulated until it is completely solidified. The inclusion injection module is used to inject and track inclusions in the discrete phase model after the shell is completely solidified, close the energy equation and solidification-melting model, open the flow equation and VOF model; The inclusion velocity setting module is used to obtain the position of inclusions in the continuous casting billet at each time step and determine whether each inclusion meets the criteria of being captured by the solidified billet shell. If so, the inclusion velocity is cleared and a new velocity is assigned to the inclusion based on its position in the continuous casting billet, so that the velocity of the inclusion in the horizontal X direction is 0 and the total velocity is equal to the casting speed. It includes: According to the position of the inclusions in the continuous casting billet, it is judged whether the inclusions are located in the vertical section, curved section or horizontal section of the continuous casting billet; If the inclusion is located in the vertical section of the continuous casting billet, the speed of the inclusion in the horizontal X and Y directions is set to 0, and the speed in the vertical Z direction is set to the pulling speed; If the inclusion is located in the bending section of the continuous casting billet, the horizontal X-direction speed of the inclusion is set to 0, and the horizontal Y-direction speed is set to , the vertical Z direction speed is set to ,in is the pulling speed, the angle between the direction of inclusion movement and the horizontal direction , and are the Z and Y coordinates of the inclusion, and are the Z and Y coordinates of the center of the curved segment respectively; If the inclusion is located in the horizontal section of the continuous casting billet, the speed of the inclusion in the horizontal X direction and the vertical Z direction is set to 0, and the speed in the horizontal Y direction is set to the pulling speed; If not, no action will be taken; The inclusion information export module is used to export inclusion information at the exit of the continuous casting slab calculation domain.

9. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method for simulating the dynamic distribution of inclusions in a continuous casting billet as claimed in any one of claims 1 to 7 are realized.

Images