Simulation method and device for cement clinker firing process of rotary kiln and electronic equipment

Abstract

The invention provides a simulation method and device for a cement clinker firing process of a rotary kiln and electronic equipment. The simulation method comprises the following steps: setting initial parameters of material particles in a simulation area in the rotary kiln; simulating a mixed flow process of the material particles to reach a steady state, and obtaining an initial state of the material particles; and simulating the reaction and fusion process of the material particles to obtain the state of the material particles at any moment. The simulation method provided by the invention is based on a discrete particle simulation method, the sintering process in the rotary kiln is simulated from the particle size, simulation can be carried out for different reaction sections according to the characteristics of the rotary kiln, the particle size reaction process of each section is finely calculated, and simulation of local details such as radial distribution in the kiln is realized.

Description

Technical Field

[0001] The present invention belongs to the technical field of industrial simulation, and in particular relates to a simulation method, device and electronic equipment for a rotary kiln cement clinker burning process. Background Art

[0002] The cement rotary kiln is a multi-physics coupled system, and its clinker production involves a variety of complex physical and chemical processes. Under high temperature, raw materials undergo complex reactions such as pyrolysis, redox reactions, and hardening. Complex multiphase flows occur between the gas, particulate, and liquid phases, and significant temperature gradients within the kiln significantly influence the properties of the clinker.

[0003] Therefore, modeling a cement rotary kiln involves complex kiln structures, different air intake methods, combustion of particulate materials (such as pulverized coal particles), and the melting, reaction, and bonding growth of cement raw materials (such as limestone and clay). The movement of gases, particles, and the kiln, as well as the movement, mass transfer, and reactions of gases and particles, vary by multiple orders of magnitude in time and space, making accurate modeling of the kiln difficult. However, detailed simulation of the clinker burning process in a cement rotary kiln is crucial for optimizing production processes and improving energy efficiency.

[0004] Based on this, it is necessary to provide a particle-scale-based simulation method for the cement clinker burning process in a rotary kiln, so as to study in detail the chemical reactions and particle fusion processes at local locations in the rotary kiln, and provide support for studying the influence of factors such as the structural parameters of the rotary kiln, air intake components and methods, traditional fuels and alternative fuels, traditional raw materials and alternative raw materials on the energy consumption, emissions, product performance, etc. of the cement rotary kiln, as well as exploring carbon reduction optimization solutions for the rotary kiln. Summary of the Invention

[0005] The purpose of the present invention is to provide a simulation method, device and electronic equipment for the cement clinker burning process in a rotary kiln, to simulate the burning process of the rotary kiln based on the particle size, and to achieve a detailed study of the burning process at a local position in the rotary kiln.

[0006] In order to achieve the purpose of the invention, the present invention adopts the following technical solutions:

[0007] In a first aspect, the present invention provides a method for simulating a cement clinker burning process in a rotary kiln, characterized in that the simulation method comprises:

[0008] Set the initial parameters of the material particles in the simulation area of ​​the rotary kiln;

[0009] Simulate the mixing flow process of material particles to reach a steady state and obtain the initialization state of the material particles;

[0010] The reaction and fusion process of material particles are simulated to obtain the state of material particles at any moment.

[0011] The simulation method provided by the present invention is based on the discrete particle simulation method, and simulates the firing process in the rotary kiln from the particle scale. According to the characteristics of the rotary kiln, simulation can be performed for different reaction sections, and the particle scale reaction process of each section can be accurately calculated to realize the simulation of local details such as the radial distribution in the kiln.

[0012] Preferably, the initial parameters of the material particles include composition, particle size, position and temperature.

[0013] Specifically, among the initial parameters of the material particles, the composition and particle size are set according to the temperature environment of the simulation area, the position is generated by uniform random numbers, and the temperature is set to the average temperature of the simulation area.

[0014] In the present invention, the rotary kiln can be decomposed into multiple temperature sections according to the axial temperature distribution characteristics of the rotary kiln, including a preheating section, a calcining section, and a cooling section. A simulation area is selected, and the change process of the particle size in each section is precisely calculated to obtain the local details of the radial distribution of each section.

[0015] The temperature environment of the simulation area of ​​the present invention can be predicted based on the macroscopic model of the rotary kiln, and can also be measured based on the actual working conditions of the rotary kiln.

