METHOD FOR PREDICTING PARAMETERS OF THE AEROIONIC MODE IN OPEN TERRAIN GROUND AREAS

Purpose. The scientific work is devoted to the development of a new method for forecasting aeroionic mode in working areas at industrial sites, taking into account ionizing sources and surrounding obstacles. Methodology. To find the concentration of positive aeroions, dust and negative aeroions, we use 3D mass transfer equations that take into account the rate of recombination of ions having different polarity and the rate of recombination of ions with dust particles. The numerical solution is based on the integration of the mass transfer equations using the finitedifference method, which turns out to be stable for any step in time. To determine the components of the air velocity vector, a three-dimensional model of the potential flow is used, where the Laplace equations for the velocity potential are the modeling equations. Findings. The mathematical method of numerical calculation of the concentration of positive, negative aeroions and dust has been developed. A feature of the method is the possibility of predicting the aeroionic mode, taking into account all physical factors that significantly affect the formation of concentration zones of aeroions in working areas at industrial sites. The method is not tied to a particular industrial site, it allows us to estimate the value of the concentration of aeroions both locally and in the entire calculated region. Originality. A method for prediction of aeroionic mode in working areas is developed based on 3D modeling of the propagation of negative, positive aeroions and dust under the influence of wind and diffusion, which allows to obtain results at each point of space or in a specific cross-section. Practical value. The proposed method of forecasting was used to solve the problem of estimating aeroionic mode in industrial zones in the open area of the industrial site of the Dnipro oil extraction plant in the presence of emission sources: positive aeroions during the operation of vehicles and respiration of workers; dust during the movement of workers and vehicles; negative aeroions of the ionizer installed in the working zone. The regularities of changes in the concentration of aeroions of various polarity and dust at a height of 1.7 m are determined, which corresponds to the position of the respiratory organs of workers. Quantitative results are needed in assessing the permissible working conditions in the workplaces of industrial sites of enterprises when creating new jobs and reengineering existing ones.


Introduction
Providing a certain degree of ionization of the air environment in the working area is one of the most important factors in maintaining a good state of health and high efficiency of staff. It has been established that a significant decrease in the content of charged particles (ions) in the air coincides with the appearance of unusual sickliness in working people, fatigue complaints, depression, nausea, insomnia, irritability, respiratory disorders, and others. At the same time, the presence of people in conditions with moderately high ionization of the atmosphere, with the vast majority of negative ions -on the contrary, there is a beneficial effect on the body. The process of ionization of air consists in the transformation of neutral atoms and air molecules into electrically charged particles (ions). Ions in the air (aeroions) can be formed as a result of natural and artificial ionization.
Natural ionization occurs as a result of the influence on the air of space and solar radiation and ionizing radiation (particles) that arise during the decay of long-lived radionuclides of the earth's crust (potassium-40, uranium-238, thorium-232, etc.). Natural ionization of the air environment is everywhere and constant in time. Artificial aerionization occurs under the influence of ionizing factors, which accompany some technological processes (X-ray and ultraviolet radiation, thermo-emission, photoelectric effect, etc.); as well as the effect of special devicesionizers. The degree of ionization of air depends on the ratio of ionization and deionization processes. The latter is due to the recombination of two ions of different polarities, adsorption of light ions on uncharged nuclei, condensation, neutralization of light and heavy ions by charges of the opposite sign, and others.
Assessment of the aeration degree of the environment is carried out by comparing the measured values with the standard ones. Artificial ionizers (high-voltage, induction, radiation, etc.) and efficient, well-organized ventilation of premises are widely used to normalize the ionic mode of the air environment, as the external clean air contains 2-5 times more ions than the air of indoor premises (50-100 ions/cm 3 ).
Many works of domestic authors are devoted to the study of aeroionic mode. Modeling of aeroionic mode in working zones under conditions of artificial ionization is presented in the work [1]. The problem of the evaluation and the relationship between aerosol pollution and the aeroionic composition of air in the working zone is considered in the works [8][9]. The works [2,11] are devoted to the modeling of the aeroionic mode in premises with artificial ionization using the means of numerical Modeling. The proposed program information system for the optimization of the microclimate in the premises is in [6]. The study of ways to improve the aeroionic mode of the working environment is presented in work [7]. The study of the microclimatic parameters influence of air exchange on the aeroionic composition of the air in the working premises and the modeling and prediction of the aeroionic composition dynamics of air are given in the works [4][5]. The methods for determining the concentration of aeroions and dust content are presented in the work [3]. The normalization of the aeroionic composition of air in office premises is considered in the work [15]. A large number of foreign works are devoted to the study of ionization processes for improving the health of employees [12][13][14]16].

