Method and device for modeling and simulation of oxygen-enriched side-blown smelting process based on key field coupling

Abstract

This invention discloses a modeling and simulation method and apparatus for oxygen-enriched side-blown smelting processes based on key field coupling, relating to the field of non-ferrous metal smelting process simulation technology. It includes: decomposing the oxygen-enriched side-blown smelting process into multiple stages; extracting features from each stage to obtain the key fields for each stage; and performing coupling calculations on the key fields of each stage to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process. This invention decomposes the complex multiphase, multi-field coupling process in side-blown smelting into multiple stages. For the key fields involved in feature extraction of different stages, it weakens or simplifies non-key fields, performs coupling calculations only on the key fields, and employs specific methods for coupling calculations on the key fields of different stages to achieve efficient and accurate analysis of the side-blown smelting process.

Description

Technical Field

[0001] This invention relates to the field of simulation technology for non-ferrous metal smelting processes, and in particular to a modeling and simulation method and apparatus for oxygen-enriched side-blown smelting processes based on key field coupling. Background Technology

[0002] The oxygen-enriched smelting process used in copper smelters can be broadly categorized into flash smelting and pool smelting. Pool smelting, based on the location of the lance, can be further divided into submerged top-blown pool smelting, oxygen bottom-blown pool smelting, and oxygen-enriched side-blown pool smelting. Oxygen-enriched side-blown smelting technology originated from the Vanyukov process in the Soviet Union and boasts advantages such as low energy consumption, strong raw material adaptability, high metal recovery rate, and low investment, making it increasingly prevalent in the copper smelting industry. In the side-blown furnace, oxygen-enriched air with a concentration of 75%–80% is blown into the lower primary tuyeres, creating vigorous agitation in the molten pool, enhancing heat and mass transfer, and strengthening the smelting kinetics. Low-pressure air is blown into the upper secondary tuyeres using lances, where the combustion of monomeric sulfur and carbon monoxide in the flue gas releases a large amount of heat, which is beneficial for subsequent acid production. After the furnace charge is fed into the side-blown furnace, the slag and matte formation reactions are completed under the agitation of the blown air and the high temperature, with slag and matte separation occurring in the lower part of the hearth. However, in actual production, the increasing complexity of copper concentrate raw material composition, the increase in feed amount after capacity expansion, and the different arrangements of spray guns will all cause changes in the flow and mixing state of the melt in the furnace, resulting in poor stirring of the molten pool and reduced smelting efficiency.

[0003] Side-blown molten pool smelting is a multiphase chemical reaction and mass and heat transfer process involving a complex gas-liquid-solid system. This includes the movement of ore particles during concentrate feeding, multiphase flow of gas, matte, and slag throughout the furnace and injection zone, component mixing, and reaction processes. The stirring and mixing of the melt in high-temperature systems is difficult to observe and measure directly. Therefore, studies on melt flow within the furnace, furnace structure optimization, and improvement of smelting efficiency often employ simulation methods.

[0004] With the rapid improvement of computing power, numerical simulation has become an indispensable research tool in many fields and has been widely used in the study of multiphase systems in metallurgy. While the research on multiphase flow calculation models suitable for side-blown smelting processes is relatively complete, numerical simulation studies on smelting reactions, component diffusion, and particle movement within side-blown furnaces are limited, failing to meet the requirements for accurate simulation calculations of the side-blown molten pool process. Since side-blown smelting is a complex mass and heat transfer process, multi-field coupled calculations are a better solution for achieving accurate analysis of the multiphase reaction flow field, thermal field, and component field within the furnace. However, solving multiple strongly coupled equations requires enormous computational resources, and the convergence of the calculation results decreases significantly. Summary of the Invention

[0005] This invention addresses the problems of complex models, difficult convergence, and long computation time encountered in modeling and simulating side-blown connected multi-field strongly coupled systems.

[0006] To solve the above-mentioned technical problems, the present invention provides the following technical solution:

[0007] On one hand, this invention provides a modeling and simulation method for an oxygen-enriched side-blown melting process based on key field coupling. This method is implemented by electronic devices and includes:

[0008] S1. The oxygen-enriched side-blown smelting process is broken down into multiple steps.

[0009] S2. Extract features from each of the multiple stages to obtain the key fields of each stage.

[0010] S3. Coupled calculations are performed on the key fields of each stage to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process.

[0011] Optionally, the multiple stages in S1 include: the penetration process of side-blown gas in liquid, the multiphase movement of gas and slag in side-blown smelting, the component mixing in the blowing process, the chemical reaction in the blowing process, the feeding of concentrate particles, and the scouring and corrosion of tuyeres bricks by gas and melt.

[0012] Optionally, the key fields for the multiphase motion of the gas matte slag in the side-blown melting process in S2 include: the multiphase flow field of the gas matte slag.

[0013] The key fields for component mixing in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0014] The key fields of chemical reaction in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0015] The key fields for feeding concentrate particles include: the multiphase flow field of gas matte slag and the particle motion field.

[0016] The key fields for the erosion and corrosion of tuyer bricks by gas and melt include: the gas-slag multiphase flow field, the temperature field, and the solid field.

