A method, system and information data processing terminal for obtaining fluid parameters in a flue

By establishing a three-dimensional model of the flue and performing iterative calculations of turbulent kinetic energy and diffusion equations, the problem of lack of prediction of emission parameter monitoring in the prior art is solved, and scientific prediction and more accurate control of flue emission data are achieved.

CN114638177BActive Publication Date: 2025-01-10TIANJIN CEMENT IND DESIGN & RES INST CO LTD
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202210293309.9
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2022-03-23
Publication Date
2025-01-10
Estimated Expiration
2042-03-23

AI Technical Summary

Technical Problem

The monitoring of emission parameters in the flue in the prior art mainly relies on the monitoring and analysis of actual numerical values, and lacks early prediction of emission working parameters, which makes it difficult to control emission time, rate and quantity.

Method used

By designing a three-dimensional model of the flue, a turbulent kinetic energy equation and a turbulent diffusion equation are established, and iterative calculations are performed in combination with boundary conditions to obtain the flue gas flow rate and other physical quantities, and then the emission data in the flue are predicted.

Benefits of technology

Scientific guidance and planning of flue emission work have been achieved, and the accuracy and efficiency of emission control have been improved.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN114638177B_ABST
    Figure CN114638177B_ABST
Patent Text Reader

Abstract

The present invention discloses a method, system and information data processing terminal for obtaining fluid parameters in a flue, belonging to the field of information processing technology, and comprising: S1, designing a three-dimensional model of the flue according to the structural parameters of the flue; S2, pre-processing the three-dimensional model to generate a calculation grid; S3, establishing a turbulent kinetic energy equation k and a turbulent diffusion equation e based on the grid and given boundary conditions; S4, obtaining the physical quantity of each grid through iterative calculation, and calculating the residual between each iteration and the previous one, and obtaining the range of stable distribution of the flue gas velocity by judging the residual; S5, first calculating the velocity of the continuous phase using the turbulent kinetic energy equation k and the turbulent diffusion equation e, or the Reynolds stress equations in three directions; then obtaining the fluid density according to the turbulent kinetic energy equation k and the turbulent diffusion equation e; calculating the kinetic energy of the particle according to the energy conservation equation; obtaining the temperature according to the energy equation, and obtaining the gas pressure according to the loss along the Bernoulli equation.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention belongs to the field of information processing technology, and in particular relates to a method and system for acquiring parameters of a fluid in a flue, and an information data processing terminal. Background Art

[0002] As we all know, flue pipe is a tubular device for exhaust gas and smoke emission. The size and shape of flue pipe is one of the important factors affecting the emission rate and emission amount of flue gas. In recent years, with the rapid development of industry, environmental pollution has become increasingly serious. In order to achieve sustainable development, humans need to control the total emission of polluted gas sources while purifying the polluted gas before emission. To this end, humans set up various sensors in the flue (such as harmful gas content monitoring, flow rate monitoring, temperature and humidity monitoring, etc.), and then perform statistical analysis on various parameters of flue gas based on the data detected by the sensors. Finally, emission rules are formulated. However, practice has found that:

[0003] The above traditional scheme is based on the monitoring and analysis of actual numerical values. However, in actual production, if the various working parameters of the flue can be predicted in advance, it will play a certain guiding role in production emissions (such as emission time, emission rate and emission amount control). Summary of the invention

[0004] Technical Purpose

[0005] The present invention provides a method and system for acquiring parameters of fluid in a flue, and an information data processing terminal; through a computer model, the emission data of the flue can be acquired in advance, thereby providing scientific guidance and planning for production emission work.

[0006] Technical Solution

[0007] The first object of the present invention is to provide a method for obtaining parameters of a fluid in a flue, comprising the following steps:

[0008] S1. Designing a three-dimensional model of the flue according to the structural parameters of the flue; the structural parameters include the inner diameter of the flue, the height of the flue, the position of the fan inlet, and the angle of the fan inlet;

