CN117972999A - Converter molten pool fire point area heat transfer behavior simulation method based on Fluent-UDF - Google Patents

Abstract

The present invention discloses a method for simulating heat transfer behavior of a converter molten pool fire point zone 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 zone 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 zone. The present invention can provide a reference for the molten pool temperature control in the converter smelting process.

Description

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

Technical Field

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 technique

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 phenomena 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.

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. As one of the main heat sources in the converter smelting process, the fire point area of the converter molten pool 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.

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 derivation 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

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.

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

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:

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 fire point area of the converter molten pool.

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 the 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 R and an intersection of the symmetry axis of the steelmaking converter geometric model and the metal liquid surface as the center is determined;

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.

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

Furthermore, a UDF program is developed, including:

Obtain the content and activity of elements undergoing oxidation reactions in the hot zone of the converter molten pool, and calculate the ratio of oxygen consumed by the reaction of unit mole of _O2 in the hot zone with the corresponding element;

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;

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

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;

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.

Furthermore, the calculation formula of the oxygen ratio is:

in, is the ratio of oxygen consumed by the chemical reaction between unit mole of O ₂ and element i in the molten metal; ΔG _i is the Gibbs free energy of the reaction between unit mole of O ₂ and element i in the molten metal; ΔG is the sum of the Gibbs free energies of the reaction between unit mole of O ₂ and each element in the molten metal;

When the temperature of the fire zone is lower than 2386K, When the temperature of the ignition zone is not less than 2386K,/> Among them, ΔG _C is the Gibbs free energy of the reaction between unit mole of O ₂ and C element in the metal liquid; ΔG _Si is the Gibbs free energy of the reaction between unit mole of O ₂ and Si element in the metal liquid; ΔG _Mn is the Gibbs free energy of the reaction between unit mole of O ₂ and Mn element in the metal liquid; ΔG _Fe is the Gibbs free energy of the reaction between unit mole of O ₂ and Fe element in the metal liquid; ΔG _P is the Gibbs free energy of the reaction between unit mole of O ₂ and P element in the metal liquid.

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:

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:

Where, q _c is the total heat release rate in the molten pool fire point area; It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and element C in the metal liquid; ΔH _C represents the enthalpy change when element C participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Si element in the metal liquid; ΔH _Si represents the enthalpy change when Si element participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Fe in the molten metal; _ΔHFe represents the enthalpy change when Fe participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Mn in the molten metal; ΔH _Mn represents the enthalpy change when Mn participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and element P in the molten metal; _ΔHP represents the enthalpy change when element P participates in the reaction; /> represents the molar flow rate of _O2 .

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.

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.

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 _ST through the DEFINE_SOURCE macro.

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.

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

The present invention uses Fluent and UDF custom functions to establish a heat transfer model of the fire point area of the converter molten pool, 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.

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

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 problem, and then adds UDF to the energy source term equation to simulate the heat release behavior of the fire point area of the molten pool, thereby avoiding the construction of a complex chemical reaction model in Fluent and greatly reducing the amount of calculation.

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

In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. 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.

1 is a schematic diagram of the execution flow of a method for simulating heat transfer behavior of a hot spot area in a converter molten pool based on Fluent-UDF according to an embodiment of the present invention;

FIG2 is a schematic diagram of a two-dimensional converter model and boundary conditions provided by an embodiment of the present invention;

FIG3 is a schematic diagram of a quasi-steady-state flow field and an impact pit morphology provided by an embodiment of the present invention;

FIG4 is a schematic diagram of a hot spot area of a converter molten pool provided by an embodiment of the present invention;

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

6 is a cloud diagram of converter molten pool velocity distribution and flow field velocity vector direction provided by an embodiment of the present invention;

FIG. 7 is a system block diagram of an electronic device provided by an embodiment of the present invention.

Detailed ways

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.

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.

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.

In addition, in the embodiments of the present invention, sometimes a subscript (such as W ₁ ) 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.

First embodiment

This embodiment provides a method for simulating heat transfer behavior of a converter molten pool fire point area based on Fluent-UDF, which 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 shown in Figure 1, including:

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;

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

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;

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;

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.

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

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