[0017] Optionally, the key fields in S3 for the penetration process of the side-blown gas in the liquid are coupled and calculated, including:

[0018] The flow region of melt and gas exceeding a preset threshold in the furnace is selected as the computational domain for geometric modeling. The mesh on the wall at the tuyeres is refined. Pressure inlet boundary conditions are used for gas inlet and pressure outlet boundary conditions are used for gas outlet.

[0019] Steady-state calculations were performed using a volumetric multiphase flow (VOF) model and a realizable k-ε turbulence model, bringing the equation residuals to converge to 10. -4 ~10-5 .

[0020] Optionally, the key fields of the multiphase motion of the gas matte slag in the side-blown smelting process in S3 are coupled and calculated, including:

[0021] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0022] The pressure-velocity coupled algorithm was used for calculation, and a first-order upwind difference scheme was employed to converge the equation residuals to 10. -3 The time step was set to 10 during the calculation. -4 .

[0023] Optionally, the key fields for component mixing in the injection process in S3 are coupled and calculated, including:

[0024] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0025] The mean field was calculated using the VOF model and the standard k-ε turbulence model. The mean field is the average value of fluid velocity, gas-matte multiphase flow field distribution, turbulent kinetic energy and turbulent dissipation rate over the statistical time after the flow field stabilizes.

[0026] The mean field is used as the transient flow field within the molten pool. A component of a predetermined concentration is added at any location within the molten pool, and the concentration-diffusion equation of the component is solved under the mean flow field. Transient calculations are employed, with a time step set to 10. -1 ~10 2 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0027] Optionally, the key fields of the chemical reaction in the injection process in S3 are coupled and calculated, including:

[0028] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0029] The mean field was calculated using the VOF model and the standard k-ε turbulence model.

[0030] Mass and heat transfer equations for the gas-liquid two-phase reaction in the oxygen-enriched side-blown smelting process were established. The mass transfer process within the molten pool was calculated using the component transport equation, and the heat transfer was calculated using the energy equation. Transient calculations were performed with a time step of 10. -3 ~10 -1 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0031] For the matte-forming and slag-forming reactions in the molten pool, the eddy dissipation model is selected, and the turbulent flow rate is used instead of the chemical reaction rate to solve for the concentration changes of the substances.

[0032] Optionally, the key fields for concentrate particle feeding in S3 are coupled and calculated, including:

[0033] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0034] The VOF model was used for calculations, and the particle motion was calculated using the Discrete Element Method-Fluid Volume Model (DEM-VOF). The momentum exchange process was achieved through the drag effect between the liquid and solid phases. Transient calculations were performed based on the steady-state flow field obtained from the VOF model. The time step for the multiphase fluid was set to 10 during the calculations. -3 ~10 -4 s, the residuals of the equation converge to 10 -3 ~10 -4 The time step for particle motion calculation is set to 10. -5 ~10 -4 s.

[0035] Optionally, the key fields for the erosion corrosion of the tuyeres bricks by gas and melt are coupled in S3, including:

[0036] A model of the wind gap area was established using the finite element method.

[0037] The key fields of each stage were loaded into the wind vent region model, and steady-state calculations were used. The equation residuals converged to 10. -4 ~10 -5 .

[0038] On the other hand, the present invention provides a modeling and simulation device for an oxygen-enriched side-blown smelting process based on key field coupling. This device is used to implement a modeling and simulation method for an oxygen-enriched side-blown smelting process based on key field coupling. The device includes:

[0039] The decomposition module is used to break down the oxygen-enriched side-blown smelting process into multiple stages.

[0040] The feature extraction module is used to extract features from each of the multiple stages to obtain the key fields of each stage.

[0041] The output module is used to perform coupled calculations on the key fields of each stage, enabling efficient and accurate analysis of the oxygen-enriched side-blown smelting process.

[0042] Optionally, multiple stages include: the penetration process of side-blown gas in the liquid, the multiphase movement of gas and slag in the side-blown smelting process, the mixing of components in the blowing process, the chemical reaction in the blowing process, the feeding of concentrate particles, and the scouring and corrosion of the tuyeres brick by gas and melt.

[0043] Optionally, the key fields for the multiphase motion of gas matte slag in the side-blown smelting process include: the multiphase flow field of gas matte slag.

[0044] The key fields for component mixing in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0045] The key fields of chemical reaction in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0046] The key fields for feeding concentrate particles include: the multiphase flow field of gas matte slag and the particle motion field.

[0047] The key fields for the erosion and corrosion of tuyer bricks by gas and melt include: the gas-slag multiphase flow field, the temperature field, and the solid field.

[0048] Optionally, the output module is further used for:

[0049] The flow region of melt and gas exceeding a preset threshold in the furnace is selected as the computational domain for geometric modeling. The mesh on the wall at the tuyeres is refined. Pressure inlet boundary conditions are used for gas inlet and pressure outlet boundary conditions are used for gas outlet.

[0050] Steady-state calculations were performed using a volumetric multiphase flow (VOF) model and a realizable k-ε turbulence model, bringing the equation residuals to converge to 10. -4 ~10 -5 .