[0009] S2, pre-processing the three-dimensional model to generate a computational grid;

[0010] S3, establishing a turbulent kinetic energy equation k and a turbulent diffusion equation e based on the grid and given boundary conditions;

[0011] The turbulent kinetic energy equation k is:

[0012] The turbulent diffusion equation e is:

[0013]

[0014] Among them, G k is the turbulent kinetic energy generated by the laminar velocity gradient; G b is the turbulent kinetic energy generated by buoyancy; Y m is the fluctuation caused by excessive diffusion in compressible turbulence; σ k , σ ε are the turbulent Prandtl numbers of the turbulent kinetic energy equation k and the turbulent diffusion equation e respectively; S k and S ε is the custom coefficient; C 1ε , C 2ε , C 3ε is a constant; a{i, j} = du i / dx j It represents 9 equations, i=1, j=1, 2, 3; i=2, j=1, 2, 3; i=3, j=1, 2, 3; that is, the derivatives of u in the x, y, z directions, the derivatives of v in the x, y, z directions, and the derivatives of w in the x, y, z directions.

[0015] S4. Obtaining the physical quantity of each grid through iterative calculation, wherein the physical quantity includes the flue gas flow rate, and calculating the residual between each iteration and the previous iteration, and obtaining the range of the stable distribution of the flue gas flow rate by judging the residual;

[0016] S5. First, use the turbulent kinetic energy equation k and the turbulent diffusion equation e, or the Reynolds stress equations in three directions to calculate the velocity of the continuous phase; then obtain the fluid density based on the turbulent kinetic energy equation k and the turbulent diffusion equation e; obtain the fluid mass based on the fluid density and volume; calculate the particle kinetic energy based on the energy conservation equation; calculate the gas buoyancy based on the external density, mass and volume; obtain the temperature based on the energy equation, and then obtain the gas density based on the ideal gas state equation; obtain the gas pressure based on the along-the-line loss of the Bernoulli equation.

[0017] Preferably, S1 specifically includes: designing a three-dimensional model using Catia according to structural parameters of the flue.

[0018] Preferably, in S2, the pre-processing is specifically as follows: based on the three-dimensional shape of the flue generated in the early stage, Gambit software is used for pre-processing calculation to generate a computational grid with a grid scale within 10 mm. The computational grid uses a hybrid grid form combining a hexahedral structured grid and a polyhedral unstructured grid to complete modeling and meshing.

[0019] Preferably, in S3, the boundary conditions include: the blowing volume of the fan, the air inlet velocity, the flue gas volume, the flue outlet pressure, the flow rate at the flue inlet, the volume composition of the flue gas, the flue gas temperature, the flue gas density, and the dynamic viscosity of the flue gas.

[0020] A second object of the present invention is to provide a system for acquiring parameters of fluid in a flue, comprising:

[0021] A modeling module, which designs a three-dimensional model of the flue according to the structural parameters of the flue; the structural parameters include the inner diameter of the flue, the height of the flue, the position of the fan inlet, and the angle of the fan inlet;

[0022] A pre-processing module performs pre-processing on the three-dimensional model to generate a computational grid;

[0023] An association module, based on the grid and given boundary conditions, establishes a turbulent kinetic energy equation k and a turbulent diffusion equation e;

[0024] The turbulent kinetic energy equation k is:

[0025] The turbulent diffusion equation e is:

[0026]

[0027] Among them, G k is the turbulent kinetic energy generated by the laminar velocity gradient; G b is the turbulent kinetic energy generated by buoyancy; Y m is the fluctuation caused by excessive diffusion in compressible turbulence; σ k , σ ε are the turbulent Prandtl numbers of the turbulent kinetic energy equation k and the turbulent diffusion equation e respectively; S k and S ε is the custom coefficient; C 1ε , C 2ε , C 3ε is a constant; a{i, j} = du i / dx j It represents 9 equations, i=1, j=1, 2, 3; i=2, j=1, 2, 3; i=3, j=1, 2, 3; that is, the derivatives of u in the x, y, z directions, the derivatives of v in the x, y, z directions, and the derivatives of w in the x, y, z directions.

