Next Article in Journal
Luenberger Disturbance Observer-Based Deadbeat Predictive Control for Interleaved Boost Converter
Next Article in Special Issue
A Study of the Critical Velocity of the Droplet Transition from the Cassie to Wenzel State on the Symmetric Pillared Surface
Previous Article in Journal
Generalization of Some Fractional Integral Operator Inequalities for Convex Functions via Unified Mittag–Leffler Function
Previous Article in Special Issue
Self-Optimizing Path Tracking Controller for Intelligent Vehicles Based on Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring of the Behaviour and State of Nanoscale Particles in a Gas Cleaning System of an Ore-Thermal Furnace

The Automation of Technological Processes and Production Department, Saint Petersburg Mining University, 2, 21st Line, 199106 St Petersburg, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(5), 923; https://doi.org/10.3390/sym14050923
Submission received: 24 February 2022 / Revised: 22 April 2022 / Accepted: 27 April 2022 / Published: 1 May 2022

Abstract

:
The aim of this paper is to define and select stable zones in the off-gas duct of an ore-thermal furnace using a mathematical model. This is needed to increase the effectiveness of exhaust gas composition control in metallurgical silicon production. Methods. The goals of this study were achieved by means of computational fluid dynamics. A model with a water-cooled furnace roof as well as a model comprising steel gas passes with a sliding shutter was developed using ANSYS Fluent software. Both models were symmetrical to ensure a uniform gas-dust distribution, which allowed us to test the adequacy of the obtained models. The models were based on the Navier–Stokes equations system as well as on a discrete phase model (DPM) that was developed using the Euler–Lagrange method. Results. As a result of the modelling, a transition flow mode (Re 0-7437) was revealed behind the sliding shutter. As such, it can be assumed that the most suitable place for measuring equipment to be installed is directly behind the closed part of the sliding shutter.

1. Introduction

A large number of fine particles with various compositions are emitted into the atmosphere during carbothermic silicon reduction. Analyses of waste gases show that most emissions are related to the consumption of carbon materials that contribute to the overall carbon balance. These are carbon-graphite electrodes that are used for heating quartz and charcoal to reduce the amount of silicon released from oxides (quartz) [1,2,3]. The main component of waste fume emissions is SiO2 microsilica (up to 85 percent), which is present in a mixture of solid carbon in various forms and states (7–8 percent) and in silicon carbide (5 percent) [4,5]. As a rule, industrial emissions are not controlled, and dust is collected from gas pass systems and from deposits on equipment and building structures.
Granulometric analyses of the fumes captured by GCS electro-filters indicate the presence of particles in microsilica fumes that are 200–250 µm in size and have an elevated carbon nanoparticle content (up to 8 percent) that can be removed with the off-gas.
Currently, carbon-free microsilica is widely used as a modifying additive to base materials in the construction industry. Thus, the use of microsilica makes it possible to produce concrete with special properties: increased durability (resistance to the action of weak acids and seawater) and increased compression strength. Production methods for silicon-carbide powder materials, such as micronized carbide (particle size < 1 m) for ceramics and nanocarbides (particle size < 1 nm) for high-quality structural ceramics and galvanics, are being intensively developed, as they create a high-value of the final product [6,7].
The main sources of exhaust gas carbon emissions are carbon electrodes, activated carbon, and carbonaceous materials.
For electrodes, the main components of carbon flow are arc heating fractures, the main quartz reduction reaction (stoichiometry), the decrease and oxidation of the lateral faces, and the destruction of the electrode soles upon contact with the charge. The key factor in electrode mass erosion is the oxygen concentration in the furnace atmosphere. As a result of the increased formation and oxidation of carbon monoxide, the oxygen concentration in the furnace atmosphere will decrease [8,9,10].
There is scientific and technical interest in modelling the distribution of waste gas concentration fields during carbothermic silicon reduction. It is necessary to ensure that the furnace’s thermal conditions and the balance of consumable carbon are controlled effectively [11].

2. Problem Overview

2.1. Causes of Carbon Nanoparticle Formation during Carbothermic Silicon Reduction in Ore-Thermal Furnaces