[0051] Optionally, the output module is further used for:

[0052] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0053] The pressure-velocity coupled algorithm was used for calculation, and a first-order upwind difference scheme was employed to converge the equation residuals to 10. -3 The time step was set to 10 during the calculation. -4 .

[0054] Optionally, the output module is further used for:

[0055] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0056] The mean field was calculated using the VOF model and the standard k-ε turbulence model. The mean field is the average value of fluid velocity, gas-matte multiphase flow field distribution, turbulent kinetic energy and turbulent dissipation rate over the statistical time after the flow field stabilizes.

[0057] The mean field is used as the transient flow field within the molten pool. A component of a predetermined concentration is added at any location within the molten pool, and the concentration-diffusion equation of the component is solved under the mean flow field. Transient calculations are employed, with a time step set to 10. -1 ~10 2 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0058] Optionally, the output module is further used for:

[0059] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0060] The mean field was calculated using the VOF model and the standard k-ε turbulence model.

[0061] Mass and heat transfer equations for the gas-liquid two-phase reaction in the oxygen-enriched side-blown smelting process were established. The mass transfer process within the molten pool was calculated using the component transport equation, and the heat transfer was calculated using the energy equation. Transient calculations were performed with a time step of 10. -3 ~10 -1 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0062] For the matte-forming and slag-forming reactions in the molten pool, the eddy dissipation model is selected, and the turbulent flow rate is used instead of the chemical reaction rate to solve for the concentration changes of the substances.

[0063] Optionally, the output module is further used for:

[0064] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0065] The VOF model was used for calculations, and the particle motion was calculated using the Discrete Element Method-Fluid Volume Model (DEM-VOF). The momentum exchange process was achieved through the drag effect between the liquid and solid phases. Transient calculations were performed based on the steady-state flow field obtained from the VOF model. The time step for the multiphase fluid was set to 10 during the calculations. -3 ~10 -4 s, the residuals of the equation converge to 10 -3 ~10 -4 The time step for particle motion calculation is set to 10. -5 ~10 -4 s.

[0066] Optionally, the output module is further used for:

[0067] A model of the wind gap area was established using the finite element method.

[0068] The key fields of each stage were loaded into the wind vent region model, and steady-state calculations were used. The equation residuals converged to 10. -4 ~10 -5 .

[0069] On the one hand, an electronic device is provided, comprising a processor and a memory, wherein the memory stores at least one instruction, which is loaded and executed by the processor to implement the above-mentioned modeling and simulation method for oxygen-enriched side-blown melting process based on key field coupling.

[0070] On the one hand, a computer-readable storage medium is provided, wherein at least one instruction is stored in the storage medium, and the at least one instruction is loaded and executed by a processor to implement the above-described modeling and simulation method for oxygen-enriched side-blown melting process based on key field coupling.

[0071] The above technical solution has at least the following advantages compared with the existing technology:

[0072] The above scheme, oxygen-enriched side-blown molten pool smelting is a multiphase and multi-field coupled system involving a complex gas-liquid-solid system. In order to obtain the in-furnace blowing and stirring and the flow and reaction of multiphase fluids, a large number of calculations are required. (1) Modeling is difficult, multiple equations need to be coupled, and the interaction relationship between multiple fields is not yet clear; (2) The amount of calculation is large, and the time step needs to be controlled when solving the above equations, which incurs a high time cost; (3) The results are divergent, and it is difficult to converge the coupled calculation of multiple equations, and the accuracy of the results needs to be verified.

[0073] The modeling and simulation method based on key field coupling proposed in this invention can solve these problems well. (1) The complex reaction process is decomposed into multiple links, and only the extracted key fields are coupled for calculation, which reduces the modeling difficulty. (2) By processing the coupling between key fields, the computational cost can be reduced and the time step can be increased. (3) After the model is simplified, the convergence of the solution process is greatly improved, and it also has good accuracy. Attached Figure Description

[0074] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0075] Figure 1 This is a schematic diagram of the modeling and simulation method for oxygen-enriched side-blown melting process based on key field coupling provided in the embodiments of the present invention;

[0076] Figure 2 This is a schematic diagram of multiple stages provided in the embodiments of the present invention;

[0077] Figure 3 This is a block diagram of the oxygen-enriched side-blown melting process modeling and simulation device based on key field coupling provided in the embodiments of the present invention;

[0078] Figure 4 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0079] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the described embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0080] like Figure 1 As shown, this embodiment of the invention provides a modeling and simulation method for an oxygen-enriched side-blown melting process based on key field coupling, which can be implemented by electronic devices. Figure 1 The flowchart shown is a modeling and simulation method for oxygen-enriched side-blown melting process based on key field coupling. The processing flow of this method may include the following steps:

[0081] S1. The oxygen-enriched side-blown smelting process is broken down into multiple steps.

[0082] like Figure 2 As shown, multiple stages may include: the penetration process of side-blown gas in the liquid, the multiphase movement of gas and slag in the side-blown smelting process, the mixing of components in the blowing process, the chemical reaction in the blowing process, the feeding of concentrate particles, and the scouring and corrosion of the tuyeres brick by gas and melt.

