2. ENVIRONMENTAL ChEMISTRy & TEChNOLOGy 2.1. Lectures
2. ENVIRONMENTAL ChEMISTRy & TEChNOLOGy 2.1. Lectures
2. ENVIRONMENTAL ChEMISTRy & TEChNOLOGy 2.1. Lectures
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Chem. Listy, 102, s265–s1311 (2008) Environmental Chemistry & Technology<br />
In our case, we have inspired in particle visualization<br />
of pollutant concentration. After many experiments of different<br />
ways of visualization, we conclude to the one showed in<br />
Fig. 3. The concentration is expressed by two basic methods<br />
– by color and by particle count. The level of the top concentration<br />
can be adjusted – in Fig. 3. the white means approximately<br />
1 % of source concentration, thus the user can see that<br />
this level of concentration will occurs in the presented scenario<br />
on the ground too. For the better perception of depth, the<br />
further particles are darkened, they do not interfere with the<br />
foreground particles.<br />
A c c u r a c y<br />
The problem of numerical calculation is the stability of<br />
the system and the accuracy of the calculation. Both depend<br />
especially on the size of calculation steps. In our case, when<br />
we have transformed the PDE (3) into the system of ODEs (7),<br />
three calculation steps exist.<br />
The first step Δx is a step along x axis and it is used as<br />
integration step in described experiment. Thus, its size primarily<br />
influences the accuracy of obtained results.<br />
Last two steps Δy and Δz discretize y and z variables and<br />
they have impact on the size of ODEs system (7). These steps<br />
subdivide the space in lateral directions and both are used<br />
for approximation of the second derivatives appearing in<br />
PDE (3). Therefore, they influence both the stability and the<br />
accuracy and they thus indirectly affect the size of Δx.<br />
Fig. 4. The absolute error of the numerical calculation. The<br />
vertical axis shows the error, the horizontal axis shows the distance<br />
from source along x axis (in wind direction)<br />
In case of our numerical solution of A-DE there exist<br />
specific properties of the calculation behavior, as you can<br />
see from our error measurement shown in Fig. 4., where the<br />
absolute error is shown. The error evolution is shown along<br />
the wind direction. More specifically the error is the mean<br />
error in all points (one plane) of the specific x distance.<br />
s389<br />
The greatest error along x axis has been measured close<br />
to the source, which is caused by definition of the initial condition<br />
(7a). Fig. 4. shows other peak of error which is around<br />
0.7 m distance from the source. That is the place where the<br />
plume reaches the ground. The boundary condition (7e) which<br />
is the approximation of boundary condition (3e) is the main<br />
reason of this error existence. It must be noted that the numerical<br />
calculations were stable in spite of measured errors.<br />
Conclusion<br />
The numerical method of A-DE solution has been proposed<br />
and implemented and the results have been presented.<br />
In addition, the accuracy of our model has been verified by<br />
comparison with the analytical solution.<br />
Presented method is relatively simple and easy to implement,<br />
therefore, it can be relatively easy extend to more<br />
general form of the A-DE which will be the main goal of the<br />
future project progress. The extensions could be the general<br />
wind direction and speed, the non-flat ground with obstacles<br />
(trees, buildings etc.), non-stationary point source/sources<br />
etc. Moreover, other atmospheric parameters such as temperature<br />
and humidity and substance chemical properties will<br />
be added to the model for even more physically and chemically<br />
correct behavior.<br />
This preliminary work will be served as a base for more<br />
sophisticated model that will be a part of the intelligent system<br />
for human protection against consequences of industrial<br />
accidents and its analysis. To do that many experiments and<br />
measurements of real substance outflows, gas dispersion etc.<br />
will have to be done.<br />
REFEREnCES<br />
1. Builtjes P. J. H.: Air Pollution Modeling and Its Application<br />
XIV (Gryning S.E., Schiermeier F.A., ed,), p. 3,<br />
Major Twentieth Century Milestones in Air Pollution<br />
Modelling and Its Application, Springer US, 2004.<br />
<strong>2.</strong> Reynolds O.: Philos. Trans. R. Soc. London, Ser. A ,<br />
1895, 123.<br />
3. Ermak D. L.: Atmos. Environ. 11, 231 (1977).<br />
4. Jacobson M. Z.: Fundamentals of Atmospheric Modeling.<br />
Cambridge University Press, new York 2005.<br />
5. Slanco P., Bobro M., Hanculak J., Geldova E.: Acta<br />
Montanistica Slovaca,, 313 (2000).<br />
6. Roussel G., Delmaire G., Ternisien E., Lherbier R.:<br />
Environmental Modelling & Software 15, 653 (2000).<br />
7. Pagendarm H. G.: Visualization and Intelligent Design<br />
in Engineering and Architecture, p. 315, Scientific Visualization<br />
in computational fluid dynamics, Computational<br />
Mechanics Publications, 1993.