[0028] The information processing module obtains the physical quantity of each grid through iterative calculation, wherein the physical quantity includes the flue gas flow rate, and calculates the residual between each iteration and the previous iteration, and obtains the range of the stable distribution of the flue gas flow rate by judging the residual;

[0029] The acquisition module first calculates the velocity of the continuous phase using the turbulent kinetic energy equation k and the turbulent diffusion equation e, or the Reynolds stress equations in three directions; then obtains the fluid density based on the turbulent kinetic energy equation k and the turbulent diffusion equation e; obtains the fluid mass based on the fluid density and volume; calculates the particle kinetic energy based on the energy conservation equation; calculates the gas buoyancy based on the external density, mass and volume; obtains the temperature based on the energy equation, and then obtains the gas density based on the ideal gas state equation; and obtains the gas pressure based on the along-the-line loss of the Bernoulli equation.

[0030] Preferably, Catia is used to design a three-dimensional model according to the structural parameters of the flue.

[0031] Preferably, in the pre-processing module, the pre-processing is specifically as follows: based on the three-dimensional shape of the flue generated in the early stage, Gambit software is used for pre-processing calculation to generate a computational grid with a grid scale within 10 mm. The computational grid uses a hybrid grid form combining a hexahedral structured grid and a polyhedral unstructured grid to complete modeling and meshing.

[0032] Preferably, in the associated module, the boundary conditions include: the blower volume, air inlet velocity, flue gas volume, flue outlet pressure, flue inlet flow, flue gas volume composition, flue gas temperature, flue gas density, and flue gas dynamic viscosity.

[0033] The third invention objective of this patent is to provide an information data processing terminal that implements the above-mentioned method for obtaining the fluid parameters in the flue.

[0034] The fourth invention objective of this patent is to provide a computer-readable storage medium, including instructions, which, when run on a computer, enables the computer to execute the above-mentioned method for obtaining the fluid parameters in the flue.

[0035] The advantages and positive effects of the present invention are as follows: the present invention can more intuitively and directly characterize the flow process of gas in the pipeline, help understand parameters such as gas velocity distribution in the pipeline, and provide a basis for arranging flue gas measurement equipment in the pipeline, thereby improving the accuracy of flue gas velocity measurement in the pipeline. Through the computer model, the emission data of the flue can be obtained in advance, thereby providing scientific guidance and planning for production emission work. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 is a flow chart of a preferred embodiment of the present invention;

[0037] Figure 2 is the three-dimensional model obtained in S1 in the preferred embodiment of the present invention;

[0038] Figure 3 is a simulation diagram obtained in S4 in a preferred embodiment of the present invention;

[0039] Figure 4 It is a system block diagram of a preferred embodiment of the present invention. DETAILED DESCRIPTION

[0040] In order to further understand the content, features and effects of the present invention, the following embodiments are given and described in detail as follows.

[0041] See also Figure 1 , a method for obtaining fluid parameters in a flue, comprising the following steps:

[0042] S1. According to the structural parameters of the flue (the structural parameters are derived from the flue drawings and on-site measured data), use Catia (other software can also be used) to design the three-dimensional model of the flue; keep it consistent with the shape of the actual chimney, and provide real data input for the subsequent flow field simulation and optimization design. The structural parameters include the inner diameter of the flue, the height of the flue, the position of the fan outlet, the angle of the fan outlet, the fan outlet diameter, the fan outlet flow rate and temperature; taking a 150m high cylindrical flue as an example, the flue inner cavity diameter is 4m. There is a fan inlet 2 inserted at a certain angle at the bottom of the flue, and the fan inlet diameter is 3.6m to 4m. The fan outlet and the entire flue are made of metal, and the inner and outer walls are relatively smooth. The wall of the flue 1 is regarded as an insulating wall, and there is no heat exchange with the external environment; the three-dimensional model is as follows Figure 2 As shown;

