Contribution of operations to the particulate generation (Ghose & Majee, 1998)

## 1. Introduction

Dust generation and its dispersion has been the major concern in ambient air quality in deep cavities such as open cast mines. Major mining activities in deep open pit coal mines range from exploration to the processing of end product that primarily contribute particulate matter (PM), dominantly PM_{10} leading to the problem of air pollution and related health hazards (Falk and Jurgelski, 1979; Pless-Mulloli et al, 2000). As a result, the simulation of dispersion characteristics within the deep open pit coal mines becomes essential to analyze the complex wind flow patterns that significantly affect the dispersion of PM_{10}. The dispersion equations developed within the deep pit boundary provide a reasonable accurate estimate of PM_{10} dispersion within the near field region of the deep open pit coal mines (Silvester et al, 2009). The fundamental equations of continuity and momentum describe the in pit dispersion mechanisms within the atmospheric boundary layer (ABL). In addition, the meteorological conditions within the deep open pit coal mine are significantly affected by temperature (stability) and roughness conditions which ultimately generate complex dispersion phenomenon including separation of atmospheric boundary layer, recirculation, resuspension and settling of PM_{10} (Bitkolov, 1969; Grainger and Merony, 1993)_{.} However, the in-field measurements of PM_{10} within the deep open pit coal mines are constrained by safety regulations, complex geometry of the pit, uncontrolled wind flows and different operation types that make it extremely difficult to carry out monitoring of the PM_{10} conventionally (Roy et al, 2011). Further, the simulation of dispersion characteristics using the fundamental governing equations may require modifications to incorporate the in pit microclimatic effects on the flow regimes (Markov et al, 1978; Aloyan et al, 1982). Therefore, it is essential to analyze and evaluate micro-climatic parameters including the wind turbulence and shear in order to simulate the dispersion of PM_{10} (Turner, 1994). This chapter presents a comprehensive description of the dispersion mechanisms in the deep open pit coal mines considering the topographic, thermal and meteorological factors.

## 2. Dispersion mechanism in deep open pit coal mine

Dispersion mechanisms consist of *diffusion* and *advection* processes. The atmospheric motions *transport* and *diffuse* the pollutants that are released from sources. Richardson (1926) has investigated the dispersion characteristics by using tracers. The fundamental scaling and power laws have been studied by Taylor (1921) and Roberts (1923). The temporal and spatial scales of motion serve to disperse pollutants in the atmosphere by mixing and thus lowering the ambient pollutant concentrations (Turner, 1994). Mikkelson (2003) and Hanna et al (1982) have discussed various aspects of dispersion and associated parameters that are primarily influenced by source dimensions, buoyancy and momentum of the release, roughness and surrounding topography, atmospheric stability and large scale differential heating. In case of PM_{10} its dispersion is additionally affected by intertial and dynamic forces acting on the particulates acting on them. Open cast coal mining involves various operations e.g overburden removal, drilling, blasting, mineral loading, haulage and unloading that generate particulates due to various mechanisms. It has been observed that out of total particulates generated, the PM_{10} constitute one-third to half (Ghose & Majee, 2007; Trivedi et al, 2009). The operations and their contribution to the generation of the particulates have been listed in table 1.

Operation | % of particulates generated |

Overburden removal | 7% |

Top soil removal | 1% |

Coal processing | 72% |

Coal extraction | 3% |

Wind erosion | 17% |

## 3. Flow dynamics in deep open pit coal mines

The flow in the deep pit coal mine is a special case of a deep cavity flow where shear stresses are dominant (Tani et al, 1961; Chang, 1970). It depends significantly upon the state of the boundary layer ahead of the cavity, pressure and forces due to the flow. Sometimes, vortices are formed in the cavities due to the deflection of part of the separated boundary layer. Cavities are generally classified according to their length to depth (L/D) and length to width (L/W) ratios, respectively. A cavity is deep, if L/D < 1 ; and shallow, if L/D > 1; two dimensional if L/W <1 ; and three dimensional if L/W>1 (Larcheveque et al, 2003). The cavity flows are further classified into *open type* (L/D <10) and *closed type* (L/D>13) according to the *aspect ratio*. In the closed type cavity, the shear layer generated at the leading edge collides with the cavity floor. The layer is reflected from the floor forming the expansion waves and the flow escapes the trailing edge. Therefore, two small separate zones are formed in the cavity. In the open type cavity, the free stream shear layer is reattached to the trailing edge and divides the flow into internal and external flow (Woo et al, 2008). Table 2 describes different types of cavities and their characteristics.

Sl no | Property /parameter | Shallow | Intermediate | Deep |

