🗊Презентация Движение атмосферы

Раздел III. Рассеяние примесей в атмосфере В настоящее время в практике инженерных расчетов, базирующейся на исследованиях макрогидродинамической обстановки, считается, что перенос загрязнений в атмосфере обуславливается двумя составляющими атмосферного движения: полем среднего ветра, переносящего загрязнение от одной точки к другой (детерминированная составляющая); турбулентными движениями, которые рассеивают загрязнения относительно некоторого центра (случайная составляющая). В связи с этим можно выделить следующие особенности процессов движения: 1. Особенности атмосферных движений вблизи земли; 2. Влияние трения атмосферы о поверхность земли; 3. Эффект дифференциального потока теплоты на поверхность; 4. Влияния поля ветра на турбулентность и перенос загрязнения. Вывод: Большая часть этих процессов, важных для принятия инженерных решений, реализуется в планетарном пограничном слое.

Движение атмосферы Справочно: Планетарный пограничный слой – часть атмосферы, расположенный между поверхностью земли и свободной атмосферой, где движения в основном не подвержены влиянию приземного трения, а также охлаждения и нагревания поверхности земли. В первых 50-100 (150) м над поверхностью земли мы имеем дело с приземным пограничным слоем, где движение ветра существенно зависит от шероховатости земной поверхности и вертикального градиента температуры, который в свою очередь зависит от температуры подстилающей земной поверхности.

Например, в случае, когда холодный воздух проходит над теплой поверхностью, в очень теплый весенний день при сильном солнечном сиянии конвективные потоки от поверхности земли могут подниматься до высоты 3000 м. Следует добавить, что конвективная деятельность внутри грозы, по разным причинам, может вообще нарушить понятие планетарного пограничного слоя. Например, в случае, когда холодный воздух проходит над теплой поверхностью, в очень теплый весенний день при сильном солнечном сиянии конвективные потоки от поверхности земли могут подниматься до высоты 3000 м. Следует добавить, что конвективная деятельность внутри грозы, по разным причинам, может вообще нарушить понятие планетарного пограничного слоя. Замечание: Большинство моделей рассеяния в атмосфере, которые используются в настоящее время, не учитывают таких интенсивных движений, поэтому далее ограничимся условиями более нормального пограничного слоя.

В тропосфере сосредоточено 4/5 массы всей атмосферы. Масса последней равна около 15.5*1015т. В этом слое сильно развиты турбулентность, и конвекция, сосредоточена преобладающая часть водяного пара (почти весь), здесь возникают облака, развиваются циклоны и антициклоны. В тропосфере сосредоточено 4/5 массы всей атмосферы. Масса последней равна около 15.5*1015т. В этом слое сильно развиты турбулентность, и конвекция, сосредоточена преобладающая часть водяного пара (почти весь), здесь возникают облака, развиваются циклоны и антициклоны. Внимание! Температура в тропосфере уменьшается с высотой в среднем на 6ОС/км. Верхний переходный слой тропосферы получил название – тропопауза. Далее следует стратосфера – слой атмосферы, расположенный между тропосферой и мезосферой, с нижней границей порядка 8-10 км над полюсами и 16-18 км над экватором.

Замечание! В средней части стратосферы наблюдается максимальная концентрация озона поглощающего большую часть УФ солнечной радиации и предохраняющего живую природу Земли от её вредного воздействия. Замечание! В средней части стратосферы наблюдается максимальная концентрация озона поглощающего большую часть УФ солнечной радиации и предохраняющего живую природу Земли от её вредного воздействия. Верхняя граница стратосферы определяется на высоте порядка 50 км. Верхний пограничный слой стратосферы на высоте 50-55 км получил название – стратопауза. Потом располагается мезосфера – слой атмосферы примерно от 50 до 80 км, где температура падает с высотой примерно от 0 ОС на нижней границе, до – 90 ОС на верхней. Замечание! Далее от высот 80-90 км идет термосфера, где температура снова растет до высот 200-300 км, где достигает значений порядка 1500 ОС, после чего остается почти постоянной до больших высот.

Особенности атмосферных движений вблизи земли очень разнообразны и в основном определяются структурным ландшафтом конкретной местности и ее метеорологическими характеристиками, включая их периодичность. Особенности атмосферных движений вблизи земли очень разнообразны и в основном определяются структурным ландшафтом конкретной местности и ее метеорологическими характеристиками, включая их периодичность. Влияние поверхностного слоя Главное воздействие поверхностного трения – его влияние на профиль ветра, т.е. на изменение скорости и направления с высотой над поверхностью земли. Рассмотрим напряжение сдвига и профиль ветра. Если анализировать ламинарный поток газа между двумя параллельными пластинами, то напряжение сдвига на единицу площади будет пропорционально градиенту скорости, т.е. (1) где τ – напряжение сдвига на единицу площади, μ - динамическая вязкость, не зависящая от w и h. Здесь w – скорость газа, h – расстояние от неподвижной плоскости до одной из движущихся. Для жидкостей иногда удобнее использовать кинематическую вязкость – ν=μ/ρ (2), где ρ – плотность жидкости. Замечание: Однако, в атмосфере практически нет ламинарных потоков. Как правило, мы имеем дело с нижней границей, которая фиксируется у поверхности земли и некоторой границей, которая обычно обусловлена существованием ветра. Соответственно мы строим профиль ветра.