[0043] S2. Based on the three-dimensional model generated in the early stage, the three-dimensional geometric shape of the flue is pre-processed to generate a computational grid with a total number of about 5 million grids;

[0044] The pre-processing is specifically as follows: based on the three-dimensional shape of the flue generated in the early stage, Gambit software is used for pre-processing calculation to generate a computational grid with a grid scale within 10 mm. The computational grid uses a hybrid grid form combining a hexahedral structured grid and a polyhedral unstructured grid to complete modeling and meshing.

[0045] S3. Based on the grid and given boundary conditions, assuming that the flue gas is an incompressible ideal flow, according to fluid theory, the airflow in the flue is a three-dimensional turbulent flow, and the turbulent kinetic energy equation k and the turbulent diffusion equation e are established;

[0046] The turbulent kinetic energy equation k is:

[0047] The turbulent diffusion equation e is:

[0048]

[0049] Among them, G k is the turbulent kinetic energy generated by the laminar velocity gradient; G b is the turbulent kinetic energy generated by buoyancy; Y mis the fluctuation caused by excessive diffusion in compressible turbulence; σ k , σ ε are the turbulent Prandtl numbers of the turbulent kinetic energy equation k and the turbulent diffusion equation e respectively; S k and S ε is the custom coefficient; C 1ε , C 2ε , C 3ε is a constant; a{i, j} = du i / dx j It represents 9 equations, i=1, j=1, 2, 3; i=2, j=1, 2, 3; i=3, j=1, 2, 3; that is, the derivatives of u in the x, y, z directions, the derivatives of v in the x, y, z directions, and the derivatives of w in the x, y, z directions.

[0050] i and j are representations of the free coordinate system. Since turbulence is the partial derivative of a three-dimensional equation, that is, the derivative of each coordinate direction relative to the three spatial coordinates, a total of 9 different partial derivatives will be generated. i represents the current direction and j represents the relative direction, with values ​​ranging from 1 to 3.

[0051] The given parameters (boundary conditions) are: the flue gas flow rate is determined by the flue gas velocity. Different blower 3 air volumes determine different flue gas velocity (flow rate). Given the working condition air volume inlet velocity of 16-28 m / s, the flue gas volume Q = 7.24 × 10 5 ~1.27×10 6 m 3 / h; given flue outlet pressure (slightly positive or slightly negative), inlet pressure is calculated based on outlet pressure; flue inlet flow rate is set; given flue gas composition (volume fraction): moisture (10%), O 2 (9%), CO 2 (22%), N 2 (59%); the fluid is incompressible; given the flue gas temperature, density and dynamic viscosity;

[0052] S4. Through iterative calculation, the physical quantities such as the flue gas velocity of each grid can be obtained, and the residual between each iteration and the previous one can be calculated. By judging the residual, the range of the stable distribution of the flue gas velocity can be obtained. Figure 3 The results show that the flue gas velocity of the profile is relatively stable (6D range) within the height range of 20-30m under different working conditions;

[0053] S5. Use the turbulent kinetic energy equation k and the turbulent diffusion equation e or the three-direction Reynolds stress equation to calculate the velocity of the continuous phase, and then calculate the density based on the turbulent kinetic energy equation k and the turbulent diffusion equation e. The mass can be calculated based on the density (m=ρV). The kinetic energy of the particle is calculated based on the law of conservation of energy. Then the gas buoyancy is calculated based on the external density, mass and volume (G=mg, m=ρV). The temperature is calculated based on the energy equation, and then the gas density is calculated based on the ideal gas state equation. The gas pressure is calculated based on the loss along the Bernoulli equation.

[0054] See also Figure 4 , a system for acquiring fluid parameters in a flue, comprising:

[0055] The modeling module uses Catia (other software can also be used) to design the three-dimensional model of the flue according to the structural parameters of the flue (the structural parameters are derived from the flue drawings and on-site measured data); make it consistent with the shape of the real chimney, and provide real data input for the subsequent flow field simulation and optimization design. The structural parameters include the inner diameter of the flue, the height of the flue, the position of the fan outlet, the angle of the fan outlet, the fan outlet diameter, the fan outlet flow rate and temperature; taking a 150m high cylindrical flue as an example, the flue inner cavity diameter is 4m. There is a fan inlet inserted at a certain angle at the bottom of the flue, and the fan inlet diameter is 3.6m to 4m. The fan outlet and the entire flue are made of metal, and the inner and outer walls are relatively smooth. The flue wall is regarded as an adiabatic wall, and there is no heat exchange with the external environment; the three-dimensional model is as follows Figure 2 As shown;

[0056] The pre-processing module performs calculation pre-processing on the three-dimensional geometric shape of the flue based on the three-dimensional model generated in the early stage, and generates a calculation grid with a total grid size of about 5 million;

[0057] The pre-processing is specifically as follows: based on the three-dimensional shape of the flue generated in the early stage, Gambit software is used for pre-processing calculation to generate a computational grid with a grid scale within 10 mm. The computational grid uses a hybrid grid form combining a hexahedral structured grid and a polyhedral unstructured grid to complete modeling and meshing.

[0058] Association module, based on the grid and given boundary conditions, assuming that the flue gas is an incompressible ideal flow, according to fluid theory, the airflow in the flue is a three-dimensional turbulent flow, and the turbulent kinetic energy equation k and the turbulent diffusion equation e are established;

[0059] The turbulent kinetic energy equation k is:

[0060] The turbulent diffusion equation e is:

[0061]

[0062] Among them, Gk is the turbulent kinetic energy generated by the laminar velocity gradient; G b is the turbulent kinetic energy generated by buoyancy; Y m is the fluctuation caused by excessive diffusion in compressible turbulence; σ k , σ ε are the turbulent Prandtl numbers of the turbulent kinetic energy equation k and the turbulent diffusion equation e respectively; S k and S ε is the custom coefficient; C 1ε , C 2ε , C 3ε is a constant; a{i, j} = du i / dx j It represents 9 equations, i=1, j=1, 2, 3; i=2, j=1, 2, 3; i=3, j=1, 2, 3; that is, the derivatives of u in the x, y, z directions, the derivatives of v in the x, y, z directions, and the derivatives of w in the x, y, z directions.

[0063] The given parameters (boundary conditions) are: the flue gas flow rate is determined by the flue gas velocity. Different fan air volumes determine different flue gas velocity (flow rate). Given the working condition, the inlet air velocity is 16-28 m / s, and the flue gas volume Q = 7.24 × 10 5 ~1.27×10 6 m 3 / h; given flue outlet pressure (slightly positive or slightly negative), inlet pressure is calculated based on outlet pressure; flue inlet flow rate is set; given flue gas composition (volume fraction): moisture (10%), O 2 (9%), CO 2 (22%), N 2 (59%); the fluid is incompressible; given the flue gas temperature, density and dynamic viscosity;

[0064] The information processing module can obtain the physical quantities such as the flue gas velocity of each grid through iterative calculation, and calculate the residual between each iteration and the previous one. By judging the residual, the range of the stable distribution of the flue gas velocity can be obtained. Figure 3 The results show that the flue gas velocity of the profile is relatively stable (6D range) within the height range of 20-30m under different working conditions;

[0065] The acquisition module calculates the velocity of the continuous phase using the turbulent kinetic energy equation k and the turbulent diffusion equation e or the three-direction Reynolds stress equation, and then calculates the density based on the turbulent kinetic energy equation k and the turbulent diffusion equation e. The mass can be calculated based on the density (m=ρV). The kinetic energy of the particle is calculated based on the law of conservation of energy. The gas buoyancy is then calculated based on the external density, mass and volume (G=mg, m=ρV). The temperature is calculated based on the energy equation, and then the gas density is calculated based on the ideal gas state equation. The gas pressure is calculated based on the loss along the Bernoulli equation.

