^{1}

^{1}

^{*}

^{2}

The dispersion is mainly governed by wind field and depends on the planetary boundary layer (PBL) dynamics. Accurate representation of the meteorological weather fields would improve the dispersion assessments. In urban areas representation of wind around the obstacles is not possible for the pollution dispersion studies using Gaussian based modeling studies. It is widely accepted that computational fluid dynamics (CFD) tools would provide reasonably good solution to produce the wind fields around the complex structures and other land scale elements. By keeping in view of the requirement for the micro-scale dispersion, a commercial CFD model PANACHE with PANEPR developed by Fluidyn is implemented to study the micro-scale dispersion of air pollution over an urban setup at Indira Gandhi Centre for Atomic Research (IGCAR), Kalpakkam a coastal station in the east coast of India under stable atmospheric conditions. Meso-scale module of the PANACHE model is integrated with the data generated at the site by IGCAR under RRE (Round Robin Exercise) program to develop the flow fields. Using this flow fields, CFD model is integrated to study the micro-scale dispersion. Various pollution dispersion scenarios are developed using hypothetical emission inventory during stably stratified conditions to understand the micro-scale dispersion over different locations of coastal urban set up in the IGCAR region of Kalpakkam.

The study of air pollution dispersion is pre-requisite for urbanization, industrial growth and expansion of coastal population [1,2]. The air pollution dispersion studies at various scales gives the information about local and long range scenarios impact [

Turbulence is classically calculated using k - ε losure methods to calculate the isotropic eddy viscosity parameter present in both the momentum and pollution transport equations, which assumes that a pollutant is diluted equally in all directions [

A meso-domain is chosen in the geographical coordinates of 12.2818˚N to 13.3200˚N and 79.5043˚E to 80.6603˚E. The size of the domain is 145 km × 110 km × 0.1 km. The meso-scale study domain is depicted in

Two CFD models namely PANACHE and PANEPR have been employed in the present study [

equations governing mass, momentum, and energy transfer (http://www.fluidyn.com/fluidyn/). It solves the mass, momentum, and energy conservation equations for both laminar and turbulent flow [19-22]. The equations for the two cases differ primarily in the form and magnitude of the transport coefficients (diffusivity, viscosity, and thermal conductivity). Eddy hypothesis is used to compute the turbulent contribution to the exchange coefficients.

The continuity equation for species m is

;

ρ_{m} = mass density of species m;_{}

ρ = total mass density;

u = fluid velocity;

D_{m }= diffusion coefficient for species m in air;

y_{m }= mass fraction of species;

m = ρ_{m}/ρ;

ρ_{m}^{s} = source term for species m due to pollutant emissions;

ρ_{m}^{d} = source term for species m due to dry deposition in canopies;

ρ_{m}^{p }= source term for species m due to droplet evaporation/condensation;

ρ_{m}^{c}^{ }= source term for species m due to chemical reactions;

Sc_{l} = Schmidt number for laminar flow;

Sc_{t} = Schmidt number for turbulent flow;_{}

µ_{l} = molecular viscosity of air;

µ_{t} = turbulent eddy viscosity.

We get the continuity equation for total fluid density

ρ^{s} = total mass source due to pollutant emission;

ρ^{d} = total mass source due to dry deposition in canopies;

ρ^{p} = total mass source due to droplet Evaporation and condensation.

The momentum equation for the fluid mixture is

P = Pressure;

σ = Newtonian viscous stress tensor

= (4)

μ, λ = First and second coefficients of viscosity;

λ = −2/3μ;

T = transpose;

I = unit dyadic;

A_{0} = 1 when the k - e or the k - L turbulence model is active, = 0 otherwise;

k = turbulence kinetic energy;

F^{s} = rate of momentum gain per unit volume due to Pollutant emissions;

F^{g} = force due gravitational acceleration;

F^{d} = force due to plant canopy effect, dry deposition, etc;

F^{p} = force due to interaction with droplets/particles.

The internal energy equation is

I = specific internal energy;

J = heat flux vector

=(6)

k = thermal conductivity;

T = Temperature;

h_{m} = specific enthalpy of species m;

Q^{s} = rate of specific internal energy gain due to pollutant emissions;

Q^{d} = rate of specific internal energy gain due to dry depositions in canopies;

Q^{p} = rate of specific internal energy gain due to interaction with particles;

Q^{h} = rate of specific internal energy gain from surface energy budget;

Q^{c} = rate of specific internal energy gain due to chemical reactions;

Dissipation = σ. u for laminar flow

= ρε for turbulent flow.

In the study domain, several in situ meteorological measurement stations, i.e. masts and micro-meteorological towers are in operation during the period of the study as a part of the RRE conducted by IGCAR. To start simulation, model requires initial value at all grid points. For this purpose the available meteorological parameters (e.g. wind velocity, direction, and temperature) at 2, 8, 16, 32, and 50 m elevations from ground were used. Geographical coordinates in terms of longitude and latitude of all stations are specified at appropriate locations in the developed terrain fields. The CFD model generates wind fields with the given boundary conditions by solving Navier Stokes equations at every time step. The hourly wind observations then ingested along with the CFD winds and with the objective analysis method, the model generates final wind vectors at every time step. The wind field was generated for the entire mesoscale domain for the period 15 to 22 September 2010. The model generated wind fileds at various heights viz. 2, 50 and 100 m at different times in the meso domain, and wind fields at 2 m height for 15 September 2010 is depicted in

Using the nesting option, a micro-scale domain is chosen incorporating urban location in the mesodomain in the IGCAR area. The domain is divided into various grids and the number of grids in the three directions (x-, yand

z-directions) is 30 × 30 × 12. The grid resolution in xdirection is around 66 m, y-direction is 66 m and in zdirection is about 8 m. The model generated mesodomain wind vectors serve as boundary conditions for the PANACHE-PANEPR CFD model. The model is run on daily basis using the wind solver module. As explained above wind fields will be generated at every time step and hourly wind fields are stored for the analysis. The model generated wind fields at various heights viz. 2, 50 and 100 m at three times in the nighttime (e.g. 0100, 0300 and 0500 LST hours) in the micro domain, and for 2 m height for 15 September 2010, is depicted in

The micro-scale dispersion of various pollutants over a coastal region of south-east coast of Tamil Nadu has been studied. The micro-scale domain chosen for the present study is depicted in ^{*} sign) and green field (named as HEP3 with R sign).

Numerical experiment has been made to study the dispersion using two arbitrary point sources with steady simulation condition. In the steady simulation solution marches with iterations (cycles) only and the simulation time remains unchanged from the initial defined value. During this type of simulation, following things will not be taken into consideration:

• Duration, Reference time and output dumping time;

• Transient weather data;

• Transient pollutant emissions;

• Time-averaged concentrations.

Two locations with the hypothetical emission inventory are considered, namely HEP1 and HEP2 (_{2}) respecttively with mass flow rate 3 kg·s^{–}^{1}, temperature of 25˚C, release duration of 180 s, and exit velocity of 1.5 ms^{–}^{1}. Similarly, the concentration at location HEP2 are 30% of CO and 70% of SO_{2} respectively with mass flow rate 2 kg·s^{–}^{1}, temperature of 25˚C, release duration of 180 s with a delay of 2 minutes from the HEP1, and exit velocity of 1 ms^{–}^{1}.

Using the dispersion solver module using the microscale wind field and given emission inventory, the dispersion at various stages of model integration are generated. Basically, these are ground level concentrations and the model implemented in PANACHE only. After 890 iterations the dispersion of CO from the two locations are shown in _{2} as well. The release of CO and SO_{2}, at HEP2 was occurred after 2 minutes from HEP1. But that are effect much in terms of dispersion. This could be attributed to the presence of dense urban structures present in HEP2 location. These structures will absorb solar radiation during daytime and release or re-radiate in the nighttime. They act like a secondary source of energy and will generate local thermal regime within the surrounding stable environment. This will generate a local turbulent mixing leading to more dispersion. Even though HEP1 is also an urban setup but with scattered structures, the intensity of reemission of radiation is less compared to HEP2 leads to comparatively less dispersion.

In this case, attempt was made to study the dispersion at three arbitrary point sources with unsteady simulation condition. Unsteady simulation solver uses an iterative method to march in time with pre-specified iterations are used to solve the governing equations before proceeding

to the next time-step. Three Point sources HEP1, HEP2 and HEP 3 are shown in _{2}) respectively with mass flow rate 2 kg·s^{–}^{1}, temperature of 25˚C, release duration of 120 s, and exit velocity of 2 ms^{–}^{1}. Concentration at HEP2 are 90% of CO and 10% of N_{2} respectively with mass flow rate 2.5 kg·s^{–}^{1}, temperature of 25˚C, and exit velocity of 2.2 ms^{–}^{1}. The concentration at HEP3 are 50% of CO and 50% of N_{2} respectively with mass flow rate 1.5 kg·s^{–}^{1}, temperature of 25˚C, release duration of 60 s with a delay of 30 s from HEP1 and HEP2, and exit velocity of 1.8 ms^{–}^{1}. For this case PANACHE is implemented with wind and dispersion solvers which will provide ground level concentrations only. The dispersion of CO and N_{2} from the three locations are shown in

In this case PANEPR module is implemented to study the dispersion scenarios of three arbitrary sources (as given above) with respect to time at different heights 8, 50 and 100 m. Similar to that of Case 2, wind solver and at HEP1 and HEP2 are at an altitude of 10 m and HEP3 is at an altitude of 5 m from the ground. The concentrations of pollutants for 15 September 2010 at 0100 hour at HEP1 are 70% of CO and 30% of SO_{2} respectively with mass flow rate 2 kg·s^{–}^{1}, temperature of 25˚C, release duration of 3600 s, and exit velocity of 1.8 ms^{–}^{1} and concentration at HEP2 are 30% of CO and 70% of SO_{2} respectively with mass flow rate 2.5 kg·s^{–}^{1}, temperature of 25˚C, release duration of 3600 s and exit velocity of 1.7 ms^{–}^{1}. The concentration at HEP3 are 50% of CO and 50% of SO_{2} respectively with mass flow rate 1.5 kg·s^{–}^{1}, temperature of 25˚C, release duration of 3600 s, and exit velocity of 2.1 ms^{–}^{1}. The dispersion scenarios generated at these arbitrary sources at different heights for CO and SO_{2} are given in Figures 7 and 8 respectively. Two specific scenarios are studied one is after 5 minutes 13 seconds and another one after 10 min 54 seconds. Interestingly, both the pollutants completely dispersed and not present at 100 m height of HEP3 location even after 5 minutes 13 seconds. This could be attributed to higher magnitude of wind of the order of 6 ms^{–}^{1} at those heights are noticed and the emission point is at 5 m height only. But at 8 and 50 m heights after the same time of integration, the dispersion is slightly more spread in HEP2 location compared to HEP1 and HEP3. After 10 minutes and 54 seconds, all the three sources plumes merged together.

The concentration and spread is more at 8 and 50 m compared to 100 m height. The main purpose of this exercise is to assess the capability of CFD model PANACHE and PANEPR to generate the scenarios for multiple sources of release of any pollutant of known emission inventory and weather information, the dispersion rate with reference to time, space (both horizontal and vertical) for the hazard impact studies. After examining all the above case studies, it is clearly seen that pollution dispersion over urban area is represented well in this CFD model, where Gaussian models does not serve well.

This is due to non representation of turbulent flows and wake flows around the obstacles in Gaussian models where as CFD models can approximate to better accuracy of these turbulent flows and hence with efficient dispersion solver, one can get reasonable picture of pollutant dispersion._{}

Micro-scale dispersion of air pollution over an urban setup in a coastal station in the east coast of India (IGCAR, Kalpakkam) has been studied using a commercial CFD model PANACHE with PANEPR developed by Fluidyn. Winds are generated in the meso-scale domain in the east coast of Tamilnadu, and later these winds are used as boundary conditions to generate the wind field over the urban set up of Kalpakkam region. Locations in the micro-domain (IGCAR, Kalpakkam) with varied land scale such as dense building areas (urban set up), scattered building areas and open fields are selected. Dispersion during stably stratified conditions, i.e. 0100, 0300 and 0500 LST hours during the study period are evaluated. At these locations, hypothetical emission inventories are provided and the CFD model with dispersion module is implemented and various scenarios of pollution dispersion are generated during stably stratified conditions. The dispersion of air pollutants are more over the urban domain due to its reemission of absorbed radiation during daytime. This is significantly changing the thermal structure of the PBL and the mixing which is different from the open fields. The model results also reveal the dispersion of pollutants from individual source can be visualized on horizontal as well as vertical scale with time. This study has given confidence that PANACHE and PANEPR can be used for understanding pollution dispersion by providing actual emission inventories and can be validated reasonably with the available air quality data. The model results suggest that this CFD model has generated reasonable dispersion scenarios over urban scale region which is having profound utility in pollution hazard planning management in major cities and important industrial locations.

We would like to express our sincere thanks to IGCAR, Kalpakkam providing data of the meteorological towers and Automatic weather stations which was collected during RRE programme, and Fluidyn, Bangalore for providing all facilities and access to use their modeling system in conducting the work. Mr. Srikanth gratefully acknowledges Indian Institute of Technology Kharagpur for providing assistantship to conduct the present study.