At present, reducing the carbon footprint in pyrometallurgical processes [6] is relevant to the resolution of global decarbonization issues. The manufacturing process in ore- thermal furnaces (OTFs) involves multi-component carbon systems (carbon electrodes and lining, activated charcoal, modifying additives) that affect the overall carbon balance during smelting. Activated carbon is consumed by the main silicon reduction reaction by means of stoichiometry, as well as by the side reactions caused by interactions with the impurities in quartz that take place according to thermodynamic conditions. For example, one of these irrevocable losses can be attributed to the transfer of active carbon into the crust and the subsequent production of silicon carbide when there are temperature disturbances of more than 1700–1900 °C [12]. Additionally, thermal disturbances in the furnace and increased carbon consumption are associated with carbon oxidation (an increase in CO and CO2 content) or with the abrupt release of soot and carbon particles into the atmosphere of the furnace and into the gaseous space [13].
Charge materials are consumed for the main silicon reduction reaction to form silicon dioxide, which may be accompanied by a transition into carbon monoxide in parallel with the formation of intermediate silicon monoxide, as well as silicon carbide during overheating, especially at the start of the smelting process and during the primary arc [14]. A separate consumption item is the unresponsive carbon that is derived from charcoal, which forms nanoparticles in the form of amorphous carbon, active pure carbon, or fullerenes.
The carbon resulting from electrode chipping and falling passes into a fume mixture together with microsilica. In this case, carbon particles without interaction are adsorbed on the highly developed surface of the silica fumes to form microsilica, the main form of silicon waste. This is an item of consumption due to the changes in the silicon balance that range from 380 to 450 kg per ton of produced silicon and requires operational control. In this instance, sampling and relevant data on the chemical composition of the microsilica captured in the GCS are needed for analytical comparison with the source quartz and impurities to ultimately determine the quality of silicon products [15,16,17].

2.2. Microsilica Monitoring and Nanoparticle Capture in the Ore-Thermal Furnace Gas Cleaning System

The main problem with capturing and controlling nanoscale particles in the OTF GSC is the timely determination of the emission composition of the gas fumes (in mass). Conventional GSC implements a detailed analysis after the end of the electric filter cycle and after the filter has been cleaned and the sediment has been weighed according to the time. As a result, there is a delay in the response when the electric and technological modes of the furnace are disrupted and when the electrodes are burned, which is indicated by an increase in the carbon content in the waste gases [18]. Monitoring the waste composition during the production of silicon and its alloys is associated with a number of industry-specific problems:
  • The temperatures of the gases in the immediate vicinity of the furnace roof are very high (600–850 °C), resulting in the need to cool and ventilate the gas-dust flow before it makes contact with the sensitive instruments when taking extractive measurements;
  • The high concentration of particulates in front of the filter can cause rapid abrasive damage to the measuring instruments in the gas streams that are in contact with the particles and can also affect the measured data [19,20];
  • The turbulent mode of gas flow has a wide range of time and space scales for the pulsations of all of the flow characteristics. This makes the gas flow faster than the laminar flow and results in intensive mass exchanges with high-impulse and energy levels between different stream regions due to the intensive mixing of the dispersed medium. This results in the substance having an uneven distribution in the gas flow and a consequent distortion of the measurement results.
Taking into account the described flow conditions in the furnace atmosphere and in the gas duct system of GCSs, the question of how to determine the temperature and velocity control points of the particles arises. For this purpose, it is necessary to define stable zones on the mass transfer interaction path.
To effectively control the electric mode of the furnace, installing a gas analyzer between the gas duct and the gas purification equipment in the stable flow region of the gas stream is recommended. This allows additional data to be obtained to create a digital database of information about the process [21].
The continuous growth of computing capacities and measuring technologies in recent years has opened up radically new possibilities for the computer modelling of large production units [22].
The novelty of the proposed solution is in justifying the location of the gas analyzer installation with the aid of the created computational fluid dynamics (CFD) model. The ANSYS software package, specifically its CFD modelling module ANSYS Fluent, was used for the mathematical model.

2.3. Main Features of Flue Gas Movement in the Ore-Thermal Furnace Gas Cleaning System

In most cases, when taking the fundamental laws of gas-fluid dynamics [4] in the furnace gas duct into account, the flue flow has a high Reynolds value and is turbulent due to different densities.
This movement of gases with dust is accompanied by an intensive mixing with velocity and pressure pulsations, and in addition to the main longitudinal movement, transverse movements and rotational movements are observed in individual flow volumes, especially in the vortex zones that are close to the filters and ejectors. As a result, the local cross-sectional velocities of the entire gas content flow from the furnace vault to the GCS during the mixing process under the influence of a temperature gradient between the layers for the entire mixing period.
It should be noted that turbulence has a continuous influence on the main flow parameters such as the concentration of components, temperature, velocity, and heat state. Thus, each turbulent vortex volume has its own substance and temperature concentration [23].
To avoid components and pulsations mixing and to thus obtain adequate data on flow parameters and the adequacy of the matrix, it is necessary to define data for the transient flow by taking the measurements into account in practice. On the other hand, the transient flow mode is characterized by the low mixing rate of the pulsation particles compared with the turbulence mode. In this instance, the key value for determining the flow of gas is the Reynolds number (Re) that characterizes the mode changes.
For this process, it has been experimentally determined that the Re for straight, smooth pipes with the most perturbed flow at the entrance of the pipe is 2300. However, at Re values above this value and up to a certain limit, a transient (mixed) flow mode is observed, and after this limit, turbulent flow is more likely [24]. In cylindrical receptacles (pipes), this transition interval can be varied significantly by reducing the initial disturbance of the medium by up to 50,000. For this process, we determined the factors influencing the turbulent state of the system in the gas pass system from the furnace:
  • The current state of the fume environment (kinetics, thermodynamics, etc.): compressibility (velocity of motion of silica fume) taking into account multi-phase flow (interphase exchange).
  • The outer limits of gas flow by zone: the movement along the side lining into the vault and output into the gas channel without the influence of the channel geometry.
  • The process limit stage and GCS outlet (GCS electric filter surface pressure): the flow area on the surface of the electric filter and the pressure.