[0066] An information data processing terminal for implementing the above-mentioned method for acquiring parameters of fluid in a flue.

[0067] A computer-readable storage medium includes instructions, which, when executed on a computer, enable the computer to execute the above-mentioned method for obtaining parameters of fluid in a flue.

[0068] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When the use is implemented in whole or in part in the form of a computer program product, the computer program product includes one or more computer instructions. When the computer program instructions 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 may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may 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 may be transmitted from one website site, computer, server or data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL) or wireless (e.g., infrared, wireless, microwave, etc.) mode) to another website site, computer, server or data center. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated. The available medium may 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 (e.g., a solid-state hard disk SolidState Disk (SSD)), etc.

[0069] The above description is only a preferred embodiment of the present invention and does not limit the present invention in any form. Any simple modification, equivalent changes and modifications made to the above embodiments based on the technical essence of the present invention are within the scope of the technical solution of the present invention.

Claims

1. A method for obtaining fluid parameters in a flue, characterized in that: The steps include: S1. Designing a three-dimensional model of the flue according to the structural parameters of the flue; the structural parameters include the inner diameter of the flue, the height of the flue, the position of the fan inlet, and the angle of the fan inlet; S2, pre-processing the three-dimensional model to generate a computational grid; S3, establishing a turbulent kinetic energy equation k and a turbulent diffusion equation e based on the grid and given boundary conditions; The turbulent kinetic energy equation k is: The turbulent diffusion equation e is: Among them, G k is the turbulent kinetic energy generated by the laminar velocity gradient; G b is the turbulent kinetic energy generated by buoyancy; Y m is the fluctuation caused by excessive diffusion in compressible turbulence; k , σ ε are the turbulent Prandtl numbers of the turbulent kinetic energy equation k and the turbulent diffusion equation e respectively; S k and S ε is the custom coefficient; C 1ε , C 2ε , C 3ε is a constant; a{i, j} = du i / dx j It represents 9 equations, i=1, j=1, 2, 3; i=2, j=1, 2, 3; i=3, j=1, 2, 3; i.e. the derivatives of u in the x, y, z directions, the derivatives of v in the x, y, z directions, and the derivatives of w in the x, y, z directions; S4. Obtaining the physical quantity of each grid through iterative calculation, wherein the physical quantity includes the flue gas flow rate, and calculating the residual between each iteration and the previous iteration, and obtaining the range of the stable distribution of the flue gas flow rate by judging the residual; S5. First, use the turbulent kinetic energy equation k and the turbulent diffusion equation e, or the Reynolds stress equations in three directions to calculate the velocity of the continuous phase; then obtain the fluid density based on the turbulent kinetic energy equation k and the turbulent diffusion equation e; obtain the fluid mass based on the fluid density and volume; calculate the particle kinetic energy based on the energy conservation equation; calculate the gas buoyancy based on the external density, mass and volume; obtain the temperature based on the energy equation, and then obtain the gas density based on the ideal gas state equation; obtain the gas pressure based on the along-the-line loss of the Bernoulli equation.

2. The method for obtaining the parameters of the fluid in the flue according to claim 1, characterized in that: S1 specifically includes: using Catia to design a three-dimensional model based on the structural parameters of the flue.

3. The method for obtaining the parameters of the fluid in the flue according to claim 1, characterized in that: In S2, the pre-processing is specifically as follows: based on the three-dimensional shape of the flue generated in the early stage, Gambit software is used for pre-processing calculation to generate a computational grid, and the computational grid adopts a hybrid grid form combining a hexahedral structured grid and a polyhedral unstructured grid to complete modeling and meshing.