1 | Recompression layer | Not Investigated | Not Investigated | Decreases with increased L/D ratio |

2 | Boundary layer and its influence | Not Investigated | Not Investigated | Bridges the cavity and has little influence over the cavity flow parameters |

3 | Redevelopment length | Redevelopment is small | Redevelopment is maximum | Redevelopment is small |

4 | Variation of coefficient of pressure | Low | Maximum | Low |

5 | Boundary layer thickening | More | Less | More |

6 | Free shear layer flow | Less strongly distributed | More strongly distributed | Less strongly distributed |

7 | Drag Coefficient | High | Medium | Low |

8 | Flow type | unsteady | Not Investigated | Dead/laminar |

9 | Flow reattachment | Occurs on the floor of the cavity | Not Investigated | Occurs on the rear shoulder of the cavity |

10 | Reynolds number | Large | Medium | Small |

11 | Velocity | Large | Medium | Small |

12 | Skin friction | Large | Medium | Small |

13 | Base pressure | First decreases and then increases | Not Investigated | Constant in the middle and slightly irregular in the front and rear portions |

14 | Three dimensionality of shear layer flow | Slightly | Strongly | Slightly |

### 3.1. The cavity regions

The cavity region contains different zones where the boundary layer undergoes various changes (Figure 1). The *Upstream zone* comprises the developing boundary layer zone which is upstream of the front corner of the notch with a finite approaching boundary layer. *Expansion zone* is a small zone just downstream of the separation corner where the pressure falls to the base pressure level. In this zone, the inner 10-15% of the approaching boundary layer adjusts rapidly to expansion at the separation corner and forms a new viscous sublayer slightly below the separation corner. *Mixing zone* is in the initial part of the free shear layer. The flow is similar to a separated jet type wake and lies between the isoenergic free stream flow and the dividing stream line. This intermediate buffer zone acts as a coupling between the outer free stream flow and the inner dissipative flow and brings about the transport of mass and momentum from the former to the latter. The pressure is nearly constant in this zone. This region is associated with large velocity gradients and substantial viscous activity. After the constant pressure mixing region, the shear layer negotiates the pressure rise which continues up to the reattachment point where a apart of the shear layer gets reversed and forms a recirculating cavity flow. The recirculating cavity flow lies between the separation point and the reattachment point, having a low Reynolds number and a constant static pressure. This zone is called as c*ompression zone*. This recirculating cavity flow has also been termed as separation wake. After the reattachment point, the shear layer negotiates the further pressure rise which continues upto the free stream value. This is called as the *recompression zone*. The flow in the ramp corner as well as in the recompression zone is assumed to be isentropic and invicid. Further, the recompression zone entrains an unsteady vertical flow, which is known as the *recompression wake*. After the recompression region the flow forms a new uniform stream with a *redeveloping zone* which becomes fully developed at some distance downstream. Various factors affecting the flow dynamics in the case of deep open pit coal mines have been discussed in subsequent sections.

#### 3.1.1. Topographic factors

The open pit coal mine topography is mostly rugged and uneven. In general, it is a few kilometers wide and more than 50 m deep. In the case of deep open pits, the depth may be as high as 300-400 meters. The complexity arises due to the presence of overburden dumps, construction of access roads and blasting activities. The topography affects the dispersion of pollutants due to changes in the mean flow (Castro & Apsley, 1997). It also affects the trajectory of the wind that channels and confine the plume dispersion and also causes flow re-circulations within the pit (Shi et al, 2000; Appleton et al, 2006).

#### 3.1.1.1. Slopes & benches

Slopes facilitate an easy penetration of the wind because the shape of the pit guides the wind flow (Grainger & Merony, 1993; Peng & Lu, 1995). Besides, the slopes play a major role in the reduction of the size of the primary eddies (Mcquain et al, 1994). The presense of slopes also result in increase in the wind speed (Figure 2). Additionally, the faces and elevation of the benches cause the local deflection of the plume affecting the direction of the dispersion (Figure 3). Moreover, the presence of corners enable the particulate matter to get confined to a particular region within the deep open pit till they reach a particular height where the effect of topography gets diminished and plume expansion is observed. This is called the *venturi effect* (Appleton et al, 2006). The sloping terrain may result in *katabatic* and *anabatic* flows (i.e drainage of air down or uphill sides in response to changing vertical temperature profiles). Further, the terrain elevation may either restrict or exacerbate aereal dispersion when the lower boundary layer heights confine the particulate plume closer to the terrain surface.

#### 3.1.1.2. The depth of the open pit coal mine