[0083] S2. Extract features from each of the multiple stages to obtain the key fields of each stage.

[0084] Among them, the key field for the multiphase motion of gas matte slag in the side-blown smelting process can include: the multiphase flow field of gas matte slag.

[0085] The key fields for component mixing in the injection process can include: the multiphase flow field of the gas-slag mixture and the component fields in each phase.

[0086] The key fields of chemical reaction in the injection process can include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0087] The key fields for feeding concentrate particles can include: the multiphase flow field of gas matte slag and the particle motion field.

[0088] The key fields for the erosion and corrosion of tuyer bricks by gas and melt can include: gas-slag multiphase flow field, temperature field and solid field.

[0089] S3. Coupled calculations are performed on the key fields of each stage to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process.

[0090] Optionally, the key fields in S3 for the penetration process of the side-blown gas in the liquid are coupled and calculated, including:

[0091] The flow region of melt and gas exceeding a preset threshold in the furnace is selected as the computational domain for geometric modeling. The mesh on the wall at the tuyeres is refined. Pressure inlet boundary conditions are used for gas inlet and pressure outlet boundary conditions are used for gas outlet.

[0092] Steady-state calculations were performed using a volumetric multiphase flow (VOF) model and a realizable k-ε turbulence model, and the equation residuals were converged to 10. -4 ~10 -5 .

[0093] In one feasible implementation, step 1 is the simulation of the flow state of the gas injected at the nozzle after entering the melt. For the penetration process of the side-blown gas in the liquid in step (1), the penetration process of the side-blown gas in the liquid is a bubbling flow pattern. Steady-state calculations are performed using the VOF (Volume of Fluid) model and the Realizable k-ε turbulence model. All residuals converge to 10. -4 ~10 -5 The main melt and gas flow regions within the furnace are selected as the computational domain. External mechanical structures and wall thickness are disregarded; geometric modeling is performed only for the fluid computational domain within the furnace. The velocity is highest near the wall, so the wall mesh at the tuyeres is appropriately refined. Temperature boundary conditions are avoided for the spray lance; pressure inlet boundary conditions are used for the gas inlet, and pressure outlet boundary conditions are used for the gas outlet. Using the above calculations, the gas penetration process in the liquid can be obtained. The main blowing area within the molten pool is divided into a high-speed blowing zone and a circulating stirring zone. The gas flow pattern exhibits a periodic change in three stages: growth, detachment, and rising and breaking up. The circulating stirring zone, located in the middle of the molten pool, is crucial for stirring and mixing the melt over a large volume.

[0094] Optionally, the key fields of the multiphase motion of the gas matte slag in the side-blown smelting process in S3 are coupled and calculated, including:

[0095] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary. Multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0096] The pressure-velocity coupled algorithm was used for calculation, and a first-order upwind difference scheme was employed to converge the equation residuals to 10. -3 The time step was set to 10 during the calculation. -4 .

[0097] In one feasible implementation, the gas-matte-slag multiphase motion in the side-blown smelting process of stage (2) involves the coupling between the gas, matte, and slag multiphase fields. The boundary conditions used in the calculation are determined based on the actual production conditions of a certain factory. Among them, the wall condition is set as a no-slip boundary, the nozzle inlet is set as a velocity inlet, the velocity of a single nozzle inlet is 100-250m / s, the turbulence intensity is 3%-5%, and the flue gas outlet is set as a pressure outlet. The gas flow penetration depth obtained in stage (1) can be used as the region boundary. In the equation solving process, the pressure-velocity coupling algorithm is selected, and the first-order upwind difference scheme is adopted. The equation residual convergence standard is 10. -3 The time step was set to 10 during the calculation. -4 The above calculations yield the results of the changes in the velocity field, turbulence field, and gas-mând slag phase field distribution within the furnace under different injection conditions with respect to injection time. To ensure the stability of the settling separation zone and the gas-liquid stirring and circulation intensity of the weak stirring zone and the circulating stirring zone, several circulating injection technologies were proposed, including parallel circulating injection (horizontally deflecting the spray guns on both sides of the furnace wall by 15° in opposite directions), staggered circulating injection (horizontally shifting the spray guns on both sides of the furnace wall by 100mm in opposite directions), and differential circulating injection (horizontally shifting the spray guns at opposite positions on both sides of the furnace wall by 100mm in opposite directions), so that the melt in the injection zone forms large or small circulations under the action of gas. Among them, the circulation generated by parallel circulating injection has the widest influence range and plays a driving role in the entire gas-liquid stirring zone.

[0098] Optionally, the key fields for component mixing in the injection process in S3 are coupled and calculated, including:

[0099] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary. Multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0100] The mean field was calculated using the VOF model and the standard k-ε turbulence model. The mean field is the average value of fluid velocity, gas-matte multiphase flow field distribution, turbulent kinetic energy and turbulent dissipation rate over the statistical time after the flow field stabilizes.

