Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace
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 a 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 problem 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) 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, F b Is the volume force on the infinitesimal;
energy conservation equation:
wherein, C p Specific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, S T -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 of each componentFraction, D m Is the turbulent diffusion coefficient;
(2) 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 according to the geometric parameters of the high-temperature carbonization furnace in the actual engineering, 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 performing 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 dynamic distribution characteristic of the oxygen concentration in the high-temperature carbonization furnace according to claim 1, wherein the method comprises the following steps: 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 Energy options in the Models options, selecting a Laminar model from the Viscous Models options, and introducing a Reynolds number for describing in order to judge the airflow motion state in the furnace cavity, 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 Specifes Model in the Model option is hooked to the Specifes Transport and is positioned in
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) setting the Fluid1 part as nitrogen and the Fluid2 part as oxygen in the Cell Zone Conditions option;
(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 distribution characteristics of the oxygen concentration of the high-temperature carbonization furnace, the conventional method is to perform a sensor test singly before equipment is started, start the equipment after the requirement of controlling the oxygen concentration is met, and cannot dynamically monitor the oxygen concentration of the equipment in the production process. The invention creatively utilizes the oxygen concentration and nitrogen concentration change curves of the monitoring points and the oxygen concentration and nitrogen concentration cloud charts 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 at different time periods and the distribution rules of different gas flow speeds and different oxygen concentrations, 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 visualized cloud pictures of the dynamic characteristics of different oxygen concentrations in the furnace chamber of the present invention.
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) 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, F b Is the volume force on the infinitesimal;
energy conservation equation:
wherein, C p Specific heat capacity, T temperature, k fluid heat transfer coefficient, S T -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, D m Is 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) 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, and simultaneously ensuring that the grid quality of the integral structure is more than 0.5 according to a judgment standard of the grid quality in the software;
(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 visual access 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; through the calculation of Reynolds number, a turbulence model is selected as a laminar model;
(4.4) the Specifes Model in the Models option is hooked to the Specifes Transport, and the Inlet Diffusion and Diffusion Energy Source are hooked to the Options of Options.
(4.5) selecting oxygen and nitrogen in Materials Fluid options section, selecting oxygen and nitrogen in Selected specifices in mix Template;
(4.6) setting the Fluid1 part as nitrogen and the Fluid2 part as oxygen in the Cell Zone Conditions option;
(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) - (5) are repeated, multiple times of simulation calculation are carried out, and 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 furnace chamber sealing effect 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 determined, the oxygen monitoring in the furnace can be better realized, the reasonable nitrogen supplementing gas flow speed in the furnace can be determined, theoretical data are provided for reducing the oxygen concentration, and the optimal scheme of the oxygen concentration distribution and the nitrogen supplementation of the muffle furnace chamber 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, it can be seen from the cloud diagrams of different oxygen concentration distributions and nitrogen concentrations on the detection surface in the furnace chamber that the nitrogen concentration distribution is opposite to the oxygen distribution, the oxygen concentration distribution is not uniform, and a large amount of oxygen is gathered at the inlet sealing position, which is caused by the adsorption and entrainment of the tows, so that the concentration of the oxygen can be significantly reduced by increasing the concentration of the nitrogen in the furnace chamber, 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 described embodiments of the present invention are only for describing the preferred embodiments of the present invention, and do not limit the concept and scope of the present invention, and the technical solutions of the present invention should be modified and improved by those skilled in the art without departing from the design concept of the present invention, and the technical contents of the present invention which are claimed are all described in the claims.