The Reynolds criteria for the different flow mode models are as follows:
  • At Re < 2 × 103
    laminar flow is observed.
  • At Re > 104
    the flow becomes turbulent, but when gases begin to exit the furnace vault and move into the gas channel, it is preserved.
  • At Re > 5 × 104
    a turbulent boundary layer begins to form during the beginning of the gas flow process due to a sharp change in the temperature.
  • At 2 × 104 < Re 104
    a transient flow and heat exchange mode are observed when the flow approaches the electric filter. Thus, given the turbulent flow conditions, the ratio of the average to maximum velocity increases with the turbulence. In the Re, this ratio is asymptotically close to 1 [25,26,27].

3. Materials and Methods

3.1. Computer Simulation of the Dispersion Fluid Dynamics

Offsetting the limit values for transient flows is a major difficulty in determining the flow parameters of a multi-component gas. This problem is particularly pronounced in the case of flow limits since the existing standard models can only estimate the average flow rate of a gas duct when averaging the Reynolds criteria values in ANSYS Fluent [28].
The exhaust ducts and flows of the carbothermic silicon reductions in the OTF were simulated to detect stable flow zones in which the control point for determining the concentration of the exhaust gases could be installed.
The classification of existing methods and approaches for the numerical modelling of turbulent flows is based on the level of detail required to refine and detect turbulent pulsations, as well as their energy spectrum and type of flow. Depending on this, methods can be divided into three groups (DNS—direct numerical simulation, LES—large eddy simulation, and RANS—Reynolds-averaged Navier–Stokes) [29].
For the current task, the best solution to account for large-scale turbulent vortices and flows in the boundary layer is a hybrid approach in which an RANS model is activated in a wall area with an LES model close to it. This is a detached eddy simulation (DES) method for simulating disconnected vortices that takes into account the reduction in turbulent vortices near a solid surface [30]. This method requires significant computational resources but ensures its adequacy compared to others.
In the considered case, taking the influence of the boundary layer in the gas flow transition zones on turbulence into account, the explicit resolution of even the largest vortices in the boundary layer requires considerable grinding in the calculation grid, resulting in an increase in the calculation time. Thus, a wall-adapted local eddy-viscosity model (WALE) [31] can solve this problem, taking the particularities of the gas flow and its thermal potential into consideration.
The WALE method requires a high-quality grid. The number of grid elements should be more than 1,000,000 at least, and boundary layers should be used for inflation. Otherwise, it is possible that uncorrected data will be obtained as a result of the numerical simulation.

3.2. Problem Statement: Modelling the Off-Gas Mixture under Ore-Thermal Furnace Conditions

The symmetry of the gas duct system ensures uniform flow separation. An increase in the waste gas concentration in one part of the duct will result in changes in the pressure drops in another part. This will consequently lead to asymmetric flows, making it impossible to define a stable flow area. Gas duct system symmetry allows us to evaluate the adequacy of the resulting models.
The following boundary conditions (Table 1, Table 2 and Table 3) were used as input components for modelling the furnace exhaust gas system: the water-cooled roof of the furnace and the steel gas passes. A sliding shutter was placed at a 45-degree angle at the fork of the two gas passes.
The average values of the fumes in the off-gas mixture are 10–12 percent.
According to the goals of this research, two models were developed:
  • A system with a water-cooled portion of the roof for gas passes.
  • A gas pass system without water cooling.

3.3. Governing Equations

The mathematical modelling of turbulence is based on the set of ratios used to describe the movement of a fluid/gas. The basic governing Navier–Stokes equation system consists of the equation of continuity (1) and the equation of motion (2):
ρ t + ρ u j x j = 0
ρ u i t + ρ u j u i x i = P x i + τ i j x j
where:
  • P —pressure of the gas mixture;
  • µ—gas mixture density;
  • t—time;
  • ui,j—velocity components in the i and j directions;
    τi,j—shear strain tensor.
