A method for simulating heat transfer behavior in the hot zone of converter molten pool based on Fluent-UDF

Abstract

The present invention discloses a method for simulating heat transfer behavior of a converter molten pool fire point area based on Fluent-UDF, which belongs to the technical field of iron and steel metallurgy. The method comprises: establishing a steelmaking converter geometric model and meshing the steelmaking converter geometric model, importing the meshed steelmaking converter geometric model into Fluent software, using Fluent software to perform simulation calculations, obtaining a quasi-steady-state flow field, and determining a designated area by the obtained quasi-steady-state flow field; developing a UDF program; the UDF program is used to track the fire point area morphology of the designated area, and to perform stable heat release in the designated area; importing the UDF program into Fluent for compilation, and adding the heat release rate to the Fluent energy source term equation for simulation calculation, thereby obtaining the temperature field distribution of the converter molten pool fire point area. The present invention can provide a reference for the molten pool temperature control in the converter smelting process.

Description

Technical Field

[0001] The invention relates to the technical field of iron and steel metallurgy, and in particular to a method for simulating heat transfer behavior of a molten pool fire point zone of a converter based on Fluent-UDF. Background Art

[0002] The converter steelmaking method is the main steelmaking method. In recent years, with the development and adoption of a series of new processes and equipment such as the converter sliding nozzle slag blocking process, converter gas analysis and blowing control, and converter automatic steelmaking, the converter steelmaking process control level has been greatly improved. However, in the actual smelting process, there are still problems such as inaccurate temperature control of the converter molten pool process, serious erosion of the converter refractory lining, and low hit rate of the chemical composition and temperature of the molten steel at the end of smelting. In essence, this is due to the lack of a deep understanding of the flow and heat transfer behavior of the multiphase fluid inside the converter molten pool. The fluctuations in the quality of the molten steel at the end of the steelmaking process caused by the above problems will cause fluctuations in the quality of the products in the subsequent production processes. This phenomenon has also attracted widespread attention.

[0003] The converter smelting process occurs in a complex high-temperature multi-component multi-phase system, involving phenomena such as the flow and heat transfer of multi-phase fluids. A deep understanding of the above phenomena is conducive to the control of the process parameters of the converter smelting process and the improvement of the converter smelting effect. However, the converter smelting environment is extremely complex, which restricts the direct analysis of the flow and heat transfer behavior of the molten steel in the converter molten pool. Therefore, researchers often study the fluid flow and heat transfer behavior inside the converter through physical simulation, numerical simulation, and theoretical derivation. The fire point area of ​​the converter molten pool is one of the main heat sources in the converter smelting process, which significantly affects the temperature of the molten steel in the molten pool. A deep understanding of the heat transfer behavior between the fire point area and the molten steel in the molten pool is crucial to the temperature control of the converter smelting process.

[0004] Due to the complexity of the heat transfer behavior of the converter molten pool fluid, there has been no report on the use of numerical simulation methods to study the heat transfer behavior between the molten pool fire point area and the molten steel in the molten pool. At present, the average temperature change law of the converter molten pool smelting process is often obtained by sampling analysis and theoretical deduction during the smelting process. The temperature distribution characteristics of the molten steel in the smelting process and the relationship between the molten pool fluid flow and heat transfer behavior are still unclear. Summary of the invention

[0005] The present invention provides a method for simulating the heat transfer behavior of the fire point zone of a converter molten pool based on Fluent-UDF. A heat transfer model of the fire point zone of the converter molten pool is established by Fluent and UDF custom function method, and the flow field and temperature field inside the converter molten pool are numerically solved to solve the technical problem of the difficulty in establishing the current heat transfer model of the converter molten pool.

[0006] In order to solve the above technical problems, the present invention provides the following technical solutions:

[0007] On the one hand, the present invention provides a method for simulating heat transfer behavior of a converter molten pool fire point area based on Fluent-UDF, and the method for simulating heat transfer behavior of a converter molten pool fire point area based on Fluent-UDF includes:

[0008] Establishing a geometric model of a steelmaking converter and meshing the geometric model of the steelmaking converter, importing the meshed geometric model of the steelmaking converter into Fluent software, performing simulation calculations using the Fluent software to obtain a quasi-steady-state flow field, and determining a designated area using the obtained quasi-steady-state flow field;

[0009] Developing a UDF (User defined functions) program; the UDF program is used to track the morphology of the fire point area in the designated area and to perform stable heat release in the designated area;

[0010] The UDF program is imported into Fluent for compilation, and the heat release rate is added to the Fluent energy source term equation for simulation calculation, so as to obtain the temperature field distribution of the fire point area of ​​the converter molten pool.

[0011] Furthermore, a geometric model of a steelmaking converter is established and meshed, the meshed geometric model of the steelmaking converter is imported into Fluent software, and simulation calculation is performed using Fluent software to obtain a quasi-steady-state flow field, and a designated area is determined by the obtained quasi-steady-state flow field, including:

[0012] Open the Fluent software, import the computational grid of the steelmaking converter computational domain, set the computational parameters, and perform preliminary simulation calculations to obtain the quasi-steady-state flow field; wherein the computational parameters include: physical property parameters, boundary conditions, solution methods, and computational residuals;

[0013] Based on the quasi-steady flow field, the radius and depth of the impact pit are obtained, and a circle with a radius R and an intersection of the symmetry axis of the steelmaking converter geometric model and the metal liquid surface as the center is determined;

[0014] The C_VOF macro is used to identify the computational grid region with a volume fraction of the metal liquid phase between 0 and 1, and the intersection of the identified computational grid region and the region enclosed by a circle with a radius of R is taken as the designated region.

[0015] Furthermore, R is a set value, which is not less than the radius and depth of the impact pit.

[0016] Furthermore, the UDF program is developed, including:

[0017] Obtain the content and activity of elements undergoing oxidation reactions in the fire point area of ​​the converter molten pool, and calculate the ratio of oxygen consumed by the reaction of unit mole of O2 in the fire point area with the corresponding element;

[0018] Based on the ratio of oxygen consumed by the reaction of unit mole of O2 in the fire point area with the corresponding element, the total heat release rate of the fire point area of ​​the molten pool is calculated;

[0019] A UDF program was developed based on the obtained total heat release rate.

[0020] Furthermore, when the temperature of the fire point zone is lower than 2386K, the elements undergoing oxidation reaction in the fire point zone of the converter molten pool include: C element, Si element, Mn element, P element and Fe element;

[0021] When the temperature of the fire point zone is not lower than 2386K, the elements undergoing oxidation reaction in the fire point zone of the converter molten pool include: C element, Si element, Mn element and Fe element.

[0022] Furthermore, the calculation formula of the oxygen ratio is:

[0023]

[0024] in, ΔG is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and element i in the molten metal; i is the Gibbs free energy of the reaction of unit mole of O2 with element i in the metal liquid; ΔG is the sum of the Gibbs free energy of the reaction of unit mole of O2 with each element in the metal liquid;

[0025] When the temperature of the fire zone is lower than 2386K, When the temperature of the fire point area is not less than 2386K, Among them, ΔG C ΔG is the Gibbs free energy per mole of O2 reacting with C in the metal liquid; Si ΔG is the Gibbs free energy per mole of O2 reacting with Si in the metal liquid; Mn ΔG is the Gibbs free energy per mole of O2 reacting with Mn in the metal liquid; Fe ΔG is the Gibbs free energy per mole of O2 reacting with Fe in the molten metal; P It is the Gibbs free energy of the reaction of unit mole of O2 with element P in the molten metal.

[0026] Furthermore, when the temperature of the fire point area is lower than 2386K, the calculation formula for the total heat release rate of the fire point area of ​​the molten pool is:

[0027]

[0028] When the temperature of the fire point area is not lower than 2386K, the calculation formula for the total heat release rate of the fire point area of ​​the molten pool is:

[0029]

[0030] Among them, q c is the total heat release rate in the fire point area of ​​the molten pool; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and element C in the molten metal; C It represents the enthalpy change when the C element participates in the reaction; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and Si element in the metal liquid; Si It represents the enthalpy change when Si element participates in the reaction; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and Fe in the molten metal; Fe It represents the enthalpy change when Fe element participates in the reaction; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and Mn element in the molten metal; Mn It represents the enthalpy change when the Mn element participates in the reaction; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and element P in the molten metal; P It represents the enthalpy change when element P participates in the reaction; Represents the molar flow rate of O2.

[0031] Furthermore, the morphology of the fire point area in the designated area is tracked, specifically: the volume fraction of the metal liquid phase in the designated area is tracked by using the C_VOF macro pre-packaged in the Fluent software to obtain the morphology change of the fire point area.

[0032] Furthermore, stable heat release is performed in the designated area, specifically: the DEFINE_SOURCE macro function pre-packaged in the Fluent software is used to perform stable heat release in the designated area at the calculated total heat release rate of the molten pool fire point area.

[0033] Furthermore, the heat release rate is added to the Fluent energy source term equation. Specifically, the calculated total heat release rate is added to the Fluent energy source term S through the DEFINE_SOURCE macro. T middle.

[0034] On the other hand, the present invention further provides an electronic device, comprising a processor and a memory; wherein the memory stores at least one instruction, and the instruction is loaded and executed by the processor to implement the above method.

[0035] In yet another aspect, the present invention further provides a computer-readable storage medium, wherein at least one instruction is stored in the storage medium, and the instruction is loaded and executed by a processor to implement the above method.