When the depth of a mining pit is increased, an open air space with characteristic properties under extreme physical conditions is created (Baklonov, 1986). Under such conditions, the depth of the open pit affects the penetration angle α, which is defined as the angle of the shear interface line defining the boundary between the penetration of the external ABL velocity profile and the induced leeward in pit recirculation zone (Figure 4). Further, the depth has a significant effect on the air flow patterns within the mine. In case of shallow pits, air flow expands gradually, without any gap in mine area. There is a little compression of flow line on outlet area due-to action of flow inertia. In case of medium depth mines, fluid vortex has shape of flattened ellipse whose height and width differ significantly.

The viscous forces, pressure gradients and turbulence friction act as driving mechanisms (Figure 5). Close to leeward side, there is detached point area of thickened boundary layer. In windward side zone, there is flow connection point zone, near which part of flow turns back in detached vortex zone and part moves up and toward exit from zone (Markov et al, 1978). Moreover, for the same incoming airflow and the same slope angle, the mechanical forcing and turbulence increase with increase in the depth, which will result in or strengthen the reverse airflow and the recirculation (Shi et al, 2000).

#### 3.1.1.3. Width and breadth of deep open pit coal mine

The width and breadth of the coal mine affects the rate of cooling in the mine, which is defined in terms of terrain *amplification factor (TAF)* (Equation 1).

where W_{1,} width of the mine at the top,

B_{1,} breadth of the mine at the top,

W_{2,} width of the mine at the bottom and

B_{2,} breadth of the mine at the bottom.

#### 3.1.1.4. Aspect ratio

The *aspect ratio* decides whether the air outside the pit intrudes into the pit or not. Recirculation zones are formed based on the aspect ratio of the pit (Figure 6). The formation of the secondary vortex has been observed as the aspect ratio increases. It is due to the fact that the penetration of the external boundary layer in to the cavity acts as a driving force for the formation of the secondary vortex.

#### 3.1.1.5. Effect of overburden dumps

Baklonov and Regina (1998) has investigated the effect of overburden dumps on dispersion and observed that the velocity increases as the wind passes through the dump from bottom to the top. Further, if the dump is located at the upwind side, it alters the wind flow and turbulence characteristics from those measured at the nearest meteorological station. The rough terrain changes the wind speed directions and turbulence characteristics and enhance mixing in the air stream (Kirchgessner et al, 1993). Huertas (2009) has observed the effect of deposition by altering the area of dumps which shows that the retention of particulates is dependent on the volume and the area of the overburden.

#### 3.1.2. Meteorological effects

#### 3.1.2.1. Wind speed & direction

It is observed that the wind speed is positively correlated and reduced to three fourth after entering into the pit (USEPA, 1985). The higher wind speeds attribute to the lesser deposition fluxes and lesser concentrations of particulates due to the elongation of the particle plumes generated inside the pit. Further, higher wind speeds cause an elongation of particulate plume and a reduction in the airborne concentration and deposition flux gradients observed near source (Appleton et al, 2006). The wind direction inside and outside the pit are not positively co- related. The wind direction displaces the air to a particular direction, thus causing thermal inhomogeneities within the open pit mine. The standard deviations of vertical and horizontal wind directions enable the prediction of the deviation of the wind from its original direction. It has been found that the wind direction changes by an angle of 60º with respect to the direction at the surface (USEPA, 1985).

#### 3.1.2.2. Stability

The stability affects the flow in the mine by influencing the vertical motion. It is observed that stable atmosphere suppresses the vertical motion of the pollutants resulting into deposition of pollutants. However, for unstable and neutral conditions, the escape fraction of the particulates has been found to be more. Under the effect of stratification, the flows in a pit are expected to be similar to the night time flows occurring in the mountains due to the accumulation of cold air and nocturnal cooling (Grainger & Merony, 1993). Further, the inversion depth and strength tend to be stronger than those found over valleys or flat terrain (Maki et al, 1986). Grainger & Merony (1993) have observed that inversion effects are more predominant in pits rather than on flat terrains. The occurrence of the inversion in the pit depends upon the Froude number (Fr), which is independent of the pit geometry. Further, it has been found that the coal pit dispersion under stably-stratified conditions is dominated by the buoyancy-inertia forces (Equation 2).

where,

If Fr > 1.6, then inversion is swept away;

If 1.3 < Fr < 1.6, the air in the basin is coupled to the air above the basin;

and if Fr < 1.3, the air is totally decoupled

Further, the dispersion is dependent on the pasquill stability classes (Table 3).

Stability class | Effects on dispersion in the coal mines |

A | (2), (3), (4) |

B | (2),(3), (4) |