To simulate the particulate matter (fume), a discrete phase model (DPM) was developed using the Euler–Lagrange method. It was constructed by solving the time-averaged Navier–Stokes equations for the liquid phase, which was treated as a continuum, while the dispersed phase was described by tracing the trajectory of a large number of particles through the calculated flow fields. The dispersed phase can exchange momentum, mass, and energy with the liquid phase.
The balance of the forces acting on a particle predicts the trajectory of the discrete phase particles by integrating the balance of the forces on the particle, which is written in the Lagrangian reference frame. This force balance equates inertia with the forces acting on the particle and can be written (based on the directions in the Cartesian coordinate system) as:
U p t = F D ( U U p ) + g x ( ρ p ρ ) ρ p + F x
where FD(UUp) is the drag force per unit mass of the particle.

4. Results and Discussion

4.1. Model 1: Combination of the Furnace’s Roof and Water-Cooled Gas Passes

As initial data, in addition to the variables specified in Table 1, Table 2 and Table 3 as well as the dimensions, the following parameters were made available: water-cooled thickness, 0.06 m; side lining thickness, 0.05 m; and thickness of the non-water-cooled part of the roof made of chamotte bricks, 0.5 m.
The grid of this model consisted of 2,037,611 elements to obtain adequate simulation accuracy and had a minimum orthogonal quality of at least 0.1 (Figure 1). The LES WALE turbulence model was used to solve the variable-density subsonic flow problem by means of a pressure–based coupled solver. The real Soave–Redlich–Kwong gas model [32,33] was used for the gas mixture.
The boundary conditions for both streams (exhaust gases and cooling water) are shown in Table 4, Table 5 and Table 6.
As a consequence of simulation, the temperature distribution and gas velocity profiles were obtained to define the boundary conditions for the second model as the area in which the intended gas analyzer could be installed.
At the outlet from the water-cooled part of the gas pass system, the gas flow rate was 5.5 m/s, which is equivalent to a mass flow rate of 15.54 kg/s. The temperature of the mixture was 450 °C.
Because the model had a horizontal symmetry, we were able to evaluate the temperature and velocity distribution using the cross-sectional area at the centre of one of the water-cooled parts of gas ducts 2 (Figure 2). This model makes it possible to determine the input parameters to assess possible variation among them. In this way, a prediction story was created for the next part of the GCS, in which the steel gas passes through a sliding shutter model. This provides stable initial characteristics for further modelling tasks.

4.2. Model: 2 Gas Passes without Water-Cooling

The simulated part of the exhaust gas system consists of gas passes with a rigid steel frame and a sliding shutter, which acts as a gas velocity regulator and is located at the gas pass junction. This part is of greatest interest for modelling and further analysis because stable zones were predicted to be here.
There were 2,739,629 grid elements in this model (Figure 3). The grid had an acceptable average orthogonal quality of 0.79.
Based on the fact that the most stable flows in the models will not be laminar ones because of the high velocities, considering the transient mode that was identified by means of ANSYS Fluent is advised. The Re ratio was calculated using the conventional equation:
Re = u d m i x ρ m i x μ m i x
where:
  • u—velocity, m/s;
  • dmix—hydraulic diameter, m;
  • ρmix—mixture density, kg/m3;
  • µmix—dynamic viscosity of the gas mixture, Pa · s;
  • dmix—hydraulic diameter, m;
  • ρmix—mixture density, kg/m3;
  • μmix—dynamic viscosity of the gas mixture, Pa · s.