[0016] Preferably, the simulation of the mixed flow process of the material particles to reach a steady state and obtain an initialization state of the material particles includes:

[0017] Establish a state model of the rotary kiln and derive the effect of the rotary kiln on material particles;

[0018] Based on the movement of material particles in the rotary kiln, the normal force, tangential force and torque acting on the material particles are calculated, and the speed, position and rotation speed of the material particles are updated;

[0019] When the inclination angle of the material particle accumulation state fluctuates within a threshold range, the movement of the material particles reaches a steady state, and the initialization state of the material particles is output.

[0020] Preferably, the state of the rotary kiln is divided into tilting and rotation. The rotation of the rotary kiln is achieved by setting the direction of the gravitational acceleration acting on the particles. The gravitational acceleration is decomposed into g according to the inclination angle of the rotary kiln. y and g z Two directions, the specific formula is as follows:

[0021] g z =cos(θ)g

[0022] g y =sin(θ)g

[0023] Where θ is the inclination angle of the cement rotary kiln.

[0024] Preferably, the movement of the material particles is divided into translation and rotation. The translation follows Newton's second law, and the rotation follows Euler's second law. The specific formula is as follows:

[0025]

[0026]

[0027] Where m and I represent the mass and moment of inertia of the particle; v and ω represent the translational and rotational velocities of the particle; t is time; g represents the acceleration due to gravity; T is the torque caused by the particle collision; μ r is the rolling friction coefficient of the particles; F n and F t represent the normal and tangential forces respectively.

[0028] Preferably, the normal force and the tangential force are calculated using the discrete element method, and the specific formula is as follows:

[0029] F n =k n δn ij -γ n v nij

[0030]

[0031] Among them, F n and F t are the collision normal force and tangential force, k n , γ n 、k t with γ t are the elastic and damping coefficients of normal and tangential action, δn ij and δt ij is the normal and tangential overlap of the particle, v nij and v tij is the normal and tangential relative velocity of the particle, t is the time, t c,0 is the time of initial collision of particles, μ is the friction coefficient of particles;

[0032] The torque is calculated as follows:

[0033] T=R×F t

[0034] Where T is the torque caused by particle collision, R is the vector from the particle mass center to the collision center, and F t is the collision tangential force;

[0035] The updating of the speed, position and rotation speed of the material particles adopts the Verlet iteration method.

[0036] Preferably, the inclination angle of the accumulation state of the material particles is calculated by drawing the accumulation state of the material particles according to the position and radius of the material particles using visualization software.

[0037] Preferably, the simulation of the reaction and fusion process of the material particles to obtain the state of the material particles at any moment includes:

[0038] The collision of material particles is determined based on their initial state. The colliding particles react according to the reaction kinetics of the rotary kiln and the reaction rate is calculated.

[0039] Calculate and update the composition, mass, density and particle size of material particles based on the reaction rate;

[0040] For two reacting material particles, when the mass of the material particle with reduced mass is less than the threshold at any time, the mass of the material particle with reduced mass is completely added to the material particle with increased mass that reacts with it, and the material particle with reduced mass is deleted from the calculation to realize the fusion process of the material particles;

[0041] The spatial distribution of the material particles at any moment is updated.

[0042] In a second aspect, the present invention provides a device for simulating a cement clinker burning process in a rotary kiln, the device comprising:

[0043] The initial parameter module is used to set the initial parameters of the material particles in the simulation area of ​​the rotary kiln;

[0044] The material particle initialization module is used to simulate the mixing flow process of material particles to reach a steady state and obtain the initialization state of the material particles;

[0045] The reaction and fusion module is used to simulate the reaction and fusion process of material particles and obtain the state of material particles at any moment.

[0046] In a third aspect, the present invention provides an electronic device, comprising:

[0047] at least one processor; and

[0048] a memory communicatively connected to the at least one processor; wherein,

[0049] The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the simulation method of the rotary kiln cement clinker burning process described in the first aspect.

[0050] Compared with the prior art, the present invention has the following beneficial effects:

[0051] The simulation method provided by the present invention is based on the discrete particle simulation method, and simulates the firing process in the rotary kiln from the particle scale. According to the characteristics of the rotary kiln, simulation can be performed for different reaction sections, and the particle scale reaction process of each section can be accurately calculated to realize the simulation of local details such as the radial distribution in the kiln. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] Figure 1 is a flow chart of a method for simulating a rotary kiln cement clinker burning process provided in Example 1;

[0053] Figure 2 1 is a schematic diagram of the rotary kiln provided in Example 1;

[0054] Figure 3 Schematic diagram of the spatial distribution of material particles in the preheating section of the rotary kiln provided in Example 2;

[0055] Figure 4 3 is a schematic structural diagram of a simulation device for a rotary kiln cement clinker burning process provided in Example 3;

[0056] Figure 5 This is a schematic diagram of the structure of an electronic device for simulating the cement clinker burning process in a rotary kiln provided in Example 4. DETAILED DESCRIPTION

[0057] The technical solution of the present invention is further described below by way of specific embodiments. It should be understood by those skilled in the art that the embodiments are merely to help understand the present invention and should not be regarded as specific limitations of the present invention.

[0058] Example 1

[0059] This embodiment provides a method for simulating the cement clinker burning process in a rotary kiln. The method is based on discrete particle simulation and can obtain details of the movement, transfer and reaction processes of particles at different positions in the cement rotary kiln, thereby realizing the modeling of the microscopic movement of cement particles in the rotary kiln. The method can be executed by a simulation device for the cement clinker burning process in a rotary kiln. The device can be implemented in the form of hardware and / or software, and the device can be configured in an electronic device with the ability to simulate the cement clinker burning process in a rotary kiln.

[0060] like Figure 1 As shown, the method provided in this embodiment includes:

[0061] S110. Setting initial parameters of material particles in a simulation area within the rotary kiln.

[0062] The initial parameters include composition, particle size, location and temperature.

[0063] The method of the present invention can simulate the firing process at a local position in a rotary kiln.

[0064] Specifically, based on the temperature measurement of the actual production process or the rotary kiln temperature prediction model, the rotary kiln can be divided into a preheating section, a calcining section, and a cooling section according to the temperature environment, and the firing process can be simulated for any temperature section.

[0065] Specifically, the composition and particle size of the material particles are set according to the temperature environment of the simulation area, the position of the material particles is generated by uniform random numbers, and the temperature of the material particles is set to the average temperature of the simulation area.

[0066] In order to reduce the computational complexity of the simulation as much as possible, the length of the simulation region is set to 50–100 times the maximum particle size.

[0067] S120, simulating the mixing flow process of the material particles to reach a steady state and obtain an initialization state of the material particles, including:

[0068] Establish a state model of the rotary kiln and derive the effect of the rotary kiln on material particles;

[0069] Based on the movement of material particles in the rotary kiln, the normal force, tangential force and torque acting on the material particles are calculated, and the speed, position and rotation speed of the material particles are updated;

[0070] When the inclination angle of the material particle accumulation state fluctuates within a threshold range, the movement of the material particles reaches a steady state, and the initialization state of the material particles is output.

[0071] Specifically, (1.1) a state model of a rotary kiln is established to derive the effect of the rotary kiln on material particles.

[0072] The rotary kiln is abstracted as a cylindrical wall described by mathematics. The cylinder is constructed by mathematical equations, and the main parameters are the position of the center of the cylinder, the positions of the upper and lower bottom surfaces, and the radius of the cylinder.

[0073] The state of the rotary kiln is tilting and rotating. The rotation of the rotary kiln is achieved by setting the direction of the gravitational acceleration acting on the particles. Assuming that the cement rotary kiln is tilted around the x-axis, the gravitational acceleration is decomposed into the y and z directions, such as Figure 2 As shown, the specific calculation formula is:

[0074] g z =cos(θ)g

[0075] g y =sin(θ)g

[0076] Where θ is the inclination angle of the cement rotary kiln.

[0077] The range of sin(θ) is 0.02-0.05.

[0078] (1.2) Based on the movement of the material particles in the rotary kiln, the normal force, tangential force and torque acting on the material particles are calculated, and the velocity, position and rotation speed of the material particles are updated.

[0079] The motion of material particles is divided into translation and rotation. The translation follows Newton's second law, and the rotation follows Euler's second law. The specific formula is as follows:

[0080]

[0081]

