
книги из ГПНТБ / Бызова, Н. Л. Рассеяние примеси в пограничном слое атмосферы
.pdfгде
' о
дисперсия пульсаций скорости, рассчитанная по укороченной эйле ровой реализации, длиной Т.
При g ^ l |
и т]С1, т. е. вблизи |
от источника короткого |
времени |
||
действия, можно получить |
|
|
|
||
4 ( 0 |
2 <Г v2 |
t'A |
/ Т \2''3 |
(2-13) |
|
= |
+ 0 - 4 5 |
[^-) |
<v*>t*. |
Если второй член в выражении (2.13) преобладает, то источник можно считать квазипостояиным, причем а\ изменяется как Г2 '3 . В противоположном случае при
источник работает как квазимгновенный.
На рис. 2.1 дана зависимость cr /oj от £ при разных значе ниях ср(г|).
О\ |
z |
г |
I I I 1 I W1 . |
I |
| I I I I |
10° |
| | |
IJ |
w- |
I |
1.1 |
|
|
i l l |
Рис. 2.1. Зависимость ат /о, от ^=t\xL. Значения ср(т]) нане
сены около кривых
В том случае, когда средняя скорость сноса отсутствует, при аппроксимации начального расстояния между частицами в одно мерном случае можно положить г = У < v2 >> х.Тогда все соотно-
шения останутся прежними |
при замене |
в них |
величины U на |
|
/ |
|
|
t V o 2 > |
|
]/ <v2^> |
, соответственно в этом случае |
f\~ |
—— , |
|
|
|
|
|
' Е |
|
2.1.2. Учет неоднородности приземного |
слоя |
||
при |
расчете диффузии |
в поперечном |
ветру |
направлении |
В пограничном слое атмосферы, и в особенности в приземной его части, турбулентность неоднородна, и характеристики диффу зии существенно зависят от высоты. В том случае, когда источник примеси и ее факел расположены «а одной высоте — около 50 м и выше, можно считать, что условия однородности приближенно соблюдаются, что позволяет непосредственно пользоваться приве денными методами расчета. Однако влияние неоднородности на формирование распределения примеси на уровне земли должно быть весьма существенным. Изложенный далее способ учета не однородности нижней части пограничного слоя атмосферы при расчете поперечной дисперсии факела от высотного источника был применен в работе Бызовой и Иванова (1967) и Бызовой (1970а).
х
Рис. 2.2. Схема движения частицы в неоднородном слое:
I — профиль среднего ветра; 2 — случайная траектория. 3 — средняя траектория
Пусть источник расположен на уровне И, а подстилающая по верхность — на уровне zo. Рассмотрим те частицы, которые после выхода из источника попадают на уровень z0 на некотором рас стоянии х от его эпицентра. Каждая из них двигалась по некото рой случайной траектории (рис. 2.2), причем скорость ее верти кального пульсационного движения ws в среднем по траектории была отлична от нуля и составила
ws = -^Z^. |
, |
(2.15) |
to |
|
|
где t0 — время пребывания ее в атмосфере. |
|
|
Поскольку за время <t0 по оси х частица |
смещается под дейст- |
41
вием средней скорости ветра, то, пренебрегая горизонтальным про дольным пульсационным движением, имеем
x = (N-z0)^, |
(2.16) |
где
н
U„--=-jj^U{z)dz -
средняя в слое скорость ветра, причем для всех траекторий, при водящих в точку (х, z0), величина ws одинакова. Заменим эти случайные траектории средней, которая получается в результате равномерного движения частиц со скоростью ws вниз при одно временном сносе в направлении х со скоростью U (z). Эта траекто рия, после исключения координаты х, описывается соотношением
z = H — wst; |
(2.18) |
Частицы, попавшие на уровень z0 на расстоянии х, очевидно, двигались вместе с теми индивидуальными турбулентными вихря ми, которые имели преимущественно нисходящую вертикальную составляющую пульсаций скорости и пересекали неоднородный слой. Будем считать, что корреляционная функция скоростей при движении вдоль этого пути меняется так, как если бы турбулент ные характеристики зависели только от высоты и не зависели от способа попадания частицы на данную высоту
|
BL |
{Х, z) = v\z) exp |
|
|
|
(2.19) |
||
Применяя формулу |
(.1.10) для всех частиц, |
попавших на уровень |
||||||
z0 на расстоянии х, поперечную дисперсию |
факела на уровне Zo |
|||||||
получим в виде |
|
|
|
|
|
|
|
|
|
o*(t0) = 2 / ( * 0 - т ) BL |
( Х , |
Z) dx |
|
(2.20) |
|||
при условии |
(2.18), |
где to определяется |
соотношениями |
(2.15) и |
||||
(2.16), так что в конечном |
счете получается зависимость |
попереч |
||||||
ной дисперсии на уровне z0 |
от координаты х. В частности, в при |
|||||||
земном слое |
атмосферы |
при безразличной |
стратификации |
<v2> |
||||
не меняется |
с высотой, |
скорость ветра подчиняется логарифмиче |
скому профилю |
(1.62), а лагранжев масштаб |
времени имеет вид |
(Монин, Яглом, |
1965) |
|
|
T £ ( Z ) = — . |
(2.21) |
AIL I
Подставив (2.21) в (2.20), после интегрирования получим
|
|
|
|
(2.22) |
|
|
V<&> |
xL |
(И), |
|
|
|
|
(2.23) |
Масштабы г н |
и х н |
определяются в этом случае высотой |
источника |
|
согласно выражениям (2.23). В пределе при малых |
С ^ ^ - г ^ С 2 , |
|||
а при больших |
а- |
х . |
|
|
р |
|
и |
|
|
Отметим, что нормирующие множители Гн — rL{H), |
хн~ |
xL{H), |
а также предельные соотношения определяются значениями тур булентных характеристик на уровне источника, а осредняется в слое только скорость ветра. Это, на первый взгляд, неожиданное обстоятельство представляется естественным, так как на началь
ном этапе |
(малые g) величина а2 |
определяется интенсивностью |
||
поперечной |
компоненты скорости |
ветра при |
условии, |
что <v2> |
не зависит |
от г, а при g -*- оо к поверхности |
попадают |
частицы с |
очень пологой траекторией, которая расположена на уровне источ
ника вплоть до тех |
значений времени диффузии, когда вступает |
в силу соотношение |
(1.14). Основное влияние неоднородности за |
ключается в том, что для однородного слоя нормированные а не |
много больше и при больших g быстрее приближаются к предель
ной |
зависимости. |
|
|
|
|
|
|
|
2.1.3. О рассеянии |
тяжелых |
частиц |
в |
турбулентной |
среде |
|
М. И. Юдиным (1959) сформулированы три механизма |
влия- |
||||||
ния |
гравитационного |
оседания |
частиц |
на |
диффузию их |
в |
турбу- |
лентной среде. Первый из них — смещение центра тяжести обла ка, частиц — учитывается соответствующим членом в полуэмпири ческом _у_равнении турб"улентной диффузии. Второй связан с тем, что из-за своей инерции частицы не полностью воспринимают вы
сокочастотные |
флюктуации среды; |
этот механизм рассматривался |
в ряде работ |
(В. С. Синельщиков, |
1967; Хинце, 1959); при диффу |
зии даже довольно тяжелых, но мелких частиц в турбулентной ат мосфере вдали от границ его действием можно пренебречь. Тре тий механизм связан с тем, что при своем падении тяжелая части ца пересекает траектории окружающих частиц воздуха, последо вательно попадая в сферу влияния разных вихрей. Расчет влияния этого эффекта был сделан для постоянно действующего источника Юдиным (1945; 1959), для облака отдельных частиц — Смитом (1959) и Смитом и Хэем (1961). Приведем некоторое обобщение результатов Юдина.
Следуя его идеям, используем для расчета выражение (1.10), где корреляционная функция скоростей должна быть взята вдоль средней траектории рассеивающихся частиц. Рассмотрим два пре дельных случая. Если скорость гравитационного оседания частиц w мала по сравнению с пульсационной скоростью в направлении диффузии Y < у 2 > , так что за время корреляции частица не ус певает существенно сместиться по направлению действия гравита ционной силы z, то ее можно считать практически невесомой и в выражении (1.10) пользоваться лагранжевой временной корреля
ционной |
функцией BL{I)- |
Если же w > Y<C'V2^> |
• частицы на |
столько |
быстро уходят по направлению действия |
гравитационной |
силы, что причиной уменьшения корреляции в этом случае будет перемещение их в пространстве по направлению z, и в выражении (1.10) следует взять эйлерову пространственную корреляционную функцию.
В промежуточном случае корреляционную функцию скоростей с некоторым приближением можно взять вдоль среднего пути па дающих частиц
BJr, T) = < V (л-0 , z, t) v (х0, z + г, t + т) > , |
(2.24) |
где г = ш т . Выражение (2.24) представляет собой смешанную про странственно-временную корреляционную функцию скоростей жидкой частицы. Зависимость ее от времени и координат известна только в инерционном интервале по осям координат, где она под чиняется (1.28) и (1.76).
Приняв Bw{r, 0)=0 при г > г я и Bw{0, -i) |
— |
0 |
при t > ^L , Юдин |
построил мажоранту и миноранту функции |
B |
W |
(г, Д , что дало ему |
возможность рассчитать зависимость условного коэффициента тур
булентной |
диффузии для тяжелых |
частиц |
KW = \lm-^— |
от |
|||||||
|
|
|
|
|
|
|
|
|
/~~ dt |
||
скорости гравитационного оседания w в двух вариантах, |
|
причем |
|||||||||
истинные |
значения |
KW |
должны находиться между |
этими |
вариан |
||||||
тами. |
|
|
|
|
|
|
|
|
|
|
|
В целях получения более наглядных результатов целесообраз |
|||||||||||
но приближенно положить |
B L ( Т ) в соответствии с |
(1.33), a B E ( Х ) |
|||||||||
также аппроксимировать |
экспоненциальной |
функцией. |
|
(Для |
|||||||
сохранения соотношения |
масштабов необходимо |
только |
|
вместо |
|||||||
гЕ ввести |
г*Е~\,ЪгЕ |
) . |
|
|
|
|
|
|
|
|
|
Для значений B W (г, т) приближенно можно принять |
|
|
|||||||||
|
BJr, |
,)= |
< |
^ |
> |
е х р |
[ - ] / ^ ) |
2 + |
^ |
|
(2.25) |
Положив г = ш т |
и проинтегрировав |
выражение |
(1.10) при условии |
||||||||
(2.25), получим |
|
2<г>2>т2 |
i_t_ |
_J_ |
|
|
|
|
|||
|
2 |
|
, |
• |
|
(2-26) |
|||||
|
|
|
|
|
|
|
1 + e |
|
44
где |
|
|
|
|
|
|
|
t, |
г. w |
|
(2.27) |
|
|
|
|
|
|
Таким образом, |
выражение для а\ |
отличается от аналогичного |
|||
выражения в случае w = Q только тем, что масштаб i L |
заменяется |
||||
на т™ . Для Кщ |
в этом случае при w -у с» имеем |
|
|
||
|
|
^2 - -> — , |
|
(2.28) |
|
что совпадает |
с |
предельной зависимостью, найденной |
Юдиным |
||
(1959) и Смитом |
(1959). |
|
|
= UvE , |
|
Если принять, согласно гипотезе |
замороженности, |
г £ |
|||
где U — средняя скорость ветра, то |
|
|
|
2.1.4. Поперечное рассеяние |
тяжелой примеси |
в неоднородном |
слое |
Объединяя схему учета неоднородности турбулентной среды со |
|
схемой учета влияния гравитационного |
оседания, будем считать, |
что частица, выйдя из источника, расположенного на уровне Н, по
падает |
на уровень z0 под совместным |
действием сноса средним |
ветром |
U(z) и нисходящего движения |
со скоростью W=w+ws , |
где o/s , как и раньше,;—средняя скорость пульсационного движе
ния частиц вдоль их траектории, a w—скорость |
гравитационного |
|||||||
оседания. Интегрирование |
корреляционной |
функции для расче |
||||||
та а2 в этом случае будем вести вдоль пути |
(2.18), где ws |
заменим |
||||||
величиной W. Выражение для корреляционной функции, объединяя |
||||||||
(2.18) |
и (2.25) |
, |
примем |
в виде |
(2.25), |
где |
величины %L (z) |
|
и r*E (z) |
будем |
считать зависящими от z, a r=w%, что определяет |
||||||
отклонение траектории тяжелой частицы от траектории |
жидкой |
|||||||
частицы |
среды, |
вызванное |
гравитационным |
оседанием. |
Выраже |
|||
ние (2.25) в этом случае принимает вид |
|
|
|
|||||
|
B(z, t) = |
< ^ > e x p |
Г |
|
|
(2.30) |
||
|
V |
> + |
|
|||||
|
|
|
|
|
|
|
||
Естественно |
предположить, что |
tL (z) и |
rE{z) |
меняются с вы |
сотой одинаковым образом. Тогда, по аналогии с (2.27), введя обо значение
где |
vi — безразмерная величина, |
не зависящая от |
z, |
видим, |
что |
учет |
влияния гравитационного |
оседания приводит |
к |
такому |
же |
уменьшению эффективного лагранжева масштаба, что и в отсутст
вии неоднородности. Пропорционально уменьшается и масштаб |
хн> |
||||
Очевидно, что формула для a^(Q , а также предельные |
соотноше- |
||||
|
|
|
* |
х |
|
«ия остаются |
без изменения и'Имеют вид (2.22), где 1 =—. Т а к ж е , |
||||
|
|
|
хн |
|
|
как и в однородном слое, |
при£->-0 |
влияние гравитационного |
осе |
||
дания исчезает, а при ш->-со характер предельной зависимости |
KW |
||||
от w тот же. |
|
|
|
|
|
2.2. |
Применение |
полуэмпирического уравнения |
|
||
|
(вертикальное рассеяние) |
|
|
||
|
2.2.1. |
Невесомая |
примесь |
|
|
Вертикальное рассеяние невесомой примеси от стационарного линейного источника, расположенного в нормальном ветру направ лении *, подчиняется уравнению
|
U { z ) d l = |
J - |
k { 2 ) ° ± . |
(2.32) |
||
|
|
dz |
dz |
dz |
|
|
Это же уравнение определяет функцию а{(х, z) |
в выражении |
(2.1) |
||||
для точечного |
источника, |
который |
считаем расположенным |
в точ |
||
ке z=H, х=0 |
и имеющим |
скорость |
испускания |
Q. Будем считать, |
что условие на подстилающей поверхности может соответствовать отсутствию взаимодействия ее с примесью
3±=о |
(2.33) |
dz |
|
или полному ее поглощению |
|
<7 = 0. |
(2.34) |
Уравнение (2.32) решалось многими авторами для различных случаев задания U(z) и k(z) (см. обзоры в книгах Монин, Яглом, 1965; МАЭ, 1968), однако для широких практических применений в настоящее время пригодны только два результата: степенная модель, в которой полагается
U(z) = u1(—Y |
(2.35) |
к(г) = ь(Л.у*, |
(2.36) |
* В дальнейшем всегда будем считать, что линейный источник |
раоположен |
нормально к направлению ветра. |
, |
и модель Юдина — Швеца, использованная Берляндом и соавтора ми (Берлянд и др., 1964б) для численного решения при некотором наборе исходных параметров и логарифмическом профиле ветра.
Степенная модель позволяет, варьируя показатели т и б, вы
яснить общую зависимость результата от характера |
вертикальных |
||||||||||
измененией U.(z) и k(z). |
Выражения |
для распределения |
концент |
||||||||
рации в этом случае приведены в книге |
Лайхтмана |
(1970). Поль |
|||||||||
зуясь ими, для 2 = 0 |
легко |
получить |
при условии |
(2.33) |
|
||||||
|
q(x, |
0) = |
|
Qe |
х |
I R \ Р |
|
|
(2.37) |
||
|
Г(1 + |
р ) с 7 я Я \ |
х |
|
|
||||||
где |
|
|
|
|
Н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2.38) |
||
|
|
|
(1 + т + Ъ)*Ь(Н) |
|
|
||||||
|
|
|
|
|
|
||||||
некоторый |
масштаб, |
|
|
|
|
|
|
|
|
||
|
|
|
b(z) |
= U{z)z |
|
|
|
(2.39) |
|||
безразмерный параметр, |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
\+т |
|
|
|
|
(2.40) |
|
|
|
|
Р = 1 + от + й |
|
|
|
|||||
|
|
|
|
|
|
|
|||||
U{H) |
средняя |
в слое от 0 до Я |
скорость |
ветра |
всоответ- |
||||||
\-\-m |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
ствии с (2.17). |
|
|
|
|
|
|
|
|
|
||
В случае условия на |
границе (2.34) для потока примеси при |
||||||||||
z=0 имеем |
(обозначения те же) |
|
|
|
|
|
|||||
|
|
dq_ |
|
|
Qe |
|
RV-f |
|
(2.41) |
||
|
|
dz |
,o |
Г ( 2 - Р ) ( 1 + / и + 8 ) * \ х |
|
|
|||||
|
|
|
|
|
Таким образом, зависимость приземной концентрации от х оп ределяется некоторой комбинацией показателей т и б и безразмер ным параметром (2.39), зависящим от высоты источника. Можно ввести также
(2.42)
безразмерный параметр на уровне измерения и\ и k\ и |
|
К |
(2.43) |
|
47
где |
К — среднее |
в слое от 0 до Я значение |
k(z). |
|
|
|
||||||||||
Между параметрами |
6 Ь |
В и Ь(Н) |
имеет |
место |
соотношение |
|||||||||||
|
|
В = ± ± ^ |
b(H) |
= bt |
^ |
f |
i |
^ |
r 5 |
, |
|
- |
( 2 .44) |
|||
|
|
|
2 — 5 v |
; |
|
2 - 8 [И J |
|
|
|
V |
' |
|||||
Приземная |
концентрация примеси |
от стационарного |
точечного |
|||||||||||||
источника на высоте Я при граничном |
условии |
(2.33) выражается |
||||||||||||||
в виде (2.1), где |
о у |
(л:) |
и |
q{(x,Q) |
|
определяются |
соотношениями |
|||||||||
(2.2) и (2.37); она имеет максимальное |
значение |
|
|
|
|
|||||||||||
|
|
|
|
|
_ |
Q(tt+ |
p)4-pJ(«+9) |
|
|
|
|
|
||||
|
|
|
|
|
9 ° - Г ( 1 + р ) 1 У „ Я а , Д « |
|
|
|
|
{ Z A b ) |
||||||
на |
расстоянии до проекции источника на поверхность земли |
|||||||||||||||
|
|
|
|
|
|
|
*о = — |
• |
|
|
|
|
|
|
(2-46) |
|
|
|
|
|
|
|
|
р - f а |
|
|
|
|
|
|
|
||
вид |
Аналогичные |
выражения в случае |
полного |
поглощения |
имеют |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
Q ( 2 - p |
+ g ) g - P 4 « g - ( 2 - P + « > |
|
|
|
|
||||||
|
|
|
Р ° |
|
|
Г ( 2 - р ) а ^ 1 + - |
|
|
|
|
|
|
||||
|
|
|
|
|
|
* 0 = |
|
|
|
|
|
|
|
|
(2.48) |
|
|
|
|
|
|
|
|
2 — р + <* |
|
|
|
|
|
|
|||
для |
В табл. 2.1 приведены значения х0 |
и ?o/Q, |
выраженные через Б, |
|||||||||||||
линейного |
источника |
и точечного |
при |
условии |
(2.33) |
и а = 1 , |
||||||||||
ау= |
Ьу. На основе этих |
формул |
расстояние |
до зоны |
максималь |
|||||||||||
ной концентрации и ее значение можно выразить в виде |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
л |
Н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q0 = D,-£- |
|
|
' |
|
|
|
(2.49) |
||
для линейного источника и |
|
|
|
|
|
|
|
|
|
|
||||||
для |
точечного. |
Крайние |
значения |
безразмерных |
множителей |
|||||||||||
А и D, зависящих от формы профилей, даны в табл. 2.2 при усло |
||||||||||||||||
вии, что т и б меняются от 0 до 1. |
|
|
|
|
|
|
|
|
||||||||
|
Из |
табл. 2.2 |
следует, |
|
что зависимость |
таких |
характеристик |
|||||||||
распределения |
приземной |
концентрации, |
как до и х0, |
от |
формы |
профилей сравнительно невелика. При согласовании средних в слое от 0 до Я значений U(z) и k(z) и при отсутствии информации о значениях т и б (в пределах изменения их от 0 до 1) в худшем случае можно ошибиться не более чем вдвое.
Таблица 2.1
Значения -JJ- и -Q- ДЛЯ линейного и точечного источника
|
|
при степенном задании профилей |
|
|
|
||||
величина |
I |
|
|
Источник |
|
|
|
||
|
линейный |
|
|
точечный |
|
||||
|
|
|
|
|
|||||
XolH |
|
[fi(l + m + o)(2-B)]-i |
Р 2B(\ |
+ m |
+ - М |
-1 |
|||
|
(2-6) |
||||||||
lolQ |
|
|
рр e—t |
B(I+p)i+P |
e-(i + p)(2—о) (1+/и) |
||||
|
T0+p)U„H |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Таблица 2.2 |
||
|
Значения множителей А и D формул |
(2.49) — (2.50) |
|||||||
|
|
при крайних значениях /га и б (0 и единица) |
|
|
|||||
Множитель |
|
|
Источник |
|
|
|
|
||
т |
линейный |
|
|
точечный |
|
||||
|
|
|
|
|
|||||
|
|
|
о = 0 |
a=i |
|
о=0 |
|
5=1 |
|
|
А |
0 |
0,50 |
0,50 |
|
0,25 |
|
0,17 |
|
|
|
1 |
0,25 |
0,33 |
|
0,12 |
|
0,13 |
|
|
D |
0 |
0,38 |
0,48 |
|
0,43 |
|
0,fi9 |
|
|
|
1 |
0,38 |
0,44 |
|
0,86 |
|
0,88 |
|
В работе Бызовой (1972) множители |
А и D были |
рассчитаны |
|||||||
по результатам численного решения уравнения |
(2.32), опублико |
||||||||
ванным Берляндом и соавторами |
(1964) |
(модель Юдина — Швеца |
|||||||
для k(z). |
Оказалось, что в целом |
значения |
этих множителей укла |
||||||
дываются |
в пределы, |
предписываемые |
степенной |
моделью при |
|||||
m = 0,15 (что соответствует аппроксимации |
логарифмического про |
||||||||
филя скорости |
ветра) |
и изменению б от 0 до 1. Из этих сравнений |
|||||||
следует, что различные способы задания |
профилей |
скорости ветра |
и коэффициента диффузии дают близкие результаты, если согла
совать их средние в слое от 0 до Я значения. |
|
2.2.2. Уравнение и граничное условие для оседающей |
примеси |
Для расчета вертикального рассеяния оседающей примеси вме
сто уравнения (2.32) имеем уравнение |
|
||
U(z)*L-w^- |
= -?-k(z)*L, |
(2.51) |
|
дх |
dz |
dz |
dz |
где со — скорость гравитационного |
оседания |
примеси. |
|
Граничное условие на |
подстилающей поверхности в общем |
4-1294 |
49 |