C | (2), (4) |

D | (2) |

E | (1), (2), (3), (5) |

F | (1), (3), (5) |

G | (1), (3), (5) |

Where,

Aids dispersion if terrain is flat and restricts dispersion if the terrain is rising

Impact by particulate deposition

Impact by airborne particulates

Higher values of Surface heat flux aids dispersion due to the generation of upward air currents

Lower values of surface heat flux restricts aerial dispersion

#### 3.1.3. Thermal effects

The combined action of the mechanical shear of the ABL across the surface opening and the thermal buoyancy forces produce the air flow by the differential heating in the pit surface (Silvester et al, 2009). Solar energy controls the physical processes of the atmosphere through the mediation of the surface. Open pit mining modifies the ground surface material composition, structure, cover, morphology, colour which affects the meteorological processes and also affects the partitioning of the incoming solar radiation into sensible and latent heat fluxes which affects the amount of solar radiation that is scattered back to the space. The airspace of a mine pit tends to warm up rapidly when insolation is intense, while cold air tends to accumulate in it when insolation is low or missing (Loska, 2007). More over the solar radiation has different effects during different seasons, which affects the total amount of air heated in the open pit mine (Figure 7). Further, the solar radiation influences many physical processes. The physical processes and their interactions form chains of interactions, which are highly dependent on the type of the terrains and surrounding conditions, like radiation (Figure 8).

#### 3.1.3.1. Thermal inhomogenities

Due to the combined effect of nonuniform horizontal temperature field and the wind direction, thermal inhomogenities are created in the deep open pit. With an active influx of solar radiation, the effect of thermal forces is manifested not only as the component of the resulting wind velocity, but also affects the stability within the deep open pit. The inhomogenities can be classified into four types

With increasing depth, the air temperature rises by a value greater than the adiabatic temperature gradient

The air temperature increases with depth

The air temperature remains constant with depth

The air temperature decreases with depth

Thus two main types of stratification exist within the deep open pit space: unstable (lst type) and stable (2nd, 3rd, and 4th types), which may be present in single form or as combinations (Bitkolov, 1969).

## 4. Modelling the flow in open pit coal mines

Dispersion uses mathematical equations describing the atmosphere along with chemical and physical processes to estimate and predict concentrations at various locations. Dispersion models for predicting particulate concentrations have been discussed by Mikkelson (2003), Reed (2005) and Holmes & Morwska (2006). The conventional models of dispersion for the plane terrain retain flow stationarity and horizontal uniformity hypothesis. As a result, they are no longer valid when the open pit terrain complexity produces mesoscale and local scale circulations. Most of the conventional models employ the hydrostatic approximation which neglects the vertical acceleration versus the pressure gradient and the gravitational terms. This implies that the vertical scale of motion is smaller than the horizontal scale. This applicability of the hydrostatic approximation depends not only on the geometry of the wind, but also on the vertical stability. Hence the non-hydrostatic multiscale approach is preferred. Penenko & Aloyan (1976) have proposed a three-dimensional non stationary model for mesoscale boundary layer of atmosphere. The model includes the effects of horizontal gradients of potential temperatures, stratification parameters and vertical and horizontal turbulence factors. Later, a two dimensional model has been proposed by improving this model to allow considering substantial two-dimensional effects of flow in the deep open pit coal mines having complex geometry (Markov et al, 1978). Further, Aloyan et al (1982) have used a fictious region method in which the equations are generated in a rectangular cartesian coordinate system. The model has helped in simulation of the winds in a quarry under various thermal conditions (stable, unstable and equilibrium stratifications) and investigating the effects of the external wind on the temperature inversion. The model describes the flow in deep open pit coal mines and helps analyzing influence of their geometric parameters by detailed consideration of characteristics of impurity sources. Later, these mathematical models describing the dynamics of the atmosphere in the deep open pits have been combined with the system for calculating the thermal and radiational balance on the surface (Baklonov, 1984). Baklonov (1986) has also established a thermodynamic model for the deep open pit ventilation systems. The model simulates the diurnal changes that take place during the differential solar heating of the surface during the day and the release of this energy during the evening. Further, Baklanov (1995) has presented the results of a series of numerical modeling case studies that consider the pollutant dispersion within naturally ventilated deep open pits, for a range of different external wind speeds and directions, internal thermal stratifications, recirculation flow regimes and the interaction of local open pit thermal circulations with the external atmospheric flow field. In a later study, a high-resolution, non-hydrostatic and three dimensional PBL model for characterizing the flow in the open mine has been developed by Shi et al (2000). The classical 2-D Navier-Stokes equation and a 0-D bond-graph model have been used to investigate the flow dynamics of the mines by discussing the interconnections, nonlinearities and turbulence. The 0-D model efficiently provides for a fast and representative global model for the mine ventilation problem (Witrant et al, 2008). Table 4 summarizes the various studies conducted on open pit mines.