[0082] Where m and I represent the mass and moment of inertia of the particle; v and ω represent the translational and rotational velocities of the particle; t is time; g represents the acceleration due to gravity; T is the torque caused by the particle collision; μ r is the rolling friction coefficient of the particles; F n and F t represent the normal and tangential forces respectively.

[0083] As mentioned above, the rotation of the rotary kiln is achieved by setting the direction of the gravitational acceleration acting on the particles. In the above calculation formula, the gravitational acceleration g is selected according to the inclination angle of the rotary kiln. z and g y Sets the gravitational acceleration experienced by particles.

[0084] F n and F t They include the normal and tangential forces generated by the contact between material particles and the contact between material particles and the wall of the rotary kiln.

[0085] The normal force and tangential force are calculated using the discrete element method. The specific formulas are as follows:

[0086] F n =k n δn ij -γ n v nij

[0087]

[0088] Among them, F n and F t are the collision normal force and tangential force, k n , γ n 、k t with γ t are the elastic and damping coefficients of normal and tangential action, δn ij and δt ijis the normal and tangential overlap of the particle, v nij and v tij is the normal and tangential relative velocity of the particle, t is the time, t c,0 is the time of initial collision of particles, μ is the friction coefficient of particles;

[0089] The torque is calculated as follows:

[0090] T=R×F t

[0091] Where T is the torque caused by particle collision, R is the vector from the particle mass center to the collision center, and F t is the collision tangential force.

[0092] According to the normal force, tangential force and torque, the Verlet iteration method is used to update the velocity, position and rotation speed of the material particles, and parallel computing and graphics processing units are used to calculate the material particles.

[0093] (1.3) When the inclination angle of the material particle accumulation state fluctuates within the threshold range, the movement of the material particles reaches a steady state, and the initialization state of the material particles is output.

[0094] By calculating the position and radius of the material particles output in (1.2), the accumulation state of the material particles is plotted using visualization software, and the inclination angle of the accumulation state of the material particles is analyzed. When the fluctuation of the inclination angle remains no more than 5% of the mean, it is considered that the movement of the material particles has reached a steady state, and the initialization state of the material particles is output.

[0095] S130: Simulate the reaction and fusion process of the material particles to obtain the state of the material particles at any moment.

[0096] Specifically, the collision of the material particles is determined based on their initial state, and the reaction proceeds according to the reaction kinetics shown in Table 1. If the two colliding particles contain the species listed in Table 1, the reaction proceeds according to the reaction types listed in Table 1. The reaction rate k is calculated based on the particle temperature and the mass of the different species, using the following formula:

[0097]

[0098] Among them A j 、E j They represent the pre-exponential factor and activation energy of the reaction respectively. The specific values ​​of each reaction are shown in Table 1.

[0099] The firing reaction kinetics of the cement rotary kiln are shown in Table 1. The particle reactions in the rotary kiln are divided into three categories: single-particle reaction, that is, the reaction that can be carried out by a single material, including the decomposition of calcium carbonate and magnesium carbonate; two-particle reaction, that is, the reaction of two materials, including the reaction of calcium oxide and silicon dioxide, the reaction of dicalcium silicate and calcium oxide, and the reaction of calcium oxide and aluminum oxide; three-particle reaction, that is, the reaction involving three materials, including the reaction of calcium oxide, tricalcium silicate and ferric oxide. The two-particle reaction and the three-particle reaction are unified and simplified into the two-particle reaction.

[0100] Table 1

[0101]

[0102]

[0103] Update the parameters of the material particles according to the reaction rate, including composition, mass, density and particle size.

[0104] The calculation formula for the component mass of material particles is as follows:

[0105]

[0106]

[0107] Where subscript k represents the kth particle, w represents other particles that collide with the current particle, subscripts i and j represent reactants and products, respectively. The reactants and products in each reaction are shown in Table 1, t represents the current moment, and t+Δt represents the next moment of the simulation, and R is the reaction rate, in kg / s.

[0108] The calculation formula of material particle mass is as follows:

[0109]

[0110] Wherein, the subscript k represents the kth particle, Nr represents the number of reactants in the particle, Np represents the number of products in the particle, and m i Represents the mass of different types of reactants in the particle, m j Indicates the mass of different types of products in the particle.