[0036] The present invention uses Fluent and UDF custom functions to establish a heat transfer model of the converter molten pool fire point area, and numerically solves the flow field and temperature field inside the converter molten pool. It can analyze the heat transfer behavior between the molten pool fire point area and the molten pool metal liquid under different working conditions and the temperature field distribution of the molten pool metal liquid, thereby providing a reference for the molten pool temperature control in the converter smelting process.

[0037] Compared with the prior art, the technical solution provided by the present invention has at least the following beneficial effects:

[0038] 1. The present invention obtains the heat release rate of the fire point area of ​​the converter molten pool at a certain moment by solving the solution, and then simulates the heat release behavior of the fire point area of ​​the molten pool by adding UDF to the energy source term equation, thereby avoiding the construction of a complex chemical reaction model in Fluent and greatly reducing the amount of calculation.

[0039] 2. Common methods such as physical simulation, sampling analysis of the converter smelting process and theoretical derivation can only obtain the influence of operating parameters on the average temperature of the converter molten pool. However, the scheme of the present invention can study the relationship between the fluid flow and heat transfer behavior of the molten pool in the converter smelting process and the influence of operating parameters on the distribution of the molten pool temperature field through numerical simulation, thereby better controlling the molten pool temperature in the converter smelting process. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0041] Figure 1 It is a schematic diagram of the execution flow of a method for simulating heat transfer behavior of a molten furnace fire point zone based on Fluent-UDF provided in an embodiment of the present invention;

[0042] Figure 2 is a schematic diagram of a two-dimensional converter model and boundary conditions provided by an embodiment of the present invention;

[0043] Figure 3 It is a schematic diagram of the quasi-steady-state flow field and impact pit morphology provided by an embodiment of the present invention;

[0044] Figure 4 Schematic diagram of the fire point area of ​​the converter molten pool provided by an embodiment of the present invention;

[0045] Figure 5 It is a cloud diagram of converter molten pool temperature distribution and flow field velocity vector direction provided by an embodiment of the present invention;

[0046] Figure 6 It is a converter molten pool velocity distribution and flow field velocity vector direction cloud map provided by an embodiment of the present invention;

[0047] Figure 7 It is a system block diagram of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0048] In order to make the objectives, technical solutions and advantages of the present invention more clear, the embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

[0049] First of all, it should be noted that in the embodiments of the present invention, words such as "exemplarily" and "for example" are used to indicate examples, illustrations or explanations. Any embodiment or design described as "exemplary" in the present invention should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of the word "exemplarily" is intended to present the concept in a concrete way. In addition, in the embodiments of the present invention, the meaning expressed by "and / or" can be both, or it can be either of the two.

[0050] Moreover, in the embodiments of the present invention, "image" and "picture" can sometimes be used interchangeably. It should be noted that when the difference between them is not emphasized, the meanings they intend to express are the same. "of", "corresponding, relevant" and "corresponding" can sometimes be used interchangeably. It should be noted that when the difference between them is not emphasized, the meanings they intend to express are the same.

[0051] In addition, in the embodiments of the present invention, sometimes a subscript (such as W1) may be mistakenly written as a non-subscript form (such as W1). When the difference is not emphasized, the meanings to be expressed are consistent.

[0052] First embodiment

[0053] This embodiment provides a method for simulating heat transfer behavior of a converter molten pool fire point area based on Fluent-UDF. The method can be implemented by an electronic device, which can be a terminal or a server. The execution flow of the method for simulating heat transfer behavior of a converter molten pool fire point area based on Fluent-UDF is as follows: Figure 1 As shown, including:

[0054] S1, establishing a geometric model of a steelmaking converter and meshing the geometric model of the steelmaking converter, importing the meshed geometric model of the steelmaking converter into Fluent software, performing simulation calculation using the Fluent software to obtain a quasi-steady-state flow field, and determining a designated area using the obtained quasi-steady-state flow field;

[0055] Specifically, in this embodiment, the implementation process of the above S1 is as follows:

[0056] S11, open the Fluent software, import the calculation grid of the steelmaking converter calculation domain, set the calculation parameters, and perform preliminary simulation calculation to obtain the quasi-steady-state flow field; wherein the calculation parameters include: physical property parameters, boundary conditions, solution method and calculation residual;

[0057] S12, based on the quasi-steady-state flow field, the radius and depth of the impact pit are obtained, and a circle with a radius R as the center and an intersection of the symmetry axis of the steelmaking converter geometric model and the metal liquid surface is determined; wherein R is a set value, and its value is not less than the radius and depth of the impact pit, so that the created circle can include the entire impact pit;

[0058] S13, using C_VOF macro to identify the calculation grid area with the volume fraction of the metal liquid phase between 0 and 1, and taking the intersection of the identified calculation grid area and the area surrounded by the circle with a radius of R as the designated area.

[0059] S2, establish the heat transfer model of the converter molten pool fire point area;

[0060] Specifically, in this embodiment, the implementation process of the above S2 is as follows:

[0061] S21, obtaining the content and activity of the element undergoing oxidation reaction in the hot spot area of ​​the converter molten pool, and calculating the ratio of oxygen consumed by the reaction of unit mole of O2 in the hot spot area with the corresponding element;

[0062] Among them, the elements are determined by thermodynamic analysis; when the temperature of the fire point zone is lower than 2386K, the elements that undergo oxidation reaction in the fire point zone of the converter molten pool include: C element, Si element, Mn element, P element and Fe element; when the temperature of the fire point zone is not lower than 2386K, the elements that undergo oxidation reaction in the fire point zone of the converter molten pool include: C element, Si element, Mn element and Fe element.

[0063] Furthermore, the calculation formula of the oxygen ratio is:

[0064]

[0065] in, ΔG is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and element i in the molten metal; i is the Gibbs free energy of the reaction of unit mole of O2 with element i in the molten metal; ΔG is the sum of the Gibbs free energies of the reaction of unit mole of O2 with each element in the molten metal.

[0066] When the temperature of the fire zone is lower than 2386K, When the temperature of the fire point area is not less than 2386K, Among them, ΔG CΔG is the Gibbs free energy per mole of O2 reacting with C in the metal liquid; Si ΔG is the Gibbs free energy per mole of O2 reacting with Si in the metal liquid; Mn ΔG is the Gibbs free energy per mole of O2 reacting with Mn in the metal liquid; Fe ΔG is the Gibbs free energy per mole of O2 reacting with Fe in the molten metal; P It is the Gibbs free energy of the reaction of unit mole of O2 with element P in the molten metal.

[0067] ΔG C The calculation formula is: In the formula, represents the standard Gibbs free energy of the reaction between element C and oxygen, and the Gibbs free energy of the other elements is ΔG i The calculation formula is: In the formula, represents the standard Gibbs free energy of the reaction of element i with oxygen, α i is the activity of element i, and its calculation formula is α i =f i ·W i , where f i is the activity coefficient of element i in the molten metal pool. When the temperature of the ignition zone is lower than 2386K, its calculation formula is: When the temperature of the fire point area is not lower than 2386K, the calculation formula is: In the formula, are the interaction coefficients of element i on carbon, silicon, manganese and phosphorus, respectively, and W i is the mass percentage of element i in the molten metal pool.

[0068] S22, based on the ratio of oxygen consumed by the reaction of unit mole of O2 in the fire point area with the corresponding element, the total heat release rate of the fire point area in the molten pool is calculated;

[0069] Among them, when the temperature of the fire point area is lower than 2386K, the calculation formula of the total heat release rate of the fire point area of ​​the molten pool is:

[0070]

[0071] When the temperature of the fire point area is not lower than 2386K, the calculation formula for the total heat release rate of the fire point area of ​​the molten pool is:

[0072]

[0073] Among them, q c is the total heat release rate of the molten pool fire point area, in J / s; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and element C in the molten metal; C It represents the enthalpy change when the C element participates in the reaction, in J / mol; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and Si element in the metal liquid; Si It represents the enthalpy change when Si element participates in the reaction, the unit is J / mol; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and Fe in the molten metal; Fe It represents the enthalpy change when Fe element participates in the reaction, the unit is J / mol; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and Mn element in the molten metal; Mn It represents the enthalpy change when Mn element participates in the reaction, in J / mol; ΔH is the ratio of oxygen consumed by the chemical reaction between unit mole of O2 and element P in the molten metal; P It represents the enthalpy change when element P participates in the reaction; It represents the molar flow rate of O2 in mol / s.

[0074] S3, develop user-defined function (UDF) program;

[0075] Among them, the UDF program can use the C_VOF macro pre-packaged in the Fluent software to track the volume fraction of the metal liquid phase in the specified area, so as to obtain the morphological changes of the fire point area, realize the tracking of the morphology of the fire point area in the specified area, and use the DEFINE_SOURCE macro function pre-packaged in the Fluent software to perform stable heat release in the specified area with the total heat release rate of the molten pool fire point area calculated by S2.

[0076] S4, importing the UDF program into Fluent for compilation, and adding the heat release rate into the Fluent energy source term equation for simulation calculation, so as to obtain the temperature field distribution of the molten pool fire point area of ​​the converter.

[0077] The heat release rate is added to the Fluent energy source term equation. Specifically, the calculated total heat release rate is added to the Fluent energy source term S through the DEFINE_SOURCE macro. T middle.

[0078] Below, this embodiment uses a practical application scenario to illustrate the application process of the method of the present invention.

[0079] This example takes a 100t four-hole top-blown oxygen lance converter in a steel plant as the research object, and constructs a two-dimensional axis rotation model of a single-hole oxygen lance converter when the oxygen lance position is 1.35m. The converter model and corresponding boundary conditions of this example are as follows: Figure 2 As shown. In this embodiment, the initial temperature of the molten pool of the converter model is 1773K, and the reaction temperature of each element in the molten pool fire point area is 2773K. Based on this, the specific implementation process of the method of the present invention is as follows:

[0080] (1) Open the Fluent software, import the computational grid of the steelmaking converter calculation domain, set the calculation parameters, including physical properties, boundary conditions, solution method, and calculation residuals.

[0081] Among them, the setting of physical property parameters: In this embodiment, only oxygen and molten steel are considered, and the basic physical property parameters of the two phases are shown in Table 1:

[0082] Table 1 Basic physical properties of the two phases

[0083] parameter Molten Steel oxygen <![CDATA[Density (Kg / m 3 )]]> 7100 Ideal gas Specific heat (J / Kg / K) 670 919.31 Thermal conductivity (W / m / K) 40 0.0246 Viscosity(Kg / m / s) 0.0065 Sutherland Temperature (K) 1773 300

[0084] Boundary condition setting: Set the oxygen lance inlet boundary type to pressure inlet, inlet pressure to 9000000Pa, inlet temperature to 300K, converter outlet boundary type to pressure outlet, outlet pressure to 101325Pa, outlet temperature to 1773K, converter wall to adiabatic wall. Solution selection: Pressure-velocity coupling uses PISO format, pressure interpolation is performed using the staggered pressure (PREssure Staggering Option, PRESTO!) method, free interface capture uses CISCAM format, and the transport equation is discretized using the second-order upwind format. When the energy residual is less than 10 -6 , the residuals of other variables are less than 10 -3 When , the simulation results are considered to have reached the convergence standard. Then, a preliminary simulation calculation is performed to obtain the quasi-steady flow field, such as Figure 3 As shown, the pit depth determined according to the quasi-steady flow field is 0.24 m and the diameter is 0.4 m.

[0085] (2) Determine the type and activity of elements undergoing oxidation reactions in the molten pool fire zone through thermodynamic analysis

[0086] There are five elements, C, Si, Mn, P, and Fe, that undergo direct oxidation reactions in the hot zone of the converter molten pool. However, when the reaction temperature is 2773K, there are four elements, C, Si, Mn, and Fe, that can undergo direct oxidation reactions in the hot zone. In this embodiment, the total heat release rate in the hot zone is calculated based on the element content of a furnace of steel at a certain smelting moment. At this time, the mass percentages of the three elements C, Si, and Mn in the molten steel are 3.32%, 0.0104%, and 0.1%, respectively. Table 2 below lists the chemical reactions of each element.

[0087] Table 2 Chemical reactions

[0088] Serial number Chemical reactions 1 <![CDATA[[C]+1 / 2O2=CO]]> 2 <![CDATA[[Si]+O2=SiO2]]> 3 <![CDATA[[Mn]+1 / 2O2=MnO]]> 4 <![CDATA[[Fe]+1 / 2O2=FeO]]>

[0089] Based on this, the activities of four elements C, Si, Mn and Fe in molten steel are solved.

[0090] The activity coefficient of a single element is shown in formula (1):

[0091]

[0092] In formula (1) are the interaction coefficients of element i on carbon, silicon and manganese respectively, and their values ​​are shown in Table 3. i is the mass percentage of element i in the molten steel. After solving the activity coefficients of all elements, the activity of each element can be obtained through formula (2).

[0093] α [i] =f i ·W i (2)

[0094] Table 3 Interaction coefficients of each element

[0095]

[0096] (3) Find the ratio of oxygen that reacts with each element per unit mole of O2

[0097] Gibbs free energy ΔG of a single element reaction i It can be obtained from equations (3) and (4): The gas partial pressure in the equations is one standard atmospheric pressure:

[0098]

[0099]

[0100] The expression of the sum of Gibbs freedom ΔG of a single element reaction is shown in formula (5):

[0101]

[0102] The ratio of oxygen that reacts with each element per unit mole of O2 The calculation formula is shown in formula (6):

[0103]

[0104] (4) Determine the total heat release rate of the above four elements in the fire point area.

[0105] The total heat release rate q in the molten pool fire zone at a certain moment c The calculation formula is shown in formula (7):

[0106]

[0107] ΔH in formula (7) i It represents the enthalpy change when element i participates in the reaction, in J / mol. It represents the molar flow rate of O2 in mol / s.

[0108] In this embodiment, the enthalpy changes of the four elements C, Si, Mn, and Fe are -125339.2, -908322.73, -385397.42, and -270887.89 J / mol, respectively, and the molar flow rate of O2 is 275 mol / s. The total heat release rate of the fire point area of ​​the four-hole oxygen gun calculated based on formula (7) is 69918223.7 J / s. It is believed that the total heat release rate of the fire point area of ​​the single-hole oxygen gun is one-fourth of that of the four-hole oxygen gun. Therefore, the total heat release rate of the fire point area of ​​the single-hole oxygen gun is q c It is 17479555.9 J / s.