[0101] The mean field is used as the transient flow field within the molten pool. A component of a predetermined concentration is added at any location within the molten pool, and the concentration-diffusion equation of the component is solved under the mean flow field. Transient calculations are employed, with a time step set to 10. -1 ~10 2 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0102] In one feasible implementation, the component mixing in the injection process of step (3) involves the coupling between the multiphase flow field of the gas matte slag and the component fields in each phase. Based on the model established in step (2), the VOF model and the standard k-ε turbulence model are used. After the flow field stabilizes, the average values ​​of fluid velocity, gas matte slag phase distribution, turbulent kinetic energy and turbulent dissipation rate under statistical time are calculated, which is the average field. This is used as the transient flow field in the molten pool, and a certain concentration of component is added at a certain position in the molten pool. The calculation of the flow field equations (continuity equation, momentum equation and turbulence equation) is turned off, and the concentration diffusion equation of the component is solved only under the average flow field. Transient calculation is used, with a time step of 10. -1 ~10 2 s, all residuals converge to 10 -4 ~10 -5 The above calculations can be used to obtain the diffusion and concentration distribution of components in each phase under different injection conditions, as well as the component mixing at different times and locations in the side-blown furnace.

[0103] Optionally, the key fields of the chemical reaction in the injection process in S3 are coupled and calculated, including:

[0104] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary. Multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0105] The mean field was calculated using the VOF model and the standard k-ε turbulence model.

[0106] Mass and heat transfer equations for the gas-liquid two-phase reaction in the oxygen-enriched side-blown smelting process were established. The mass transfer process within the molten pool was calculated using the component transport equation, and the heat transfer was calculated using the energy equation. Transient calculations were performed with a time step of 10. -3 ~10 -1 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0107] For the matte-forming and slag-forming reactions in the molten pool, the eddy dissipation model is selected, and the turbulent flow rate is used instead of the chemical reaction rate to solve for the concentration changes of the substances.

[0108] In one feasible implementation, the chemical reaction of the injection process in step (4) includes an exothermic heterogeneous reaction and a homogeneous reaction for matte and slag formation, involving the coupling between the multiphase flow field of the gas-matte-slag mixture and the component field. Based on the model established in step (2), the mean field is first obtained through statistics. The mass and heat transfer equations for the gas-liquid two-phase reaction in the smelting process are established, and the mass transfer process in the molten pool is calculated through the component transport equations; considering that the reaction is an exothermic process, the heat transfer is calculated by the energy equation. Transient calculation is adopted, with a time step of 10. -3 ~10 -1 s, all residuals converge to 10 -4 ~10 -5 The matte-forming and slag-forming reactions within the molten pool occur within the slag phase. These reactions are rapid at high temperatures. A vortex dissipation model is chosen, and turbulent flow rates are used to replace chemical reaction rates to solve for the concentration changes of substances. Using the above calculations, the concentration and temperature fields within a side-blown furnace under different injection conditions can be simulated, revealing the concentration changes of reactants and products during the reaction process.

[0109] Optionally, the key fields for concentrate particle feeding in S3 are coupled and calculated, including:

[0110] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary. Multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0111] The VOF model was used for calculations, and the discrete element method-fluid volume model (DEM-VOF) was employed to calculate particle motion. The momentum exchange process was achieved through the drag force between the liquid and solid phases. Transient calculations were performed based on the steady-state flow field obtained from the VOF model. The time step for the multiphase fluid was set to 10 during the calculations. -3 ~10 -4 s, the residuals of the equation converge to 10 -3 ~10 -4 The time step for particle motion calculation is set to 10. -5 ~10 -4 s.

[0112] In one feasible implementation, the concentrate particles are fed in step (5), which involves the coupling between the multiphase flow field of the gas matte slag and the particle motion field. Based on the model established in step (2), the multiphase flow field is calculated using the VOF model; the particle motion is bidirectionally coupled using DEM-VOF (Discrete Element Method-Volume of Fluid Model), and the momentum exchange process is achieved through the drag effect between the liquid and solid phases. Transient calculations are performed based on the stable flow field calculated by the VOF model. The time step of the multiphase fluid is 10. -3 ~10 -4 s, all residuals converge to 10 -3 ~10 -4 The time step for particle calculations is 10. -5 ~10 -4 The Fluent time step is set to an integer multiple of the EDEM time step to facilitate data exchange between the two. Using the above calculations, the motion and reaction of concentrate particles under the influence of the gas-matte multiphase flow field can be obtained under different injection conditions.

[0113] Optionally, the key fields for the erosion corrosion of the tuyeres bricks by gas and melt are coupled in S3, including:

[0114] A model of the wind gap area was established using the finite element method.

[0115] The key fields of each stage were loaded into the wind vent region model, and steady-state calculations were used. The equation residuals converged to 10. -4 ~10 -5 .

[0116] In one feasible implementation, the erosion corrosion of the tuyer brick by gas and melt in step (6) involves coupled calculations of multiphase flow field, temperature field, and solid field. After the jet gas exits the tuyer, it expands, forming turbulence and eddies in the circulating stirring zone, which drive the high-temperature melt to impact the tuyer and the nearby furnace lining. A tuyer region model is established using the finite element method, and the key fields mentioned above are loaded into the finite element model. Steady-state calculation is used, and all residuals converge to 10. -4 ~10 -5 The above calculations can be used to obtain the pressure and stress distribution of the tuyeres bricks under different injection conditions, and to analyze the damage to the furnace lining.