S21, obtaining the content and activity of the elements 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;

Wherein, 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.

Furthermore, the calculation formula of the oxygen ratio is:

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

When the temperature of the fire zone is lower than 2386K, When the temperature of the ignition zone is not less than 2386K,/> Among them, ΔG _C is the Gibbs free energy of the reaction between unit mole of O ₂ and C element in the metal liquid; ΔG _Si is the Gibbs free energy of the reaction between unit mole of O ₂ and Si element in the metal liquid; ΔG _Mn is the Gibbs free energy of the reaction between unit mole of O ₂ and Mn element in the metal liquid; ΔG _Fe is the Gibbs free energy of the reaction between unit mole of O ₂ and Fe element in the metal liquid; ΔG _P is the Gibbs free energy of the reaction between unit mole of O ₂ and P element in the metal liquid.

The calculation formula of ΔG _C is: In the formula, /> represents the standard Gibbs free energy of the reaction between element C and oxygen. The Gibbs free energy ΔG _i of other elements is calculated as follows: In the formula, /> It represents the standard Gibbs free energy of the reaction between element i and oxygen, α _i is the activity of element i, and its calculation formula is α _i = _fi · _Wi , where _fi is the activity coefficient of element i in the molten metal pool. When the temperature of the ignition point zone is lower than 2386K, its calculation formula is:/> When the temperature of the fire point area is not less than 2386K, the calculation formula is:/> In the formula, are the interaction coefficients of element i on carbon, silicon, manganese and phosphorus respectively, and Wi _is the mass percentage of element i in the molten metal pool.

S22, based on the obtained ratio of oxygen consumed by the reaction of unit mole of _O2 in the fire point area with the corresponding element, calculate the total heat release rate of the fire point area in the molten pool;

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:

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:

Where, q _c is the total heat release rate of the molten pool fire point area, in J/s; It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and element C in the molten metal; ΔH _C represents the enthalpy change when element C participates in the reaction, and its unit is J/mol; It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Si element in the metal liquid; ΔH _Si represents the enthalpy change when Si element participates in the reaction, and the unit is J/mol; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Fe in the molten metal; _ΔHFe represents the enthalpy change when Fe participates in the reaction, and its unit is J/mol;/> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Mn in the molten metal; ΔH _Mn represents the enthalpy change when Mn participates in the reaction, and its unit is J/mol; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and element P in the molten metal; _ΔHP represents the enthalpy change when element P participates in the reaction; /> It represents the molar flow rate of _O2 in mol/s.

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

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.

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.

Among them, 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 _ST through the DEFINE_SOURCE macro.

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

This embodiment 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 of this embodiment and the corresponding boundary conditions are shown in Figure 2. The initial temperature of the molten pool of the converter model in this embodiment 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:

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

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:

Table 1 Basic physical properties of the two phases

parameter Molten Steel oxygen Density (Kg/m ³ ) 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

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 method 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 second-order upwind format is used to discretize the transport equation. When the energy residual is less than 10 ^-6 and the residuals of other variables are less than 10 ^-3 , the simulation results are considered to have reached the convergence standard. Subsequently, a preliminary simulation calculation was performed to obtain a quasi-steady-state flow field, as shown in Figure 3. The pit depth determined by the quasi-steady-state flow field is 0.24m and the diameter is 0.4m.

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

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 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.

Table 2 Chemical reactions

Serial number chemical reaction 1 [C] + 1/2O ₂ = CO 2 [Si] + O ₂ = SiO ₂ 3 [Mn] + 1/2O ₂ = MnO 4 [Fe] + 1/2O ₂ = FeO

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

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

In formula (1) are the interaction coefficients of element i on carbon, silicon and manganese respectively, and their values are shown in Table 3. Wi _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 by formula (2).

α _[i] = _fi · _Wi (2)

Table 3 Interaction coefficients of each element

(3) Find the ratio of oxygen that reacts with each element per unit mole of O ₂

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

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

The ratio of oxygen that reacts with each element per mole of O ₂ The calculation formula is shown in formula (6):

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

The calculation formula of the total heat release rate _qc in the molten pool fire point area at a certain moment is shown in formula (7):