[0111] The calculation formula of the density of material particles is as follows:

[0112]

[0113] Wherein, the subscript k represents the kth particle, Nr represents the number of reactants in the particle, Np represents the number of products in the particle, and Y i represents the mass fraction of different types of reactants in the particle, Yj It represents the mass fraction of different types of products in the particle, and ρ is the true density of different types of substances.

[0114] The calculation formula of the particle size of the material particles is as follows:

[0115]

[0116] Wherein, the subscript k represents the kth particle, m represents the total mass of the particles, and ρ is the density of the particles.

[0117] During the reaction process, the particles fuse. This is achieved by adding the mass of the particle with the reduced mass to the particle with the increased mass, while removing the particle with the reduced mass from the calculation.

[0118] The formula for calculating the component mass of the material particles with increased mass is as follows:

[0119]

[0120]

[0121] Among them, the subscript l represents large particles, the subscript s represents small particles, i represents the reactant component, and j represents the product component.

[0122] Small particles are directly deleted from the simulation system without calculating the composition of each substance.

[0123] When the reaction and fusion process of the material particles at any moment is simulated, the position, particle size, and mass of different materials in the particles at the current moment are obtained, and the spatial distribution of the material particles in the rotary kiln can be output.

[0124] The particle size and composition of the material particles can be compared with the measured values ​​of actual production to verify the accuracy of the model.

[0125] Example 2

[0126] This embodiment provides a simulation method for the cement clinker burning process in a rotary kiln, taking the preheating section of the rotary kiln as the research object to further explain the composition of the specific technical solution.

[0127] The method provided in this embodiment includes:

[0128] S110. Setting initial parameters of material particles in a simulation area within the rotary kiln.

[0129] According to the temperature measurement of the actual production process, the rotary kiln is divided into preheating section, calcining section and cooling section, and the specific position of each section is determined according to the temperature.

[0130] In this embodiment, the length of the cement rotary kiln is 75 meters, the preheating section is from the kiln tail to the kiln temperature of about 1200 ° C, the length is about 20 meters, the calcination section is the middle transition zone, the temperature is about 1200-1300 ° C, the length is about 25 meters, and the cooling section is the subsequent about 30 meters, the temperature is lower than 1300 ° C.

[0131] For the preheating section, the initial temperature is approximately 800°C. The material particles consist of limestone, iron ore, and bauxite, and the particle size distribution is set between 1 and 100 μm based on actual manufacturer measurements. Assuming a uniform mixture of material particles, the particles are evenly spaced according to their size distribution. The positions of the material particles are generated using a random uniform number, and the temperature of the material particles is set to the average preheating section temperature of 1000°C.

[0132] According to the diameter of the rotary kiln, the diameter of the cylindrical simulation area is set to 4.7 meters.

[0133] The length of the simulation region was set to 50 times the maximum particle size.

[0134] S120, simulating the mixing flow process of the material particles to reach a steady state, and obtaining an initialization state of the material particles.

[0135] The multi-scale particle calculation software DEMms, a self-developed software of the Institute of Process Engineering of the Chinese Academy of Sciences, is used to simulate the particle falling and mixing flow process.

[0136] (1.1) Establish a state model of the rotary kiln and derive the effect of the rotary kiln on material particles.

[0137] The rotary kiln is in a state of tilting and rotating, and the rotation of the rotary kiln is achieved by setting the direction of gravity acceleration acting on the particles.

[0138] The sine value of the rotary kiln's inclination angle is selected as 0.03, and the gravitational acceleration in the two decomposition directions is calculated as follows:

[0139] g z =(1-0.03 2 ) 1 / 2 ×9.81m / s 2 =9.80m / s 2

[0140] g y =0.03×9.81m / s 2 =0.29m / s 2

[0141] (1.2) Calculate the normal force, tangential force, and torque acting on the material particles, and update the velocity, position, and rotational speed of the material particles.

[0142] The movement of material particles is divided into translation and rotation. The translation follows Newton's second law, and the rotation follows Euler's second law.

[0143] As mentioned above, the rotation of the rotary kiln is achieved by setting the direction of the gravitational acceleration acting on the particles. In the calculation of the movement of the material particles, the gravitational acceleration g is calculated according to the obtained acceleration g z and g y Sets the gravitational acceleration experienced by particles.