[0117] Furthermore, by performing key field coupling calculations in the above-mentioned steps, relatively accurate simulation results for different field information can be obtained relatively efficiently.

[0118] In this embodiment of the invention, oxygen-enriched side-blown molten pool smelting is a multiphase and multi-field coupled system involving a complex gas-liquid-solid system. In order to obtain the in-furnace blowing and stirring and the flow and reaction of multiphase fluids, a large number of calculations are required. (1) Modeling is difficult, multiple equations need to be coupled, and the interaction relationship between multiple fields is not yet clear; (2) The amount of calculation is large, and the time step needs to be controlled when solving the above equations, which incurs a high time cost; (3) The results are divergent, and it is difficult to converge the coupled calculation of multiple equations, and the accuracy of the results needs to be verified.

[0119] The modeling and simulation method based on key field coupling proposed in this invention can solve these problems well. (1) The complex reaction process is decomposed into multiple links, and only the extracted key fields are coupled for calculation, which reduces the modeling difficulty. (2) By processing the coupling between key fields, the computational cost can be reduced and the time step can be increased. (3) After the model is simplified, the convergence of the solution process is greatly improved, and it also has good accuracy.

[0120] like Figure 3 As shown, this embodiment of the invention provides a modeling and simulation device 300 for an oxygen-enriched side-blown smelting process based on key field coupling. This device 300 is used to implement a modeling and simulation method for an oxygen-enriched side-blown smelting process based on key field coupling. The device 300 includes:

[0121] The decomposition module 310 is used to decompose the oxygen-enriched side-blown smelting process into multiple stages.

[0122] The feature extraction module 320 is used to extract features from each of the multiple stages to obtain the key fields of each stage.

[0123] Output module 330 is used to perform coupled calculations on the key fields of each stage, so as to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process.

[0124] Optionally, multiple stages include: the penetration process of side-blown gas in the liquid, the multiphase movement of gas and slag in the side-blown smelting process, the mixing of components in the blowing process, the chemical reaction in the blowing process, the feeding of concentrate particles, and the scouring and corrosion of the tuyeres brick by gas and melt.

[0125] Optionally, the key fields for the multiphase motion of gas matte slag in the side-blown smelting process include: the multiphase flow field of gas matte slag.

[0126] The key fields for component mixing in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0127] The key fields of chemical reaction in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase.

[0128] The key fields for feeding concentrate particles include: the multiphase flow field of gas matte slag and the particle motion field.

[0129] The key fields for the erosion and corrosion of tuyer bricks by gas and melt include: the gas-slag multiphase flow field, the temperature field, and the solid field.

[0130] Optionally, the output module 330 is further used for:

[0131] The flow region of melt and gas exceeding a preset threshold in the furnace is selected as the computational domain for geometric modeling. The mesh on the wall at the tuyeres is refined. Pressure inlet boundary conditions are used for gas inlet and pressure outlet boundary conditions are used for gas outlet.

[0132] Steady-state calculations were performed using the volumetric multiphase flow (VOF) model and the realizable k-ε turbulence model, and the equation residuals were converged to 10. -4 ~10 -5 .

[0133] Optionally, the output module 330 is further used for:

[0134] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0135] The pressure-velocity coupled algorithm was used for calculation, and a first-order upwind difference scheme was employed to converge the equation residuals to 10. -3 The time step was set to 10 during the calculation. -4 .

[0136] Optionally, the output module 330 is further used for:

[0137] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0138] The mean field was calculated using the VOF model and the standard k-ε turbulence model. The mean field is the average value of fluid velocity, gas-matte multiphase flow field distribution, turbulent kinetic energy and turbulent dissipation rate over the statistical time after the flow field stabilizes.

[0139] The mean field is used as the transient flow field within the molten pool. A component of a predetermined concentration is added at any location within the molten pool, and the concentration-diffusion equation of the component is solved under the mean flow field. Transient calculations are employed, with a time step set to 10.-1 ~10 2 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0140] Optionally, the output module 330 is further used for:

[0141] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0142] The mean field was calculated using the VOF model and the standard k-ε turbulence model.

[0143] Mass and heat transfer equations for the gas-liquid two-phase reaction in the oxygen-enriched side-blown smelting process were established. The mass transfer process within the molten pool was calculated using the component transport equation, and the heat transfer was calculated using the energy equation. Transient calculations were performed with a time step of 10. -3 ~10 -1 s, the residuals of the equation converge to 10 -4 ~10 -5 .

[0144] For the matte-forming and slag-forming reactions in the molten pool, the eddy dissipation model is selected, and the turbulent flow rate is used instead of the chemical reaction rate to solve for the concentration changes of the substances.

[0145] Optionally, the output module 330 is further used for:

[0146] The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The wall condition is set as a no-slip boundary, multiple spray gun inlets are set as velocity inlets, the velocity of a single spray gun inlet is set to 100-250 m / s, the turbulence intensity of a single spray gun inlet is set to 3%-5%, and the flue gas outlet is set as a pressure outlet.