### 4.1. The model and its features

The governing equations which serve the basis of modeling the deep open pit mines have been described (Markov et al, 1978). Equations 3-6 include non stationary hydrodynamic equations of Navier-Stokes, continuity equation for incompressible fluid and equation for carryover of impurity concentration in plane cartesian system of coordinates.

where,

u, w - projection of velocity vectors on horizontal and vertical directions (m/sec);

p - deviation of pressure from hydrostatic (Kpa);

c - specific volumetric concentration of impurities (g/m^{3});

k - turbulent kinematic viscosity factor (kg/m/s);

L_{B -} characteristic longitudinal dimension of upper section of the mine(m);

Re_{T} =

Pr_{d}= k/k_{e,} Prandtl’s diffusion number, which characterizes ratio between turbulence viscosity factor k and diffusion factor k_{e}.

For the values LB/

This model overcomes the limitations of the single dimensional models, which does not allow considering substantial two-dimensional effects of flow in open pit coal mines having complex geometry of surface. This approach helps in solving wide range of problems concerning natural air exchange. The model describes dispersion in deep open coal mines and helps analyzing influence of their geometric parameters on natural air exchange, by considering pollutant characteristics. Thus, the model with sufficient practical accuracy describes the natural ventilation process and can be useful for calculation of flows in deep open pit coal mines.

### 4.2. The model limitations

The present model has no provision for considering the effects of settling of the particles. Moreover, it does not consider: stability (temperature effect)/ buoyancy effects that govern the dispersion of PM within the mine boundary. The model can be modified to incorporate the effects of the buoyancy interms of the temperature difference (Equation 5). Due to the temperature difference existing between the two vertical layers of the atmosphere within the mine, the buoyant force may rise the particle vertically (Grainger & Merony, 1993). Moreover, the temperature equation (Equation 6) can be modified to incorporate the effect of the vertical temperature gradient and the existing temperature gradient. These limitations are presently being addressed as Ph.D problem in Civil Engineering Department at IIT Delhi (Chinthala, 2010).

AQ model/ Technique | Type | Reference |

3-D Non Hydrostatic Model | Numerical | Gresho et al (1976) |

Non Stationary N-S model | Numerical | Markov et al (1978) |

3-D Non Stationary Model | Numerical | Aloyan et al (1982) |

Cole and Fabrick | Box model | Cole & Fabrick (1984) |

2-D Hydrostatic Model | Numerical | Herwehe (1984) |

ISC3 | Gaussian | USEPA(1995c) |

Wind Tunnel | Physical | Peng & Lu (1995) |

K-ε hydrodynamic model | Numerical | Fomin (1996) |

Wind Tunnel | Physical | Shi et al (2000) |

Wind Tunnel | Physical | Grainger & Merony (2003) |

Fabrick | Gaussian | Reed (2005) |

Gelekin | Gaussian | Reed (2005) |

Herwere | Eulerian | Reed (2005) |

Kalgoorlie | Unknown | Reed (2005) |

PerieraSoares&Branquinho | Gaussian | Reed (2005) |

Winges | Gaussian | Reed (2005) |

2-D Non Hydrostatic model | Numerical | Kharytonov et al (2005) |

ADMS | Gaussian | Appleton et al (2006) |

RANS | Numerical | Bodnor et al (2008) |

CFD | Numerical | Silvester et al (2009) |

k-ε model | Numerical | Kakosimos et al (2011) |

## 5. Conclusions

Simulating the dispersion in deep open pit coal mines is an effective tool to give advanced warning of potential emission problems and providing the basis for future planning applications. The meteorological conditions within an open pit coal mine are significantly affected by temperature (stability) and roughness conditions that generate complex fluid flows phenomenon e.g. separation of atmospheric boundary layer, recirculation, resuspension and settling of particulates_{.} The complex geometry, uncontrolled wind flows and different operations taken place within the deep open pit coal mines, makes it extremely difficult to measure the particulates. Therefore, an appropriate mathematical model is needed which may predict the complex processes in the near field regions of the deep open pit coal mines effectively. The fundamental understanding of the microclimate that is created in the deep open pit mine due to influence of the external ABL and the contributions of the mechanical and thermal forces present, can be achieved by developing a dynamic numerical based inverse hill model to investigate the complex wind flows and dispersion mechanisms.