Computational & Technology Resources
an online resource for computational,
engineering & technology publications
Civil-Comp Proceedings
ISSN 1759-3433
CCP: 80
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY
Edited by: B.H.V. Topping and C.A. Mota Soares
Paper 2

Air Pollution Modelling using Finite Differences in a Terrain Conformal Coordinate System

N. Sanín and G. Montero

Institute of Intelligent Systems and Numerical Applications in Engineering, University of Las Palmas de Gran Canaria, Spain

Full Bibliographic Reference for this paper
, "Air Pollution Modelling using Finite Differences in a Terrain Conformal Coordinate System", in B.H.V. Topping, C.A. Mota Soares, (Editors), "Proceedings of the Fourth International Conference on Engineering Computational Technology", Civil-Comp Press, Stirlingshire, UK, Paper 2, 2004. doi:10.4203/ccp.80.2
Keywords: wind modelling, mass consistent models, air pollution model, eulerian model, finite differences, accurate time-stepping.

Summary
In this paper a 3D model for pollutant transport in the atmosphere is proposed. We consider a set of coupled convection-diffusion-reaction equations. The convective phenomenon is mainly produced by a wind field [1,2,3] which is obtained from a 3D mass consistent model. In particular, the modelling of oxidation and hydrolysis of sulphur and nitrogen oxides released to the surface layer is carried out by using a linear module of chemical reactions [4]. Dry deposition process is represented by the so-called deposition velocity and it is introduced as a boundary condition. Wet deposition is included in the source term of the governing equations using the washout coefficient [5].

To obtain a numerical solution, first, the problem is transformed using a conformal coordinate system. This allows us to work with a simpler domain in order to build a mesh that provides high consistency finite difference schemes. Then, the convection-diffusion-reaction equations are solved using a high order time discretization which is obtained following the technique of Lax and Wendroff. The model have been tested with a numerical experiment in the Island of La Palma (at the Canary Islands).

Though, in general, the variable vertical spacing, used in our finite difference discretization leads to schemes of first consistence order, some strategies here proposed provides second order schemes. So, the proposed formulation for the convection-diffusion-reaction problem have interesting properties of consistence and stability. The model does not only allow to generate wind maps from the measurements obtained in few stations, but also to obtain the history of a pollution episode for the considered species.

References
1
Montero, G., Montenegro, R., Escobar, J.M. "A 3-D Model for Wind Field Adjustment". J. Wind. Eng. and Ind. Aerd., v. 76-46, pp. 249-261, 1998. doi:10.1016/S0167-6105(98)00022-1
2
Montero, G., Sanín, N. "Modelling of Wind Field Adjustment Using Finite differences in a Terrain Conformal Coordinate System". J. Wind. Eng. and Ind. Aerd., v. 89, pp. 471-488, 2001. doi:10.1016/S0167-6105(00)00075-1
3
Winter, G., Montero, G., Ferragut, L., Montenegro, R. "Adaptive strategies using standard and mixed finite elements for wind field adjustment". Solar Energy, v. 54, 1, pp. 49-56, 1995. doi:10.1016/0038-092X(94)00100-R
4
Winter, G., Betancor, J., Montero, G. "3D Simulation in the lower troposphere: Wind field adjustment to observational data and dispersion of air pollutants from combustion of sulfur-containing fuel", in J.I.Díaz(ed), Ocean circulation and pollution control - A mathematical and numerical investigation, Springer Berlin-Heidelberg, pp. 29-51, 2004.
5
Zannetti, P. Air Pollution Modeling. Comput. Mech. Publ., Boston, 1990.

purchase the full-text of this paper (price £20)

go to the previous paper
go to the next paper
return to the table of contents
return to the book description
purchase this book (price £95 +P&P)