[0147] The VOF model was used for calculations, and the discrete element method-fluid volume model (DEM-VOF) was employed to calculate particle motion. The momentum exchange process was achieved through the drag force between the liquid and solid phases. Transient calculations were performed based on the steady-state flow field obtained from the VOF model. The time step for the multiphase fluid was set to 10 during the calculations. -3 ~10 -4 s, the residuals of the equation converge to 10 -3 ~10 -4 The time step for particle motion calculation is set to 10. -5 ~10 -4s.

[0148] Optionally, the output module 330 is further used for:

[0149] A model of the wind gap area was established using the finite element method.

[0150] The key fields of each stage were loaded into the wind vent region model, and steady-state calculations were used. The equation residuals converged to 10. -4 ~10 -5 .

[0151] In this embodiment of the invention, oxygen-enriched side-blown molten pool smelting is a multiphase and multi-field coupled system involving a complex gas-liquid-solid system. In order to obtain the in-furnace blowing and stirring and the flow and reaction of multiphase fluids, a large number of calculations are required. (1) Modeling is difficult, multiple equations need to be coupled, and the interaction relationship between multiple fields is not yet clear; (2) The amount of calculation is large, and the time step needs to be controlled when solving the above equations, which incurs a high time cost; (3) The results are divergent, and it is difficult to converge the coupled calculation of multiple equations, and the accuracy of the results needs to be verified.

[0152] The modeling and simulation method based on key field coupling proposed in this invention can solve these problems well. (1) The complex reaction process is decomposed into multiple links, and only the extracted key fields are coupled for calculation, which reduces the modeling difficulty. (2) By processing the coupling between key fields, the computational cost can be reduced and the time step can be increased. (3) After the model is simplified, the convergence of the solution process is greatly improved, and it also has good accuracy.

[0153] Figure 4 This is a schematic diagram of the structure of an electronic device 400 provided in an embodiment of the present invention. The electronic device 400 can vary considerably due to different configurations or performance. It may include one or more central processing units (CPUs) 401 and one or more memories 402. The memory 402 stores at least one instruction, which is loaded and executed by the processor 401 to implement the following modeling and simulation method for oxygen-enriched side-blown melting process based on key field coupling:

[0154] S1. The oxygen-enriched side-blown smelting process is broken down into multiple steps.

[0155] S2. Extract features from each of the multiple stages to obtain the key fields of each stage.

[0156] S3. Coupled calculations are performed on the key fields of each stage to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process.

[0157] In an exemplary embodiment, a computer-readable storage medium is also provided, such as a memory including instructions that can be executed by a processor in a terminal to complete the above-described modeling and simulation method for oxygen-enriched side-blown melting process based on key field coupling. For example, the computer-readable storage medium may be a ROM, random access memory (RAM), CD-ROM, magnetic tape, floppy disk, or optical data storage device.

[0158] Those skilled in the art will understand that all or part of the steps of the above embodiments can be implemented by hardware or by a program instructing related hardware. The program can be stored in a computer-readable storage medium, such as a read-only memory, a disk, or an optical disk.

[0159] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A modeling and simulation method for oxygen-enriched side-blown smelting process based on key field coupling, characterized in that, The method includes: S1. The oxygen-enriched side-blown smelting process is broken down into multiple steps; S2. Extract features from each of the multiple steps to obtain the key fields of each step; S3. Perform coupled calculations on the key fields of each step to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process; The multiple steps in S1 include: the penetration process of side-blown gas in liquid, the multiphase movement of gas and slag in side-blown smelting, the component mixing in the blowing process, the chemical reaction in the blowing process, the feeding of concentrate particles, and the scouring and corrosion of tuyeres bricks by gas and melt. The key fields for the multiphase motion of gas matte slag in the side-blown smelting process in S2 include: the multiphase flow field of gas matte slag; The key fields for component mixing during the injection process include: the multiphase flow field of the gas-molten metal and the component fields in each phase; The key fields of chemical reaction in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase; The key fields for concentrate particle feeding include: the multiphase flow field of gas matte slag and the particle motion field; The key fields for the erosion and corrosion of tuyer bricks by gas and melt include: gas-molten slag multiphase flow field, temperature field, and solid field; The coupling calculation of the key fields of each link in S3 includes: The flow region of melt and gas exceeding a preset threshold within the furnace is selected as the computational domain for geometric modeling. The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The tuyeres region model is established using the finite element method.

2. The method according to claim 1, characterized in that, The key field coupling calculation in S3 for the penetration process of the side-blown gas in the liquid includes: The wall mesh at the air outlet is refined, and the gas inlet adopts the pressure inlet boundary condition, while the gas outlet adopts the pressure outlet boundary condition. The fluid volume multiphase flow (VOF) model and Realizable were adopted. k-ε Steady-state calculations were performed using the turbulence model, and the equation residuals converged to 10. -4 ~10 -5 .