Вывод: Мы вводим в рассмотрение понятие - профиль ветра, вид которого зависит от скорости ветра. Вывод: Мы вводим в рассмотрение понятие - профиль ветра, вид которого зависит от скорости ветра. Известно множество различных профилей ветра, тем не менее основной теоретический профиль обычно задается выражением: (3) где К - постоянная Кармана и (4) - "скорость" трения. Замечание: В пределах слоя воздуха вблизи земли эффект вязкости преобладает. По мере удаления от поверхности, напряжение сдвига становиться частично турбулентным, частично вязким и, наконец, достаточно далеко от земли напряжение сдвига определяется главным образом турбулентными движениями. Слой, в котором вязкое напряжение преобладает, очень тонкий. Внимание! С увеличением скорости ветра над некоторой заданной поверхностью или при заметном увеличении шероховатости подстилающей поверхности при постоянной скорости ветра может наступить такое состояние, когда эффект воздействия сил давления, связанных с вихреобразованием, которое возникает от элементов шероховатости (инерционные силы), станет преобладать над чисто вязкостным напряжением (вязкостные силы). Вывод: Поток в этом случае, как иногда говорят, становится "аэродинамически шероховатым" (потенциальным) и вязкость больше не влияет на профиль.

В таком случае, как оказалось, профиль скорости зависит от параметра шероховатости в виде: В таком случае, как оказалось, профиль скорости зависит от параметра шероховатости в виде: (5) где z0 - параметр шероховатости. Замечание: Вообще говоря, порядок величины параметра меньше, чем действительная высота элемента шероховатости. Для гладкой плоской поверхности, снега, и даже скошенной травы, его типичные значения составляют менее 1 см, для травы высотой 4-5см они превышают 1см и для травы высотой 60-70 см увеличиваются до 10 см. На больших площадях, где имеются кусты, деревья, заборы, живые изгороди и тому подобные препятствия, этот параметр может увеличиваться до 2 м. Другой способ представления профиля скорости ветра в виде функции от высоты, который применяется в моделях и связывает скорость приземного ветра со скоростью ветра на высоте трубы, использует соотношение: (6) здесь r - соответствует некоторой заданной высоте, где известна средняя скорость ветра , а z - интересующей нас высоте, на которой мы хотим найти скорость ветра w(z), где α - табличный коэффициент. Замечание: В соответствии со степенным законом средняя скорость ветра на некоторой высоте z равна средней скорости ветра на заданной высоте умноженной на отношение этих высот, возведенная в некоторую степень α.

Замечание: Градиентный ветер на верхней границе пограничного слоя со снижением высоты постепенно подвергается воздействию приземного трения, вызывающего в северном полушарии поворот направления ветра вниз. Замечание: Градиентный ветер на верхней границе пограничного слоя со снижением высоты постепенно подвергается воздействию приземного трения, вызывающего в северном полушарии поворот направления ветра вниз. На правом рисунке были показаны векторы и их проекции на подстилающую поверхность, демонстрирующие поворот ветра в этом случае. Соответственно проекция вектора ветра на поверхность земли меняется. Теоретически спираль Экмана дает разницу примерно в 45О между направлением ветра у поверхности и ветром на верхней границе пограничного слоя. Замечание: Однако, экспериментальные наблюдения показали, что изменение направления часто оказываются меньше 6О и обычно не превышают 30О. Вывод: Спираль Экмана теоретически важна, но на практике имеет существенно меньшее значение.

Влияние температуры подстилающей поверхности Влияние температуры подстилающей поверхности Тепловая конвекция или конвекция некоторых свойств материи означает перенос этого качества от одного места к другому посредством движения воздуха. Конвекция теплоты рассматривается применительно к процессам загрязнения воздуха как механизм вертикального переноса. Свободная конвекция имеет место тогда, когда движение целиком или частично обуславливается только собственной плавучестью, возникающей при нагревании воздуха, непосредственно прилегающего к земной поверхности. Вынужденная конвекция – движение, которое вызывается другими причинами. При полностью вынужденной конвекции плавучесть не влияет на движение и коэффициент теплопередачи. Она возникает только тогда, когда земля намного теплее проходящего над ней воздуха. Согласно климатологическим данным такие явления могут возникать только в 10-15% времени. Замечание: Однако за эти промежутки времени турбулентная структура атмосферы становится заметно отличающейся от обычной, и рассеяние примеси, выброшенной в атмосферу, сильно увеличивается в вертикальном направлении.

