1 / 36

# TOXIC RELEASE & DISPERSION MODELS - PowerPoint PPT Presentation

TOXIC RELEASE & DISPERSION MODELS. Prepared by Associate Prof. Dr. Mohamad Wijayanuddin Ali Chemical Engineering Department Universiti Teknologi Malaysia. The applicable conditions are - - Constant mass release rate, Q m = constant, - No wind, < u j > = 0,

Related searches for TOXIC RELEASE & DISPERSION MODELS

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

## PowerPoint Slideshow about 'TOXIC RELEASE & DISPERSION MODELS' - flavio

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

DISPERSION MODELS

Prepared by

Associate Prof. Dr. Mohamad Wijayanuddin Ali

Chemical Engineering Department

Universiti Teknologi Malaysia

- Constant mass release rate, Qm = constant,

- No wind, <uj> = 0,

- Steady state, <C>/t = 0, and

- Constant eddy diffusivity, Kj = K* in all directions.

For this case, Equation 9 reduces to the form,

(10)

Case 1 : Steady-state, Continuous Point

Release with No Wind

Equation 10 is more tractable by defining a radius as r² = x² + y² + z². Transforming Equation 10 in terms of r yields

(11)

For a continuous, steady state release, the concentration flux at any point, r, from the origin must equal the release rate, Qm (with units of mass/time). This is represented mathematically by the following flux boundary condition.

(12)

The remaining boundary condition is

(13)

(14)

Solving Equation 14 for <C> yields,

(15)

It is easy to verify by substitution that Equation 15 is also a solution to Equation 11 and thus a solution to this case. Equation 15 is transformed to rectangular coordinates to yield,

(16)

- Puff release, instantaneous release of a fixed mass of material, Qm* (with units of mass),

- No wind, <uj> = 0, and

- Constant eddy diffusivity, Kj = K*, in all directions.

Equation 9 reduces, for this case, to

(17)

Case 2 : Puff with No Wind

(18)

The solution to Equation 17 in spherical coordinates is

(19)

and in rectangular coordinates is

(20)

- Constant mass release rate, Qm = constant,

- No wind, <uj> = 0, and

- Constant eddy diffusivity, Kj = K* in all directions.

For this case, Equation 9 reduces to Equation 17 with initial condition, Equation 18, and boundary condition, Equation 13. The solution is found by integrating the instantaneous solution, Equation 19 or 20 with respect to time. The result in spherical coordinates is

(21)

Case 3 : Non Steady-state, Continuous

Point Release with No Wind

(22)

As t, Equations 21 and 22 reduce to the corresponding steady state solutions, Equations 15 and 16.

- Continuous release, Qm = constant,

- Wind blowing in x direction only, <uj> = <ux> = u = constant, and

- Constant eddy diffusivity, Kj= K* in all directions.

For this case, Equation 9 reduces to

(23)

Point Source Release with Wind

Equation 23 is solved together with boundary conditions, Equation12 and 13. The solution for the average concentration at any point is

(24)

If a slender plume is assumed (the plume is long and slender and is not far removed from the x-axis),

(25)

and, using , Equation 24 is simplified to

(26)

Along the centreline of this plume, y = z = 0 and

(27)

This is the same as Case 2, but with eddy diffusivity a function of direction. The applicable conditions are -

- Puff release, Qm* = constant,

- No wind, <uj> = 0, and

- Each coordinate direction has a different, but constant eddy diffusivity, Kx, Ky and Kz.

Equation 9 reduces to the following equation for this case.

(28)

The solution is

(29)

Case 5 : Puff with No Wind.

Eddy Diffusivity a Function of Direction

Case 6 : Steady-state, Continuous Point Source function of direction. The applicable conditions are -

Release with Wind.

Eddy Diffusivity a Function of Direction

This is the same as Case 4, but with eddy diffusivity a function of direction. The applicable conditions are -

- Puff release, Qm* = constant,

- Steady state, <C>/t = o,

- Wind blowing in x direction only, <uj> = <ux> = u = constant,

- Each coordinate direction has a different, but constant eddy diffusivity, Kx, Ky and Kz, and

- Slender plume approximation, Equation 25.

Equation 9 reduces to the following equation for this case. function of direction. The applicable conditions are -

(30)

The solution is

(31)

Along the centreline of this plume, y = z = 0 and the average concentration is given by

(32)

Case 7 : Puff with Wind function of direction. The applicable conditions are -

This is the same as Case 5, but with wind. Figure 8 shows the geometry. The applicable conditions are -

- Puff release, Qm* = constant,

- Wind blowing in x direction only, <uj> = <ux> = u = constant, and

- Each coordinate direction has a different, but constant eddy diffusivity, Kx, Ky and Kz,.

The solution to this problem is found by a simple transformation of coordinates. The solution to Case 5 represents a puff fixed around the release point.

If the puff moves with the wind along the function of direction. The applicable conditions are -x-axis, the solution to this case is found by replacing the existing coordinate x by a new coordinate system, x - ut, that moves with the wind velocity. The variable t is the time since the release of the puff, and u is the wind velocity. The solution is simply Equation 29, transformed into this new coordinate system.

(33)

This is the same as Case 5, but with the source on the ground. The ground represents an impervious boundary. As a result, the concentration is twice the concentration as for Case 5. The solution is 2 times Equation 29.

(34)

Case 8 : Puff with No Wind

with Source on Ground