Sutherland’s formula was used to determine the dynamic viscosity:
μ t = μ 0 273 + C T + C T 3 2 273
where:
µt—dynamic viscosity of the gas at the temperature, Pa · s;
µ0—dynamic viscosity of the gas at 0 °C, Pa · s;
T—the absolute temperature of the gas, K;
C—Sutherland’s constant.
Thus, the dynamic viscosity of the gas mixture can be found using the following equation:
M m i x μ m i x = a 1 M 1 μ 1 + a 2 M 2 μ 2 + + a n M n μ n
where Mmix, M1, M2, and Mn are the molecular masses of the gas mixtures and their components; a1, a2, and an are the volumetric fractions of the components in the gas mixtures; and µmix, µ1, µ2, and µn are the dynamic viscosities of the gas mixtures and their components, P a s [34,35].
The density of the furnace gas mixture is a very important indicator for determining the concentration field profile that is repeated on the GCS electric filters:
ρ m i x = y 1 ρ 1 + y 2 ρ 2 + + y n ρ n
where y1, y2, and yn are the volumetric fractions, and ρ1, ρ2, and ρn are the densities of the components, kg/m3.
ρ 0 273 273 + t
The kinematic viscosity was calculated using the following formula:
υ = μ m i x ρ m i x
As a result of calculating the required auxiliary parameters in ANSYS CFD-Post, the contours of the main off-gas parameters, such as velocity and kinematic viscosity, were obtained. These parameters have a direct impact on the Re as a key criterion for the flow conditions.
Any obstacle in the gas flow path changes the dynamic flow characteristics. As such, it is necessary to know what happens before and after the gas flow encounters the sliding shutter to evaluate any changes in the cross-section of the gas path. Knowing this, we can determine the control actions for the system.
Figure 4 illustrates the variations in the dynamic characteristics of the gas-dust flowover and the volume of the gas duct. A transient flow mode can be observed in the volume-rendering region behind the closed part of the sliding shutter (shown by the pointer) with the smallest Re. Taking into account the initial conditions of developing the best possible environment (laminar flow) for measurements to be taken in because of the high velocities, the transition mode revealed here can be considered to be an appropriate result of the modelling.
The flow velocity is the most significant parameter affecting both the flow mode and Re. The linear dependence of the Re on the flow velocity, as shown in Figure 5, allows us to estimate the flow mode using the volumetric velocity profile alone for future calculations. Based on this chart, the appropriate velocity should be less than 1 m/s to reach the transient Re.
Histograms, such as those in Figure 6, are needed to confirm the existence of a transitional mode and to estimate the proportion of its volumetric distribution in particular figures. A small fraction of transient off-gas flow (red column on the histogram (0-7437) Re) can be observed. These values are appropriate for minimizing turbulent pulsations in this case.
Thus, despite the wide variety of turbulence models in ANSYS Fluent, in cases where there is strong turbulence, considering transient flow modes is advised.

5. Conclusions

The silicon smelting process in OTFs is only controlled in the furnace and at the inlets and outlets of the gas duct. Thus, everything that goes into the gas ducts is a “black box”. Emissions can only be controlled when the electrostatic precipitators are clogged. Modelling allows us to look inside the process and to react to the changes that take place in a short period of time.
Something else that is important for process control is minimizing the effects of turbulent pulsations on measurements in the gas duct. The main innovation of this paper is the determination of stable zones by means of CFD modelling to avoid measurement errors. All of the methods mentioned above can profoundly increase the transparency and controllability of the silicon production process.
When producing silicon in an OTF when the exhaust gas was at a volume flow rate of 250,000 m3/h (under normal conditions), strong turbulence developed in the gas ducts, preventing the adequate measurement of the concentrations of each gas mixture component. In order to control the velocity at the fork of the two gas passes, a sliding shutter was installed in one of them at a 45° angle. This resulted in a 4–5-fold reduction in the velocity at the exit of the gas duct. A transient flow mode (Re < 10,000) was formed directly behind the closed part of the shutter, allowing the concentrations of the flue gas components to be measured with the required accuracy and to be controlled using an additional parameter.
In modern conditions, the production of silicon from quartz raw materials in OTFs needs to address energy efficiency issues by taking into account the distribution of gas-fume streams. As a result of the simulation presented here:
  • The contours of the main parameters defining the flow mode in the exhaust gas transfer line, namely the kinematic viscosity and velocity, were obtained.
  • The flow mode was determined by calculating the Reynolds criterion along the exhaust gas transfer line from the OTF to the GCS.
  • It was revealed that the most suitable place for the installation of measuring equipment is directly behind the closed part of the sliding shutter. In this area, there is a transient flow mode with the lowest velocity and lowest Reynolds criterion value. In this location, the flow is influenced by turbulent forces at least, allowing the concentrations of the flow components to be measured with the required accuracy.

Author Contributions

The conceptualization and methodology of this research have been developed by V.B.; the model development and validation have been performed by O.M. As for writing, the original text has been written by O.M.; editing of the article, as well as supervision, has been done by V.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian Science Foundation grant no. 22-29-00397.

Institutional Review Board Statement

The study was conducted according to the guidelines of Saint Petersburg Mining University and approved by the University Review Board (or Ethics Committee) of Saint Petersburg Mining University.

Data Availability Statement