[0144] The discrete element method is used to calculate the collision force and torque between particles.

[0145] The Verlet iteration method is used to update the position, velocity and rotation speed of the particles, and parallel computing and graphics processing units (GPUs) are used to achieve high-speed calculations for a large number of particles.

[0146] (1.3) Output the initialization state of the material particles.

[0147] By calculating the position and radius of the material particles output in (1.2), the accumulation state of the material particles is drawn using the visualization software Paraview-5.8.1, and then the inclination angle of the accumulation state of the material particles is analyzed. When the fluctuation of the inclination angle remains no more than 5% of the mean, it is considered that the movement of the material particles has reached a steady state, and the initialization state of the material particles is output, such as Figure 3 As shown in the 0s state.

[0148] S130: Simulate the reaction and fusion process of the material particles to obtain the state of the material particles at any moment.

[0149] The collision conditions of the material particles are determined according to the initialization state of the material particles, and the reaction simulation is performed according to the reaction kinetics shown in Table 1 to calculate the reaction rate.

[0150] The component mass, particle mass, density and particle size of the material particles are updated according to the reaction rate calculation.

[0151] During the reaction of material particles, the fusion process of material particles is simulated: for two reacting material particles, when the mass of the material particle with reduced mass is less than one-fifth of its initial mass at the current moment, the mass of the material particle with reduced mass is completely added to the material particle with increased mass that reacts with it, and the material particle with reduced mass is deleted from the calculation.

[0152] When the reaction and fusion process of the material particles at a certain moment is simulated, the position, particle size, and mass of different materials in the particles at the current moment are obtained, and the spatial distribution of the material particles in the rotary kiln is output, such as Figure 3 As shown in FIG, the spatial distribution results of the material particles in the rotary kiln at different times (5s, 10s, 15s, 20s and 25s) during the rotary kiln firing process were obtained respectively.

[0153] Example 3

[0154] This embodiment provides a simulation device for the cement clinker burning process in a rotary kiln. Figure 4 As shown, the device includes: an initial parameter module 110, a material particle initialization module 120 and a reaction and fusion module 130, wherein:

[0155] The initial parameter module 110 is used to set the initial parameters of the material particles in the simulation area of ​​the rotary kiln;

[0156] The material particle initialization module 120 is used to simulate the mixing flow process of the material particles to reach a steady state and obtain the initialization state of the material particles;

[0157] The reaction and fusion module 130 is used to simulate the reaction and fusion process of material particles to obtain the state of the material particles at any time.

[0158] In the initial parameter module 110, the initial parameters of the material particles in the simulation area of ​​the rotary kiln are set, including:

[0159] Set the composition and particle size of the material particles according to the temperature environment of the simulation area;

[0160] Generate the position of material particles by uniform random numbers;

[0161] The temperature of the material particles is set to the average temperature of the simulation area.

[0162] In the material particle initialization module 120, the mixed flow process of the simulated material particles reaches a steady state, and the initialization state of the material particles is obtained, including:

[0163] Establish a state model of the rotary kiln and derive the effect of the rotary kiln on material particles;

[0164] Based on the movement of material particles in the rotary kiln, the normal force, tangential force and torque acting on the material particles are calculated, and the speed, position and rotation speed of the material particles are updated;

[0165] When the inclination angle of the material particle accumulation state fluctuates within a threshold range, the movement of the material particles reaches a steady state, and the initialization state of the material particles is output.

[0166] In the reaction and fusion module 130, the reaction and fusion process of the material particles is simulated to obtain the state of the material particles at any time, including:

[0167] The collision of material particles is determined based on their initial state. The colliding particles react according to the reaction kinetics of the rotary kiln and the reaction rate is calculated.

[0168] Calculate and update the composition, mass, density and particle size of material particles based on the reaction rate;

[0169] For two reacting material particles, when the mass of the material particle with reduced mass is less than the threshold at any time, the mass of the material particle with reduced mass is completely added to the material particle with increased mass that reacts with it, and the material particle with reduced mass is deleted from the calculation to realize the fusion process of the material particles;

[0170] The spatial distribution of the material particles at any moment is updated.

[0171] The device provided in this embodiment can execute the simulation method of the rotary kiln cement clinker burning process provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0172] Example 4