Замечание: Обычно при конвективных условиях вблизи приподнятых источников наблюдаются очень высокие концентрации. Струя, как говорят, петляет или имеет волнообразную форму, если смотреть на неё сбоку. Замечание: Обычно при конвективных условиях вблизи приподнятых источников наблюдаются очень высокие концентрации. Струя, как говорят, петляет или имеет волнообразную форму, если смотреть на неё сбоку. Вывод: Вертикальная температурная структура вблизи земли зависит от нагрева и в свою очередь влияет на плотностную структуру атмосферы. Это определяет вид струи. Различают струи: 1. Веерообразную (устойчивые условия) 2. Надымляющую или задымляющую (устойчивые условия наверху и нейтральные внизу) 3. Волнообразную (неустойчивые условия) 4. Конусообразную (нейтральные условия) 5. Приподнятую (устойчивые условия внизу, нейтральные наверху

Различные типы дымовых струй, наблюдаемых в атмосфере

Локальные (местные) ветровые потоки Локальные (местные) ветровые потоки Существуют определенные локальные эффекты, которые влияют как на траектории, так и на турбулентность поля ветра. Один из таких локальных потоков, который может оказывать влияние на траекторию веществ – долинный поток.

Другой тип возмущения локального ветрового потока – морской бриз, который возникает вследствие разницы температур между сушей и водной поверхностью, которая приводит к возникновению циркуляции воздуха вблизи береговой черты. Другой тип возмущения локального ветрового потока – морской бриз, который возникает вследствие разницы температур между сушей и водной поверхностью, которая приводит к возникновению циркуляции воздуха вблизи береговой черты. Внимание! Примесь, перемещающаяся внутри этой циркуляции, будет возвращаться, и увеличивать загрязнение. При этом в действительности бризовая циркуляция намного сложнее. Воздушный поток, непрерывно движущийся параллельно поверхности земли, увеличивает диффузию в результате взаимодействия между воздухом, который дует над поверхностью воды и который подвергается воздействию суши. Это вызывает увеличение перемешивания и приводит к высоким приземным концентрациям.

Еще один третий важный фактор, влияющий на атмосферное загрязнение – влияние города. Еще один третий важный фактор, влияющий на атмосферное загрязнение – влияние города. Город действует не только как источник загрязнений, но и выступает в роли "острова" тепла, т.е. как область, имеющая более высокую температуру, чем окружающая местность. Это способствует возникновению зимой и летом восходящего движения воздуха у поверхности земли. В результате, возникает перенос загрязнений внутрь города, слабое восходящее движение воздуха в центре и нисходящее движение загрязнений на его окраине. Это движение в верхней части общего потока может способствовать появлению более высоких концентраций на уровне земли, чем можно было ожидать в другом случае. Вывод: Движение атмосферы, вызванное локальными ветровыми потоками, во многом определяется не только топографическими препятствиями, но и зонами циркуляции, которые могут приводить к резкому увеличению загрязненности атмосферного воздуха.

Турбулентность и флуктуации Турбулентность и флуктуации Турбулентность – та характеристика атмосферного движения, которая ответственна за рассеяние примеси, попавшей в воздушный поток. "Механическая" турбулентность вызывается элементами шероховатости подстилающей поверхности, – будь то короткая трава, которая дает очень маленький эффект, или строения, которые очень заметно влияют на турбулентность. Замечание: Турбулентность постоянно генерируется и разрушается в атмосфере, а её интенсивность, обычно, циклически меняется со временем суток.

Вывод: Исследования поведения ветра позволили установить, что флуктуации направления ветра, в первом приближении, подчиняются закону нормального распределения и, следовательно, могут быть охарактеризованы при помощи величины стандартного отклонения (М=0). Вывод: Исследования поведения ветра позволили установить, что флуктуации направления ветра, в первом приближении, подчиняются закону нормального распределения и, следовательно, могут быть охарактеризованы при помощи величины стандартного отклонения (М=0). Это стандартное отклонение связано с рассеянием вещества в диффундирующей струе и является основой для гауссовых моделей рассеяния. Аналогичным образом, т.е. близко к нормальному распределению, может быть описано рассеяние поперечных и вертикальных флуктуаций и также охарактеризовано стандартным поперечным и вертикальным отклонением. Замечание: Однако экспериментальные данные по этому вопросу применительно к вертикальным отклонениям носят ограниченный характер.

Обобщение особенностей Обобщение особенностей Таким образом, можно считать, что вещество (загрязнение, примесь) выброшенные в атмосферу: переносится потоком ветра (детерминированная составляющая) и рассеивается турбулентными вихрями (случайная составляющая). Скорость ветра увеличивается с увеличением расстояния от поверхности земли до верхней границы планетарного пограничного слоя. Направление ветра гораздо меньше меняется с высотой, чем это предсказывает теория. Ветровой поток зависит также от разности температур. Береговой и морской бризы влияют на направление переноса примеси. Направление потока может каналироваться или изменяться топографией. Препятствия, образуемые рельефом местности, влияют на турбулентность и рассеяние.

Выделяют два вида турбулентности: Выделяют два вида турбулентности: механическая, вызванная элементами шероховатости, вихри которой не столь велики конвективная, которая возникает при нагревании поверхности, где вихри способны достигать очень больших размеров (торнадо). Механическая турбулентность в сочетании с высокими скоростями относительно однородна, причем максимальные размеры вихрей невелики. Существенное влияние плавучести указывает на значительно более медленное изменение турбулентности и более крупные размеры вихрей. Размеры вихрей определяют максимальный размер облака рассеивающейся примеси, чем больше будет облако, тем меньше будет концентрация примеси в его центре. Вывод:

Лекция 4. Теории рассеяния. Градиентная и статистическая теории Замечание: Проблема турбулентной диффузии в атмосфере еще не сформулирована в том смысле, что пока нет возможности объяснить все важные аспекты этой проблемы при помощи одной единой модели. В настоящее время обычно пользуются двумя альтернативными подходами: 1. Теория градиентного переноса (К-теория) 2. Статистическая теория. Первый подход основан на применении дифференциального уравнения (системы уравнений), описывающего адвективное и турбулентное движение в атмосфере и их влияние на рассеяние атмосферной примеси. С точки зрения градиентного переноса диффузия в фиксированной точке пропорциональна градиенту концентраций. При этом свойства атмосферы рассматриваются в закрепленной (неподвижной) системе координат. Справочно: Адвекция – доставка, горизонтальный перенос воздуха и (вместе с ним) тех или иных его свойств, например, теплоты, влаги, запыленности. Второй подход использует основы статистической теории и теории вероятностей. Статистические диффузионные теории рассматривают происходящие в атмосфере движения как движения центра облака рассеивающего вещества. Оба этих подхода рассматриваются как теории непрерывных движений (непрерывных процессов) и теории дискретных движений (дискретных процессов) в зависимости от того, как моделируется движение частиц, как непрерывное или как дискретные шаги. Замечание: Обе теории представлены для атмосферы и между ними прослеживается тесная взаимосвязь!!!

Математическая схема этой теории взята из физики тепло- и электропроводности. Математическая схема этой теории взята из физики тепло- и электропроводности. Она базируется на известном уравнении: (7) С граничными (краевыми) условиями для точечного источника: при t→∞ при t→0 для всех х, кроме х=0 где , так что Здесь Q - мощность источника, т.е. полный выброс вещества, представляемого через концентрацию В трехмерной пространственной форме для некоторой точки на поверхности земли уравнение записывается в виде: (8)

В этом случае нетрудно увидеть, что степень изменения концентрации q со временем связана с градиентами по всем трем направлениям, а также с коэффициентами диффузии и их изменением по тем же направлениям. В этом случае нетрудно увидеть, что степень изменения концентрации q со временем связана с градиентами по всем трем направлениям, а также с коэффициентами диффузии и их изменением по тем же направлениям. Краевые условия при этом имеют тот же вид, что и для точечного источника. Замечание: Представленное уравнение является классическим уравнением теплопроводности и его решения для различных случаев, полученные различными способами, хорошо известны. Фундаментальное решение этого уравнения известно как функция Гаусса: (9) где и Здесь х - расстояние от центра диффундирующего облака. Коэффициенты a и b могут быть получены, как из условий нормировки, так и из уравнения неразрывности, которое дает полную мощность выброса.

Соответственно, для мгновенного точечного источника, производящего выброс мощностью Q в момент времени t=0 можно записать: Соответственно, для мгновенного точечного источника, производящего выброс мощностью Q в момент времени t=0 можно записать: (10) Представленное решение для одномерного случая может быть записано и для трехмерного случая турбулентности в виде: (11) Поскольку дифференцирование выполняется для коэффициента K в интересующем нас направлении, т.е. по координате (x, y или z) соответствующее интегрирование этого уравнения позволяет рассмотреть другие источники, помимо мгновенного точечного источника. Например, интегрирование по пространству дает уравнение для объемного мгновенного источника, такого как взрыв. Соответственно, интегрирование по времени дает решение для непрерывного действующего точечного источника. Замечание: Одна из трудностей K-теории состоит в предположении постоянства турбулентной диффузии в пограничном слое, в котором существуют явно выраженный сдвиг среднего ветра и большая изменчивость вертикальных температурных градиентов вследствие различного нагрева подстилающей поверхности.

Соответственно, нет единства по этому вопросу. Разные авторы по разному определяют коэффициент диффузии K, как функцию устойчивости и высоты над поверхностью. Соответственно, нет единства по этому вопросу. Разные авторы по разному определяют коэффициент диффузии K, как функцию устойчивости и высоты над поверхностью. Например, коэффициент вертикальной диффузии может быть представлен в виде степенной функции, подобной той, которая используется для представления скорости ветра в зависимости от высоты над подстилающей поверхностью. Некоторые авторы вводили изменяющиеся значения коэффициентов боковой и вертикальной диффузии для объяснения изменения направления ветра с высотой в планетарном пограничном слое. Однако получение аналитических решений в этих случаях крайне затруднено. Замечание: Другим ограничением K-теории является предположение о том, что диффузия, пропорциональна градиенту. В общем случае это допущение до сих пор в строгом виде не доказано. Однако, его правомерность в тех или иных случаях теории рассеяния подтверждена экспериментально. Вывод: В настоящее время K-теория используется, главным образом, для таких территорий как города, где обычно наблюдаются слабые градиенты и имеются много поверхностных (площадных) и линейных источников.

Наиболее часто используемую форму математического описания переноса и рассеяния примеси в атмосфере, записанную для простоты в одномерном виде можно записать так: Наиболее часто используемую форму математического описания переноса и рассеяния примеси в атмосфере, записанную для простоты в одномерном виде можно записать так: (12) где с - обозначает концентрацию интересующего нас вещества; w - средняя скорость ветра в направлении x; K- коэффициент вихревой или турбулентной диффузии; x - расстояние до источника и t - время. Первый член в правой части представляет собой перенос средним ветром, а второй диффузию за счет турбулентности. Замечание: Если первый член отсутствует (или исключен) мы имеем дело с описанием рассеяния примеси относительно центра покоящегося (или движущего) облака (клубы дыма). Мы можем получить гауссово уравнение для концентраций из уравнений К-теории. Однако при этом надо сделать и проверить следующие основные допущения: 1. Решение не зависит от времени. 2.Скорость ветра не меняется ни поперек общего потока, ни вертикально. 3. Коэффициенты диффузии не зависят от местоположения (не являются функцией координат), хотя в реальной атмосфере это далеко не так. 4. Диффузия в направлении оси основного движения мала по сравнению со средним потоком переноса примеси в этом направлении.

Измерить коэффициенты рассеяния К-теории в атмосфере трудно, поэтому были разработаны модели, исходно использующие нормальное распределение, как средство представления рассеяния газа или мелких частиц в атмосфере. Измерить коэффициенты рассеяния К-теории в атмосфере трудно, поэтому были разработаны модели, исходно использующие нормальное распределение, как средство представления рассеяния газа или мелких частиц в атмосфере. В системе координат, которая используется для выполнения расчетов по гауссовым моделям рассеяния, предполагается, что ось х совпадает с горизонтальным направлением распространения струи вдоль среднего ветра, ось у перпендикулярна оси х в горизонтальной плоскости, а ось z направлена вертикально, обычно, образуя правую тройку векторов.

Основное уравнение рассеяния в гауссовой теории имеет вид: Основное уравнение рассеяния в гауссовой теории имеет вид: (13) где C - концентрация в некоторой точке с координатами (x,y,z) при данном выбросе мощностью Q на высоте H, представленная гауссовым распределением. Высота H соответствует положению осевой линии струи для приподнятого источника и равняется сумме геометрической высоты трубы h и подъема струи ΔH; Q - постоянная мощность выброса; w - действующая на струю средняя скорость ветра. При этом обычно делаются следующие основные допущения: 1. Рассеяние струи в горизонтальной и вертикальной плоскости описывается гауссовым распределением со стандартными отклонениями распределения концентрации σY и σz по осям y и z соответственно. 2. Нет осаждения или взаимодействия струи с подстилающей поверхностью, т.е. имеет место отражение струи от поверхности земли. Для описания этого эффекта экспоненциальный член, содержащий z и H - удваивается. 3. Турбулентность всюду одинакова, а σY и σz - функции расстояния от источника. 4. Скорость ветра постоянна во всем слое распространения струи. 5. Направление ветра не меняется при движении потока.

Замечание: В этой модели обычно в качестве единиц измерения используются: для концентрации C – кг/м3 (г/м3); для мощности выброса Q - кг/с (г/с); для w - м/с; σY и σz , H и все координаты – м. Замечание: В этой модели обычно в качестве единиц измерения используются: для концентрации C – кг/м3 (г/м3); для мощности выброса Q - кг/с (г/с); для w - м/с; σY и σz , H и все координаты – м. Существует много уравнений, которые могут быть использованы для отдельных случаев рассеяния, которые обычно заимствуются из теории броуновского движения и зачастую много проще основного уравнения. Пример: Когда концентрации рассчитываются только на поверхности земли, и предполагается полное отражение примеси от поверхности, можно пользоваться уравнением: (14) Здесь исключены второй экспоненциальный член и координата z в экспонентах, т.к. на уровне земли z=0. Пример: Если нас интересуют только концентрации вдоль осевой линии струи, где они максимальны, то можно использовать следующую форму записи: (15)

В этом уравнении экспоненциальный член, содержащий y и σz, опущен, т.к. осевой линии соответствует случай y=0. В этом уравнении экспоненциальный член, содержащий y и σz, опущен, т.к. осевой линии соответствует случай y=0. Если источник наземный (не приподнятый) может использоваться следующая форма: (16) Здесь концентрация вдоль осевой линии струи на уровне земли определяется как функция мощности выброса, скорости ветра и рассеяния примеси вдоль направления ветра. Вывод: Эти упрощенные гауссовы модели и целый ряд их модификаций нашли широкое применение у нас в стране и за рубежом при решении многих задач охраны окружающей среды. Замечание: В частности, было произведено включение этих уравнений и других модификаций в модели качества воздуха, которые используются агентством охраны окружающей среды (АООС).

Коэффициенты К-теории. Коэффициенты, которые используются в различных моделях этой теории обычно считают зависимыми от высоты над поверхностью земли и числа Ричардсона или градиента температуры. Коэффициенты К-теории. Коэффициенты, которые используются в различных моделях этой теории обычно считают зависимыми от высоты над поверхностью земли и числа Ричардсона или градиента температуры. Число Ричардсона Ri - является функцией термической устойчивости или неустойчивости атмосферы и определяется соотношением: (17) где g - ускорение свободного падения; T - средняя температура воздуха в слое ∂z; - градиент потенциальной температуры; ∂w/ ∂z - вертикальный градиент скорости ветра. Справочно: В общем случае, число Ричардсона – критерий подобия в гидродинамике, равный отношению потенциальной энергии тела, погруженного в жидкость к его кинетической энергии. Под «телом» понимается рассматриваемая жидкость или газ, при этом число Ричардсона определяется следующим образом: где ρ – плотность тела, Δρ - разность плотностей тела и среды, g – ускорение свободного падения, L – характерная длина, w – характерная скорость.

Для расчета коэффициентов горизонтальной диффузии Лью и Пандольфо предложили следующие выражения: Для расчета коэффициентов горизонтальной диффузии Лью и Пандольфо предложили следующие выражения: для для (18) для где k - постоянная Кармана (~0,4); z - высота; z0 - параметр шероховатости, обычно малый по сравнению с z; α- постоянная Обухова; h - постоянная Пристли. Слейд предложил более упрощенную форму: (19) (20) (21) где σΘ - стандартное отклонение угловых колебаний флюгера, u, v, w - ортогональные компоненты ветра (u - запад-восток, v - север-юг, w - вертикаль).

Замечание: Для вертикальных коэффициентов методы определения столь же разнообразны. В моделях, развитых Пандельфо и Лью, для вертикальной дисперсии используются те же коэффициенты, что и для горизонтальной. Предполагается, что струя рассеивается в виде конуса и в целом симметрична. Замечание: Для вертикальных коэффициентов методы определения столь же разнообразны. В моделях, развитых Пандельфо и Лью, для вертикальной дисперсии используются те же коэффициенты, что и для горизонтальной. Предполагается, что струя рассеивается в виде конуса и в целом симметрична. Коэффициент вертикальной диффузии по ассоциации с уравнением (19) обычно записывают в виде: (22) где z - высота над поверхностью земли, H - высота слоя турбулентного перемешивания, а (23) т.е. рассчитывается из вертикального профиля температуры и ее изменения во времени.

Коэффициент боковой (поперечной) диффузии (СКО) σY представляет собой горизонтальное расширение струи в направлении, перпендикулярном направлению движения. Коэффициент боковой (поперечной) диффузии (СКО) σY представляет собой горизонтальное расширение струи в направлении, перпендикулярном направлению движения. Коэффициент вертикальной дисперсии (СКО) σZ представляет собой расширение струи по вертикали. Справочно: Дисперсионные кривые, изображенные на представленных рисунках, основаны на экспериментальных данных, собранных в США и Великобритании, при этом все эти данные относятся к источникам, выброс из которых производится вблизи поверхности земли, как правило, на высоте не более 10 м и для периодов времени порядка 10 мин. Вывод: Следовательно, строго говоря, коэффициенты дисперсии (СКО), представленные на рисунках, могут быть использованы только для источников, расположенных вблизи земли, и для 10-мин периодов осреднения по времени. Данные, использованные Пасквилом, действительны до расстояния 1 км от источника и для территорий, имеющих климат, подобный климату Южной Англии. Кривые А, B, C, D, E ,F на представленных рисунках соответствуют различным классам устойчивости атмосферы.

Атмосферная устойчивость – это то состояние турбулентности в атмосфере, которое характеризует поле ветра, переносящее и рассеивающее примесь. Атмосферная устойчивость – это то состояние турбулентности в атмосфере, которое характеризует поле ветра, переносящее и рассеивающее примесь. Замечание: Классификация устойчивости была применена Пасквилом для того, чтобы разделить периоды слабой турбулентности, представленные кривой F, и периоды сильной турбулентности, например, в солнечный полдень с легким ветром – кривая А. Паксвил определил свои классы устойчивости в зависимости от поступающей солнечной радиации, эти категории представлены в таблице, где D – нейтральный класс, включая сплошную облачность и днем и ночью.

Пример: Пример: В дневное время, когда приходящая солнечная радиация велика, а скорости ветра малы, наблюдается категория А. При сильной солнечной радиации и скоростях ветра, больших 6 м/с, наблюдается категория С. Т.к. сильные ветры препятствуют развитию термической конвекции в такой степени, какая имеет место при слабых ветрах. То же самое наблюдается в ночное время, если небо менее чем наполовину покрыто облаками, в то время как слабые ветры и относительно безоблачное небо ночью соответствует случаю большей устойчивости. Нейтральный класс D соответствует также случаю сплошной облачности и днем и ночью. Замечание: Все эти таблицы приводятся в руководствах по моделям качества воздуха и проверены АООС в 1980г. Имеются также другие методы определения категорий устойчивости.

Ниже приведена классификация устойчивости по Паксвиллу-Гиффорду как функция вертикальных флуктуаций. Ниже приведена классификация устойчивости по Паксвиллу-Гиффорду как функция вертикальных флуктуаций. Эти значения можно привести к шероховатости подстилающей поверхности путем умножения на (z0/15 cм)0,2, где z0 – средний параметр шероховатости подстилающей поверхности в радиусе 1-3 км от источника Замечание: Данный способ может быть использован, если флюгер дает данные изменения ветра, по которым может быть рассчитано стандартное отклонение вертикальной составляющей направлений ветра для периодов осреднения от 30 мин до 1 часа.

Классы устойчивости по Паксвиллу-Гиффорду могут быть представлены как функции горизонтальных флуктуаций ветра . (Таблица 2) Классы устойчивости по Паксвиллу-Гиффорду могут быть представлены как функции горизонтальных флуктуаций ветра . (Таблица 2) При этом соответствующую категорию рассеяния характеризует диапазон изменения стандартного отклонения осредненного за период 30 мин - 1 час.

Гауссовы коэффициенты дисперсии, как правило, определяются, исходя из данных полученных за 10-мин период. В тоже время, в большинстве моделей качество воздуха рассматривается как среднее значение концентрации загрязнений за период времени, приблизительно равный 1ч и более. Способ, позволяющий получить значение концентраций за период больше 10 мин, состоит в использовании соотношения (27): Гауссовы коэффициенты дисперсии, как правило, определяются, исходя из данных полученных за 10-мин период. В тоже время, в большинстве моделей качество воздуха рассматривается как среднее значение концентрации загрязнений за период времени, приблизительно равный 1ч и более. Способ, позволяющий получить значение концентраций за период больше 10 мин, состоит в использовании соотношения (27): (27) где средняя концентрация за период порядка 60 мин приравнивается к средней концентрации за 10-мин период, и умножается на соотношение выбранных периодов. Величина Р изменяется в пределах 0,17...0,20. Замечание: При использовании моделей Агентства по охране среды США (АООС) обычно не обращают внимания на разницу средних значений концентраций рассматриваемых периодов и часто в качестве осредненных за часовой период коэффициентов принимают коэффициенты дисперсии, осредненные за 10 мин. Это делается с целью получения более осторожных оценок концентрации. Фактически здесь используют метод предельных оценок!!!

В то время как боковое (поперечное) расширение облака увеличивается со временем, вертикальное расширение облака растет существенно меньше. Следовательно, значения вертикального коэффициента дисперсии, полученные за 10-мин период, могут также применяться и для периода один час. В то время как боковое (поперечное) расширение облака увеличивается со временем, вертикальное расширение облака растет существенно меньше. Следовательно, значения вертикального коэффициента дисперсии, полученные за 10-мин период, могут также применяться и для периода один час. Средние значения концентраций за периоды 3 и 24 ч обычно получают путем проецирования и сложения проекций, средних за 1 ч струй при соответствующих направлениях и скоростях ветра за каждый час. Для больших периодов (таких, как месяц, сезон и год) при вычислении концентраций используется соотношение типа (28): (28) Здесь Q - средняя мощность источника за данный период времени; σ2S - соответствует σZ для расстояния в каждом классе устойчивости; f(θ,S,N) - повторяемость направления ветра, класса устойчивости, и градации скорости ветра. В знаменателе, uN - соответствует скорости ветра в классе устойчивости N и в секторе 2πx/16 (обычно предполагается, что окружность 360о делится по направлениям на 16 секторов). В экспоненциальном члене высота Hu представляет собой эффективную высоту струи при скорости ветра u. Члены внутри фигурных скобок суммируются по всем градациям скорости ветра и классам устойчивости.

В теории рассеяния имеются определяющие процессы, которые требуют более детального рассмотрения в связи с гауссовыми моделями дисперсии. В теории рассеяния имеются определяющие процессы, которые требуют более детального рассмотрения в связи с гауссовыми моделями дисперсии. Подъем струи; Химические взаимодействия; Вымывание примеси; Сухое осаждение примеси; Условия замыкания (ограничение переноса), когда существует верхняя "крышка", которая вместе с поверхностью земли ограничивает вертикальное рассеяние загрязнений; Задымление, которое наступает, когда сплошной устойчивый слой воздуха (например, ночью) быстро разрушается снизу в результате нагревания поверхности земли; Затененность зданиями.

В моделях АООС и в большинстве других моделей используются уравнения для подъема струи, полученные Бриггсом. В моделях АООС и в большинстве других моделей используются уравнения для подъема струи, полученные Бриггсом. Эти выражения основываются на теории вовлечения воздуха в струю и плавучести струи, причем используются коэффициенты, определенные эмпирически из наблюдений для широкого класса источников. Уравнения Бриггса имеют вид: где Δh - подъем струи, м; ; g - ускорение свободного падения, м/с2; Vs - средняя скорость выхода газов из трубы, м/с; Rs - внутренний радиус трубы, м; Ts - средняя температура газов в струе, К; Tα - температура окружающего воздуха, К; U - скорость ветра на высоте трубы, м/с; ρ - расстояние от источника до рецептора, м. Δh - представляет собой подъем струи в метрах, который следует добавить к физической высоте трубы, чтобы определить эффективную высоту источника. F - вертикальный поток за счет плавучести (архимедовых сил), который определяется в соответствии с приведенным выражением. Расстояние от трубы по направлению ветра до точки, где прекращается подъем струи, определяется как 3,5X*

Справочно: Имеются и другие уравнения для определения подъема струи, однако они используются не часто по сравнению с уравнениями Бриггса. Обобщение уравнений для подъема струи выполнил Мозес. Он сопоставил порядка 20 уравнений и разработал систему вероятностных оценок действительного подъема струи. При этом расчетные значения подъема струи по Бриггсу, оказались близки к средним величинам, а возможный разброс значений находился в пределах 50%. Справочно: Имеются и другие уравнения для определения подъема струи, однако они используются не часто по сравнению с уравнениями Бриггса. Обобщение уравнений для подъема струи выполнил Мозес. Он сопоставил порядка 20 уравнений и разработал систему вероятностных оценок действительного подъема струи. При этом расчетные значения подъема струи по Бриггсу, оказались близки к средним величинам, а возможный разброс значений находился в пределах 50%.

Падающие капли промывают объем воздуха, через который они проходят. Частицы примеси, находящиеся в этом объеме и следующие по линиям тока воздушного потока, в силу инерции, электрического притяжения и молекулярной диффузии могут частично смещаться с линии тока и пересекать путь падающих капель. Падающие капли промывают объем воздуха, через который они проходят. Частицы примеси, находящиеся в этом объеме и следующие по линиям тока воздушного потока, в силу инерции, электрического притяжения и молекулярной диффузии могут частично смещаться с линии тока и пересекать путь падающих капель. По Слейду, предполагается, что процесс вымывания описывается экспоненциальным законом: (32) где С0 - концентрация в начале вымывания; Λ - коэффициент вымывания, с-1, t - время. В особом случае горизонтально однородных осадков над непрерывной струей максимальное вымывание, возможное на некотором расстоянии, когда дождь начался одновременно с выбросом или раньше его, записывается в виде: (33) где wmax - максимальная скорость вымывания, г/(м2с). Скорость вымывания для осадков, начинающихся непосредственно над струей, в гауссовой модели записывается в виде: (34) где Λ - функция размера частиц, размера капель и интенсивности осадков, причем Λ/u - может изменяться от 10-3 до 10-5 (м-1).

В проблемах загрязнения воздуха наиболее важны два типа осаждения. Если частицы достаточно велики, осадки образуются за счет гравитационного осаждения. В проблемах загрязнения воздуха наиболее важны два типа осаждения. Если частицы достаточно велики, осадки образуются за счет гравитационного осаждения. Для газов и мелких частиц осаждение на поверхность происходит в результате их непосредственного контакта с этой поверхностью. Гравитационное поле земли играет важную роль в осаждении аэрозоля на поверхность земли. Скорость осаждения частиц зависит от баланса между аэродинамической силой трения и силой притяжения земли. Для гладких сферических частиц при условии пренебрежения эффектом скольжения условия баланса может быть записано следующим образом: При ламинарном движении Re = 10-4...1 эти уравнения сводятся к известной формуле Стокса: (35) где В – коэффициент трения.

При ламинарном движении Re = 10-4...1 эти уравнения сводятся к известной формуле Стокса: При ламинарном движении Re = 10-4...1 эти уравнения сводятся к известной формуле Стокса: для гладких сферических частиц: (36) Замечание: Влияние формы частиц на скорость седиментации проявляется в уменьшении скорости седиментации до 2/3 ее значения по отношению к скорости гладких сферических частиц. Для гладких частиц овальной формы изменение скорости седиментации описывается коэффициентом 0,5...1,04. Струя, состоящая из частиц, имеющих некоторую осредненную скорость седиментации, может быть представлена отклонением по высоте ее осевой линии, которое пропорционально скорости осаждения, и тогда общее уравнение для концентраций согласно Оверкампу может быть записано в виде:

(37) где обычные гауссовы члены определяются, как показано ранее и α0(xg) - фактор отражения. Для газов и частиц, которые полностью отражаются от земли, он равен единице. Фактор отражения вводится только у поверхности земли и имеет вид: (38) где uD - скорость осаждения или эффективного захвата газа поверхностью. В точке с координатой xg, α(x,z)=α0(xg) и на расстоянии xg может быть получено методом итерации соотношение: (39) Замечание: Модель Оверкампа используется в комплексной модели АООС для промышленных источников. Особенностью этой модели и модели Слейда является то, что примесь, по предположению, немедленно перераспределяется после осаждения с нижней границы облака. Хорст предложил модель поверхностного истощения, чтобы скорректировать этот эффект, однако она не нашла широкого применения из-за больших затрат машинного времени.

Ограничение переноса – это условия, происходит быстрое перемещение струи под слоем инверсии, в то время как устойчивый инверсионный слой воздуха блокирует рассеяние примеси вверх. Обычно для высоких приподнятых струй эти условия не играют существенной роли, но они могут привести к быстро растущим концентрациям для низких струй в утреннее время, когда инверсия медленно разрушается при нагревании и примесь перемешивается между землей и слоем приподнятой инверсии. Ограничение переноса – это условия, происходит быстрое перемещение струи под слоем инверсии, в то время как устойчивый инверсионный слой воздуха блокирует рассеяние примеси вверх. Обычно для высоких приподнятых струй эти условия не играют существенной роли, но они могут привести к быстро растущим концентрациям для низких струй в утреннее время, когда инверсия медленно разрушается при нагревании и примесь перемешивается между землей и слоем приподнятой инверсии. Один из способов определения концентраций от струи, которая распространяется в таких условиях основан на уравнении: (40) Здесь L - это высота слоя инверсии, ограничивающего рассеяние. Знак суммирования вводится в связи с учетом эффекта отражения, суммирование производится по N от 1 до I. В большинстве случаев достаточно принять I равным 3 или 4, чтобы получить устойчивое решение для концентраций.

Задымление наступает в том случае, когда струя от приподнятого источника распространяется в атмосфере быстро нагревающейся снизу и струя перед этим слабо рассеивающаяся в устойчиво стратифицированном слое атмосферы, начинает интенсивно переноситься к земле в результате возросшего вертикального перемешивания. Обычно такие условия существуют в течение короткого периода времени. Выражение для концентрации при этом имеет вид: Задымление наступает в том случае, когда струя от приподнятого источника распространяется в атмосфере быстро нагревающейся снизу и струя перед этим слабо рассеивающаяся в устойчиво стратифицированном слое атмосферы, начинает интенсивно переноситься к земле в результате возросшего вертикального перемешивания. Обычно такие условия существуют в течение короткого периода времени. Выражение для концентрации при этом имеет вид: (41) где Концентрации при задымлении оценивают интегрально. Значения для интегралов в квадратных скобках могут быть найдены в статистических таблицах для нормального распределения и обычно принимаются во внимание для той части струи, которая распространяется в направлении ветра. Когда инверсия ограничивает эффективную высоту, предполагается, что половина струи перемещается вдоль направления ветра, другая половина остается в устойчивом воздухе наверху. Уравнение может быть преобразовано, если задымление близко к максимальному в виде: (42) где (43)

Для коэффициентов бокового рассеяния добавляется в индекс F. Это делается для преодоления некоторых трудностей, возникающих при оценке примерного значения горизонтального рассеяния, т.к. при перемешивании устойчивой струи по вертикали возникает некоторое дополнительное горизонтальное рассеяние. Для коэффициентов бокового рассеяния добавляется в индекс F. Это делается для преодоления некоторых трудностей, возникающих при оценке примерного значения горизонтального рассеяния, т.к. при перемешивании устойчивой струи по вертикали возникает некоторое дополнительное горизонтальное рассеяние. Если пренебречь этим рассеянием и пользоваться значением σY для устойчивых условий, то вполне вероятно, что расчетные значения концентраций окажутся выше, чем реальные. Тейлор предложил аппроксимацию, которая увеличивает коэффициент бокового рассеяния на величину, равную 1/8 высоты струи: (44) Замечание: Расчеты условий, когда происходит задымление, делаются редко, за исключением специальных экстремальных случаев и обычно не включаются в модели качества воздуха АООС, однако само ограничение рассеяния, в эти модели включается.

Одна из проблем, которая часто встречается в рассеянии примеси от низких источников, состоит в том, что начальное рассеяние существенно увеличивается в тени строений, окружающих источник выброса. Одна из проблем, которая часто встречается в рассеянии примеси от низких источников, состоит в том, что начальное рассеяние существенно увеличивается в тени строений, окружающих источник выброса. Имеется ряд способов учета влияния аэродинамической тени зданий, однако АООС предпочитает способ разработанный Хьюбертом, основанный на исследованиях в аэродинамической трубе. КМПИ-модель АООС включает поправку на эффект аэродинамической тени; программы этой модели учитывают эффект аэродинамической тени зданий, вводя модифицированное значение σZ, когда отношение высоты трубы к высоте зданий больше, чем 1,2, но меньше, чем 2,5; если отношение высоты струи к высоте здания меньше или равно 1,2 вводятся модифицированные значения σY и σZ. При этом здания квалифицируются как низкие (hu ≥hb) или высокие (hu ≤ hb) . Замечание: Ширина здания hu принимается равной диаметру круга, площадь которого примерно равна площади, занимаемой зданием. включает общую процедуру преобразования σY и σZ на расстояниях, больших 3hb для низких зданий и 3hu - для высоких. Внимание! Воздушный поток в аэродинамической тени здания обычно становится более турбулентным и завихренным. КМПИ-модель непригодна для оценки концентраций внутри таких зон.

Модифицированное выражение для σZ в случае низких зданий имеет вид: Модифицированное выражение для σZ в случае низких зданий имеет вид: (45) где высота здания hb принимается в метрах. Для высоких зданий, вместо hb в уравнение (45) обычно подставляется hw и модифицированное выражение принимает вид: (46) где ширина здания hw принимается в метрах. Важно отметить, что σ’Z не может быть меньше, чем в случае точечного источника (см. уравнение (26)). Для расстояний в направлении ветра больших 10hb (низкие здания) или 10hW (высокие здания), вводится вертикальное виртуальное расстояние xz, которое прибавляется к действительному расстоянию по направлению ветра x, чтобы учесть прирост начального подъема струи вызванный аэродинамической тенью здания.

Расстояние xz до виртуального источника А, который выбрасывает облако примеси, эквивалентное по размерам облаку, выбрасываемому в аэродинамическую тень здания. Расстояние xz до виртуального источника А, который выбрасывает облако примеси, эквивалентное по размерам облаку, выбрасываемому в аэродинамическую тень здания. Другими словами, это расстояние до воображаемого точечного источника, на котором размеры рассеивающего облака становятся сравнимыми с размерами аэродинамической тени здания. Таким образом, xz для низкого здания может быть получено в виде: (47) где a=p и b=q - постоянные, зависящие от устойчивости представлены в таблице 6.

Замечание: Для зданий с отношением ширины к высоте hw/hb, большим 5, в настоящее время нет удовлетворительных данных для построения выражения для σY . Замечание: Для зданий с отношением ширины к высоте hw/hb, большим 5, в настоящее время нет удовлетворительных данных для построения выражения для σY . Для здания, у которого ширина значительно больше высоты и труба расположена приблизительно в его центре (т.е. удалена от какого-либо края), следует принимать во внимание только масштаб его высоты. Модифицированное выражение для σY в случае низкого здания, аналогично соотношению с поперечным виртуальным расстоянием хY : (50) Модифицированное выражение для σY в случае высокого здания: (51) (52)