In this work, the task is to construct a mathematical model and calculate with its help concentration and temperature fields, as well as the velocity distribution of gas sweats movement during the combustion of the heating gas in the vertical of the coke oven Structures of TDCs. Is being studied the effect of the excess air factor and the nature of the gas dynamic in the heating section of the coke oven to the concentration of the combustion products. As a modeling area, two heating piers (vertical) were considered, located near the middle of the coke oven of the PVR system with a volume of 41.6 m3. We considered a three-dimensional problem with a partition of the modeling region along the x axis into 24 cells; along the y axis, by 13 cells; on the z axis, 24 cells. Vertical dimensions: along the axis x0.48 m; on the y axis, 0.40 m; on the z-axis - 7,2 m. Simultaneous solution of three-dimensional equations of motion, heat energy transfer, mass transfer and kinetic interactions of gas flow components is carried out. The methods of mathematical modeling are used to estimate the velocity fields, traces and temperatures in the heating section of the coke oven at an excess air factor in the range 1.2-1.6. The possibility of overheating of the upper part of the vertical due to the appearance of a turbulent vortex is shown. Extreme behavior of the nitrogen monoxide content is obtained depending on the excess air factor. It is shown that the presence of asymmetry of the boundary velocities at the inlet and outlet from the vertical leads to a decrease in the recirculation and an increase in the concentration of nitrogen monoxides in the flue gases.
Optimization of coke oven operating modes is actual task. Simulation of the combustion process of gaseous fuels will make it possible to trace the distribution of the temperature, concentration, velocity fields along the height of the combustion zone and after it. This will allow us to analyze the effect of gas dynamics and combustion on the heating of various zones Refractory masonry, the dependence of the composition of flue gases, especially such harmful constituents such as formaldehyde and nitrogen oxides, on the nature of motion of gas flows. To solve this type of problem it is necessary simultaneous consideration of gas-dynamic, concentration and thermal problems.
The task is to construct a mathematical model and calculate with its help the concentration and temperature fields, and also distribution of speeds of movement of gas sweats at combustion of a heating gas in a vertical of coke oven of system PVR. The influence coefficient of excess air and the nature of gas dynamics in the heating space of the coke oven on the concentration of combustion products.
As the modeling area, two heating piers (vertical) were considered, located near the middle of the coke
furnaces of the PVR system with a volume of 41.6 m3. We considered a three-dimensional problem with
partitioning the modeling domain along the x-axis into 24 cells; on the y-axis - by 13 cells; on the z axis, 24 cells.
The vertical dimensions are chosen according to [
In solving the gasdynamic problem, the following equations were used [
- is the density.
Integration of equations (1) was carried out within the framework of the method of finite differences using
the "chess" grid [
The equations for the transfer of the concentrations of the components of the gas mixture in the form [
constant; chemical reaction by - concentration of component in mole fractions; the order of the component
the Schmidt number (it was chosen equal to 0.9 [
When integrating Eq. (3), boundary conditions of the first kind were applied at the entrance to the vertical and on the laying of the vertical; At the exit from the vertical, we used boundary conditions of the third kind in the form
The heat energy transfer equation is chosen in the form [
When integrating equation (4), the boundary conditions were applied I-th kind on the basis of experimental values of temperature measurements in the conditions of Coke-chemical production of JSC "NTMK". The temperature of the gas flow at the entrance to the vertical is of the order of 1000 ° C, the temperature on the masonry about 12500C.
In modeling the combustion process, a coke oven gas was used as a fuel, the interaction of its
components with air was carried out using model constructions detailed in [
The results of the calculations in the vertical section (the middle along the y axis) heating partition with = 1, 2 are shown in Fig. 1 and 2. In this case, the same boundary conditions are set for the velocity at the input and output from the vertical (1 m / s). There is a fairly intensive recirculation movement of gas flows. There is a turbulent vortex in the upper right side of the vertical (Fig. 1 (a)). In the same place is the zone of elevated temperature (Figure 1 (b)), which can affect negatively for the life of refractories of this part of the masonry. Formation of nitrogen oxides and formaldehyde occurs throughout the vertical region (Figure 2 (a, b)). (a)
(b) in the excess air factor. The content of nitrogen oxides passes through a maximum at about = 1.4 - 1.45. This may be due to the onset with an increase in NO due to an increase in the N2 concentration in the input stream, and then to a decrease in the formation of NO, due to a decrease in the temperature in the vertical due to a decrease in the proportion of coke oven gas in the composition of the heating mixture.
1. Methods of mathematical modeling have been used to estimate the velocity, concentration and temperature fields in the coke furnace heating section with an excess air factor in the range 1.2-1.6.
2. The possibility of overheating of the upper part of the vertical due to the appearance of a turbulent vortex is shown.
3. Extreme behavior of the content of nitrogen mono-oxides is obtained depending on the excess air factor. 4. It is shown that the presence of asymmetry of the boundary velocities at the inlet and outlet from the vertical leads to a decrease in the recirculation and an increase in the concentration of nitrogen monoxides in the flue gases.
The implementation of the model was carried out using the author's software package.
Limit velocity at 1,29 m/s vertical output