CN110147594A - Probe into the analogy method of FGD by spraying mist tower interior flow field distribution - Google Patents
Probe into the analogy method of FGD by spraying mist tower interior flow field distribution Download PDFInfo
- Publication number
- CN110147594A CN110147594A CN201910384251.7A CN201910384251A CN110147594A CN 110147594 A CN110147594 A CN 110147594A CN 201910384251 A CN201910384251 A CN 201910384251A CN 110147594 A CN110147594 A CN 110147594A
- Authority
- CN
- China
- Prior art keywords
- fgd
- spraying mist
- tower
- mist tower
- flow field
- 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
- 239000003595 mist Substances 0.000 title claims abstract description 48
- 238000005507 spraying Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000009826 distribution Methods 0.000 title claims abstract description 28
- 239000000523 sample Substances 0.000 title description 7
- 230000003009 desulfurizing effect Effects 0.000 claims abstract description 21
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims abstract description 17
- 239000003546 flue gas Substances 0.000 claims abstract description 16
- 239000012530 fluid Substances 0.000 claims abstract description 14
- 239000002245 particle Substances 0.000 claims abstract description 11
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000010586 diagram Methods 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 10
- 239000000779 smoke Substances 0.000 claims description 5
- 238000007789 sealing Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 11
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 238000005457 optimization Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 4
- 239000000571 coke Substances 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 14
- 238000006477 desulfuration reaction Methods 0.000 description 9
- 239000012071 phase Substances 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 239000007789 gas Substances 0.000 description 6
- 230000023556 desulfurization Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000002002 slurry Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 239000007921 spray Substances 0.000 description 3
- RAHZWNYVWXNFOC-UHFFFAOYSA-N sulfur dioxide Inorganic materials O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 239000005864 Sulphur Substances 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- UBAZGMLMVVQSCD-UHFFFAOYSA-N carbon dioxide;molecular oxygen Chemical compound O=O.O=C=O UBAZGMLMVVQSCD-UHFFFAOYSA-N 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000032258 transport Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/14—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by absorption
- B01D53/1456—Removing acid components
- B01D53/1481—Removing sulfur dioxide or sulfur trioxide
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/14—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by absorption
- B01D53/18—Absorbing units; Liquid distributors therefor
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/14—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols by absorption
- B01D53/18—Absorbing units; Liquid distributors therefor
- B01D53/185—Liquid distributors
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/34—Chemical or biological purification of waste gases
- B01D53/46—Removing components of defined structure
- B01D53/48—Sulfur compounds
- B01D53/50—Sulfur oxides
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/34—Chemical or biological purification of waste gases
- B01D53/74—General processes for purification of waste gases; Apparatus or devices specially adapted therefor
- B01D53/77—Liquid phase processes
- B01D53/78—Liquid phase processes with gas-liquid contact
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D2258/00—Sources of waste gases
- B01D2258/02—Other waste gases
- B01D2258/0283—Flue gases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Environmental & Geological Engineering (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of analogy methods for probing into the distribution of FGD by spraying mist tower interior flow field, the division of determination and computational domain grid file first to the foundation for the FGD by spraying mist tower physical model probed into, solution domain: selection solver, computational domain grid file is imported to selected solver, it solves the problems, such as to select corresponding mathematical model as needed, given primary condition and boundary condition, carry out Physical Quantity Calculation in required FGD by spraying mist tower；Solving result is finally shown as cloud atlas and scatterplot diagram form, the clear Flow Field Distribution situation that must be understood in desulfurizing tower.This analogy method is simulated calculating to flow process of the coke oven flue gas in FGD by spraying mist tower, the intracorporal Flow Field Distribution of FGD by spraying mist tower tower is accurately, easily obtained, such as distribution of VELOCITY DISTRIBUTION, pressure and discrete phase distribution of particles, fluid distrbution state under the influence of different factors can industrial design to FGD by spraying mist tower and optimization industrial process conditions provide reference guide to make enterprise preferably reduce production cost.
Description
Technical field
The present invention relates to a kind of desulfurization technology of coke oven tail gas, in particular to one kind probes into FGD by spraying mist tower interior flow field point
The analogy method of cloth.
Background technique
Semidry method FGD by spraying mist technology is the more mature skill that can be applied with market few in number in existing desulfurization technology
One of art, a common feature of Summary of Semi-dry Flue Gas Desulfurization are the slurry drop and cigarette that desulfurizer slurry is atomized into certain partial size
Gas conducted heat, mass transfer and physics and chemical reaction, by SO2It is fixed in the form of solid-state sulphur in slurry drop, slurry drop is heated dry
It is dry at after powdered, enter deduster with flue gas.The Flow Field Distribution situation clearly grasped in desulfurizing tower first is to understand liquid phase
The most important thing of mass transfer, heat transfer theory and theory analysis, but due to tower body equipment is huge, biggish flue gas flow and higher
The influence of many factors such as flue-gas temperature causes the Flow Field Distribution situation wanted in clear grasp desulfurizing tower extremely difficult, theoretical
Analysis is it is found that be extremely complex Three dimensional Turbulent flow field in tower, therefore need by hydrodynamic method and fluent software
The physical fields such as velocity field, pressure field inside desulfurizing tower are carried out calculating to be further analysis desulfuration efficiency and flow field change
Rule lay the foundation.
A large amount of related data display technologies are mature, being widely used in a whole set of desulfurizer device of comercial operation, not only equipment is huge
Greatly, facility occupied area is wide and pipeline arrangement complex operations are not easy.Interior flow field research is carried out not using these huge tower bodies
But it brings various inconvenient for the test job of engineer and also results in serious waste of resources, research cost increases, flue gas
The appearance for the problems such as parameter is not easy to control, therefore, there is an urgent need to a kind of new technical solutions qualitative can estimate operating parameter
The reasonability of change.
Computational fluid dynamics (Computational Fluid Dynamics, abbreviation CFD) is hydromechanical one
Subdiscipline, it is using electronic computer as tool, using the mathematical method of various discretizations, all kinds of problems of Fluid Mechanics into
Row numerical experiment, computer mould, which fit analysis and research, transports a kind of engineering of being optimal of result to solve various practical problems
Calculation method.With the rapid development of computer technology, CFD technology has become the important means of research Three-dimensional Flow, has both at home and abroad
A variety of method for numerical simulation can simulate the flow field of experimental provision, compare from calculated result and experimental result, numerical value
No matter advantageous on time and precision simulation is, and this method has network analysis quantification, performance analysis mobilism, quality point
Analyse the main features such as reliability, design result optimization, design process high efficiency.
Therefore, based on the above thinking analysis, the fluent software simulation computing platform based on CFD technology is in heavy construction
Application just come into being.After establishing a kind of analogy method of FGD by spraying mist tower interior flow field analysis on the platform, do not need
Actually do any experiment can be obtained Flow Field Distribution situation ancillary works Shi Shixian inside tower body to FGD by spraying mist tower structure and
The optimization of Process operating parameters etc. finally improves desulfuration efficiency and saves operating cost.
Find that few similar patents described the interior flow field point of this kind of spray column by the retrieval to existing patent
Cloth situation, especially for the FGD by spraying mist tower handled for coke oven tail gas.It is inhaled since the sweetening process is not only related to physics
Receipts further relate to chemical reaction, significantly affect so Flow Field Distribution (variables such as speed, vacuum degree) can all have desulfuration efficiency,
So the fundamental physical quantity how conveniently to obtain in FGD by spraying mist tower on each position is particularly important.
Summary of the invention
The present invention be directed to the problems that the analysis of the interior flow field of desulfurizing tower is difficult, propose one kind and probe into FGD by spraying mist tower
The analogy method of portion's Flow Field Distribution, finite volume method is applied in the analysis of this simulation process.It, can be effective based on this method
Ground removes the dependence in conventional method for experimental facilities, makes sightless flow field change visualization in tower, and pass through change phase
Parameter is closed to reach optimal design effect, not only accelerates housing structure and its Process operating parameters optimization progress, Er Qieke
To reduce operating cost.
The technical solution of the present invention is as follows: a kind of analogy method for probing into the distribution of FGD by spraying mist tower interior flow field, specifically includes
Following steps:
1) foundation, the determination in solution domain and the division of computational domain grid file of FGD by spraying mist tower physical model:
Zero influenced according to the structure for the FGD by spraying mist tower probed into and equipment size and FGD by spraying mist tower internal convection field
Part carries out simplification and establishes FGD by spraying mist tower physical model；
The FGD by spraying mist tower physical model of foundation is imported into Meshing software and determines simulation computational domain, the simulation is asked
Solution domain is determined as smoke entrance and tower body outer surface is the sealing fluid domain on boundary；
It is divided on the fluid domain using structured grid；
To boundary part names and derived grid file, into solution；
2) solution procedure:
Solver is selected, grid file in step 1) is imported to selected solver, solves the problems, such as selection pair as needed
Mathematical model is answered, primary condition and boundary condition is given, carries out Physical Quantity Calculation in required FGD by spraying mist tower；
3) step 2) solving result is shown as cloud atlas and scatterplot diagram form, the clear Flow Field Distribution that must be understood in desulfurizing tower
Situation.
Step 2) the specific steps:
2.1) after grid file, to be imported to selected solver, further check to grid file first ensures nothing
Negative volume occurs and modifying dimension scale keeps grid consistent with the Unit Scale in computational domain；
2.2) it, selects mathematical model and primary condition is set；
2.3), the inlet flue gas property according to used in the FGD by spraying mist tower probed into setting simulation used in flue gas group be grouped as and
Temperature, vacuum degree and initial velocity parameters in series；
2.4), according to probe into FGD by spraying mist tower actual conditions setting solve calculate control parameter, setting solve format, from
It dissipates format, physical time scale, particle under-relaxation factor and the condition of convergence and activates monitor；
2.5) it, initializes flow field and completes iterative solution and calculate, obtain the fundamental physical quantity in required FGD by spraying mist tower.
The mathematical model is by energy model, Standard k-epsilon turbulence model, component shipping model and discrete
Phase model composition.
The beneficial effects of the present invention are: the present invention probes into the analogy method of FGD by spraying mist tower interior flow field distribution, not only
It provides fundamental basis for industrial operation and optimizes reference, the design and exploitation of desulfurizing tower can also be promoted；And the method for the present invention
The disadvantages of engineering test investment is big, the period is long is overcome, so that the acquisition of data can be with the limitation in ablation experiment place, to industry
Operation instruct and then saves production and operating cost, therefore has a good application prospect and economic value.
Detailed description of the invention
Fig. 1 is the three-dimensional physical model figure of inventive desulfurization tower；
Fig. 2 is the speed cloud atlas under present invention specific implementation operating condition；
Fig. 3 is the vacuum degree cloud atlas under present invention specific implementation operating condition；
Fig. 4 is the drop of the invention being embodied under operating condition in dwell time in the tower figure；
Fig. 5 is the drop being embodied under operating condition temperature variation in tower of the invention.
Specific embodiment
Invention provides the completely new visual angle of one kind for the interior flow field analysis of desulfurizing tower and is to reduce desulfuration efficiency to make to pollute
Object reaches discharge standard and lays a good foundation, thus for it is industrial reduce operating cost and optimization desulfurizing tower provide reference frame and
Guidance.
The present invention is illustrated below in conjunction with drawings and examples.
Probe into the analogy method of FGD by spraying mist tower interior flow field distribution, comprising the following steps:
1, pretreatment process, the pretreatment process mainly include the foundation of FGD by spraying mist tower physical model, solve domain
Determining and computational domain grid file division；
The foundation of the FGD by spraying mist tower physical model utilizes the completion of DesignModeler software, computational domain grid file
Division utilize Meshing Software on Drawing, the specific steps are as follows:
1) part that the structure and equipment size of the FGD by spraying mist tower, probed into and FGD by spraying mist tower internal convection field influence
Its geometrical model is determined after reasonably simplifying, and FGD by spraying mist tower object is made according to existing size using DesignModeler software
Manage model；
2) FGD by spraying mist tower physical model, is imported into Meshing software and determines simulation computational domain, the analog approach domain
It is determined as smoke entrance and tower body outer surface is the sealing fluid domain on boundary；
3) (structured grid includes hexahedron structure grid), is divided using structured grid on the fluid domain,
Consider that mesh quality and calculator memory situation determine hexahedral mesh number simultaneously；
4), to boundary part names and derived grid file, into solution procedure.
2, solution procedure, the solution procedure mainly include mathematical model, primary condition and the boundary that setting needs to solve
Condition and calculating solve；
Fluent software can be used as solver in the solution procedure, the specific steps are as follows:
1) it, will have been exported in step 1 the 4) step after the grid file kept imports fluent software, first to grid text
Part, which carries out further inspection, which to be ensured to occur and modify dimension scale without negative volume, keeps grid consistent with the Unit Scale in computational domain；
2) it, selects mathematical model and primary condition is set, mathematical model includes activation energy equation, standard k- ε turbulent flow mould
Type；
3), the group that flue gas used is simulated in the setting of the inlet flue gas property according to used in the FGD by spraying mist tower probed into is grouped as and warm
The series of parameters such as degree, vacuum degree and initial velocity；
4) it, is solved according to the FGD by spraying mist tower actual conditions setting probed into and calculates control parameter, setting solves format, discrete
Format, physical time scale, particle under-relaxation factor and the condition of convergence simultaneously activate monitor；
5) it, initializes flow field and completes iterative solution and calculate, obtain the fundamental physical quantity in required FGD by spraying mist tower.
3, solving result is predominantly shown as cloud atlas and scatterplot diagram form by last handling process, the last handling process；Institute
State last handling process method particularly includes: post-processed using CFD Post software and the data that simulation obtains are exported and shown
It can effectively observe with analysis mode calculating for figure or curve negotiating post-processing as a result, clear must understand in desulfurizing tower
Flow Field Distribution situation.
In order to preferably more intuitively recognize the Flow Field Distribution situation inside desulfurizing tower, need according to certain factory's actual production
Data carry out suitably simplified and analog study as initial condition and to simulated object on the basis of simulating purpose.Master of the present invention
The three-dimensional physical model of desulfurizing tower is established using DesignModeler software, mainly by gas inlet section, exhanst gas outlet section,
The part such as middle part tower body, bottom cone and nozzle forms, and following table is the major parameter of desulfurizing tower.Fig. 1 is the three-dimensional article of desulfurizing tower
Manage illustraton of model (using tower bottom as origin, tower vertical direction is that 3D model is established in z-axis direction).
Because Structure of Flue Gas Desulfuration Absorbing in the present invention is more regular simple, so using Meshing software to physical model into
The operation of row grid dividing, is divided into 40000 after comprehensively considering the factors such as mesh quality, number of grid and Computing ability
The higher hexahedron structure grid of quality simultaneously exports preservation.
Derived grid file is imported fluent software setting mathematical model parameter and is solved.Primary study of the present invention is steady
The Flow Field Distribution situation of liquid phase under state, so selecting Pressure-based and Steady type in this specific embodiment
Fluent solver.According to selection energy model, Standard k-epsilon turbulence model, group after the interpretation to relevant issues
Divide basis of the mathematical models such as shipping model and Discrete Phase Model as subsequent solution.Wherein Standard k-epsilon turbulent flow
Model is proposed by Launder and Spalding, because of its own stability, economical and higher computational accuracy for having
Making application range in turbulence model most extensively, is also most a model known to people.Standard k-ε model is by solving turbulent flow
Kinetic energy equation (k) and turbulence vortex Dissipative Equation (ε) obtain the solution of k and ε, then calculate turbulent viscosity with the value of k and ε again, finally
Assume to obtain eddy stress solution by Boussinesq, wherein k- ε equation is as follows:
Wherein: ρ is turbulent flow density, kgm-3；T indicates the time；uiFor i-th of speed of fluid, ms-1；xiFor fluid
I-th of locality, xiAnd xjAll refer to space a direction, works as xiAnd xjPhase time equivalence table space same direction, works as xiWith
xjUnequal epoch table space different directions；K is Turbulent Kinetic, m2·s-1；μ is turbulent viscosity, Pas；ε is Turbulent Kinetic consumption
The rate of dissipating, m2·s-1；CμFor empirical；GkThe tubulence energy generated for laminar velocity gradient, kgm-1·s-1；GbFor buoyancy production
Raw tubulence energy, kgm-1·s-1；YMFor the fluctuation that transition diffusion generates in compressible turbulent flow, kgm-1·s-3；C1ε、C2ε、
C3εIt is constant, C1ε=1.44, C2ε=1.92, C3ε=0.09, σε=1.3, σk=1.0.
The theoretical thought of Discrete Phase Model be by means of Lagrangian method based on studying single particle movement process
Record the particle during the motion physical quantity change with time rule and comprehensive all particles movements with constitute entirely from
The track of dephasing, the model can be indicated by following equation:
Wherein: gxIt is the acceleration of gravity on the direction x；FxFor additional acceleration item；FD(u-up) it is per unit mass particle
Resistance；U is fluid phase velocity, ms-1；upFor particle speed, ms-1；μ is turbulent viscosity, Pas；ρ is turbulent flow density,
kg·m-3；ρpFor the density of particle, kgm-3；dpFor particle diameter, m；CDIt is drag coefficient；ReFor Reynolds (Reynold's) criterion.
The setting of boundary condition
(1) smoke inlet condition:
Flue gas used is made of vapor, oxygen, carbon dioxide, sulfur dioxide and nitrogen in the present invention, wherein each component
Mass fraction be respectively 0.13,0.11,0.1,0.005,0.655.Smoke inlet speed chooses 10m/s, flue-gas temperature 165
DEG C, less turbulence 5%, hydraulic radius 3.43m, entrance is escape to the effect of discrete phase.
(2) exhanst gas outlet condition:
Back pressure is atmospheric pressure, and flue-gas temperature is 155 DEG C, less turbulence 5%, hydraulic radius 3.35m, outlet section
Effect to discrete phase is escape.
(3) wall condition:
This specific embodiment is set static wall surface for wall surface and is thought fluid on wall surface using Stationary Wall
Without sliding, fixed boundary temperature is 165 DEG C, ignores wall surface to the radiation event of external environment.
(4) discrete phase parameter:
Using the vertex of bottom cone as coordinate origin, Solid-cone class is selected according to the manufacturing parameter present invention of certain factory
It the nozzle of type and is placed at (0,0,30) coordinate, spray droplet size distribution meets uniform distribution, and spray liquid temperature is 40 DEG C
It is sprayed by Z axis negative direction (i.e. identical as flow of flue gas direction) with the speed of 20m/s and jet angle is 45 °, flow 4kg/s,
Spout radius is 0.02 and average grain diameter is 0.002m.
(5) control parameter is solved:
Usually can be used in the solver four kinds of i.e. SIMPLE, SIMPLEC of pressure-velocity form correlation,
PISO, coupled select SIMPLE as method for solving because the present invention is permanent stream calculation；In addition the present invention will fit
When reducing relaxation factor to guarantee to calculate and faster and better reach convergence, relaxation factor is set as 0.7 in the case, the condition of convergence
It is set as 1e-5, initial method is set as Hybrid initialization.
Using CFD Post the poster processing soft to above-mentioned calculated result carry out processing and will simulation obtain data output and
The form of figure or scatter plot is shown as intuitively to observe and result that analysis mode calculates and clearly to understand desulfurizing tower
Interior Flow Field Distribution situation.
Specific in conjunction with following example, the present invention will be further described: excellent to simulate inside the desulfurizing tower of certain steel enterprise
Change analysis, we can see that the VELOCITY DISTRIBUTION situation in desulfurizing tower inside center section from Fig. 2, on the section XZ only in entrance and
Exit velocity is larger, and the low regime of 1m/s is in the middle part of remaining most of region especially tower body.We can be with from Fig. 3
Find out the vacuum degree distribution situation on the section desulfurizing tower inside center XZ, vacuum degree remains basically stable in tower, only in entrance and exit
Etc. sub-fractions region in changing by a small margin.Fig. 4 be in the case of desulfurizing tower drop temperature is 40 DEG C its in the tower internal stops time
Figure, can intuitively find out drop in dwell time in the tower at 10 seconds or so.Fig. 5 is that desulfurizing tower drop temperature is 60 DEG C of situations
Its lower temperature variation in tower, it can be seen that drop temperature at nozzle is in minimum temperature, with the increasing of residence time
Add, drop temperature constantly increases until stablizing at 61 DEG C or so.
The above is only the citing of embodiments of the present invention, it is noted that for the ordinary skill of the art
For personnel, without departing from the technical principles of the invention, several improvements and modifications can also be made, these improve and become
Type also should be regarded as protection scope of the present invention.
Claims (3)
1. a kind of analogy method for probing into the distribution of FGD by spraying mist tower interior flow field, which is characterized in that specifically comprise the following steps:
1) foundation, the determination in solution domain and the division of computational domain grid file of FGD by spraying mist tower physical model:
The part influenced according to the structure for the FGD by spraying mist tower probed into and equipment size and FGD by spraying mist tower internal convection field into
FGD by spraying mist tower physical model is established in row simplification；
The FGD by spraying mist tower physical model of foundation is imported into Meshing software and determines simulation computational domain, the analog approach domain
It is determined as smoke entrance and tower body outer surface is the sealing fluid domain on boundary；
It is divided on the fluid domain using structured grid；
To boundary part names and derived grid file, into solution；
2) solution procedure:
Solver is selected, grid file in step 1) is imported to selected solver, solves the problems, such as to select corresponding number as needed
Model is learned, primary condition and boundary condition is given, carries out Physical Quantity Calculation in required FGD by spraying mist tower；
3) step 2) solving result is shown as cloud atlas and scatterplot diagram form, the clear Flow Field Distribution feelings that must be understood in desulfurizing tower
Condition.
2. probing into the analogy method of FGD by spraying mist tower interior flow field distribution according to claim 1, which is characterized in that the step
Rapid 2) specific steps:
2.1) after grid file, to be imported to selected solver, further check to grid file first ensures without negative body
Product occurs and modifying dimension scale keeps grid consistent with the Unit Scale in computational domain；
2.2) it, selects mathematical model and primary condition is set；
2.3), the group that flue gas used is simulated in the setting of the inlet flue gas property according to used in the FGD by spraying mist tower probed into is grouped as and warm
Degree, vacuum degree and initial velocity parameters in series；
2.4) it, is solved according to the FGD by spraying mist tower actual conditions setting probed into and calculates control parameter, setting solves format, discrete lattice
Formula, physical time scale, particle under-relaxation factor and the condition of convergence simultaneously activate monitor；
2.5) it, initializes flow field and completes iterative solution and calculate, obtain the fundamental physical quantity in required FGD by spraying mist tower.
3. probing into the analogy method of FGD by spraying mist tower interior flow field distribution according to claim 2, which is characterized in that the number
Model is learned to be made of energy model, Standard k-epsilon turbulence model, component shipping model and Discrete Phase Model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910384251.7A CN110147594A (en) | 2019-05-09 | 2019-05-09 | Probe into the analogy method of FGD by spraying mist tower interior flow field distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910384251.7A CN110147594A (en) | 2019-05-09 | 2019-05-09 | Probe into the analogy method of FGD by spraying mist tower interior flow field distribution |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110147594A true CN110147594A (en) | 2019-08-20 |
Family
ID=67594132
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910384251.7A Pending CN110147594A (en) | 2019-05-09 | 2019-05-09 | Probe into the analogy method of FGD by spraying mist tower interior flow field distribution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110147594A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110929460A (en) * | 2019-11-29 | 2020-03-27 | 福建龙净环保股份有限公司 | Flow field optimization method and system for flue evaporation process with zero discharge of desulfurization wastewater |
| CN113868876A (en) * | 2021-09-30 | 2021-12-31 | 湖南工商大学 | Integrated design method of flue gas pollutant control equipment |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
-
2019
- 2019-05-09 CN CN201910384251.7A patent/CN110147594A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Non-Patent Citations (1)
| Title |
|---|
| 房慧德: "旋转喷雾干燥法脱硫塔优化模拟研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110929460A (en) * | 2019-11-29 | 2020-03-27 | 福建龙净环保股份有限公司 | Flow field optimization method and system for flue evaporation process with zero discharge of desulfurization wastewater |
| CN110929460B (en) * | 2019-11-29 | 2023-07-14 | 福建龙净环保股份有限公司 | Flow field optimization method and system for desulfurization wastewater zero-emission flue evaporation process |
| CN113868876A (en) * | 2021-09-30 | 2021-12-31 | 湖南工商大学 | Integrated design method of flue gas pollutant control equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chauhan et al. | The turbulent/non-turbulent interface and entrainment in a boundary layer | |
| Zeinali et al. | Turbulent flow over partially superhydrophobic underwater structures: The case of flow over sphere and step | |
| Cattafesta III et al. | Three-dimensional boundary-layer transition on a swept wing at Mach 3.5 | |
| CN102880734A (en) | Airplane tail jet flow atmospheric diffusion modeling method based on CFD (computational fluid dynamics) | |
| CN110489777A (en) | A kind of array hole electric spark ultrasonic Compound Machining flow field simulation method | |
| CN110147594A (en) | Probe into the analogy method of FGD by spraying mist tower interior flow field distribution | |
| CN112069689A (en) | Simulation method and system for fuel atomization characteristic of aircraft engine | |
| Drikakis et al. | Azure: An advanced cfd software suite based on high-resolution and high-order methods | |
| Ibrahim et al. | Nonlinear breakup model for a liquid sheet emanating from a pressure-swirl atomizer | |
| Haider et al. | Mathematical analysis of flow passing through a rectangular nozzle | |
| Bondarev et al. | Parametric study of the accuracy of OpenFOAM solvers for the oblique shock wave problem | |
| CN109583131B (en) | Optimization design method for surface microstructure size parameter with drag reduction effect | |
| CN115310325A (en) | Multi-disciplinary coupling analysis framework and method for air refueling of hose taper sleeve | |
| Sakri et al. | Numerical simulation of liquids draining from a tank using OpenFOAM | |
| Ghahremanian et al. | Numerical and experimental verification of initial, transitional and turbulent regions of free turbulent round jet | |
| Njobuenwu et al. | Coupled RANS–LPT modelling of dilute, particle-laden flow in a duct with a 90 bend | |
| Kim | Direct numerical simulation of two-phase flow with application to air layer drag reduction | |
| Chernyshev et al. | Research of the flow-over of the mechanical core of four pipes by a wind current, using Control Volume Method | |
| CN106547968A (en) | The nearly wall escape phenomenon quantization method in absorption tower | |
| Gutović et al. | CFD Analysis of Pressure Losses in Flat-Oval Duct Fittings. | |
| Namazian et al. | Numerical simulation of particle-gas flow through a fixed pipe, using one-way and two-way coupling methods | |
| Kim et al. | Accuracy improvement of the bleed boundary condition with the effects of porosity variations and expansion waves | |
| Yeh | Numerical investigation of liquid jet emanating from plain-orifice atomizers with chamfered or rounded orifice inlets | |
| Ghods | Detailed numerical simulation of liquid jet in crossflow atomization with high density ratios | |
| Liu et al. | Research on the Similarity Scale of Flood Discharge Atomization Based on Water‒Air Two-Phase Flow. Water 2023, 15, 442 |
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 |
Application publication date: 20190820 |
|
| RJ01 | Rejection of invention patent application after publication |