3. The method according to claim 1, characterized in that, The key field coupling calculation in S3 for the multiphase motion of gas-molten matte slag in the side-blown smelting process includes: Set the wall conditions to a no-slip boundary, set multiple spray gun inlets as velocity inlets, set the velocity of a single spray gun inlet to 100-250 m / s, set the turbulence intensity of a single spray gun inlet to 3%-5%, and set the flue gas outlet as a pressure outlet. The pressure-velocity coupled algorithm was used for calculation, and a first-order upwind difference scheme was employed to converge the equation residuals to 10. -3 The time step was set to 10 during the calculation. -4 .

4. The method according to claim 1, characterized in that, The coupling calculation of the key field for component mixing in the jetting process in S3 includes: Set the wall conditions to a no-slip boundary, set multiple spray gun inlets as velocity inlets, set the velocity of a single spray gun inlet to 100-250 m / s, set the turbulence intensity of a single spray gun inlet to 3%-5%, and set the flue gas outlet as a pressure outlet. Using VOF model and standards k-ε The mean field is calculated using a turbulence model; wherein, the mean field is the average value of fluid velocity, gas-matte multiphase flow field distribution, turbulent kinetic energy, and turbulent dissipation rate over a statistical time period after the flow field stabilizes; The mean field is used as the transient flow field within the molten pool. A component of a predetermined concentration is added at any location within the molten pool. The concentration-diffusion equation of the component is solved under the mean flow field using transient calculations, with a time step set to 10. -1 ~10 2 s, the residuals of the equation converge to 10 -4 ~10 -5 .

5. The method according to claim 1, characterized in that, The coupling calculation of key fields in the injection process chemical reaction in S3 includes: Set the wall conditions to a no-slip boundary, set multiple spray gun inlets as velocity inlets, set the velocity of a single spray gun inlet to 100-250 m / s, set the turbulence intensity of a single spray gun inlet to 3%-5%, and set the flue gas outlet as a pressure outlet. Using VOF model and standards k-ε Turbulence model for calculating mean field; Mass and heat transfer equations for the gas-liquid two-phase reaction in the oxygen-enriched side-blown smelting process were established. The mass transfer process within the molten pool was calculated using the component transport equation, and the heat transfer was calculated using the energy equation. Transient calculations were performed with a time step of 10. -3 ~10 -1 s, the residuals of the equation converge to 10 -4 ~10 -5 ; For the matte-forming and slag-forming reactions in the molten pool, the eddy dissipation model is selected, and the turbulent flow rate is used instead of the chemical reaction rate to solve for the concentration changes of the substances.

6. The method according to claim 1, characterized in that, The coupling calculation of the key field for concentrate particle feeding in S3 includes: Set the wall conditions to a no-slip boundary, set multiple spray gun inlets as velocity inlets, set the velocity of a single spray gun inlet to 100-250 m / s, set the turbulence intensity of a single spray gun inlet to 3%-5%, and set the flue gas outlet as a pressure outlet. The VOF model was used for calculations, and the particle motion was calculated using the Discrete Element Method-Fluid Volume Model (DEM-VOF). The momentum exchange process was achieved through the drag effect between the liquid and solid phases. Transient calculations were performed based on the steady-state flow field obtained from the VOF model. The time step for the multiphase fluid was set to 10 during the calculations. -3 ~10 -4 s, the residuals of the equation converge to 10 -3 ~10 -4 The time step for particle motion calculation is set to 10. -5 ~10 -4 s.

7. The method according to claim 1, characterized in that, The key field coupling calculations in S3 for the erosion corrosion of the tuyeres bricks by gas and melt include: The key fields of each stage are loaded into the wind vent region model, and steady-state calculations are used. The equation residuals converge to 10. -4 ~10 -5 .

8. A modeling and simulation device for oxygen-enriched side-blown smelting process based on key field coupling, characterized in that, The device includes: The decomposition module is used to break down the oxygen-enriched side-blown smelting process into multiple stages. The feature extraction module is used to extract features from each of the multiple steps to obtain the key fields of each step. The output module is used to perform coupled calculations on the key fields of each step, so as to achieve efficient and accurate analysis of the oxygen-enriched side-blown smelting process. The multiple stages in S1 include: the penetration process of side-blown gas in liquid, the multiphase movement of gas and slag in side-blown smelting, the component mixing in the blowing process, the chemical reaction in the blowing process, the feeding of concentrate particles, and the scouring and corrosion of tuyeres bricks by gas and melt. The key fields for the multiphase motion of gas matte slag in the side-blown melting process in S2 include: the multiphase flow field of gas matte slag; The key fields for component mixing during the injection process include: the multiphase flow field of the gas-molten metal and the component fields in each phase; The key fields of chemical reaction in the injection process include: the multiphase flow field of the gas-molten laitance and the component fields in each phase; The key fields for concentrate particle feeding include: the multiphase flow field of gas matte slag and the particle motion field; The key fields for the erosion and corrosion of tuyer bricks by gas and melt include: gas-molten slag multiphase flow field, temperature field, and solid field; S3 involves coupling calculations of the key fields for each stage, including: The flow region of melt and gas exceeding a preset threshold within the furnace is selected as the computational domain for geometric modeling. The gas flow penetration depth obtained through the penetration process of side-blown gas in the liquid is used as the region boundary for geometric modeling. The tuyeres region model is established using the finite element method.

Images