The results was obtained during the authours investigation and were not published in other papers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nedosekin, A.O.; Rejshahrit, E.I.; Kozlovskij, A.N. Strategic approach to assessing economic sustainability objects of mineral resources sector of Russia. J. Min. Inst. 2019, 237, 354–360. [Google Scholar] [CrossRef]
  2. Litvinenko, V.; Tcvetkov, P.; Molodstov, K.V. The social and market mechanism of sustainable development of public companies in the mineral resource sector. Eurasian Min. 2019, 1, 36–41. [Google Scholar] [CrossRef]
  3. Litvinenko, V. Digital Economy as a Factor in the Technological Development of the Mineral Sector. Nat. Resour. Res. 2019, 29, 1521–1541. [Google Scholar] [CrossRef]
  4. Cheremisina, O.V.; Al-Salim, S.Z. Modern Methods of Ananlytical Control of Industrial Gases. J. Min. Inst. 2017, 228, 726–730. [Google Scholar]
  5. Ivanchik, N.N.; Balanovsky, A.E.; Kondratyev, V.V.; Tyutrin, A.A. Research of silicon waste processing products as ultradispersed activating fluxes of arc welding. J. Sib. Fed. Univ. 2018, 11, 155–167. [Google Scholar] [CrossRef] [Green Version]
  6. Nemchinova, N.V.; Tyutrin, A.A.; Sokolnikova, Y.V.; Fereferova, T.T. Analytical Investigations of Silicon Silicon Production Raw Materials and Products. J. Sib. Fed. Univ. 2017, 10, 37–48. [Google Scholar] [CrossRef]
  7. Nemchinova, N.V.; Mineev, G.G.; Tyutrin, A.A.; Yakovleva, A.A. Utilization of Dust from Silicon Production. Steel Transl. 2017, 47, 763–766. [Google Scholar] [CrossRef]
  8. Chigondo, F. From Metallurgical-Grade to Solar-Grade Silicon: An Overview. Silicon 2018, 10, 789–798. [Google Scholar] [CrossRef]
  9. Liu, Y.; Kong, J.; Zhuang, Y.; Xing, P.; Yin, H.; Luo, X. Recycling high purity silicon from solar grade silicon cutting slurry waste by carbothermic reduction in the electric arc furnace. J. Clean. Prod. 2019, 224, 709–718. [Google Scholar] [CrossRef]
  10. Karlina, A.I.; Kondrat’ev, V.V.; Balanovsky, A.E.; Kolosov, A.D.; Ivanchik, N.N. Results of modification of cast iron by carbon nanostructures of gas cleaning dust of silicon production. Advances in Engineering Research. In Proceedings of the International Conference on Aviamechanical Engineering and Transport (AviaENT 2018), Irkutsk, Russia, 21–26 May 2018; pp. 169–173. [Google Scholar]
  11. Vasilyeva, N.V.; Boikov, A.V.; Erokhina, O.O.; Trifonov, A.Y. Automated digitization of radial charts. J. Min. Inst. 2021, 247, 82–87. [Google Scholar] [CrossRef]
  12. Gembitskaya, I.M.; Gvozdetskaya, M.V. Transformation of grains of technological raw materials in the process of obtaining fine powders. J. Min. Inst. 2021, 249, 401–407. [Google Scholar] [CrossRef]
  13. Bartie, N.; Cobos-Becerra, L.; Froehling, M.; Reuter, M. The resources, exergetic and environmental footprint of the silicon photovoltaic circular economy: Assessment and opportunities. Resour. Conserv. Recycl. 2021, 169, 24. [Google Scholar] [CrossRef]
  14. Yolkin, K.S.; Yolkin, D.K.; Kolosov, A.D.; Ivanov, N.A.; Shtayger, M.G. Technologies, which allow to reduce an impact of metal silicon production on the environment. IOP Conf. Ser. Mater. Sci. Eng. 2018, 411, 012028. [Google Scholar] [CrossRef]
  15. Saevarsdottir, G.; Magnusson, T.; Kvande, H. Reducing the Carbon Footprint: Primary Production of Aluminum and Silicon with Changing Energy Systems. J. Sustain. Metall. 2021, 7, 848–857. [Google Scholar] [CrossRef]
  16. Legemza, J.; Findorák, F.; Bul’ko, B.; Briančin, J. New Approach in Research of Quartzes and Quartzites for Ferroalloys and Silicon Production. Metals 2021, 11, 23. [Google Scholar] [CrossRef]
  17. Rozhikhina, I.D.; Nokhrina, O.I.; Yolkin, K.S.; Golodova, V.A. Ferroalloy production: State and development trends in the world and Russia. IOP Conf. Ser. Mater. Sci. Eng. 2020, 866, 012004. [Google Scholar] [CrossRef]
  18. Asaniv, D.A.; Zapasnyi, V.V.; Ermekova, A.T.; Maratova, G.; Ivanov, A.A.; Cherepanov, N.I. Current Status of Dust Collection Systems in Aksu Ferroalloy Plant Smelting Shop 1 and Functional Improvement to These Systems. Metallurgist 2018, 62, 391–400. [Google Scholar] [CrossRef]
  19. Danilov, A.; Smirnov, Y.; Korelsky, D.S. Effective methods for reclamation of area sources of dust emission. J. Ecol. Eng. 2017, 18, 1–7. [Google Scholar] [CrossRef]
  20. Kero, I.; Dahl, S.; Tranell, G. Airborne Emissions from Si/FeSi Production. J. Met. 2017, 69, 365–380. [Google Scholar]
  21. Pashkevich, M.A.; Petrova, T.A. Development of an operational environmental monitoring system for hazardous industrial facilities of Gazprom Dobycha Urengoy. J. Phys. Conf. Ser. 2019, 1384, 012040. [Google Scholar] [CrossRef]
  22. Koteleva, N.; Kuznetsov, V.; Vasileva, N.V. Simulator for Educating the Digital Technologies Skills in Industry. Part One. Dynamic Simulation of Technological Processes. Appl. Sci. 2021, 11, 10885. [Google Scholar] [CrossRef]
  23. Qili, W.; Binbin, J.; Mingquan, Y.; Min, H.; Xiaochuan, L.; Komarneni, S. Numerical simulation of the flow and erosion behavior of exhaust gas and particles in polysilicon reduction furnace. Sci. Rep. 2020, 10, 1909. [Google Scholar] [CrossRef] [PubMed]
  24. Vasilyeva, M.A.; Vöth, S. Multiphysical Model of heterogeneous flow moving along a channel of variable cross-section. J. Min. Inst. 2017, 227, 558–562. [Google Scholar]
  25. Beloglazov, I.; Morenov, V.A.; Leusheva, E.L. Flow modeling of high-viscosity fluids in pipeline infrastructure of oil and gas enterprises. Egypt. J. Pet. 2021, 30, 43–51. [Google Scholar] [CrossRef]
  26. Belolgazov, I.I.; Petrov, P.A.; Martynov, S.A. Application of Production Processes Control Algorithm Using Adaptive Control System. In Proceedings of the International Russian Automation Conference (RusAutoCon), Sochi, Russia, 9–16 September 2018; pp. 1–4. [Google Scholar]
  27. Hamzehloo, A.; Lusher, D.J.; Laizet, S.; Sandham, N.D. Direct numerical simulation of compressible turbulence in a counter-flow channel configuration. Phys. Rev. Fluids 2021, 6, 094603. [Google Scholar] [CrossRef]
  28. Sandham, N.D.; Johnstone, R.; Jacobs, T.C. Jacobs Surface-sampled simulations of turbulent flow at high Reynolds number. Int. J. Numer. Methods Fluids 2017, 85, 25. [Google Scholar] [CrossRef] [Green Version]
  29. Modeling Transient Compressible Flow. ANSYS Fluent Tutorial Guide; Release 18.0; ANSYS, Inc.: Canonsburg, PA, USA, 2017; pp. 267–310. [Google Scholar]
  30. Zhao, Y.; Akolekara, H.; Weatheritta, J.; Michelassib, V.; Sandberg, R. RANS Turbulence Model Development using CFD-Driven Machine. J. Comput. Phys. 2020, 411, 109413. [Google Scholar] [CrossRef] [Green Version]
  31. Mishra, P.; Aharwal, K.R. A review on selection of turbulence model for CFD analysis of air flow within a cold storage. In Proceedings of the 2nd International Conference on Advances in Mechanical Engineering (ICAME 2018), Kattankulathur, India, 22–24 March 2018; Volume 402, pp. 1–9. [Google Scholar]
  32. Kuhnen, J.; Song, B.; Scarselli, D.; Budanur, N.B.; Riedl, M.; Willis, A.; Avila, M.; Hof, B. Destabilizing turbulence in pipe. Nat. Phys. 2018, 14, 386–390. [Google Scholar] [CrossRef]
  33. Salama, A. Velocity Profile Representation for Fully Developed Turbulent Flows in Pipes: A Modified Power Law. Fluids 2021, 6, 369. [Google Scholar] [CrossRef]
  34. Duraisamy, K.; Iaccarino, G.; Xiao, H. Turbulence modeling in the age of data. Annu. Rev. Fluid Mech. 2019, 51, 357–377. [Google Scholar] [CrossRef] [Green Version]
  35. Tania, S.; Ferreira, K.; Arts, T.; Croner, E. On the Influence of High Turbulence on the Convective Heat Flux on the High-Pressure Turbine Vane LS89. Int. J. Turbomach. Propuls. Power 2019, 4, 37. [Google Scholar]