ΔH _i in formula (7) represents the enthalpy change when element i participates in the reaction, and its unit is J/mol. It represents the molar flow rate of _O2 in mol/s.

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 considered 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 _qc of the fire point area of the single-hole oxygen gun is 17479555.9 J/s.

(5) Based on the calculated total heat release rate, a UDF program is developed, which 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 Figure 4.

(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.

(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.

(8) Import the UDF program into Fluent for compilation, add the heat release rate _qc to the Fluent energy source item _ST through the DEFINE_SOURCE macro, and continue the simulation calculation to obtain the temperature field distribution, flow field distribution and flow field velocity vector direction diagram of the converter molten pool, as shown in Figures 5 and 6. Comparing Figures 5 and 6, it can be seen that the temperature of the area with lower flow rate in the converter molten pool is also relatively low. Therefore, it can be seen that the temperature field distribution of the converter molten pool is closely related to the flow field distribution; analyzing Figure 5, it can be seen that the direction of the temperature gradient of the converter molten pool under 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 area, 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 molten pool temperature field distribution and the relative direction of the molten pool temperature gradient and flow field velocity vector can also be studied, so as to summarize the influence of different operating parameters on the heat transfer behavior of the molten pool, and the molten pool temperature in the converter smelting process can be better controlled in actual production.

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.

Second embodiment

This embodiment provides an electronic device, as shown in FIG7 , 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, the instruction is loaded and executed by the processor to implement the method of the first embodiment above. In addition, the electronic device may also include a transceiver, the processor and the transceiver may be connected via a communication bus, and the transceiver is used to communicate with other devices.

Next, the components of the electronic device are specifically described in conjunction with FIG. 7 :

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.

In a specific implementation, as an embodiment, the processor may include one or more CPUs, such as CPU0 and CPU1 shown in FIG. 7 . Of course, this is only an exemplary description.

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, which will not be repeated here.

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 disc, laser disc, optical disc, digital versatile disc, Blu-ray disc, 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 be coupled to the processor through the interface circuit (not shown in FIG. 7 ) of the electronic device, and the embodiment of the present invention does not specifically limit this.

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

In addition, it should be noted that the structure of the electronic device shown in FIG7 does not constitute a limitation on the device, and the actual device may include more or fewer components than shown, 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 can refer to the technical effects described in the first embodiment, so they will not be described here.

Third embodiment

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.

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 device. 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, computer, server or data center to another website, 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, data center, etc. that contains one or more available media sets. The available medium can be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium. The semiconductor medium may be a solid state drive.

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 box in the flowchart and/or block diagram, as well as the combination of the processes and/or boxes 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 a device for implementing the functions specified in one process or multiple processes in the flowchart and/or one box or multiple boxes in the block diagram.

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 product including an instruction device that implements the functions specified in one or more processes of the flowchart and/or one or more blocks of the block diagram. These computer program instructions may also be loaded onto a computer or other programmable data processing terminal device, so that a series of operation steps are performed 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 steps for implementing the functions specified in one or more processes of the flowchart and/or one or more blocks of the block diagram.

It should also be noted that, in this article, relational terms such as first and second 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 a 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 further 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 a description of the association relationship of the associated objects, indicating that there can be three relationships, for example, A and/or B can represent: A exists alone, A and B exist at the same time, and B exists alone, where A and B can be singular or plural. In addition, the character "/" in this article generally indicates that the associated objects before and after are in an "or" relationship, but it may also represent an "and/or" relationship, which can be understood by referring to the context before and after. “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 mean: a, b, c, a-b, a-c, b-c, or a-b-c, where a, b, and c can be single or plural.

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.

Those of ordinary skill 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.

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 mutual coupling or direct coupling or communication connection 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.

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.

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, once the basic creative concept of the present invention is known, a number of improvements and modifications can be made by those skilled in the art 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 (10)

1. A method for simulating heat transfer behavior in a converter molten pool fire point area based on Fluent-UDF, characterized in that the method for simulating heat transfer behavior in a converter molten pool fire point area based on Fluent-UDF 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 fire point area of the converter molten pool. 2. The method for simulating heat transfer behavior of a converter molten pool fire point zone based on Fluent-UDF according to claim 1 is characterized in that a geometric model of a steelmaking converter is established and the geometric model of the steelmaking converter is meshed, the meshed steelmaking converter geometric model 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 R and an intersection of the symmetry axis of the steelmaking converter geometric model and the metal liquid surface as the center is determined; 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. 3. The method for simulating heat transfer behavior of the fire point zone of the converter molten pool based on Fluent-UDF as described in claim 2 is characterized in that R is a set value, and its value is not less than the radius and depth of the impact pit. 4. The method for simulating heat transfer behavior of a converter molten pool fire point zone based on Fluent-UDF according to claim 1, characterized in that developing a UDF program comprises: Obtain the content and activity of elements undergoing oxidation reactions in the hot zone of the converter molten pool, and calculate the ratio of oxygen consumed by the reaction of unit mole of _O2 in the hot zone with the corresponding element; 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; A UDF program was developed based on the obtained total heat release rate. 5. The method for simulating heat transfer behavior of the hot spot zone of the converter molten pool based on Fluent-UDF as claimed in claim 4, characterized in that when the temperature of the hot spot zone is lower than 2386K, the elements undergoing oxidation reaction in the hot spot 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 undergoing oxidation reaction in the fire point zone of the converter molten pool include: C, Si, Mn and Fe. 6. The method for simulating heat transfer behavior of the hot spot zone of the converter molten pool based on Fluent-UDF according to claim 5, characterized in that the calculation formula of the oxygen ratio is: in, is the ratio of oxygen consumed by the chemical reaction between unit mole of O ₂ and element i in the molten metal; ΔG _i is the Gibbs free energy of the reaction between unit mole of O ₂ and element i in the molten metal; ΔG is the sum of the Gibbs free energies of the reaction between unit mole of O ₂ and each element in the molten metal; When the temperature of the fire zone is lower than 2386K, When the temperature of the ignition zone is not less than 2386K,/> Among them, ΔG _C is the Gibbs free energy of the reaction between unit mole of O ₂ and C element in the metal liquid; ΔG _Si is the Gibbs free energy of the reaction between unit mole of O ₂ and Si element in the metal liquid; ΔG _Mn is the Gibbs free energy of the reaction between unit mole of O ₂ and Mn element in the metal liquid; ΔG _Fe is the Gibbs free energy of the reaction between unit mole of O ₂ and Fe element in the metal liquid; ΔG _P is the Gibbs free energy of the reaction between unit mole of O ₂ and P element in the metal liquid. 7. The method for simulating heat transfer behavior of the hot spot zone of a converter molten pool based on Fluent-UDF as claimed in claim 6, characterized in that when the temperature of the hot spot zone is lower than 2386K, the calculation formula for the total heat release rate of the hot spot zone 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: Where, q _c is the total heat release rate in the molten pool fire point area; It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and element C in the metal liquid; ΔH _C represents the enthalpy change when element C participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Si element in the metal liquid; ΔH _Si represents the enthalpy change when Si element participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Fe in the molten metal; _ΔHFe represents the enthalpy change when Fe participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and Mn in the molten metal; ΔH _Mn represents the enthalpy change when Mn participates in the reaction; /> It is the ratio of oxygen consumed by the chemical reaction between unit mole of _O2 and element P in the molten metal; _ΔHP represents the enthalpy change when element P participates in the reaction; _SO2 represents the molar flow rate of _O2 . 8. The method for simulating heat transfer behavior of the fire point zone of the converter molten pool based on Fluent-UDF as described in claim 1 is characterized in that the morphology of the fire point zone in the designated area is tracked, specifically: the C_VOF macro pre-packaged in the Fluent software is used to track the volume fraction of the metal liquid phase in the designated area to obtain the morphological changes of the fire point zone. 9. The method for simulating heat transfer behavior of the converter molten pool fire point zone based on Fluent-UDF as described in claim 7 is characterized in that 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 zone. 10. The method for simulating heat transfer behavior of the hot spot zone of the converter molten pool based on Fluent-UDF as described in claim 7 is characterized in that 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 _ST through the DEFINE_SOURCE macro.