[0109] (5) Based on the calculated total heat release rate, a UDF program is developed. The UDF program can track the morphology of the fire point area within a specified range, and use the DEFINE_SOURCE macro function pre-packaged in the Fluent software to stably release heat in the tracking area. The schematic diagram of the morphology of the fire point area specified in this embodiment is shown in FIG. Figure 4 shown.

[0110] (6) According to the radius and depth of the impact pit, the radius of the circle is determined to be 0.4 m, and the center of the circle is the intersection of the symmetry axis of the converter model and the steel liquid surface.

[0111] (7) The C_VOF macro is used to identify the computational grid area with a liquid steel volume fraction between 0 and 1, and the intersection of the identified area and the circular domain with a radius of 0.4 m is used as the tracking range specified in the UDF.

[0112] (8) Import the UDF program into Fluent for compilation, and use the DEFINE_SOURCE macro to convert the heat release rate q c Add to Fluent energy source term S TBy continuing the simulation calculation, the temperature field distribution, flow field distribution and flow field velocity vector direction diagram of the converter molten pool can be obtained, such as Figure 5 and Figure 6 As shown. Figure 5 and Figure 6 It can be seen that the temperature in the converter molten pool with lower flow velocity is also lower. Therefore, it can be seen that the temperature field distribution of the converter molten pool is closely related to the flow field distribution. Figure 5 It can be seen that the direction of the converter molten pool temperature gradient in this embodiment is opposite to the direction of the flow field velocity vector, and the flow of molten steel will hinder the diffusion of the high temperature zone, causing the molten pool to heat up slowly. Through the method of the present invention, the influence of parameters such as oxygen lance position and oxygen lance operating pressure on the distribution of the molten pool temperature field and the relative direction of the molten pool temperature gradient and the flow field velocity vector can also be studied, thereby summarizing the influence of different operating parameters on the heat transfer behavior of the molten pool, and better controlling the molten pool temperature in the converter smelting process in actual production.

[0113] In summary, this embodiment provides a method for simulating the heat transfer behavior of the fire point zone of the converter melt pool based on Fluent-UDF. The method for simulating the heat transfer behavior of the fire point zone of the converter melt pool utilizes Fluent and UDF custom functions to establish a heat transfer model of the fire point zone of the converter melt pool, and numerically solves the flow field and temperature field inside the converter melt pool. It can analyze the heat transfer behavior between the melt pool fire point zone and the melt pool metal liquid under different working conditions and the temperature field distribution of the melt pool metal liquid, thereby providing a reference for the melt pool temperature control in the converter smelting process.

[0114] Second embodiment

[0115] This embodiment provides an electronic device, such as Figure 7 As shown, the electronic device includes: a processor and a memory; wherein the processor and the memory can be connected via a communication bus; the memory stores at least one instruction, and the instruction is loaded and executed by the processor to implement the method of the first embodiment. In addition, the electronic device may also include a transceiver, the processor and the transceiver can be connected via a communication bus, and the transceiver is used to communicate with other devices.

[0116] Next, combine Figure 7 The following is a detailed introduction to the various components of the electronic device:

[0117] Among them, the processor is the control center of the electronic device, and the electronic device may include multiple processors, each of which may be a single-core processor (single-CPU) or a multi-core processor (multi-CPU). The processor here may be a processor or a general term for multiple processing elements. For example, the processor is one or more central processing units (CPUs), or other general-purpose processors, application specific integrated circuits (ASICs), or one or more integrated circuits configured to implement an embodiment of the present invention, such as one or more microprocessors (digital signal processors, DSPs), or one or more field programmable gate arrays (field programmable gate arrays, FPGAs), or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc. The processor may execute various functions of the electronic device by running or executing software programs stored in the memory and calling data stored in the memory.

[0118] In a specific implementation, as an embodiment, the processor may include one or more CPUs, such as Figure 7 The CPU0 and CPU1 shown in the figure are, of course, only exemplary.

[0119] The memory is used to store the software program for executing the solution of the present invention, and the execution is controlled by the processor. The specific implementation method can refer to the above method embodiment and will not be repeated here.