This is the same as Case 6, but with the release source on the ground, as shown in Figure 9. The ground represents an impervious boundary. As a result, the concentration is twice the concentration as for Case 6. The solution is 2 times Equation 31.

(35)

with Source on Ground

Figure 9 the ground, as shown in Figure 9. The ground represents an impervious boundary. As a result, the concentration is twice the concentration as for Case 6. The solution is 2 times Equation 31.Steady-state plume with source at ground level. The concentration is twice the concentration of a plume without the ground.

For this case the ground acts as an impervious boundary at a distance H from the source. The solution is

(36)

Case 10 : Continuous, Steady-state Source.

Source at Height Ht, above the Ground

Cases 1 through 10 above all depend on the specification of a value for the eddy diffusivity, Kj. In general, Kj changes with position, time, wind velocity, and prevailing weather conditions. While the eddy diffusivity approach is useful theoretically, it is not convenient experimentally and does not provide a useful framework for correlation.

Sutton solved this difficulty by proposing the following definition for a dispersion coefficient.

(37)

with similar relations given for sy and sz. The dispersion coefficients, sx, sy, and sz represent the standard deviations of the concentration in the downwind, crosswind.

Pasquill-Gifford Model

with similar relations given for a value for the eddy diffusivity, sy and sz. The dispersion coefficients, sx, sy, and sz represent the standard deviations of the concentration in the downwind, crosswind, and vertical (x,y,z) directions, respectively. Values for the dispersion coefficients are much easier to obtain experimentally than eddy diffusivities.

The dispersion coefficients are a function of atmospheric conditions and the distance downwind from the release. The atmospheric conditions are classified according to 6 different stability classes shown in Table 2. The stability classes depend on wind speed and quantity of sunlight. During the day, increased wind speed results in greater atmospheric stability, while at night the reverse is true. This is due to a change in vertical temperature profiles from day to night.

The dispersion coefficients, a value for the eddy diffusivity, sy and sz for a continuous source were developed by Gifford and given in Figures 10 and 11, with the corresponding correlation given in Table 3. Values for sx are not provided since it is reasonable to assume sx = sy. The dispersion coefficients sy and sz for a puff release are given in Figures 12 and 13. The puff dispersion coefficients are based on limited data (shown in Table 3) and should not be considered precise.

The equations for Cases 1 through 10 were rederived by Pasquill using relations of the form of Equation 37. These equations, along with the correlation for the dispersion coefficients are known as the Pasquill-Gifford model.

Stability class for puff model : a value for the eddy diffusivity,

A,B : unstable

C,D : neutral

E,F : stable

Table 2 Atmospheric Stability Classes for Use with the Pasquill-Gifford Dispersion Model

Figure 10 a value for the eddy diffusivity, Horizontal dispersion coefficient for Pasquill-Gifford plume model. The dispersion coefficient is a function of distance downwind and the atmospheric stability class.

Figure 11 a value for the eddy diffusivity, Vertical dispersion coefficient for Pasquill-Gifford plume model. The dispersion coefficient is a function of distance downwind and the atmospheric stability class.

Figure 12 a value for the eddy diffusivity, Horizontal dispersion coefficient for puff model. This data is based only on the data points shown and should not be considered reliable at other distances.

Figure 13 a value for the eddy diffusivity, Vertical dispersion coefficient for puff model. This data is based only on the data points shown and should not be considered reliable at other distances.

Table 3 a value for the eddy diffusivity, Equations and data for Pasquill-Gifford Dispersion Coefficients

This case is identical to Case 7. The solution has a form similar to Equation 33.

(38)

The ground level concentration is given at z = 0.

(39)

Case 1 : Puff. Instantaneous Point Source at Ground Level,

Coordinates Fixed at Release Point.

Constant Wind in x Direction Only with Constant Velocity u

The ground level concentration along the similar to Equation 33.x-axis is given at y = z= 0.

(40)

The centre of the cloud is found at coordinates (ut,0,0). The concentration at the centre of this moving cloud is given by

(41)

The total integrated dose, Dtidreceived by an individual standing at fixed coordinates (x,y,z) is the time integral of the concentration.

(42)

The total integrated dose at ground level is found by integrating Equation 39 according to Equation 42. The result is -

(43)

The total integrated dose along the x-axis on the ground is

(44)

Frequently the cloud boundary defined by a fixed concentration is required. The line connecting points of equal concentration around the cloud boundary is called an isopleth.

For a specified concentration, < integrating Equation 39 according to Equation 42. The result is -C>*, the isopleths at ground level are determined by dividing the equation for the centreline concentration, Equation 40, by the equation for the general ground level concentration, Equation 39. This equation is solved directly for y.

(45)

The procedure is

1. Specify <C>*, u, and t.

2. Determine the concentrations, <C> (x,0,0,t), along the x-axis using Equation40. Define the boundary of the cloud along the x-axis.

3. Set <C> (x,y,0,t) = <C>* in Equation 45 and determine the values of y at each centreline point determined in step 2.

The procedure is repeated for each value of t required.

This case is identical to Case 9. The solution has a form similar to Equation 35.

(46)

The ground level concentration is given at z = 0.

(47)

Case 2 : Plume. Continuous, Steady-state,

Source at Ground Level, Wind Moving

in x Direction at Constant Velocity u

The concentration along the centreline of the plume directly downwind is given at y = z= 0.

(48)

The isopleths are found using a procedure identical to the isopleth procedure used for Case 1.

For continuous ground level releases the maximum concentration occurs at the release point.