1 / 36

TOXIC RELEASE & DISPERSION MODELS

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,

flavio
Télécharger la présentation

TOXIC RELEASE & DISPERSION MODELS

An Image/Link below is provided (as is) to download presentation 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. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. TOXIC RELEASE & DISPERSION MODELS Prepared by Associate Prof. Dr. Mohamad Wijayanuddin Ali Chemical Engineering Department Universiti Teknologi Malaysia

  2. The applicable conditions are - - 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

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

  4. Equation 12 is separated and integrated between any point r and r =. (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)

  5. The applicable conditions are - - 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

  6. The initial condition required to solve Equation 17 is (18) The solution to Equation 17 in spherical coordinates is (19) and in rectangular coordinates is (20)

  7. The applicable conditions are - 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

  8. and in rectangular coordinates is (22) As t, Equations 21 and 22 reduce to the corresponding steady state solutions, Equations 15 and 16.

  9. This case is shown in Figure 7. The applicable conditions are - 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) Case 4 : Steady-state, Continuous Point Source Release with Wind

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

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

  12. Case 6 : Steady-state, Continuous Point Source 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.

  13. Equation 9 reduces to the following equation for this case. (30) The solution is (31) Along the centreline of this plume, y = z = 0 and the average concentration is given by (32)

  14. Case 7 : Puff with Wind 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.

  15. If the puff moves with the wind along the 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)

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

  17. 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) Case 9 : Steady-state Plume with Source on Ground

  18. Figure 9 Steady-state plume with source at ground level. The concentration is twice the concentration of a plume without the ground.

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

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

  21. 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, 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.

  22. The dispersion coefficients, 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.

  23. Stability class for puff model : A,B : unstable C,D : neutral E,F : stable Table 2 Atmospheric Stability Classes for Use with the Pasquill-Gifford Dispersion Model

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

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

  26. Figure 12 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.

  27. Figure 13 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.

  28. Table 3 Equations and data for Pasquill-Gifford Dispersion Coefficients

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

  30. The ground level concentration along the 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)

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

  32. For a specified concentration, <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.

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

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

More Related