4. The method for obtaining the parameters of the fluid in the flue according to claim 1, characterized in that: In S3, the boundary conditions include: the blowing volume of the fan, the air volume inlet velocity, the smoke volume, the flue outlet pressure, the flow rate at the flue inlet, the volume composition of the smoke, the smoke temperature, the smoke density, and the dynamic viscosity of the smoke.

5. A system for acquiring fluid parameters in a flue, characterized in that: At least: Modeling module, which designs the three-dimensional model of the flue according to the structural parameters of the flue; The structural parameters include the inner diameter of the flue, the height of the flue, the position of the fan inlet, and the angle of the fan inlet; A pre-processing module performs pre-processing on the three-dimensional model to generate a computational grid; An association module, based on the grid and given boundary conditions, establishes a turbulent kinetic energy equation k and a turbulent diffusion equation e; The turbulent kinetic energy equation k is: The turbulent diffusion equation e is: Among them, G k is the turbulent kinetic energy generated by the laminar velocity gradient; G b is the turbulent kinetic energy generated by buoyancy; Y m is the fluctuation caused by excessive diffusion in compressible turbulence; k , σ ε are the turbulent Prandtl numbers of the turbulent kinetic energy equation k and the turbulent diffusion equation e respectively; S k and S ε is the custom coefficient; C 1ε , C 2ε , C 3ε is a constant; a{i, j} = du i / dx j It represents 9 equations, i=1, j=1, 2, 3; i=2, j=1, 2, 3; i=3, j=1, 2, 3; i.e. the derivatives of u in the x, y, z directions, the derivatives of v in the x, y, z directions, and the derivatives of w in the x, y, z directions; The information processing module obtains the physical quantity of each grid through iterative calculation, wherein the physical quantity includes the flue gas flow rate, and calculates the residual between each iteration and the previous iteration, and obtains the range of the stable distribution of the flue gas flow rate by judging the residual; The acquisition module first calculates the velocity of the continuous phase using the turbulent kinetic energy equation k and the turbulent diffusion equation e, or the Reynolds stress equations in three directions; then obtains the fluid density based on the turbulent kinetic energy equation k and the turbulent diffusion equation e; obtains the fluid mass based on the fluid density and volume; calculates the particle kinetic energy based on the energy conservation equation; calculates the gas buoyancy based on the external density, mass and volume; obtains the temperature based on the energy equation, and then obtains the gas density based on the ideal gas state equation; and obtains the gas pressure based on the along-the-line loss of the Bernoulli equation.

6. The system for acquiring parameters of fluid in a flue according to claim 5, characterized in that: In the modeling module: Catia is used to design a three-dimensional model based on the structural parameters of the flue.

7. The system for acquiring parameters of fluid in a flue according to claim 5, characterized in that: In the pre-processing module, the pre-processing is specifically as follows: based on the three-dimensional shape of the flue generated in the early stage, Gambit software is used for pre-processing calculation to generate a computational grid with a grid scale within 10 mm. The computational grid uses a hybrid grid form that combines a hexahedral structured grid and a polyhedral unstructured grid to complete modeling and meshing.

8. The system for acquiring parameters of fluid in a flue according to claim 5, characterized in that: In the associated module, the boundary conditions include: the blower volume, air inlet velocity, flue gas volume, flue outlet pressure, flue inlet flow, flue gas volume composition, flue gas temperature, flue gas density, and flue gas dynamic viscosity.

9. An information data processing terminal for implementing the method for acquiring parameters of fluid in a flue as claimed in any one of claims 1 to 4.

10. A computer-readable storage medium comprising instructions, which, when executed on a computer, enables the computer to execute the method for obtaining parameters of a fluid in a flue as claimed in any one of claims 1 to 4.

Citation Information

Patent Citations

  • Liquid metal flow heat transfer calculation method

    CN113095004A

  • Aerodynamic derivatives calculation method for flight vehicle

    US20200410147A1