Volume 1 - The Atmospheric Studies Group at TRC

Volume 1 - The Atmospheric Studies Group at TRC Volume 1 - The Atmospheric Studies Group at TRC

29.01.2015 Views

3.2 CALPUFF Platform Downwash The existing CALPUFF Huber-Snyder/Schulman-Scire (HS/SS) downwash modules are used to model downwash effects due to elevated structures by using that portion of the structure that is “solid”. A platform height (z plat ) is added to the list of variables that describe the effective building dimensions The effective building width H w and building height H b are prepared for each of 36 wind directions (10- degree intervals) by neglecting the gap below z plat . That is, the structure is defined as if the solid portion rests on the “ground”, and the EPA Building Profile Input Program (BPIP) can be used to develop the direction-specific effective height and width. Any point-source emission released on or near the structure is prescribed using the full release height h s above the “ground”, not the height above the platform deck. The full release elevation above ground is adjusted by subtracting the platform height prior to any tests that define the downwash potential (e.g., the 2.5 building-height rule for GEP), and any downwash plume enhancements that depends on the effective stack height. This adjusted stack height is not used as the physical release height in any other calculations. The Huber-Snyder (Huber and Snyder, 1976; Huber, 1977) technique is used for h − z > H + T L (3-36) s plat b bd b where L b is the lesser of the effective building height and width, and T bd has a default value of 0.5. A negative value of T bd indicates the Huber-Snyder method is used for all stacks, and a value of 1.5 results in the Schulman-Scire (Scire and Schulman, 1980; Schulman and Hanna, 1986) method always being used. If T Tb is set equal to 0.5 (its default value), the CALPUFF treatment will be equivalent to that in ISC3. When the Huber-Snyder technique is used, the first step is to compute the effective plume height H e due to momentum rise at a downwind distance of 2 H b . This rise uses the wind speed at the full stack height, h s . If (H e –z plat ) exceeds H b + 1.5 L b , building downwash effects are assumed to be negligible. Otherwise, buildinginduced enhancement of the plume dispersion coefficients is evaluated. For adjusted stack heights h s –z plat less than 1.2H b , both σ y and σ z are enhanced. Only σ z is enhanced for adjusted stack heights above 1.2 H b (but below H b + 1.5 L b ). Enhancements to σ y and σ z are not functions of h s or h s -z plat . When the Schulman-Scire technique is used, a linear decay factor is applied to the building-induced enhancement of the vertical dispersion coefficient, and plume rise is Final Report Vol.1 16

adjusted to account for the effect of downwash. The plume rise equations are not functions of h s or h s -z plat . The linear decay factor for σ z is determined as: σ ′ z = A σ ′ z (3-37) where σ ′ z is determined from the HS downwash equations, and ⎧1 ⎪ H A = ⎨ ⎪ ⎪ ⎩0 b − ( H e 2L b − z plat ) + 1 H b < ( H H b e H − z e + 2L plat b ≤ ( H ) ≤ H < ( H b e − z b plat ) + 2L − z plat b ) (3-38) AERMOD Turbulence Profile Option The turbulence velocity for horizontal fluctuations is computed from contributions from shear and buoyancy. The total turbulence velocity is obtained by summing the component variances: σ = σ + σ (3-39) v 2 vs 2 vb The shear component is modeled with a variance that is a maximum at the ground (σ vs 2 =3.6u * 2 ), and decreases linearly through the depth of the mechanically mixed layer to a residual value of 0.25 m 2 /s 2 , if this residual value is less than the value at the ground. If the residual value is larger, the value at the ground is used for all heights. The buoyancy component is constant (σ vb 2 =0.35w * 2 ) up to z i , and decreases linearly to the residual value at 1.2z i . Again, the value at the ground is used at all heights if it is less than the residual value. The turbulence velocity for vertical fluctuations is also computed from shear and buoyancy contributions. The buoyancy component of the variance is computed as: 2 wb 2 * 3 σ = 1.6w ( z z ) 2 / ( z ≤ 0.1z ) (3-40) 2 wb i i 2 σ = 0.35w ( z = 0.1z to z ) (3-41) 2 wb * i i 2 −6( z−z ) zi σ = 0.35w e i ( z > z ) (3-42) * i The shear component has both a surface-driven contribution (from u * ) and a residual contribution from turbulence aloft that is assumed to have an intensity of order 2% of the wind speed at z i : 2 ws 2 * 2 σ = 1.3u (1 − z z ) + (0.02 u z / z ) ( z ≤ z ) (3-43) i zi i i Final Report Vol.1 17

adjusted to account for the effect of downwash. <strong>The</strong> plume rise equ<strong>at</strong>ions are not<br />

functions of h s or h s -z pl<strong>at</strong> . <strong>The</strong> linear decay factor for σ z is determined as:<br />

σ ′ z<br />

= A σ ′ z<br />

(3-37)<br />

where σ ′<br />

z<br />

is determined from the HS downwash equ<strong>at</strong>ions, and<br />

⎧1<br />

⎪<br />

H<br />

A = ⎨<br />

⎪<br />

⎪<br />

⎩0<br />

b<br />

− ( H<br />

e<br />

2L<br />

b<br />

− z<br />

pl<strong>at</strong><br />

)<br />

+ 1<br />

H<br />

b<br />

< ( H<br />

H<br />

b<br />

e<br />

H<br />

− z<br />

e<br />

+ 2L<br />

pl<strong>at</strong><br />

b<br />

≤ ( H<br />

) ≤ H<br />

< ( H<br />

b<br />

e<br />

− z<br />

b<br />

pl<strong>at</strong><br />

)<br />

+ 2L<br />

− z<br />

pl<strong>at</strong><br />

b<br />

)<br />

(3-38)<br />

AERMOD Turbulence Profile Option<br />

<strong>The</strong> turbulence velocity for horizontal fluctu<strong>at</strong>ions is computed from contributions<br />

from shear and buoyancy. <strong>The</strong> total turbulence velocity is obtained by summing the<br />

component variances:<br />

σ = σ + σ<br />

(3-39)<br />

v<br />

2<br />

vs<br />

2<br />

vb<br />

<strong>The</strong> shear component is modeled with a variance th<strong>at</strong> is a maximum <strong>at</strong> the ground<br />

(σ vs 2 =3.6u * 2 ), and decreases linearly through the depth of the mechanically mixed<br />

layer to a residual value of 0.25 m 2 /s 2 , if this residual value is less than the value <strong>at</strong><br />

the ground. If the residual value is larger, the value <strong>at</strong> the ground is used for all<br />

heights. <strong>The</strong> buoyancy component is constant (σ vb 2 =0.35w * 2 ) up to z i , and decreases<br />

linearly to the residual value <strong>at</strong> 1.2z i . Again, the value <strong>at</strong> the ground is used <strong>at</strong> all<br />

heights if it is less than the residual value.<br />

<strong>The</strong> turbulence velocity for vertical fluctu<strong>at</strong>ions is also computed from shear and<br />

buoyancy contributions. <strong>The</strong> buoyancy component of the variance is computed as:<br />

2<br />

wb<br />

2<br />

*<br />

3<br />

σ = 1.6w<br />

( z z )<br />

2 / ( z ≤ 0.1z<br />

)<br />

(3-40)<br />

2<br />

wb<br />

i<br />

i<br />

2<br />

σ = 0.35w ( z = 0.1z<br />

to z )<br />

(3-41)<br />

2<br />

wb<br />

*<br />

i<br />

i<br />

2 −6(<br />

z−z<br />

) zi<br />

σ = 0.35w<br />

e<br />

i ( z > z )<br />

(3-42)<br />

*<br />

i<br />

<strong>The</strong> shear component has both a surface-driven contribution (from u * ) and a residual<br />

contribution from turbulence aloft th<strong>at</strong> is assumed to have an intensity of order 2% of<br />

the wind speed <strong>at</strong> z i :<br />

2<br />

ws<br />

2<br />

*<br />

2<br />

σ = 1.3u<br />

(1 − z z ) + (0.02 u z / z ) ( z ≤ z )<br />

(3-43)<br />

i<br />

zi<br />

i<br />

i<br />

Final Report Vol.1 17

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!