Figure 1. Furnace roof model. (1, 5—exhaust gas outputs; 2, 6—cooling water outputs; 3, 7—concrete refractory; 4—chamotte refractory; 8—exhaust gas input; 9—cooling water input.)
Figure 1. Furnace roof model. (1, 5—exhaust gas outputs; 2, 6—cooling water outputs; 3, 7—concrete refractory; 4—chamotte refractory; 8—exhaust gas input; 9—cooling water input.)
Symmetry 14 00923 g001
Figure 2. Sectional view of temperature and velocity distribution contours of model 1.
Figure 2. Sectional view of temperature and velocity distribution contours of model 1.
Symmetry 14 00923 g002
Figure 3. Pass model. (1, 2—exhaust gas inputs; 3—sliding shutter; 4—exhaust gas output.)
Figure 3. Pass model. (1, 2—exhaust gas inputs; 3—sliding shutter; 4—exhaust gas output.)
Symmetry 14 00923 g003
Figure 4. Velocity and Re distribution contours in the gas passes.
Figure 4. Velocity and Re distribution contours in the gas passes.
Symmetry 14 00923 g004
Figure 5. Re/velocity dependence chart.
Figure 5. Re/velocity dependence chart.
Symmetry 14 00923 g005
Figure 6. Histograms showing the variations in the Re of the exhaust gas mixture.
Figure 6. Histograms showing the variations in the Re of the exhaust gas mixture.
Symmetry 14 00923 g006
Table 1. Initial model data.
Table 1. Initial model data.
ParameterUnitValue
Volume (gas)Nm3/h *1000
Temperature (gas)°C500
Pressure (water)kPa250,000
Volume (water)Nm3/h *250,000
Temperature (water)°C35–45
*—m3/h under normal conditions (0 °C, 100 kPa).
Table 2. Average composition of the waste gases in percentages [6].
Table 2. Average composition of the waste gases in percentages [6].
Mixture ComponentPercent
CO88.6
CO24.81
CH41.42
N22.5
H22.67
Table 3. Average composition of the fumes in percentages [7].
Table 3. Average composition of the fumes in percentages [7].
Mixture ComponentPercent
SiO285.41
Al2O30.46
Fe2O30.30
CaO1.50
MgO1.24
C6.09
Na2O0.08
SO30.16
P2O50.12
K2O0.31
TiO20.02
SiC5.03
Table 4. Model 1 boundary conditions for exhaust gas mixture.
Table 4. Model 1 boundary conditions for exhaust gas mixture.
ParameterInputOutput
Type of boundary conditionsMass-flow inletPressure outlet
Hydraulic diameter, m10.343
Mass flow rate, kg/s29.16-
Gauge pressure, Pa-0
Temperature, °C500-
Re60,889.4236,962.6
Turbulence intensity, percentage10.343
Table 5. Model 1 boundary conditions for cooling water.
Table 5. Model 1 boundary conditions for cooling water.
ParameterInputOutput
Type of boundary conditionsMass-flow inletPressure outlet
Hydraulic diameter, m0.060.06
Mass flow rate, kg/s278-
Gauge pressure, Pa-0
Temperature, °C35-
Re150,075150,075
Turbulence intensity, percent0.060.06
Table 6. Model 2 boundary conditions.
Table 6. Model 2 boundary conditions.
ParameterInputOutput
Type of boundary conditionsMass-flow inletPressure outlet
Hydraulic diameter, m32.7
Mass flow rate, kg/s15.54-
Gauge pressure, Pa-0
Temperature, °C4300
Re117,255.6266,330
Turbulence intensity, percent3.723.35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bazhin, V.; Masko, O. Monitoring of the Behaviour and State of Nanoscale Particles in a Gas Cleaning System of an Ore-Thermal Furnace. Symmetry 2022, 14, 923. https://doi.org/10.3390/sym14050923

AMA Style

Bazhin V, Masko O. Monitoring of the Behaviour and State of Nanoscale Particles in a Gas Cleaning System of an Ore-Thermal Furnace. Symmetry. 2022; 14(5):923. https://doi.org/10.3390/sym14050923

Chicago/Turabian Style

Bazhin, Vladimir, and Olga Masko. 2022. "Monitoring of the Behaviour and State of Nanoscale Particles in a Gas Cleaning System of an Ore-Thermal Furnace" Symmetry 14, no. 5: 923. https://doi.org/10.3390/sym14050923

APA Style

Bazhin, V., & Masko, O. (2022). Monitoring of the Behaviour and State of Nanoscale Particles in a Gas Cleaning System of an Ore-Thermal Furnace. Symmetry, 14(5), 923. https://doi.org/10.3390/sym14050923

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop