Open access peer-reviewed chapter

Particle Dispersion Within a Deep Open Cast Coal Mine

By Sumanth Chinthala and Mukesh Khare

Submitted: October 8th 2010Reviewed: March 27th 2011Published: July 5th 2011

DOI: 10.5772/16326

Downloaded: 4192

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 PM10 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 PM10. The dispersion equations developed within the deep pit boundary provide a reasonable accurate estimate of PM10 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 PM10 (Bitkolov, 1969; Grainger and Merony, 1993). However, the in-field measurements of PM10 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 PM10 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 PM10 (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 diffusionand advectionprocesses. The atmospheric motions transportand diffusethe 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 PM10 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 PM10 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
Overburden removal7%
Top soil removal1%
Coal processing72%
Coal extraction3%
Wind erosion17%

Table 1.

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


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 noProperty /parameterShallowIntermediateDeep
1Recompression layerNot InvestigatedNot InvestigatedDecreases with increased L/D ratio
2Boundary layer and its influenceNot InvestigatedNot InvestigatedBridges the cavity and has little influence over the cavity flow parameters
3Redevelopment lengthRedevelopment is smallRedevelopment is maximumRedevelopment is small
4Variation of coefficient of pressureLowMaximumLow
5Boundary layer thickeningMore
6Free shear layer flowLess strongly distributed
More strongly distributedLess strongly distributed
7Drag Coefficient
8Flow typeunsteadyNot InvestigatedDead/laminar
9Flow reattachmentOccurs on the floor of the cavityNot InvestigatedOccurs on the rear shoulder of the cavity
10Reynolds numberLargeMediumSmall
12Skin frictionLargeMediumSmall
13Base pressureFirst decreases and then increasesNot InvestigatedConstant in the middle and slightly irregular in the front and rear portions

14Three dimensionality of shear layer flowSlightlyStronglySlightly

Table 2.

Types of cavities and their characteristics (Chowdhary, 1977)

3.1. The cavity regions

The cavity region contains different zones where the boundary layer undergoes various changes (Figure 1). The Upstream zonecomprises the developing boundary layer zone which is upstream of the front corner of the notch with a finite approaching boundary layer. Expansion zoneis 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 zoneis 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 compression 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 zonewhich 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.

Figure 1.

Cavity flow field operating in open configuration (solid line) and closed configuration (dotted line).Source:Lazar et al, 2008

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). 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 katabaticand anabaticflows (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.

Figure 2.

Wind speed due to the presence of slopes (Ngo and Lechford, 2008) 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.

Figure 3.

Effect of local deflection on the dispersion of particulate plume (Appleton et al, 2006)

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).

Figure 4.

Penetration angle α in deep open pit mines (Peng and Lu, 1995) 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).

TAF= 3W1B1/ [W1B1+ W2B2+ (W1B1W2B2)1/2]E1

where W1, width of the mine at the top,

B1, breadth of the mine at the top,

W2, width of the mine at the bottom and

B2, breadth of the mine at the bottom. Aspect ratio

The aspect ratiodecides 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.

Figure 5.

Vector fields of speed of steady-state flow processes in the deep open pit mines (Markov et al, 1978). 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.

Figure 6.

Effect of aspect ratio on the flow (Kang & Sung, 2009)

3.1.2. Meteorological effects 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). 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).

Fr=Uh [gh(θh θ0)θ]1/2E2


θh, is the air temperature at the ridge topθ0, is the air temperature at the valley floor

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 classEffects on dispersion in the coal mines
A(2), (3), (4)
B(2),(3), (4)
C(2), (4)
E(1), (2), (3), (5)
F(1), (3), (5)
G(1), (3), (5)

Table 3.

Factors affecting dispersion for different stability classes


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

  2. Impact by particulate deposition

  3. Impact by airborne particulates

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

  5. 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). 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

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

  2. The air temperature increases with depth

  3. The air temperature remains constant with depth

  4. 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).

Figure 7.

Effect of solar radiation during winter and summer in the mine pit

Figure 8.

Chains of interactions of physical processes (Zilitinkevich, 1990)


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.



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/m3);

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

LB - characteristic longitudinal dimension of upper section of the mine(m);

ReT =LBu/kReynold’s turbulence number;

Prd= k/ke, Prandtl’s diffusion number, which characterizes ratio between turbulence viscosity factor k and diffusion factor ke.

For the values LB/u,u,ρu2andk, the subscript characterizes the parameters of incident flow.

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/ TechniqueTypeReference
3-D Non Hydrostatic ModelNumericalGresho et al (1976)
Non Stationary N-S model
Markov et al (1978)
3-D Non Stationary Model

Aloyan et al (1982)
Cole and Fabrick

Box model
Cole & Fabrick (1984)
2-D Hydrostatic ModelNumericalHerwehe (1984)
Wind TunnelPhysicalPeng & Lu (1995)
K-ε hydrodynamic model
Fomin (1996)
Wind TunnelPhysicalShi et al (2000)
Wind TunnelPhysicalGrainger & Merony (2003)
FabrickGaussianReed (2005)
GelekinGaussianReed (2005)
HerwereEulerianReed (2005)
KalgoorlieUnknownReed (2005)
PerieraSoares&BranquinhoGaussianReed (2005)
Reed (2005)
2-D Non Hydrostatic modelNumericalKharytonov et al (2005)
ADMSGaussianAppleton et al (2006)
RANSNumericalBodnor et al (2008)
Silvester et al (2009)
k-ε modelNumericalKakosimos et al (2011)

Table 4.

Modelling studies on open pit mines


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.

© 2011 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike-3.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited and derivative works building on this content are distributed under the same license.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sumanth Chinthala and Mukesh Khare (July 5th 2011). Particle Dispersion Within a Deep Open Cast Coal Mine, Air Quality - Models and Applications, Dragana Popovi?, IntechOpen, DOI: 10.5772/16326. Available from:

chapter statistics

4192total chapter downloads

2Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Mathematical Modeling of Air Pollutants: An Application to Indian Urban City

By P. Goyal and Anikender Kumar

Related Book

First chapter

Observational Study of Black Carbon in the North Suburb of Nanjing, China

By Lili Tang, Shengjie Niu, Mingliang Yan, Xuwen Li, Xiangzhi Zhang, Yuan Zhu, Honglei Shen, Minjun Xu and Lei Tang

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us