[0120] Optionally, the memory may be a read-only memory (ROM) or other types of static storage devices that can store static information and instructions, a random access memory (RAM) or other types of dynamic storage devices that can store information and instructions, or an electrically erasable programmable read-only memory (EEPROM), a compact disc read-only memory (CD-ROM) or other optical disc storage, optical disc storage (including compressed optical discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto. The memory may be integrated with the processor or exist independently and accessed through the interface circuit ( Figure 7 (not shown) is coupled to the processor, which is not specifically limited in this embodiment of the present invention.

[0121] The transceiver may include a receiver and a transmitter ( Figure 7 The receiver is used to implement the receiving function, and the transmitter is used to implement the sending function. The transceiver can be integrated with the processor or exist independently and communicate with the electronic device through the interface circuit ( Figure 7 (not shown) is coupled to the processor, which is not specifically limited in this embodiment of the present invention.

[0122] In addition, it should be noted that Figure 7 The structure of the electronic device shown in the figure does not constitute a limitation on the device, and the actual device may include more or fewer components than shown in the figure, or combine certain components, or arrange the components differently. In addition, the technical effects achieved by the electronic device when executing the method of the first embodiment above can refer to the technical effects described in the first embodiment above, so they are not repeated here.

[0123] Third embodiment

[0124] This embodiment provides a computer-readable storage medium, which stores at least one instruction, and the instruction is loaded and executed by a processor to implement the method of the first embodiment. The computer-readable storage medium may be a ROM, a random access memory, a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, etc. The instructions stored therein may be loaded by a processor in a terminal to execute the method.

[0125] In addition, it should be noted that the present invention can be provided as a method, an apparatus or a computer program product. Therefore, the embodiment of the present invention can be in the form of a full or partial hardware embodiment, a full or partial software embodiment or an embodiment combining software and hardware. Moreover, when implemented using software, the embodiment of the present invention can be in the form of a computer program product implemented on one or more computer-usable storage media containing computer-usable program codes. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, the process or function described in the embodiment of the present invention is generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium, for example, the computer instructions can be transmitted from one website site, computer, server or data center to another website site, computer, server or data center by wired (e.g., infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device such as a server or data center containing one or more available media sets. The available medium may be a magnetic medium (eg, a floppy disk, a hard disk, a magnetic tape), an optical medium (eg, a DVD), or a semiconductor medium. The semiconductor medium may be a solid state hard disk.

[0126] The embodiments of the present invention are described with reference to the flowcharts and / or block diagrams of the methods, terminal devices (systems), and computer program products according to the embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the 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, an embedded processor, or other programmable data processing terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal 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.

[0127] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device, so that a series of operation steps are executed on the computer or other programmable terminal device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable terminal device provide for implementing the process in the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0128] It should also be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. The terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that the process, method, article or terminal device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal device. In the absence of more restrictions, the elements defined by the sentence "including one..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements. In addition, the term "and / or" is only an association relationship describing the associated objects, indicating that there can be three relationships, for example, A and / or B, which can represent: A exists alone, A and B exist at the same time, and B exists alone, wherein A and B can be singular or plural. In addition, the character " / " in this article generally indicates that the objects before and after are in an "or" relationship, but it may also indicate an "and / or" relationship. Please refer to the context for specific understanding. "At least one" means one or more, and "plurality" means two or more. "At least one of the following" or similar expressions refers to any combination of these items, including any combination of single or plural items. For example, at least one of a, b or c can be represented by: a, b, c, ab, ac, bc or abc, where a, b, c can be single or plural.

[0129] In addition, it can be understood that in various embodiments of the present invention, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0130] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0131] In several embodiments provided by the present invention, it should be understood that the disclosed equipment, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic, for example, the division of functional modules / units is only a logical function division, and there may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point, the coupling or direct coupling or communication connection between each other shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms. The unit described as a separate component may or may not be physically separated, and the component displayed as a unit may or may not be a physical unit, that is, it may be located in one place, or it may be distributed on multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit.

[0132] If the method is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0133] Finally, it should be noted that the above is a preferred embodiment of the present invention. It should be pointed out that although the preferred embodiment of the present invention has been described, for those skilled in the art, once the basic creative concept of the present invention is known, several improvements and modifications can be made without departing from the principles of the present invention, and these improvements and modifications should also be regarded as the protection scope of the present invention. Therefore, the attached claims are intended to be interpreted as including the preferred embodiment and all changes and modifications that fall within the scope of the embodiments of the present invention.

Claims

1. A method for simulating heat transfer behavior in the hot zone of a converter molten pool based on Fluent-UDF. It is characterized in that The Fluent-UDF-based method for simulating heat transfer behavior in the hot zone of the converter molten pool includes: Establishing a geometric model of a steelmaking converter and meshing the geometric model of the steelmaking converter, importing the meshed geometric model of the steelmaking converter into Fluent software, performing simulation calculations using the Fluent software to obtain a quasi-steady-state flow field, and determining a designated area using the obtained quasi-steady-state flow field; Developing a UDF (User defined functions) program; the UDF program is used to track the morphology of the fire point area in the designated area and to perform stable heat release in the designated area; The UDF program is imported into Fluent for compilation, and the heat release rate is added to the Fluent energy source term equation for simulation calculation, so as to obtain the temperature field distribution of the molten pool fire point area of ​​the converter; A geometric model of a steelmaking converter is established and meshed, the meshed geometric model of the steelmaking converter is imported into Fluent software, and simulation calculation is performed using Fluent software to obtain a quasi-steady-state flow field, and a designated area is determined by the obtained quasi-steady-state flow field, including: Open the Fluent software, import the computational grid of the steelmaking converter computational domain, set the computational parameters, and perform preliminary simulation calculations to obtain the quasi-steady-state flow field; wherein the computational parameters include: physical property parameters, boundary conditions, solution methods, and computational residuals; Based on the quasi-steady flow field, the radius and depth of the impact pit are obtained, and a circle with a radius of R and an intersection of the symmetry axis of the steelmaking converter geometric model and the metal liquid surface as the center is determined; wherein R is a set value, and its value is not less than the radius and depth of the impact pit, so that the created circle can include the entire impact pit; The C_VOF macro is used to identify the computational grid region with a volume fraction of the metal liquid phase between 0 and 1, and the intersection of the identified computational grid region and the region enclosed by a circle with a radius of R is taken as the designated region.

2. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 1, It is characterized in that Develop UDF programs, including: Obtain the element content and activity of the oxidation reaction in the hot zone of the converter molten pool, and calculate the unit mole of O 2 The proportion of oxygen consumed in reacting with the corresponding element; Based on the obtained flash point area unit mole O 2 The proportion of oxygen consumed in the reaction with the corresponding element is used to calculate the total heat release rate in the fire point area of ​​the molten pool; A UDF program was developed based on the obtained total heat release rate.

3. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 2, It is characterized in that When the temperature of the fire point area is lower than 2386K, the elements that undergo oxidation reaction in the fire point area of ​​the converter molten pool include: C, Si, Mn, P and Fe. When the temperature of the fire point zone is not lower than 2386K, the elements undergoing oxidation reaction in the fire point zone of the converter molten pool include: C element, Si element, Mn element and Fe element.

4. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 3, It is characterized in that The calculation formula of the oxygen ratio is: in, The unit is mole of O 2 The ratio of oxygen consumed in chemical reaction with element i in the molten metal; ΔG i The unit is mole of O 2 Gibbs free energy of reaction with element i in the metal liquid; ΔG is the unit mole of O 2 The sum of the Gibbs free energies of the reactions with each element in the molten metal; When the temperature of the fire zone is lower than 2386K, When the temperature of the fire point area is not less than 2386K, Among them, ΔG C The unit is mole of O 2 Gibbs free energy of reaction with C element in molten metal; ΔG Si The unit is mole of O 2 Gibbs free energy when reacting with Si in the metal liquid; ΔG Mn The unit is mole of O 2 Gibbs free energy of reaction with Mn in the molten metal; ΔG Fe The unit is mole of O 2 Gibbs free energy of reaction with Fe in molten metal; ΔG P The unit is mole of O 2 Gibbs free energy of reaction with P element in molten metal.

5. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 4, It is characterized in that When the temperature of the fire point area is lower than 2386K, the calculation formula for the total heat release rate of the fire point area of ​​the molten pool is: When the temperature of the fire point area is not lower than 2386K, the calculation formula for the total heat release rate of the fire point area of ​​the molten pool is: Among them, q c is the total heat release rate in the fire point area of ​​the molten pool; The unit is mole of O 2 The ratio of oxygen consumed by chemical reaction with C element in the molten metal; ΔH C It represents the enthalpy change when the C element participates in the reaction; The unit is mole of O 2 The ratio of oxygen consumed by chemical reaction with Si element in the metal liquid; ΔH Si It represents the enthalpy change when Si element participates in the reaction; The unit is mole of O 2 The ratio of oxygen consumed by chemical reaction with Fe in the molten metal; ΔH Fe It represents the enthalpy change when Fe element participates in the reaction; The unit is mole of O 2 The ratio of oxygen consumed by chemical reaction with Mn in the molten metal; ΔH Mn It represents the enthalpy change when the Mn element participates in the reaction; The unit is mole of O 2 The ratio of oxygen consumed by chemical reaction with P element in the molten metal; ΔH P It represents the enthalpy change when element P participates in the reaction; Indicates O 2 The molar flow rate.

6. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 1, It is characterized in that Tracking the morphology of the fire point area in the designated area, specifically: using the C_VOF macro pre-packaged in the Fluent software to track the volume fraction of the metal liquid phase in the designated area to obtain the morphology change of the fire point area.

7. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 5, It is characterized in that The heat release is carried out stably in the designated area, specifically, the DEFINE_SOURCE macro function prepackaged in the Fluent software is used to carry out the heat release stably in the designated area at the calculated total heat release rate of the molten pool fire point area.

8. The method for simulating heat transfer behavior of the hot spot area of ​​the converter molten pool based on Fluent-UDF as claimed in claim 5, It is characterized in that Add the heat release rate to the Fluent energy source term equation. Specifically, add the calculated total heat release rate to the Fluent energy source term S through the DEFINE_SOURCE macro. T middle.

Images