CN110705160A - Airflow simulation calculation method of desulfurization and dust removal integrated desulfurization tower - Google Patents
Airflow simulation calculation method of desulfurization and dust removal integrated desulfurization tower Download PDFInfo
- Publication number
- CN110705160A CN110705160A CN201910920944.3A CN201910920944A CN110705160A CN 110705160 A CN110705160 A CN 110705160A CN 201910920944 A CN201910920944 A CN 201910920944A CN 110705160 A CN110705160 A CN 110705160A
- Authority
- CN
- China
- Prior art keywords
- desulfurization
- simulation
- dust removal
- setting
- removal integrated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Treating Waste Gases (AREA)
Abstract
The invention provides an airflow simulation calculation method of a desulfurization and dust removal integrated desulfurization tower, which comprises the following steps of: according to the structure and the size of the desulfurization and dust removal integrated desulfurization tower, a geometric modeling tool is adopted to establish a numerical simulation geometric model; inputting a geometric model in a grid dividing tool, and dividing a simulation calculation grid; reading in a simulation calculation grid in computational fluid dynamics software, setting boundary conditions of a flue gas inlet and a flue gas outlet, and setting wall boundary conditions; setting a desulfurization unit area and a dust removal unit area in computational fluid dynamics software, wherein the desulfurization unit area is simulated by adopting a porous medium model, and the dust removal unit area is set as a gas passage model; and setting various simulation calculation parameters in computational fluid dynamics software, initializing a calculation domain, and starting airflow simulation calculation until a result is obtained. The invention greatly reduces the computational complexity of the simulation and improves the simulation efficiency on the premise of ensuring the simulation result to be true and credible.
Description
Technical Field
The invention relates to airflow simulation calculation of a desulfurizing tower, in particular to an airflow simulation calculation method of a desulfurizing and dedusting integrated desulfurizing tower.
Background
The wet desulphurization technology has the characteristics of high desulphurization efficiency, no wastewater discharge, secondary utilization of desulfurized products and the like, and is widely applied to flue gas treatment in industries such as electric power, steel and the like. The traditional wet desulphurization tower contains a large amount of particles in the tail gas after wet desulphurization, and a large amount of aerosol is formed, so that the phenomena of white smoke tailing, chimney rain and the like are caused, and the secondary pollution of the surrounding environment is caused. To the above-mentioned problem that traditional wet flue gas desulfurization tower exists, patent CN201721872916.1 provides a integrative desulfurizing tower of desulfurization dust removal, sets up the dust removal unit above the desulfurizing unit, effectively catches the particulate matter after the desulfurization, suppresses the production of aerosol.
The appearance of desulfurization dust removal integral type desulfurizing tower has shown the dust removal effect that has improved wet flue gas desulfurization tower, but in this desulfurizing tower, the series connection of a plurality of functional unit causes the mobile complexity of flue gas to rise, in order to make the flue gas flow in the desulfurizing tower even, better with the desulfurizer contact reaction, more fully utilize the dust removal module, reduce pressure loss in the tower simultaneously, need carry out numerical simulation to the gas flow-mass transfer process of desulfurization dust removal integral type desulfurizing tower to optimize the interior air current of tower.
Disclosure of Invention
The invention aims to provide a gas flow simulation calculation method of a desulfurization and dust removal integrated desulfurization tower, which aims to perform numerical simulation on a gas flow-mass transfer process of the desulfurization and dust removal integrated desulfurization tower so as to optimize gas flow in the tower.
The invention is realized by the following steps:
the invention provides an airflow simulation calculation method of a desulfurization and dust removal integrated desulfurization tower, which comprises the following steps of:
according to the structure and the size of the desulfurization and dust removal integrated desulfurization tower, a geometric modeling tool is adopted to establish a numerical simulation geometric model;
inputting a geometric model in a grid dividing tool, and dividing a simulation calculation grid;
reading in a simulation calculation grid in computational fluid dynamics software, setting boundary conditions of a flue gas inlet and a flue gas outlet, and setting wall boundary conditions;
setting a desulfurization unit area and a dust removal unit area in computational fluid dynamics software, wherein the desulfurization unit area is simulated by adopting a porous medium model, and the dust removal unit area is set as a gas passage model;
and setting various simulation calculation parameters in computational fluid dynamics software, initializing a calculation domain, and starting airflow simulation calculation until a result is obtained.
Further, a porous medium multiphase flow model in ANSYS Fluent software was used to simulate the desulfurization unit and set the porosity, viscous drag coefficient, and inertial loss coefficient.
Further, the porosity is calculated by the following formula:
wherein m is the flow of the desulfurizing agent and the unit is kg/s; v. of0The speed of the liquid drop sprayed by the sprayer is in m/s; s is the cross-sectional area of the desulfurization unit and is expressed in m;
after obtaining the porosity, the viscous resistance coefficient and the inertia loss coefficient C2Calculated using the formula of eurgrun, namely:
wherein D ispThe average diameter of the droplets sprayed by the sprayer is m.
Further, the geometric modeling tool uses one of Gambit, Solidworks, and SpaceClaim.
Further, the mesh partitioning tool uses one of Gambit, Workbench shifting, and Fluent shifting.
Further, the computational fluid dynamics software employs ANSYS Fluent; the flue gas Inlet boundary condition adopts Velocity Inlet, the Outlet boundary condition adopts Pressure Outlet, and the wall surface boundary condition adopts No Slip.
Further, the setting of each simulation calculation parameter in the computational fluid dynamics software specifically includes: setting gravity acceleration, flue gas temperature, density and inlet flow rate; setting a turbulence model; setting the fluid as smoke, and setting the density, viscosity and inlet speed of the smoke; setting a difference format and setting a relaxation factor.
Further, the turbulence model selects a k-epsilon model and a standard wall function.
Further, the difference format is set to a first-order center difference format.
Further, the relaxation factor is set to 0.7-0.8.
Compared with the prior art, the invention has the following beneficial effects:
according to the airflow simulation calculation method of the desulfurization and dust removal integrated desulfurization tower, provided by the invention, the desulfurization unit area is set to be the porous medium model, and the movement of particulate matters which do not affect the airflow much in the dust removal unit area is neglected, so that the calculation complexity of simulation is greatly reduced and the simulation efficiency is improved on the premise of ensuring the simulation result to be real and credible. Aiming at the porous medium model of the desulfurization unit, the invention provides a formula for calculating the porosity according to different desulfurizer flow rates, and further calculates the viscous resistance coefficient and the inertia loss coefficient of the porous medium model under the current working condition by using an Eurgunn formula, thereby changing the current situation of using empirical parameters and obviously improving the simulation precision.
Drawings
FIG. 1 is a geometric simulation model of the gas flow of the desulfurization and dust removal integrated desulfurization tower provided by the embodiment of the invention;
FIG. 2 is a longitudinal center section showing a simulation result of the gas flow of the desulfurization and dust removal integrated desulfurization tower provided in the embodiment of the present invention;
fig. 3 is a cross section of an inlet of a dust removal unit, which is a simulation result of airflow of the desulfurization and dust removal integrated desulfurization tower provided by the embodiment of the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1 to fig. 3, an embodiment of the present invention provides a method for calculating an airflow simulation of a desulfurization and dust removal integrated desulfurization tower, including the following steps:
according to the structure and the size of the desulfurization and dust removal integrated desulfurization tower, a geometric modeling tool is adopted to establish a numerical simulation geometric model; specifically, the geometric modeling tool uses one of Gambit, Solidworks, and SpaceClaim.
Inputting a geometric model in a grid dividing tool, and dividing a simulation calculation grid; specifically, the mesh partitioning tool uses one of Gambit, workbench shifting, and Fluent shifting.
Reading in a simulation calculation grid in computational fluid dynamics software, setting boundary conditions of a flue gas inlet and a flue gas outlet, and setting wall boundary conditions; specifically, the computational fluid dynamics software employs ANSYS Fluent; the flue gas Inlet boundary condition adopts Velocity Inlet, the Outlet boundary condition adopts Pressure Outlet, and the wall surface boundary condition adopts No Slip.
Setting a desulfurization unit area and a dust removal unit area in computational fluid dynamics software, wherein the desulfurization unit area is simulated by adopting a porous medium model, and the dust removal unit area is set as a gas passage model;
and setting various simulation calculation parameters in computational fluid dynamics software, initializing a calculation domain, and starting airflow simulation calculation until a result is obtained.
The principle of the desulfurization unit of the desulfurization and dust removal integrated desulfurization tower is to use a liquid desulfurizing agent such as Ca (OH)2And the ammonia water and other solutions are uniformly sprayed on the cross section of at least one desulfurizing tower through a sprayer, small droplets of a desulfurizing agent solution are uniformly distributed in a certain area, and the area is a desulfurizing unit. When the flue gas passes through the desulfurization unit, the flue gas passes through the gaps among the liquid drops, and SO in the flue gas2Performing gas-liquid two-phase chemical reaction with the desulfurizer to remove SO2. If models of the droplets are respectively established in numerical simulation software and the interaction force between the flue gas and each droplet is calculated, the required grid number and the calculation amount greatly exceed the calculation capability of the current computer. In the numerical simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower, the effect of liquid drops on the gas flow is more concerned, and the effect of the gas flow on the liquid drops can be ignored. Therefore, the liquid drops in the desulfurization unit are regarded as a static solid structure, the flue gas flows from the gaps among the liquid drops, and the porous medium model is adopted to simulate the problem, namely the porous medium model is adopted to simulate the desulfurization unit of the desulfurization and dust removal integrated desulfurization tower, so that the truth and credibility of a simulation result can be ensured, and the calculation complexity of the simulation is greatly reduced.
The principle of the dust removal unit is that a corona wire (cathode) and a particulate matter trapping polar plate (anode) are arranged above the desulfurization unit in the desulfurization tower, and a strong electric field is formed between the corona wire and the particulate matter trapping polar plate to enable the particulate matter to directionally move after being electrified. The particle trapping plate is in a tubular structure which is arranged vertically and horizontally, and the gas flows in and out in the vertical direction. Because the dust removal unit area only has influence on the movement of the particulate matter and can neglect the influence on the airflow, the invention sets the dust removal unit as the gas passage model in the fluid mechanics software, greatly reduces the calculation complexity of the simulation and improves the simulation efficiency on the premise of ensuring the simulation result to be true and credible.
The porous medium model is simulated by adding a source term in a momentum equation, wherein the source term is formed by two parts, namely, the acting force of a solid structure in a flow area is regarded as the distributed resistance added on a fluid: a viscous loss term and an inertial loss term. In the preferred embodiment, a porous medium multiphase flow model in ANSYS Fluent software is used to simulate the desulfurization unit, and porosity, viscous drag coefficient and inertial loss coefficient are set.
For desulfurization units, empirical porosity and drag coefficient parameters are currently used primarily. However, under different desulfurizing agent flow rates, the use of empirical parameters will cause great calculation errors, and the porosity and the resistance coefficient need to be calculated according to actual working conditions.
In the preferred embodiment, the porosity is calculated by the following formula:
wherein m is the flow of the desulfurizing agent and the unit is kg/s; v. of0The speed of the liquid drop sprayed by the sprayer is in m/s; s is the cross-sectional area of the desulfurization unit and is expressed in m;
after obtaining the porosity, the viscous resistance coefficient and the inertia loss coefficient C2Calculated using the formula of eurgrun, namely:
wherein D ispThe average diameter of the droplets sprayed by the sprayer is m.
Through the mode, the viscous resistance coefficient and the inertia loss coefficient of different porous medium models can be determined according to different desulfurizer flow rates, and the simulation precision is remarkably improved.
In this embodiment, the setting of each simulation calculation parameter in the computational fluid dynamics software specifically includes: setting gravity acceleration, flue gas temperature, density and inlet flow rate; setting a turbulence model, specifically selecting a k-epsilon model and a standard wall function; setting the fluid as smoke, and setting the density, viscosity and inlet speed of the smoke; setting a difference format, specifically a first-order central difference format, and setting a relaxation factor, specifically 0.7-0.8.
After the setting of each parameter is finished, initializing a calculation domain in computational fluid dynamics software, and then starting calculation until a result of the airflow simulation calculation is obtained. According to the simulation calculation result, the flue gas flow field in the tower can be analyzed.
According to the airflow simulation calculation method of the desulfurization and dust removal integrated desulfurization tower provided by the embodiment of the invention, the desulfurization unit area is set to be the porous medium model, and the movement of particulate matters which do not affect the airflow much in the dust removal unit area is neglected, so that the calculation complexity of simulation is greatly reduced and the simulation efficiency is improved on the premise of ensuring the simulation result to be true and credible. Aiming at the porous medium model of the desulfurization unit, the invention provides a formula for calculating the porosity according to different desulfurizer flow rates, and further calculates the viscous resistance coefficient and the inertia loss coefficient of the porous medium model under the current working condition by using an Eurgunn formula, thereby changing the current situation of using empirical parameters and obviously improving the simulation precision.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (10)
1. A method for simulating and calculating airflow of a desulfurization and dust removal integrated desulfurization tower is characterized by comprising the following steps of:
according to the structure and the size of the desulfurization and dust removal integrated desulfurization tower, a geometric modeling tool is adopted to establish a numerical simulation geometric model;
inputting a geometric model in a grid dividing tool, and dividing a simulation calculation grid;
reading in a simulation calculation grid in computational fluid dynamics software, setting boundary conditions of a flue gas inlet and a flue gas outlet, and setting wall boundary conditions;
setting a desulfurization unit area and a dust removal unit area in computational fluid dynamics software, wherein the desulfurization unit area is simulated by adopting a porous medium model, and the dust removal unit area is set as a gas passage model;
and setting various simulation calculation parameters in computational fluid dynamics software, initializing a calculation domain, and starting airflow simulation calculation until a result is obtained.
2. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 1, wherein: and (3) simulating a desulfurization unit by using a porous medium multiphase flow model in ANSYS Fluent software, and setting porosity, viscous resistance coefficient and inertia loss coefficient.
3. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 2, wherein: the porosity is calculated by the following formula:
wherein m is the flow of the desulfurizing agent and the unit is kg/s; v. of0The speed of the liquid drop sprayed by the sprayer is in m/s; s is the cross-sectional area of the desulfurization unit and is expressed in m;
after obtaining the porosity, the viscous resistance coefficient and the inertia loss coefficient C2Calculated using the formula of eurgrun, namely:
wherein D ispThe average diameter of the droplets sprayed by the sprayer is m.
4. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 1, wherein: the geometric modeling tool uses one of Gambit, Solidworks and SpaceClaim.
5. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 1, wherein: the grid dividing tool uses one of Gambit, Workbench Meshing and Fluent Meshing.
6. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 1, wherein: the computational fluid dynamics software adopts ANSYS Fluent; the flue gas Inlet boundary condition adopts Velocity Inlet, the Outlet boundary condition adopts Pressure Outlet, and the wall surface boundary condition adopts No Slip.
7. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 1, wherein the setting of the simulation calculation parameters in the computational fluid dynamics software specifically comprises: setting gravity acceleration, flue gas temperature, density and inlet flow rate; setting a turbulence model; setting the fluid as smoke, and setting the density, viscosity and inlet speed of the smoke; setting a difference format and setting a relaxation factor.
8. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 7, wherein: the turbulence model selects a k-epsilon model and a standard wall function.
9. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 7, wherein: and setting the difference format as a first-order central difference format.
10. The method for calculating the simulation of the gas flow of the desulfurization and dust removal integrated desulfurization tower as recited in claim 7, wherein: the relaxation factor is set to 0.7-0.8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910920944.3A CN110705160A (en) | 2019-09-27 | 2019-09-27 | Airflow simulation calculation method of desulfurization and dust removal integrated desulfurization tower |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910920944.3A CN110705160A (en) | 2019-09-27 | 2019-09-27 | Airflow simulation calculation method of desulfurization and dust removal integrated desulfurization tower |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110705160A true CN110705160A (en) | 2020-01-17 |
Family
ID=69197639
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910920944.3A Pending CN110705160A (en) | 2019-09-27 | 2019-09-27 | Airflow simulation calculation method of desulfurization and dust removal integrated desulfurization tower |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110705160A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112966452A (en) * | 2021-03-11 | 2021-06-15 | 机械工业第九设计研究院有限公司 | Method for simulating airflow speed and pressure intensity of coating spray booth |
| CN114491751A (en) * | 2022-01-20 | 2022-05-13 | 三一重型装备有限公司 | Flow field analysis system, flow field analysis method, and computer-readable storage medium |
| CN117634340A (en) * | 2023-11-20 | 2024-03-01 | 北京科技大学 | Determination method for desulfurization effect of bottom argon blowing ladle |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0923978A1 (en) * | 1997-12-22 | 1999-06-23 | Mitsubishi Heavy Industries, Ltd. | Flue gas treating process |
| CN104636552A (en) * | 2015-02-06 | 2015-05-20 | 江苏迪萨机械有限公司 | Numerical simulation computing method based on ANSYS CFX software and used for desulphurization process in amidogen wet desulphurization spray single tower |
| CN106000007A (en) * | 2016-07-06 | 2016-10-12 | 福建龙净环保股份有限公司 | Simulation system and method for wet desulphurization flow fields |
| CN207755804U (en) * | 2017-12-27 | 2018-08-24 | 武汉都市环保工程技术股份有限公司 | A kind of desulfurizing tower |
-
2019
- 2019-09-27 CN CN201910920944.3A patent/CN110705160A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0923978A1 (en) * | 1997-12-22 | 1999-06-23 | Mitsubishi Heavy Industries, Ltd. | Flue gas treating process |
| CN104636552A (en) * | 2015-02-06 | 2015-05-20 | 江苏迪萨机械有限公司 | Numerical simulation computing method based on ANSYS CFX software and used for desulphurization process in amidogen wet desulphurization spray single tower |
| CN106000007A (en) * | 2016-07-06 | 2016-10-12 | 福建龙净环保股份有限公司 | Simulation system and method for wet desulphurization flow fields |
| CN207755804U (en) * | 2017-12-27 | 2018-08-24 | 武汉都市环保工程技术股份有限公司 | A kind of desulfurizing tower |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112966452A (en) * | 2021-03-11 | 2021-06-15 | 机械工业第九设计研究院有限公司 | Method for simulating airflow speed and pressure intensity of coating spray booth |
| CN114491751A (en) * | 2022-01-20 | 2022-05-13 | 三一重型装备有限公司 | Flow field analysis system, flow field analysis method, and computer-readable storage medium |
| CN117634340A (en) * | 2023-11-20 | 2024-03-01 | 北京科技大学 | Determination method for desulfurization effect of bottom argon blowing ladle |
| CN117634340B (en) * | 2023-11-20 | 2024-05-24 | 北京科技大学 | Determination method for desulfurization effect of bottom argon blowing ladle |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110705160A (en) | Airflow simulation calculation method of desulfurization and dust removal integrated desulfurization tower | |
| Chen et al. | Experimental and numerical study on effects of deflectors on flow field distribution and desulfurization efficiency in spray towers | |
| CN104636552A (en) | Numerical simulation computing method based on ANSYS CFX software and used for desulphurization process in amidogen wet desulphurization spray single tower | |
| Narimani et al. | Optimization of vane mist eliminators | |
| Wang et al. | Numerical simulation research of flow field in ammonia-based wet flue gas desulfurization tower | |
| CN112069689B (en) | Simulation method and system for fuel atomization characteristic of aircraft engine | |
| CN105138718B (en) | A kind of projectional technique and its subsidiary engine load adjustment method of combination UDF Novel fires desulfurizing tower of electric power plant desulfuration efficiency | |
| CN105413442A (en) | Vertical desulfurization tower | |
| CN113094837A (en) | Wind resistance design method of horizontal axis wind turbine blade under strong wind action | |
| CN203507792U (en) | Synchronous desulfurization and denitrification tower | |
| CN113360979B (en) | Intelligent control system for vortex suppression and dehumidification integrated flow of steel box girder and application method | |
| Che et al. | Effect of adjusted mesoscale drag model on flue gas desulfurization in powder-particle spouted beds | |
| Guan et al. | Numerical study on the penetration of droplets in a zigzag demister | |
| CN116440689B (en) | High-efficiency injection flow equalization desulfurization spray tower | |
| Sayem et al. | Effects of baffles on flow distribution in an electrostatic precipitator (ESP) of a coal based power plant | |
| Liang et al. | Study on the performance of a vortex mist eliminator for improving the separation efficiency of fine particles | |
| Huang et al. | Particle-Scavenging prediction in sieve plate scrubber via dimension reduction in computational fluid dynamics | |
| CN108197386A (en) | Manifold clarifier structural optimization method based on CFD emulation | |
| CN107321157A (en) | A kind of horizontal wet desulphurization device | |
| CN107256283B (en) | High-precision analysis method for oil mixing characteristics in fuel tank of aircraft | |
| Kaesemann et al. | The Meaning of Droplet‐Droplet Interaction for the Wet Flue‐Gas Cleaning Process | |
| CN104239730B (en) | Non-Gaussian turbulent flow analogy method based on Lagrange random particulate matter model | |
| Xie | Hydrodynamics and mass transfer of rotating packed beds for CO2 capture | |
| Valger | Modeling solid particle transport and air flow around obstacle | |
| CN213995360U (en) | Structure of bottle-shaped adsorber in asphalt flue gas purification device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200117 |