[0173] This embodiment provides an electronic device for implementing a method for simulating a cement rotary kiln clinker burning process, such as Figure 5 As shown, the electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital assistants, cellular phones, smart phones, wearable devices (such as helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present application described and / or required herein.

[0174] like Figure 5As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12, a random access memory (RAM) 13, etc., which is communicatively connected to the at least one processor 11. The memory stores a computer program that can be executed by the at least one processor. The processor 11 can perform various appropriate actions and processes according to the computer program stored in the read-only memory (ROM) 12 or the computer program loaded from the storage unit 18 into the random access memory (RAM) 13. Various programs and data required for the operation of the electronic device 10 can also be stored in the RAM 13. The processor 11, ROM 12, and RAM 13 are connected to each other via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0175] Multiple components in the electronic device 10 are connected to the I / O interface 15, including an input unit 16, such as a keyboard, a mouse, etc.; an output unit 17, such as various types of displays, speakers, etc.; a storage unit 18, such as a second storage area, an optical disk, etc.; and a communication unit 19, such as a network card, a modem, a wireless communication transceiver, etc. The communication unit 19 allows the electronic device 10 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks.

[0176] The processor 11 can be any general-purpose and / or specialized processing component with processing and computing capabilities. Some examples of the processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various specialized artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The processor 11 executes the various methods and processes described above, such as the method for simulating the rotary kiln cement clinker burning process.

[0177] In some embodiments, the simulation method of the rotary kiln cement clinker burning process can be implemented as a computer program, which is tangibly contained in a computer-readable storage medium, such as a storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on the electronic device 10 via the ROM 12 and / or the communication unit 19. When the computer program is loaded into the RAM 13 and executed by the processor 11, one or more steps of the simulation method of the rotary kiln cement clinker burning process described above can be performed. Alternatively, in other embodiments, the processor 11 can be configured to execute the simulation method of the rotary kiln cement clinker burning process by any other appropriate means (for example, by means of firmware).

[0178] Various embodiments of the systems and techniques described herein can be implemented in digital electronic circuit systems, integrated circuit systems, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), application specific standard products (ASSPs), system-on-chip systems (SOCs), programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments can include being implemented in one or more computer programs that are executable and / or interpreted on a programmable system that includes at least one programmable processor, which can be a special purpose or general purpose programmable processor that can receive data and instructions from a storage system, at least one input device, and at least one output device, and transmit data and instructions to the storage system, the at least one input device, and the at least one output device.

[0179] Computer programs for implementing the methods of the present application can be written in any combination of one or more programming languages. These computer programs can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable target-determining device, so that when the computer program is executed by the processor, the functions / operations specified in the flowcharts and / or block diagrams are implemented. The computer program can be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0180] In the context of the present application, a computer-readable storage medium can be a tangible medium that can contain or store a computer program for use by an instruction execution system, device or equipment or used in combination with an instruction execution system, device or equipment. A computer-readable storage medium can include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared or semiconductor systems, devices or equipment, or any suitable combination of the foregoing. Alternatively, a computer-readable storage medium can be a machine-readable signal medium. A more specific example of a machine-readable storage medium can include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0181] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0182] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), a blockchain network, and the Internet.

[0183] A computing system may include clients and servers. The clients and servers are typically remote from each other and typically interact via a communication network. This client-server relationship arises through computer programs running on the respective computers, creating a client-server relationship. The server may be a cloud server, also known as a cloud computing server or cloud host. This server is a hosting product within the cloud computing service ecosystem that addresses the management difficulties and limited scalability of traditional physical hosting and VPS services.

[0184] It should be understood that the various forms of the processes shown above can be used to reorder, add, or delete steps. For example, the steps described in this application can be performed in parallel, sequentially, or in a different order, as long as the desired information of the technical solution of this application can be achieved. This document is not limited here.

[0185] The applicant declares that the above is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Those skilled in the art should understand that any changes or substitutions that can be easily thought of by those skilled in the art within the technical scope disclosed by the present invention fall within the scope of protection and disclosure of the present invention.

Claims

