CN112435716A - Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace - Google Patents
Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace Download PDFInfo
- Publication number
- CN112435716A CN112435716A CN202011275672.5A CN202011275672A CN112435716A CN 112435716 A CN112435716 A CN 112435716A CN 202011275672 A CN202011275672 A CN 202011275672A CN 112435716 A CN112435716 A CN 112435716A
- Authority
- CN
- China
- Prior art keywords
- furnace
- oxygen
- setting
- temperature carbonization
- inlet
- 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.)
- Granted
Links
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 title claims abstract description 106
- 239000001301 oxygen Substances 0.000 title claims abstract description 106
- 229910052760 oxygen Inorganic materials 0.000 title claims abstract description 106
- 238000003763 carbonization Methods 0.000 title claims abstract description 54
- 238000004088 simulation Methods 0.000 title claims abstract description 54
- 238000009826 distribution Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000000007 visual effect Effects 0.000 title claims abstract description 8
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 95
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 39
- 238000012544 monitoring process Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 239000007789 gas Substances 0.000 claims abstract description 9
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 238000007789 sealing Methods 0.000 claims abstract description 5
- 229910001873 dinitrogen Inorganic materials 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 13
- 239000012530 fluid Substances 0.000 claims description 13
- 238000009792 diffusion process Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 239000013589 supplement Substances 0.000 claims description 7
- 230000000903 blocking effect Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 238000004134 energy conservation Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 238000012800 visualization Methods 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000009413 insulation Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims description 2
- 229910001882 dioxygen Inorganic materials 0.000 claims description 2
- 229920000049 Carbon (fiber) Polymers 0.000 abstract description 7
- 239000004917 carbon fiber Substances 0.000 abstract description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 7
- 238000007380 fibre production Methods 0.000 abstract description 5
- 230000001502 supplementing effect Effects 0.000 abstract description 4
- 238000013461 design Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 229910002804 graphite Inorganic materials 0.000 description 1
- 239000010439 graphite Substances 0.000 description 1
- 239000011261 inert gas Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006068 polycondensation reaction Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000000197 pyrolysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 230000009469 supplementation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002912 waste gas Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- D—TEXTILES; PAPER
- D01—NATURAL OR MAN-MADE THREADS OR FIBRES; SPINNING
- D01F—CHEMICAL FEATURES IN THE MANUFACTURE OF ARTIFICIAL FILAMENTS, THREADS, FIBRES, BRISTLES OR RIBBONS; APPARATUS SPECIALLY ADAPTED FOR THE MANUFACTURE OF CARBON FILAMENTS
- D01F9/00—Artificial filaments or the like of other substances; Manufacture thereof; Apparatus specially adapted for the manufacture of carbon filaments
- D01F9/08—Artificial filaments or the like of other substances; Manufacture thereof; Apparatus specially adapted for the manufacture of carbon filaments of inorganic material
- D01F9/12—Carbon filaments; Apparatus specially adapted for the manufacture thereof
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/80—Data visualisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Crystallography & Structural Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Databases & Information Systems (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Textile Engineering (AREA)
- Furnace Details (AREA)
- Muffle Furnaces And Rotary Kilns (AREA)
Abstract
A visual simulation method for dynamic distribution characteristics of oxygen concentration in a high-temperature carbonization furnace relates to the technical field of a real-time monitoring method for the oxygen concentration distribution in the high-temperature carbonization furnace used in carbon fiber production. The method comprises the following steps: (1) constructing a three-dimensional mathematical model for calculating the total flow field of the high-temperature carbonization furnace; (2) establishing a three-dimensional simulation model of a muffle cavity and an inlet and outlet sealing cavity of the high-temperature carbonization furnace by adopting three-dimensional CAD software SOLIDWORKS, and setting related parameters; (3) carrying out mesh division on the three-dimensional simulation model; (4) importing the three-dimensional simulation model subjected to grid division into an Ansys Fluent module and setting; (5) setting a detection point and a detection surface, and carrying out simulation operation to obtain a simulation result; (5) and (5) under the same setting condition, setting different parameters and repeating the steps (2) to (5) to perform simulation calculation for multiple times. The method can visually judge the oxygen concentration distribution state in the high-temperature carbonization furnace, can better realize the monitoring of the oxygen in the furnace and determine the reasonable flow velocity of the nitrogen supplementing gas in the furnace, thereby providing theoretical data for reducing the oxygen concentration.
Description
Technical Field
The invention relates to the technical field of a method for monitoring the oxygen concentration distribution in a high-temperature carbonization furnace used in carbon fiber production in real time.
Background
Carbon fiber production belongs to the high energy consumption industry, and wherein high temperature carbonization stove is one of the power consumption big household in the carbon fiber production equipment, simultaneously, high temperature carbonization stove also is carbon fiber production's key equipment, mainly used carries out high temperature carbonization to the preoxidized fiber, makes it turn into the carbon fiber that carbon element content is greater than 90%. The high-temperature carbonization furnace is an integration of high-temperature technology and high-temperature equipment, and the use temperature is generally 1000-1600 ℃. The carbonization is carried out in inert gas for pyrolysis and polycondensation reaction, oxygen cannot exist in the carbonization furnace, the muffle furnace made of graphite can generate chemical reaction with the oxygen, the service life is shortened, and the performance of the carbon fiber is reduced at the same time, but trace oxygen still enters the carbonization furnace in practice, mainly because air adsorbed among monofilaments or entrained air is brought into the carbonization furnace along with the operation of the tows, and the large tows are more serious; the conventional method is to charge nitrogen at the starting stage of equipment, produce after detecting that the oxygen content in the furnace meets the requirements, and the oxygen content is dynamically changed in the normal working process of the equipment and is not usually monitored in real time, so when the oxygen content exceeds the standard, the equipment cannot give an alarm, so that chemical reaction occurs in the furnace, the service life of the furnace body is shortened for a long time, and the product quality is reduced. Therefore, a reasonable design method needs to be selected, so that the furnace chamber can reach the oxygen content meeting the specification, the reaction in the furnace is reduced, and the service life of the furnace body is prolonged.
Disclosure of Invention
In summary, the invention aims to solve the technical deficiencies that the oxygen concentration in the existing high-temperature carbonization furnace is not monitored in real time, when the oxygen content exceeds the standard, the equipment cannot give an alarm, so that chemical reaction occurs in the furnace, the service life of the furnace body is shortened for a long time, and the product quality is reduced, and provides a visual simulation method for the dynamic distribution characteristic of the oxygen concentration in the high-temperature carbonization furnace; the reasonable airflow speed of the nitrogen supplement gas in the furnace is determined by simulating the concentration of the oxygen in the furnace in the design process, and the dynamic distribution characteristic of the concentration of the oxygen in the furnace is visualized. The oxygen concentration distribution state in the high-temperature carbonization furnace can be visually judged, oxygen monitoring in the furnace can be better realized, and reasonable nitrogen gas supplementing air flow velocity in the furnace can be determined, so that theoretical data are provided for reducing the oxygen concentration.
In order to solve the technical problems provided by the invention, the technical scheme is as follows:
a visual simulation method for the dynamic distribution characteristics of oxygen concentration in a high-temperature carbonization furnace is characterized in that
The method comprises the following steps:
(1) and constructing a three-dimensional mathematical model for calculating the total flow field of the high-temperature carbonization furnace, wherein the three-dimensional continuity equation, the momentum equation and the energy conservation equation are respectively shown in formulas (1), (2) and (3):
three-dimensional continuity equation:
where ρ -fluid density; t-time; v-velocity vector, where u, V, w are the components of V in the three x, y and z directions;
Navier-Stokes equation for momentum equation:
whereinMu is dynamic viscosity, FbIs the volume force on the infinitesimal;
energy conservation equation:
wherein, CpSpecific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, ST-a viscous dissipation term;
oxygen enters the furnace chamber and then is mixed with nitrogen, the mass fraction of each local phase needs to be calculated by using a component transport equation, and a gas diffusion component equation can be obtained according to a component mass conservation law, as shown in a formula (4):
wherein ω is the mass fraction of each component, DmIs the turbulent diffusion coefficient;
(2) according to the geometric parameters of the high-temperature carbonization furnace in the actual engineering, establishing a three-dimensional simulation model of the muffle cavity of the high-temperature carbonization furnace and the inlet and outlet seal cavity by using SOLIDWORKS, and setting related parameters of the three-dimensional simulation model of the muffle cavity of the high-temperature carbonization furnace and the inlet and outlet seal cavity;
(3) transmitting the established three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace into a Blocking module of ICEM software, carrying out grid division on the three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace in an O-Block mode in the Blocking module, wherein a grid division strategy adopts a BiGeometric mode, a control ratio factor is a default value of 1.2, the grid quality of the whole structure is ensured to be more than 0.5 according to a judgment standard of the grid quality in the software, and meanwhile, the inlet and outlet and wall boundary names of the three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace are defined and comprise a seal cavity air inlet, an inlet and outlet and a furnace cavity wall boundary name;
(4) importing the three-dimensional simulation model divided into the grids in the step (3) into a FLUENT module of ANSYS software, and setting the FLUENT module;
(5) setting a detection point and a detection surface in a FLUENT module in ANSYS software, and carrying out simulation operation to obtain a simulation result which is used as an index for judging the dynamic distribution characteristic of the oxygen concentration of the muffle cavity of the high-temperature carbonization furnace;
(6) under the same setting condition, setting different parameters of the three-dimensional simulation model and repeating the steps (2) to (5) to perform multiple times of simulation calculation, and determining the nitrogen gas flow velocity for controlling the oxygen content by using an oxygen content change curve of a monitoring point in the furnace chamber, a nitrogen content change curve, an oxygen content cloud chart of a monitoring surface and a nitrogen content cloud chart of the monitoring surface as indexes for evaluating the sealing effect of the furnace chamber and the oxygen content distribution so as to visually judge the oxygen concentration distribution state in the high-temperature carbonization furnace, realize the oxygen monitoring in the furnace and determine the reasonable nitrogen gas supplement flow velocity in the furnace, provide theoretical data for reducing the oxygen concentration, and further obtain the optimal scheme of the oxygen gas concentration distribution and the nitrogen gas supplement in the muffle chamber.
The technical scheme for further limiting the invention comprises the following steps:
the parameters set in the step (2) at least comprise: the geometric shape and the geometric dimension of the muffle cavity, the geometric shape and the geometric dimension of the inlet and outlet seal, the inlet dimension of the inlet and outlet seal nitrogen gas pipe and the outlet dimension of the nitrogen gas pipe.
3. The method for simulating the visualization of the dynamic distribution characteristics of the oxygen concentration in the high-temperature carbonization furnace according to claim 1, wherein: in step (4), the process of setting the Ansys Fluent module is as follows:
(4.1) importing a User-Defined temperature parameter compiled according to the equipment operation process parameter in a User Defined option;
(4.2) in the General option, setting the y-direction gradient accumulation as a preset value according to requirements, and setting the time option as Transient heat transfer;
(4.3) selecting an Energy Equation from the Models options, selecting a laminar model from the Viscous Models options, and introducing a Reynolds number for judging the motion state of the airflow in the furnace cavity for description, wherein the Reynolds number has a calculation formula as follows:
wherein v, rho and mu are respectively the flow velocity, density and viscosity coefficient of the fluid, and d is the characteristic length; selecting a turbulence model as a laminar model through the calculation of Reynolds number;
(4.4) the Species Model in the Models option is checked against the Species Transport and
the Options check for Inlet Diffusion and Diffusion Energy Source.
(4.5) selecting oxygen and nitrogen in Materials Fluid options section, selecting oxygen and nitrogen in Selected specifices in mix Template;
(4.6) in the Cell Zone Conditions option, part of Fluid1 is set to nitrogen and part of Fluid2 is set to oxygen;
(4.7) setting an inlet Boundary condition as Pressure-inlet in a Boundary Conditions option, setting a Velocity map as a preset value according to actual parameters, setting a Thermal option as UDF tm-inlet, setting an outlet Boundary condition as Pressure-outlet, setting a furnace chamber on one side as a convection heat exchange surface, defining a comprehensive temperature value of furnace wall air in each hour by UDF, setting a convection heat exchange coefficient as a preset value according to an actual process, and setting other wall surfaces as heat insulation wall surfaces;
and (4.8) after selecting the Check case, calculating based on the three-dimensional mathematical model in the step (1).
In the step (5), a central point of the three-dimensional simulation model in the muffle furnace is set as a monitoring point in a FLUENT module in ANSYS software, and a monitoring surface is a Z-direction plane passing through the central point.
In the step (5), the simulation result includes: the oxygen and nitrogen concentration change curves of the monitoring points, the oxygen and nitrogen concentration cloud pictures of the monitoring surface and the oxygen concentration cloud pictures of the furnace chamber.
Compared with the prior art, the invention has the advantages that:
(1) regarding the problem of the oxygen concentration distribution characteristic of the high-temperature carbonization furnace, the conventional method is to perform a sensor test singly before the equipment is started, and start the equipment after the oxygen concentration control requirement is met, so that the oxygen concentration of the equipment in the production process cannot be dynamically monitored. The invention creatively utilizes the oxygen concentration and nitrogen concentration change curves of the monitoring points and the oxygen concentration and nitrogen concentration cloud pictures of the detection surface to judge the oxygen concentration distribution characteristic capability of the hearth, and can better optimize the internal airflow parameters during operation.
(2) By a numerical simulation method, a muffle furnace model, fluid parameters, different oxygen concentrations and other factors and a macroscopic performance relation model are established, and the influence of different gas flow velocities on the dynamic distribution characteristics of the oxygen concentrations is theoretically predicted, so that the rationality of design can be judged, and the design cost is reduced.
(3) The method can obtain the distribution characteristics of the oxygen concentration in the muffle furnace, the flow speed of different gases and the distribution rule of different oxygen concentrations at different time periods, thereby providing reference for the design of the high-temperature carbonization furnace.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional model created in the simulation method of the present invention.
Fig. 2 and 3 are the oxygen and nitrogen concentration change curves of the monitoring points of the invention.
Fig. 4-9 are schematic diagrams of different oxygen and nitrogen concentration clouds of the monitoring surface of the invention.
Fig. 10-13 are schematic diagrams of the dynamic characteristic visualization cloud charts of different oxygen concentrations in the furnace chamber.
Detailed Description
The invention will be further described with reference to the accompanying drawings and preferred embodiments of the invention.
A visual simulation method for the dynamic distribution characteristic of oxygen concentration in a high-temperature carbonization furnace comprises the following steps:
(1) and constructing a three-dimensional mathematical model for calculating the total flow field of the high-temperature carbonization furnace, wherein the three-dimensional continuity equation, the momentum equation and the energy conservation equation are respectively shown in formulas (1), (2) and (3):
three-dimensional continuity equation:
where ρ -fluid density; t-time; v-velocity vector, where u, V, w are the components of V in the three x, y and z directions;
Navier-Stokes equation for momentum equation:
wherein μ is dynamic viscosity, FbIs the volume force on the infinitesimal;
energy conservation equation:
wherein, CpSpecific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, ST-a viscous dissipation term;
oxygen enters the furnace chamber and then is mixed with nitrogen, the mass fraction of each local phase needs to be calculated by using a component transport equation, and a gas diffusion component equation can be obtained according to a component mass conservation law, as shown in a formula (4):
wherein ω is the mass fraction of each component, DmIs the turbulent diffusion coefficient.
(2) As shown in fig. 1, according to the geometric parameters of the high-temperature carbonization furnace in the actual engineering, adopting SOLIDWORKS software in three-dimensional auxiliary software to establish three-dimensional simulation models of a muffle cavity and an inlet and outlet seal cavity of the high-temperature carbonization furnace, and setting related parameters of the three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace; the specifically set parameters at least include: the geometric shape and the geometric dimension of the muffle cavity, the geometric shape and the geometric dimension of the inlet and outlet seal, the inlet dimension of the inlet and outlet seal nitrogen gas pipe and the outlet dimension of the nitrogen gas pipe.
(3) The method comprises the steps of transmitting established three-dimensional simulation models of a muffle cavity and an inlet and outlet seal cavity of the high-temperature carbonization furnace to a Blocking module of ICEM software, carrying out grid division on the three-dimensional simulation models in the Blocking module in an O-Block mode, adopting a BiGeometric mode for a grid division strategy, controlling a ratio factor to be a default value of 1.2, simultaneously ensuring that the grid quality of an integral structure is more than 0.5 according to a judgment standard of the grid quality in the software, and defining the names of an inlet, an outlet and a wall boundary of the three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace for the convenience of later-stage setting of calculation conditions, wherein the names of the inlet, the outlet and the wall boundary of the furnace cavity comprise;
(4) importing the three-dimensional simulation model divided into the grids in the step (3) into a FLUENT module of ANSYS software, and setting the FLUENT module; the specific process of setting the Ansys Fluent module is as follows:
(4.1) importing a User-Defined temperature parameter compiled according to the equipment operation process parameter in a User Defined option;
(4.2) in the General option, setting the y-direction gradient accumulation as a preset value according to requirements, and setting the time option as Transient heat transfer;
(4.3) selecting an Energy Equation from the Models options, selecting a laminar model from the Viscous Models options, and introducing a Reynolds number for judging the motion state of the airflow in the furnace cavity for description, wherein the Reynolds number has a calculation formula as follows:
wherein v, rho and mu are respectively the flow velocity, density and viscosity coefficient of the fluid, and d is the characteristic length; selecting a turbulence model as a laminar model through the calculation of Reynolds number;
(4.4) the Species Model in the Models option checks the Species Transport and the Inlet Diffusion and Diffusion Energy Source in the Options.
(4.5) selecting oxygen and nitrogen in Materials Fluid options section, selecting oxygen and nitrogen in Selected specifices in mix Template;
(4.6) in the Cell Zone Conditions option, part of Fluid1 is set to nitrogen and part of Fluid2 is set to oxygen;
(4.7) setting an inlet Boundary condition as Pressure-inlet in a Boundary Conditions option, setting a Velocity map as a preset value according to actual parameters, setting a Thermal option as UDF tm-inlet, setting an outlet Boundary condition as Pressure-outlet, setting a furnace chamber on one side as a convection heat exchange surface, defining a comprehensive temperature value of furnace wall air in each hour by UDF, setting a convection heat exchange coefficient as a preset value according to an actual process, and setting other wall surfaces as heat insulation wall surfaces;
and (4.8) after selecting the Check case, calculating based on the three-dimensional mathematical model in the step (1).
(5) Setting a detection point and a detection surface in a FLUENT module in ANSYS software, and carrying out simulation operation to obtain a simulation result which is used as an index for judging the dynamic distribution characteristic of the oxygen concentration of the muffle cavity of the high-temperature carbonization furnace; a central point of a three-dimensional simulation model in a muffle furnace is set as a monitoring point in a FLUENT module in ANSYS software, and a monitoring surface is a Z-direction plane passing through the central point. The simulation result comprises: the oxygen and nitrogen concentration change curves of the monitoring points, the oxygen and nitrogen concentration cloud pictures of the monitoring surface and the oxygen concentration cloud pictures of the furnace chamber.
(6) Under the same setting condition, different parameters are set by the three-dimensional simulation model, the steps (2) to (5) are repeated, multiple times of simulation calculation are carried out, the oxygen content change curve of a monitoring point in the furnace chamber, the nitrogen content change curve, the oxygen content cloud chart of a monitoring surface and the nitrogen content cloud chart of the monitoring surface are used as indexes for evaluating the sealing effect of the furnace chamber and the oxygen content distribution, so that the nitrogen gas flow speed for controlling the oxygen content is determined, the oxygen concentration distribution state in the high-temperature carbonization furnace is visually and intuitively judged, the oxygen monitoring in the furnace can be better realized, the reasonable nitrogen gas supplementing flow speed in the furnace is determined, theoretical data are provided for reducing the oxygen concentration, and the optimal scheme of the muffle oxygen furnace chamber gas concentration distribution and nitrogen gas supplementation is obtained.
The invention can obtain the change rule of the oxygen concentration dynamic distribution in the furnace cavity under different flow rates by modifying the flow rate in the inlet boundary velocity-inlet. Comparing the curves of the oxygen and nitrogen concentration with time at the monitoring point shows that the trends of the curves of the two gas concentrations increase oppositely with time, as shown in fig. 2 and 3, therefore, the increase of the nitrogen gas flow speed has obvious effect on accelerating the reduction of the oxygen concentration.
As shown in fig. 4-9, as can be seen from the cloud diagrams of different oxygen concentration distributions and nitrogen concentrations on the detection surface in the furnace chamber, the nitrogen concentration distribution is opposite to the oxygen distribution, the oxygen concentration distribution is not uniform, a large amount of oxygen is gathered at the inlet seal, which is caused by the adsorption and entrainment of the tows, the increase of the nitrogen concentration in the furnace chamber can significantly reduce the oxygen concentration, and the method has obvious effect; the schematic diagrams of the dynamic change of the oxygen concentration of the furnace chamber from different time periods are shown in fig. 10-13, and show that the oxygen concentration in the furnace chamber is gradually reduced along with the continuous supplement of nitrogen, a large amount of oxygen in the furnace chamber is discharged from a waste gas port, the oxygen concentration in the furnace chamber approaches to 0, the oxygen is gathered at an inlet sealing part by supplementing the nitrogen, the effect of isolating the oxygen by using the nitrogen is achieved, meanwhile, the visual simulation of the dynamic characteristic of the oxygen distribution in the furnace chamber is realized, and the optimal scheme of the concentration distribution of the oxygen in the muffle furnace chamber and the nitrogen supplement is obtained by verifying the simulation result of ANSYS, multiple simulations and comparing and analyzing results.
The embodiments of the present invention are described only for the preferred embodiments of the present invention, and not for the limitation of the concept and scope of the present invention, and various modifications and improvements made to the technical solution of the present invention by those skilled in the art without departing from the design concept of the present invention shall fall into the protection scope of the present invention, and the technical content of the present invention which is claimed is fully set forth in the claims.
Claims (5)
1. A visual simulation method for the dynamic distribution characteristic of oxygen concentration in a high-temperature carbonization furnace is characterized by comprising the following steps:
(1) and constructing a three-dimensional mathematical model for calculating the total flow field of the high-temperature carbonization furnace, wherein the three-dimensional continuity equation, the momentum equation and the energy conservation equation are respectively shown in formulas (1), (2) and (3):
three-dimensional continuity equation:
where ρ -fluid density; t-time; v-velocity vector, where u, V, w are the components of V in the three x, y and z directions;
Navier-Stokes equation for momentum equation:
wherein μ is dynamic viscosity, FbIs the volume force on the infinitesimal;
energy conservation equation:
wherein, CpSpecific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, ST-a viscous dissipation term;
oxygen enters the furnace chamber and then is mixed with nitrogen, the mass fraction of each local phase needs to be calculated by using a component transport equation, and a gas diffusion component equation can be obtained according to a component mass conservation law, as shown in a formula (4):
wherein ω is the mass fraction of each component, DmIs turbulent flowA diffusion coefficient;
(2) according to the geometric parameters of the high-temperature carbonization furnace in the actual engineering, establishing a three-dimensional simulation model of the muffle cavity of the high-temperature carbonization furnace and the inlet and outlet seal cavity by using SOLIDWORKS, and setting related parameters of the three-dimensional simulation model of the muffle cavity of the high-temperature carbonization furnace and the inlet and outlet seal cavity;
(3) transmitting the established three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace to a Blocking module of ICEM software, carrying out grid division on the three-dimensional simulation models in the Blocking module in an O-Block mode, wherein a grid division strategy adopts a BiGeometric mode, a control ratio factor is a default value of 1.2, the grid quality of the whole structure is ensured to be more than 0.5 according to a judgment standard of the grid quality in the software, and meanwhile, the inlet and outlet and wall boundary names of the three-dimensional simulation models of the muffle cavity and the inlet and outlet seal cavity of the high-temperature carbonization furnace are defined, wherein the inlet and outlet and wall boundary names comprise a seal cavity inlet, a furnace cavity inlet and outlet and a furnace cavity boundary name;
(4) importing the three-dimensional simulation model divided into the grids in the step (3) into a FLUENT module of ANSYS software, and setting the FLUENT module;
(5) setting a detection point and a detection surface in a FLUENT module in ANSYS software, and carrying out simulation operation to obtain a simulation result which is used as an index for judging the dynamic distribution characteristic of the oxygen concentration of the muffle cavity of the high-temperature carbonization furnace;
(6) under the same setting condition, setting different parameters of the three-dimensional simulation model, repeating the steps (2) to (5), carrying out multiple times of simulation calculation, and determining the nitrogen gas flow velocity for controlling the oxygen content by using an oxygen content change curve of a monitoring point in the furnace chamber, a nitrogen content change curve, an oxygen content cloud chart of a monitoring surface and a nitrogen content cloud chart of the monitoring surface as indexes for evaluating the sealing effect of the furnace chamber and the oxygen content distribution so as to visually judge the oxygen concentration distribution state in the high-temperature carbonization furnace, realize the oxygen monitoring in the furnace and determine the reasonable nitrogen gas supplement flow velocity in the furnace, provide theoretical data for reducing the oxygen concentration, and further obtain the optimal scheme of the oxygen gas concentration distribution and the nitrogen gas supplement in the muffle furnace chamber.
2. The method for simulating the visualization of the dynamic distribution characteristics of the oxygen concentration in the high-temperature carbonization furnace according to claim 1, wherein: the parameters set in the step (2) at least comprise: the geometric shape and the geometric dimension of the muffle cavity, the geometric shape and the geometric dimension of the inlet and outlet seal, the inlet dimension of the inlet and outlet seal nitrogen gas pipe and the outlet dimension of the nitrogen gas pipe.
3. The method for simulating the visualization of the dynamic distribution characteristics of the oxygen concentration in the high-temperature carbonization furnace according to claim 1, wherein: in step (4), the process of setting the Ansys Fluent module is as follows:
(4.1) importing a User-Defined temperature parameter compiled according to the equipment operation process parameter in a User Defined option;
(4.2) in the General option, setting the y-direction gradient accumulation as a preset value according to requirements, and setting the time option as Transient heat transfer;
(4.3) selecting an Energy Equation from the Models options, selecting a laminar model from the Viscous Models options, and introducing a Reynolds number for judging the motion state of the airflow in the furnace cavity for description, wherein the Reynolds number has a calculation formula as follows:
wherein v, rho and mu are respectively the flow velocity, density and viscosity coefficient of the fluid, and d is the characteristic length; selecting a turbulence model as a laminar model through the calculation of Reynolds number;
(4.4) the Species Model in the Models option checks the Species Transport and the Inlet Diffusion and Diffusion Energy Source in the Options.
(4.5) selecting oxygen and nitrogen in Materials Fluid options section, selecting oxygen and nitrogen in Selected specifices in mix Template;
(4.6) in the Cell Zone Conditions option, part of Fluid1 is set to nitrogen and part of Fluid2 is set to oxygen;
(4.7) setting an inlet Boundary condition as Pressure-inlet in a Boundary Conditions option, setting a Velocity map as a preset value according to actual parameters, setting a Thermal option as UDF tm-inlet, setting an outlet Boundary condition as Pressure-outlet, setting a furnace chamber on one side as a convection heat exchange surface, defining a comprehensive temperature value of furnace wall air in each hour by UDF, setting a convection heat exchange coefficient as a preset value according to an actual process, and setting other wall surfaces as heat insulation wall surfaces;
and (4.8) after selecting the Check case, calculating based on the three-dimensional mathematical model in the step (1).
4. The method for simulating the visualization of the dynamic distribution characteristics of the oxygen concentration in the high-temperature carbonization furnace according to claim 1, wherein: in the step (5), a central point of the three-dimensional simulation model in the muffle furnace is set as a monitoring point in a FLUENT module in ANSYS software, and a monitoring surface is a Z-direction plane passing through the central point.
5. The method for simulating the visualization of the dynamic distribution characteristics of the oxygen concentration in the high-temperature carbonization furnace according to claim 1, wherein: in the step (5), the simulation result includes: the oxygen and nitrogen concentration change curves of the monitoring points, the oxygen and nitrogen concentration cloud pictures of the monitoring surface and the oxygen concentration cloud pictures of the furnace chamber.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011275672.5A CN112435716B (en) | 2020-11-16 | 2020-11-16 | Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011275672.5A CN112435716B (en) | 2020-11-16 | 2020-11-16 | Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112435716A true CN112435716A (en) | 2021-03-02 |
| CN112435716B CN112435716B (en) | 2023-04-07 |
Family
ID=74701198
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011275672.5A Active CN112435716B (en) | 2020-11-16 | 2020-11-16 | Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112435716B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113484470A (en) * | 2021-07-01 | 2021-10-08 | 浙江始祖鸟环境工程有限公司 | High-sensitivity hydrogen sulfide gas detector |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140260313A1 (en) * | 2013-03-12 | 2014-09-18 | The University Of Sydney | Micro-mixer/combustor |
| CN104361195A (en) * | 2014-09-15 | 2015-02-18 | 燕山大学 | Three-dimensional flow thermal coupling modeling method for cement grate cooler |
| US20150152344A1 (en) * | 2012-06-08 | 2015-06-04 | How Kiap Gueh | Melt gasifier system |
| CN104707480A (en) * | 2015-02-12 | 2015-06-17 | 杭州天蓝环保设备有限公司 | CFD-based SNCR-SCR denitration process ammonia supplementing design method |
| CN106503380A (en) * | 2016-10-28 | 2017-03-15 | 中国科学院自动化研究所 | Coking nitrogen oxides in effluent concentration prediction method and forecasting system |
| CN107679314A (en) * | 2017-09-27 | 2018-02-09 | 福建中烟工业有限责任公司 | A kind of method and system for the luid mechanics computation model for establishing simulation cigarette burning |
| CN108119761A (en) * | 2017-12-14 | 2018-06-05 | 中国石油天然气股份有限公司 | A kind of gas-liquid two-phase flows down natural gas line extent of corrosion Forecasting Methodology |
| CN109858192A (en) * | 2019-03-15 | 2019-06-07 | 中国人民解放军陆军装甲兵学院 | Numerical simulation calculation method based on FLUENT MHD module |
| US20200082036A1 (en) * | 2018-09-12 | 2020-03-12 | AGC Inc. | Simulation method, simulation method by mbd program, numerical analysis apparatus, numerical analysis system for mbd, numerical analysis program, and mbd program |
| CN111400958A (en) * | 2020-04-15 | 2020-07-10 | 天华化工机械及自动化研究设计院有限公司 | ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace |
| CN111400934A (en) * | 2020-04-15 | 2020-07-10 | 天华化工机械及自动化研究设计院有限公司 | Method for simulating oxygen content distribution of low-temperature carbonization furnace based on WORKBENCH |
| CN111651908A (en) * | 2020-07-17 | 2020-09-11 | 天华化工机械及自动化研究设计院有限公司 | ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace |
-
2020
- 2020-11-16 CN CN202011275672.5A patent/CN112435716B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150152344A1 (en) * | 2012-06-08 | 2015-06-04 | How Kiap Gueh | Melt gasifier system |
| US20140260313A1 (en) * | 2013-03-12 | 2014-09-18 | The University Of Sydney | Micro-mixer/combustor |
| CN104361195A (en) * | 2014-09-15 | 2015-02-18 | 燕山大学 | Three-dimensional flow thermal coupling modeling method for cement grate cooler |
| CN104707480A (en) * | 2015-02-12 | 2015-06-17 | 杭州天蓝环保设备有限公司 | CFD-based SNCR-SCR denitration process ammonia supplementing design method |
| CN106503380A (en) * | 2016-10-28 | 2017-03-15 | 中国科学院自动化研究所 | Coking nitrogen oxides in effluent concentration prediction method and forecasting system |
| CN107679314A (en) * | 2017-09-27 | 2018-02-09 | 福建中烟工业有限责任公司 | A kind of method and system for the luid mechanics computation model for establishing simulation cigarette burning |
| CN108119761A (en) * | 2017-12-14 | 2018-06-05 | 中国石油天然气股份有限公司 | A kind of gas-liquid two-phase flows down natural gas line extent of corrosion Forecasting Methodology |
| US20200082036A1 (en) * | 2018-09-12 | 2020-03-12 | AGC Inc. | Simulation method, simulation method by mbd program, numerical analysis apparatus, numerical analysis system for mbd, numerical analysis program, and mbd program |
| CN109858192A (en) * | 2019-03-15 | 2019-06-07 | 中国人民解放军陆军装甲兵学院 | Numerical simulation calculation method based on FLUENT MHD module |
| CN111400958A (en) * | 2020-04-15 | 2020-07-10 | 天华化工机械及自动化研究设计院有限公司 | ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace |
| CN111400934A (en) * | 2020-04-15 | 2020-07-10 | 天华化工机械及自动化研究设计院有限公司 | Method for simulating oxygen content distribution of low-temperature carbonization furnace based on WORKBENCH |
| CN111651908A (en) * | 2020-07-17 | 2020-09-11 | 天华化工机械及自动化研究设计院有限公司 | ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace |
Non-Patent Citations (2)
| Title |
|---|
| 于宽: "大型碳纤维预氧化装备温度控制特性研究及应用", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》 * |
| 孙汉生: "循环流化床锅炉降氧抑氮燃烧模式的模拟与试验研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113484470A (en) * | 2021-07-01 | 2021-10-08 | 浙江始祖鸟环境工程有限公司 | High-sensitivity hydrogen sulfide gas detector |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112435716B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111400958B (en) | ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace | |
| CN111651908B (en) | ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace | |
| WO2022099715A1 (en) | Simulation method for visualizing dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace | |
| CN112446177B (en) | Simulation method for heat insulation performance of external heat insulation material of high-temperature carbonization furnace | |
| WO2022099713A1 (en) | Three-dimensional simulation method for tow heating process in low temperature carbonization furnace based on overset model | |
| CN111400934B (en) | Method for simulating oxygen content distribution of low-temperature carbonization furnace based on WORKBENCH | |
| CN111680432B (en) | Low-temperature carbonization furnace multi-coupling field stress distribution simulation method based on WORKBENCH | |
| WO2022011723A1 (en) | Ansys-based multi-field coupling stress distribution simulation method for high temperature carbonization furnace | |
| WO2021207952A1 (en) | High-temperature carbonization furnace oxygen content distribution simulation method based on ansys | |
| WO2022099714A1 (en) | Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace | |
| CN112528572A (en) | Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model | |
| CN112435716B (en) | Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace | |
| CN107023966A (en) | A kind of subway station air conditioning cooling water water outlet temperature setting value optimization method | |
| CN112270109B (en) | Method for simulating heating performance of graphite rod in high-temperature carbonization furnace | |
| CN107506516A (en) | A kind of communications equipment room flow field model is established and analysis method and system | |
| WO2021207953A1 (en) | Workbench-based low-temperature carbonization furnace oxygen content distribution simulation method | |
| CN114330050A (en) | Temperature and humidity numerical simulation method for underground workshop ventilation gallery of pumped storage power station | |
| CN113627020A (en) | Printing workshop VOCs diffusion model and intelligent monitoring system | |
| CN113536640B (en) | Optimization design method for internal flow channel structure of air distributor based on orthogonal test | |
| WO2022011724A1 (en) | Workbench-based method for simulating multi-coupling field stress distribution of low-temperature carbonization furnace | |
| WO2022099712A1 (en) | Simulation method for heating performance of graphite rod in high-temperature carbonization furnace | |
| CN110210137A (en) | A kind of drag computation method of SCR denitration | |
| CN111625906A (en) | Test and simulation combined type selection method for fan in power battery air cooling system | |
| WO2022099716A1 (en) | Method for simulating thermal insulation property of external thermal insulation material of high-temperature carbonization furnace | |
| CN114186353A (en) | Numerical simulation-based large-space air conditioner cross influence factor calculation method |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240130 Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Patentee after: CHINA PETROLEUM & CHEMICAL Corp. Country or region after: China Patentee after: Sinopec Shangai Petrochemical Co.,Ltd. Address before: 730060 Gansu city of Lanzhou province Xigu District Heshui Road No. 3 Patentee before: TIANHUA INSTITUTE OF CHEMICAL MACHINERY AND AUTOMATION Co.,Ltd. Country or region before: China |