11.10.2013 Views

2. ENVIRONMENTAL ChEMISTRy & TEChNOLOGy 2.1. Lectures

2. ENVIRONMENTAL ChEMISTRy & TEChNOLOGy 2.1. Lectures

2. ENVIRONMENTAL ChEMISTRy & TEChNOLOGy 2.1. Lectures

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!