1. A method for simulating a rotary kiln cement clinker burning process, characterized in that: The simulation method comprises: Set the initial parameters of the material particles in the simulation area of ​​the rotary kiln; Simulate the mixing flow process of material particles to reach a steady state and obtain the initialization state of the material particles; The reaction and fusion process of material particles are simulated to obtain the state of material particles at any moment.

2. The simulation method according to claim 1, characterized in that The initial parameters of the material particles include the composition, particle size, position and temperature of the material particles.

3. The simulation method according to claim 1 or 2, characterized in that: The simulated mixing flow process of the material particles reaches a steady state and obtains an initialization state of the material particles, including: Establish a state model of the rotary kiln and derive the effect of the rotary kiln on material particles; Based on the movement of material particles in the rotary kiln, the normal force, tangential force and torque acting on the material particles are calculated, and the speed, position and rotation speed of the material particles are updated; When the inclination fluctuation of the accumulation state of the material particles remains within the threshold range, the movement of the material particles reaches a steady state, and the initialization state of the material particles is output.

4. The simulation method according to claim 3, characterized in that The states of the rotary kiln are divided into tilting and rotation. The rotation of the rotary kiln is achieved by setting the direction of the gravitational acceleration acting on the particles. The gravitational acceleration is decomposed into g according to the inclination angle of the rotary kiln. y and g z Two directions, the specific formula is as follows: g z =cos(θ)g g y =sin(θ)g Where θ is the inclination angle of the cement rotary kiln.

5. The simulation method according to claim 3 or 4, characterized in that: The motion of the material particles is divided into translation and rotation. The translation follows Newton's second law, and the rotation follows Euler's second law. The specific formula is as follows: Where m and I represent the mass and moment of inertia of the particle; v and ω represent the translational and rotational velocities of the particle; t is time; g represents the acceleration due to gravity; T is the torque caused by the particle collision; μ r is the rolling friction coefficient of the particles; F n and F t represent the normal and tangential forces respectively.

6. The simulation method according to any one of claims 3 to 5, characterized in that: The normal force and tangential force are calculated using the discrete element method, and the specific formula is as follows: F n =k n δn ij -γ n v nij F t =k t δt ij -c t v tij , max(|F t |)=|μF n | Among them, F n and F t are the collision normal force and tangential force, k n , γ n 、k t with γ t are the elastic and damping coefficients of normal and tangential action, δn ij and δt ij is the normal and tangential overlap of the particle, v nij and v tij is the normal and tangential relative velocity of the particle, t is the time, t c,0 is the time of initial collision of particles, μ is the friction coefficient of particles; The torque is calculated as follows: T=R×F t Where T is the torque caused by particle collision, R is the vector from the particle mass center to the collision center, and F t is the collision tangential force; The updating of the speed, position and rotation speed of the material particles adopts the Verlet iteration method.

7. The simulation method according to any one of claims 3 to 6, characterized in that: The inclination angle of the material particle accumulation state is calculated by drawing the material particle accumulation state according to the position and radius of the material particles using visualization software.

8. The simulation method according to any one of claims 1 to 7, characterized in that: The simulation of the reaction and fusion process of the material particles to obtain the state of the material particles at any moment includes: The collision of material particles is determined according to the initial state of the material particles. The colliding material particles react according to the reaction kinetics of the rotary kiln and the reaction rate is calculated. Calculate and update the composition, mass, density and particle size of material particles based on the reaction rate; For two reacting material particles, when the mass of the material particle with reduced mass is less than the threshold at any time, the mass of the material particle with reduced mass is completely added to the material particle with increased mass that reacts with it, and the material particle with reduced mass is deleted from the calculation to realize the fusion process of the material particles; The spatial distribution of the material particles at any moment is updated.

9. A device for simulating the cement clinker burning process in a rotary kiln, characterized in that: The device comprises: The initial parameter module is used to set the initial parameters of the material particles in the simulation area of ​​the rotary kiln; The material particle initialization module is used to simulate the mixing flow process of material particles to reach a steady state and obtain the initialization state of the material particles; The reaction and fusion module is used to simulate the reaction and fusion process of material particles and obtain the state of material particles at any moment.

10. An electronic device, characterized in that: The electronic device comprises: at least one processor; and a memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, and the computer program is executed by the at least one processor so that the at least one processor can execute the simulation method of the rotary kiln cement clinker burning process as described in any one of claims 1-8.

Images