Purpose
The purpose of this study is to develop a new method for forecasting aeroionic mode in working areas at industrial sites, taking into account ionizing sources and surrounding obstacles (buildings).

Methodology
It is known that the formation of the concentration field of aeroions is influenced by a complex of physical factors: airflow, diffusion, electric field effect and others. In addition, there is the interaction of ions of different polarity and their interaction with particles of dust. To account for these processes during the Modeling of aeroions dispersion, the 3D equation of mass transfer is used. It is written for each component (positive aeroions, dust, negative aeroions). The Modeling equations have the form: to describe the process of dispersion of positive ions for modeling the transfer of dust to describe the process of dissipating negative ions The denotations of physical parameters in these equations are as follows: For numerical implementation such a boundary condition is recorded: where ( 1, ) the values of ions in the last boundary cell, ( , ) p A i jthe value of ions in the penultimate difference cell.
On the ground surface and all solid surfaces, the boundary condition for non-permeability is set: This condition physically means that the flow of ions at these boundaries is zero.
Using the Modeling equations (1) -(3) speed components ( , , ) u x y z , ( , , ) v x y z , ( , , ) w x y z must satisfy the equation of continuity: To fulfill this condition one should solve an aerodynamic problem for determining the velocity field. To determine the components of the air velocity vector, a three-dimensional model of the potential flow is used, where the Laplace equation for the velocity potential serves as the simulating equation: where Pspeed potential.
Components of the velocity vector of the air are related to the velocity potential by the following relationships: For the Laplace equation, the following boundary conditions are set: on firm boundaries the condition of non- on the boundary of the flow outlet from the calculated region the Dirichlet boundary condition is set P const  ; at the boundary of the air flow inlet into the calculated area, the boundary condition is set When simulating the propagation of aeroions, it is necessary to take into account the influence of the electric field on this process. Since the aeroions have a charge, they generate an electric field E. For the Modeling of the electric field, the following equation is used [2,12]: where 0 dielectric permeability; e qvolumetric charge density.
From equation (5.7) we can go to the scalar potential, if we take into account such dependence: Then we obtain the Poisson equation: p A x y zconcentration of positive aeroions; scalar potential; e elementary charge. On the basis of this equation electric field Modeling is carried out.
As the boundary condition, the condition of electrical insulation is used: where nunit vector of the external normal to the solid surface.
For numerical integration of the mass transfer equations, a difference scheme is used based on the geometric and physical splitting of the simulating equations (1) - (3). In each of the four fractional steps, only one direction of perturbation transfer is taken into account; the fifth split step takes into account the change in the concentration of aeroions due to the action of the emission source. The calculation is performed on a rectangular difference grid. Concentration of quantities p A , n A , d A is determined in the centers of the difference cells, whereas the components of the velocity vector of the air are set at the boundaries of the difference cells [2,[11][12].

Findings
To estimate the aeroionic mode in the working zones, based on the developed method, it is necessary to determine the aeroions emission intensity from the specific sources. As already noted, one of these sources is the worker, who exhales positive aeroions. Physical experiment was conducted to obtain scientifically substantiated information regarding the magnitude of such emission.
For the experiment, a group of people of different ages was selected: 16 people in each age group. The mouth area was measured S=9.5 cm 2 , the rate of air flow when exhaled was on average  Table 1.
Experiment to measure the concentration of aeroions from human exhalation was carried out using the KT-401 AIRION TESTER. The data obtained on the basis of the conducted physical experiment, were used during computational experiments. During the experiments it was taken into account that the source of positive aeroion emission during exhalation had a periodic action, namely, exhalation for 1 s in quantity Q , the inhalationemission is absent and after 2 s again exhalation with the same emission. That is, the source of periodic emission () Qt is known. Verification of the developed forecast method was carried out. A test task was solved to calculate the emission of positive aeroions with exhaust gases of vehicles.
During the physical experiment, the positive aeroions at different distances from the Daewoo Matiz exhaust pipe were measured using the KT-401 AIRION TESTER. The velocity of wind and gas flow at the exhaust pipe outlet was measured using an anemometer PM 6252 B Digital Anemometer. The speed of the air flow was Vair=1.2-1.4 m/s, and the speed of gas flow Vgas=1.7-1.9 m/s, diameter of the exhaust pipe hole d=32 mm.
The method of conducting the experiment was as follows: the concentration of positive aeroions was measured at a height of 320 mm from the ground surface at different distances from the exhaust pipe. To obtain data on the intensity of emissions, the calculation was made according to the ratio haust pipe, Sarea of the outlet cross section of the exhaust pipe hole, Cthe concentration of positive aeroions in the outlet cross section. Immediately during the experiment parameters gas V and C were measured at a distance of 0.5 cm from the outlet cross section of the exhaust pipe hole. The obtained data were used as input information for numerical calculation on the basis of the developed method. Using this method, the impurity concentration values at the points of measurements were also calculated. Data analysis shows the adequacy of the numerical calculation method. The corresponding deviation is due to the flow turbulization caused by the influence of the car body on the aerodynamics of the flow.
The source of positive aeroions on the street is primarily people, cars. During calculations, a complex of factors was taken into account, which had been overlooked in scientific publications: aeroion emission from cars; aeroion emission from humans; the influence of various obstacles in the street on the aeroionic mode formation; the influence of the uneven wind velocity field on the formation of aeroionic mode; the influence of the uneven field of atmospheric diffusion on the formation of the aeroionic mode; modeling of aeroions dispersion in a three-dimensional formulation (3D modeling).
Taking into account the listed factors extremely complicates the solving of these problems: The task solution is based on three stages: 1. Modeling the formation of a field of concentration of positive, negative aeroions and dust under the influence of wind and atmospheric diffusion.
2. Modeling the transport of aeroions under the action of an electric field.
3. Change the field of concentration of aeroion, dust due to their interaction.
The developed method of forecasting was used to solve the task of estimating the level of pollution of working areas in the open air in the presence of a source of emissions: positive aeroionsemissions from vehicles and people; dusttraffic; sources of negative aeroionsan ionizer installed in the working zone and the screen (fence), which changes the aerodynamics of the aeroions and dust. As a research object, the Dnipro Oil Extraction Plant was chosen, namely, the territory where the cargo scales ( Fig. 1) operate, where the lorries with seeds arrive. This area is a source of dust, positive aeroions, and workers responsible for unloading cars are constantly under the influence of these negative factors.
Scheme of the calculated area is presented in Fig. 2 There are the results of the computational experiment for calculating the field of positive, negative aeroion and dust concentration in the working area of the industrial site (Fig. 3, Fig. 5, Fig. 7). The speed of the air flow was set at 5 m/s, the size of the calculated area was 8x4 m, the intensity of the sources of aeroions was 1000 ions/cm 3 Fig. 3 it is evident that at the distance of 3 m from the beginning of the calculated area there are no negative aeroions, since there is no source of emission on this site, and on another section a complex picture of the distribution, which is caused by the movement of the air flow and the action of atmospheric diffusion, is observed.
The value of the concentration at the height H = 1.7 m slowly increases under the influence of the emission source and diffusion (Fig. 4).  The dust concentration at the height of H = 1.7 m slowly increases under the influence of the emission source (motor vehicles), then it is reduced due to the screen action and decreases at insignificant distance, but under the influence of the source of emissions (workers) increases again and decreases due to changing the flow direction near the awning of the temporary building and diffusion (Fig. 6). Fig. 7 shows the distribution of positive aeroions. Apparently, in the working zone there is a dynamics of growth of positive aeroions, which is caused by the transport of positive aeroions from vehicles (source of emissionscars) and the emission of positive aeroions from people in the working zone, that is, the interaction of sources of aeroion emissions takes place.
As the dust concentration values, the concentration of positive aeroions slowly increases at the height H = 1.7 m under the action of the emission source of positive aeroions (vehicles), then it is reduced by the screen action and decreases at insignificant distance, but under the influence of the source of emissions (workers) increases again and decreases due to changing the flow direction near the awning of the temporary building and diffusion (Fig. 8).
The optimal level of positive aeroions 1500-3000 ions/cm 3 , the maximum allowable level is 50000 ions/cm 3 , in the working area is observed exceeding the permissible level by 10%.