УДК 551.511.32: 532.517.4
DOI: 10.33764/2618-981 Х-2019-4-1-90-98
МОДЕЛИРОВАНИЕ ПРОНИКАЮЩЕЙ ТУРБУЛЕНТНОЙ КОНВЕКЦИИ И РАССЕЯНИЯ ЗАГРЯЗНЕНИЙ В ТЕРМИЧЕСКИ УСТОЙЧИВО СТРАТИФИЦИРОВАННОЙ СРЕДЕ
Людмила Ивановна Курбацкая
Институт вычислительной математики и математической геофизики СО РАН, 630090, Россия, г. Новосибирск, проспект академика Лаврентьева, 6, старший научный сотрудник, тел. (383)330-61-52, e-mail: [email protected]
Явление городского острова тепла и связанная с ним циркуляция, обусловленная энергией, генерируемой антропогенными источниками, наиболее интенсивна в ночное время при ясном небе и слабом окружающем ветре. Фундаментальными характеристиками ночного острова тепла являются его высота перемешивания, интенсивность, распределения кинетической энергии и температуры, а также, индуцируемая островом тепла циркуляция. При моделировании городского острова тепла необходимо аккуратное описание турбулентной конвекции в случае слабого окружающего ветра, когда турбулентное движение в центре города обусловлено, главным образом, архимедовой силой плавучести. Основная цель состояла в развитии модели турбулентного переноса для описания структуры городского острова тепла и рассеяния примесей в условиях устойчиво стратифицированной окружающей среды. Сформулирована четырехпараметрическая модель турбулентного переноса для описания эволюции полностью турбулентного плума над протяженным поверхностным источником тепла и распределения примеси от поверхностных источников.
Ключевые слова: турбулентность, планетарный пограничный слой, городской остров тепла, крупномасштабная циркуляция, рассеяние примесей, численное моделирование.
MODELING OF PENETRATING TURBULENT CONVECTION AND SCATTERING OF POLLUTION IN A THERMALLY SUSTAINABLE STRATIFIED ENVIRONMENT
Lyudmila I. Kurbatskaya
Institute Computational Mathematics and Mathematical Geophysics SB RAS, 6, Prospect Аkademik Lavrentiev St., Novosibirsk, 630090, Russia, Senior Researcher, phone: (383)330-61-52, e-mail: [email protected]
The urban-heat-island phenomenon and its associated circulation, driven by the energy generated by anthropogenic sources are found to be most intense at nighttime under clear skies and weak ambient wind. A thermal plume is generated by an underlying heat island in the form of an area source. The fundamental characteristics of most concern for a nocturnal heat island are its mixing height, intensity, temperature and turbulent kinetic energy distributions, and heat-island-induced circulation. To model the urban-heat-island structure it is necessary to describe correctly the turbulent penetrative convection in the case of near-calm conditions when the ambient wind speed approaches zero and the turbulent motion at the urban center is dominated by the thermal buoyancy force. In most existent mathematical models the turbulence of an urban-heat-island is parameterized. The principal aim of this investigation is the development of turbulent transport model for the simulation of turbulent urban-heat-island structure and pollutant dispersion in the stably stratified environment. The four-parametric turbulent model describing an evolution of fully turbulent low-aspect-ratio buoyant plume with no initial momentum and dispersion of contaminant from a ground-level source under calm and stable stratified conditions is formulated and evaluated.
Key words: turbulence, stable stratification, atmospheric boundary layer, pollution dispersion, modeling.
Введение
При мезомасштабном моделировании структуры полей ветра и температуры город может рассматриваться в виде нагретого "пятна" внутри окружающих его окрестностей [1], а эффекты городской шероховатости на турбулентный перенос импульса, тепла и массы могут быть учтены в параметризованном виде [2, 3]. Поскольку на границе раздела "город-окрестности" развивается внутренний пограничный слой, имеющий как механическое, так и термическое происхождение, желательно сохранить анизотропию турбулентного переноса (в вертикальном и горизонтальном направлениях), которая не может быть адекватно воспроизведена в традиционно используемых одно- и двухпараметрических моделях турбулентного атмосферного пограничного слоя (АПС) (см., например, [2, 3]). Улучшенная мезомасштабная модель турбулентного АПС с анизотропными выражениями для турбулентных потоков импульса и тепла [4, 5] и параметризацией эффектов городской шероховатости на турбулентный перенос позволяет получить [5] картину полей ветра и температуры над урбанизированной поверхностью, которая качественно согласуется с данными наблюдений [6], зафиксировавшими повышение скорости ветра над городом. Кроме того, в силу горизонтальной термической неоднородности, создаваемой городом и его окрестностями, модель позволяет исследовать эффект влияния продольной турбулентной диффузии тепла на такие интегральные характеристики городского пограничного слоя, как его высота [1]. Выражения для турбулентных потоков импульса и тепла модели включают три параметра: кинетическую энергию турбулентности (КЭТ) E = ( 1 / 2)(uiui ), скорость ее диссипации s
и дисперсию турбулентных флуктуаций температуры (б2 ), которые определяются из решения дифференциальных уравнений переноса трехпараметрической модели турбулентности [7, 8].
Определяющие уравнения эйлеровой диффузионной модели
Для описания рассеяния пассивной примеси в атмосферном пограничном
слое (АПС) базовая трехпараметрическая E — s—<6 > - модель турбулентности [9] должна быть дополнена уравнениями для осредненной концентрации C(xi, t), вектора турбулентного потока примеси f =< ще > и корреляции между флуктуациями концентрации и температуры < ев >. Поскольку результаты моделирования переноса активной примеси (тепла) [8] удовлетворительно согласуются с данными измерений [9], для турбулентного потока пассивного скаляра (массы вещества) также использована полностью явная анизотропная алгебраическая модель, полученная [10] в результате упрощения дифференциального уравнения переноса для потока < ще >
D < uic > д I E д < uic >
-i---1 a\s ~ < ukui >-—
Dt dxk I s dxl
dC
dU;
< uiuj >--< ujc >
dxj dxj
S dU-
-aXC e < uic > +a2C < ujc >~dXL + a2cSiß <с0> Siß <с0> С1)
в приближении локально равновесной турбулентности до алгебраического уравнения. В этом приближении левая часть уравнения (1) полагается равной нулю, и уравнение записывается в виде алгебраического уравнения для вектора турбулентного потока концентрации:
E 2 дС E
Aj ■ fj = - - (bj + - ES,) — + (1 - «2C )ßgSi 3 -< C0>, s 3 д x s
(2)
j
где Ъу = <и1и]-)-(2 / 3 - тензор анизотропии. Прогностическое уравнение для корреляции <с0)
D< с0> Dt
дС д© s
- Dlff{ с0>=-< uj 0>---< ujc^ff- ~а3С^ < C0>
j
д©
j
s
E
(3)
также упрощается в приближении слабо-равновесной турбулентности до алгебраического уравнения:
< с0> =
1 E
а3С s
д© дС
< ыс>— + <иj 0>-
д x,
v J
д x
(4)
J J
Подставляя (4) в (2), получаем алгебраическое уравнение для потока fj:
2 , (1 -а2СЬ „ E „ , дС
Bj • fj = - - {(bj + ^ E SJ) + ^ ~2С/ ßgSi з - < UJ 0 >} ,
s 3 a^r s дx.
(5)
3С
где
Bij = а1С Sij +
E ди ( E f (1 - а2С) „ „ д©
s д x,
■ +
v ь /
а
-ßg
3С
д xJ
(6)
Выражения для вертикального <wc) и горизонтального <ис) турбулентных потоков концентрации выводятся с помощью символьной алгебры из (5) и записываются в следующем виде:
-< wс> = — D
E
а1С
< w2 > + аХх < w0>
Л дС 1 ( E
— + — а1С v —
д z D s
<uw) + a X1 <u0> — ,(7) д x
дС
-(ис) = — В
— 8
(и2)Х2 +——— [{им) + а Х1(м0 8 — 2
—С
—.X
1
+—
В
. Е — — I 2\
(uw)X2 +---1(м? ) + аА1
8 — г ^
—С
—2
В (7)-(8) введены обозначения: а = (1 -а 2С) / а 3С, = р g (Е/ 8)
Е —Ж
Xз = ^(— / 8) , X2 = а1С +--+ а X
—0
8 — г — г
С учетом выражений (7)-(8) уравнение для средней концентрации
—С +А С. — + А С ■ Ж = -А( мс)-А( ис)
—I — х — г — г — х
(9)
в пренебрежении (при больших числах Рейнольдса) молекулярной диффузией приобретает замкнутый вид. Откалиброванные численные значения констант диффузионной модели (7) - (9) равны: а1С = а3С = 3,28, а2С = 0,4.
Результаты численного моделирования дисперсии пассивной примеси над городом
Тестовые испытания модели проведены для критической метеорологической ситуации, возникающей в атмосферном пограничном слое в условиях слабого окружающего ветра и устойчивой стратификации атмосферы. Такая ситуация типична для формирования над городом турбулентной циркуляции воздуха - явления, называемого городским островом тепла. Было рассмотрено два сценария. Первый сценарий строился для лабораторного эксперимента [9], в котором проникающая турбулентная конвекция индуцируется поверхностным постоянным источником тепла в виде круглой пластины заданного диаметра. Этот источник тепла моделирует прототип городского острова тепла с малым относительным удлинением (вертикальный линейный масштаб много меньше горизонтального масштаба). Уравнения гидродинамики, описывающие циркуляцию над городским островом тепла малого относительного удлинения, могут быть записаны без учета силы Кориолиса и радиации в цилиндрической системе координат.
Кроме того, может быть принято гидростатическое приближение, а эффекты плавучести учтены в приближении Буссинеска [7, 8]. Необходимые метеорологические параметры такие, как средний ветер, температура, турбулентные величины полей скорости и температуры вычислялись с помощью трехпара-метрической модели турбулентного переноса. Распределения этих величин в турбулентном термическом факеле, полученные для цилиндрической геомет-
рии турбулентной циркуляции над городским островом тепла, приведены в [7], и в настоящих тестах они используются в виде "входной информации" для диффузионных вычислений. Уравнения (1-9) модели записываются в цилиндрических координатах для искомых средней концентрации С(г,г^) и вторых моментов поля концентрации < игс >, < игс > и < св > [7].
Этот сценарий преследовал основную цель - выяснить роль эффектов плавучести на распределение средней концентрации примеси в перемешанном и инверсионном слоях.
Результаты моделирования структуры турбулентной циркуляции над городским островом тепла (различных величин турбулентных полей скорости и температуры) получены ранее, их можно найти в [1,10].
Следует при этом отметить, что поскольку нет данных измерений по рассеянию примеси от поверхностного источника над городским островом тепла в рассматриваемый критический метеорологический период, нет возможности подвергнуть результаты численной реализации диффузионных моделей прямой количественной проверке путем сопоставления с данными измерений. О степени их достоверности можно судить по косвенным признакам. Во-первых, аналогичная по своей сути модель переноса активной примеси (тепла) дает результаты [1], вполне удовлетворительно согласующиеся с данными прямых инструментальных измерений [11]. Во-вторых, точность численного решения проверена на последовательно измельченных сетках (25 х116 узлов по горизонтали и вертикали, соответственно, и сетке 50 х 232 узла).
Непрерывный источник примеси заданной производительности Q реализован при численном моделировании в двух модификациях: а) - линейный размер источника совпадает с диаметром В городского острова тепла (нагреваемой пластины радиуса В /2) и Ь) - линейный размер источника составляет одну пятую часть радиуса. Начало источника совпадает с осью расчетной цилиндрической области (г / В = 0). На источнике задается постоянный вертикальный поток примеси: -Вт ■ (—С / — г) = Нс, где Вт = V /8с - коэффициент молекулярной диффузии, V - коэффициент кинематической вязкости, 8с - молекулярное число Шмидта, величина потока примеси Нс = Q/ (А * В), Д = 0,5 или 0,1, соответственно. Величина Q задавалась из условия, чтобы число Рейнольдса Яе = Q / V обеспечивало поступление примеси из источника без начального импульса и, таким образом, лимитировалось скоростью натекающего на источник внешнего потока.
В начальный момент времени окружающая среда находится в покое, начальные поля концентрации С и корреляции < с в > равны нулю. На рис. 1 нанесены: (а) - контуры концентрации для источника примеси, протяженность которого совпадает с протяженностью источника тепла (диаметром нагреваемого диска В), (Ь) - контуры функции тока, показывающие формирование в пределах перемешанного слоя (0 < г / < 1) двух интенсивных вихрей, которые и формируют мощное восходящее движение в центральной части острова
тепла. Рис. 1, а показывает, что под влиянием интенсивного восходящего движения примесь, не накапливаясь вблизи поверхности, переносится в перемешанный слой и далее в слой инверсии, рассеивается в горизонтальном направлении, а также проникает за пределы инверсионного слоя (г / «1). На рис. 2
показаны: (а) - контуры концентрации для источника примеси, протяженность которого равна одной пятой протяженности источника тепла, (Ь) - контуры функции тока. В случае источника меньшей протяженности картина распространения примеси качественно похожа на картину, приведенную на рис. 1, а), однако в количественном распределении концентрации по высоте перемешанного слоя имеются различия.
Рис. 1. Контуры концентрации (а) и функции тока (Ь) для дисперсии примеси от непрерывного источника протяженности 0,5 В
N N
1.С
е.:
-15 -* .3 -0.5 0.С 0 5 1,0 15
. <ь>
••- ... -......_. П'";™!
- -- -...._.............._ ^ - - ■
■■ ■■■■■ •,
■ 7;. с. 17.
■ ' ■"'-'- — ' -
■0.5 о.: с
г/О
Рис. 2. Контуры концентрации (а) и функции тока (Ь) для дисперсии примеси от непрерывного источника протяженности 0,1В
Так, например, в сечении г/ = 0,4 (рис. 3, а, Ь) максимальная концентрация для источника меньшей протяженности, примерно, в пять раз меньше, чем для источника, размеры которого совпадают с протяженностью острова тепла (при одинаковой производительности Q от источника большей протяженности эмиссия примеси в ед. времени также, примерно, больше в пять раз). Проникновение примеси за пределы инверсионного слоя зафиксировано в лабораторных экспериментах по распространению пассивной примеси от непрерывного и мгновенного источников в конвективном пограничном слое [11].
N N
0,60,40,20,0
0,12-
..0,034
4 0,040,00-
(а)
2/2,
.....□.....0.4
-■Л.....1,0
\Ь -V- 1 А
-Í3
0,0
0,5
1,0
1,5
2,0
Ib)
i z; z.
-Л- 1,0
Y; -V.....1,4 í
0,0
0,5
1,0
1,5
2,0
r/D
Рис. 3. Радиальные профили безразмерной концентрации C / Cmax для непрерывного источника пассивной примеси:
а) источник протяженности 0,1D ; b) источник протяженности 0,5D
Второй сценарий рассматривался для двумерного случая. Горизонтальная протяженность области интегрирования равнялась 100 км. Разрешение по горизонтали составляло 1 км, по вертикали - 10 м в пределах первых 50 м от подстилающей поверхности, с последующим растяжением сетки в вертикальном направлении вплоть до высоты 1000 м (выше, до 3000 м, шаг сетки постоянен). Топография поверхности плоская, с урбанизированной поверхностью (городом) протяженностью 10 км, расположенной в центре вычислительной области с абсциссой от 45 км до 55 км (на рис. 4 отмечен темной полоской). Метеорологические начальные условия определялись заданием геострофического ветра (V = 3 м/с) в направлении с запада на восток и атмосферной термической стратификацией, равной 3,5 К/км для потенциальной температуры. Данные измерений температуры поверхности Земли аппроксимировались выражением 0ё (х,0,г) = 68т(%г / 43200), где г - текущее время в секундах. Это единственное нестационарное граничное условие задачи, которое моделирует 24-часовой цикл нагревания солнцем земной поверхности. Остров тепла задавался в виде контраста температуры по отношению к температуре поверхности по той же формуле, но с амплитудой, увеличенной на 4 градуса.
На рис. 4 показаны распределения концентрации на нижнем расчетом уровне вблизи урбанизированной поверхности (города) в 13 часов дня и в 7 часов утра. Можно видеть, что к 07 часам утра происходит накопление концентрации трассера (символы ■ на рис. 4) в пределах города.
Рис. 4. Концентрация пассивного трассера C на самом нижнем
расчетном уровне на 07 часов полных первых суток моделирования (■) и на 13 часов второго дня моделирования (•), как функция горизонтального расстояния
Заключение
Сформулирована эйлеровая дисперсионная модель для описания распространения пассивной примеси в городском атмосферном пограничном слое. Полностью явные алгебраические выражения для компонент вектора турбулентного потока концентрации физически корректно учитывают анизотропию турбулентного переноса, и эффекты плавучести на перенос примеси. Параметризация городской шероховатости, учитывающая как механические факторы (лобовое сопротивление зданий городской застройки), так и термические факторы (городской остров тепла), показывает, что концентрация пассивного трассера вблизи поверхности выше в ночные и ранние утренние часы по сравнению с концентрацией трассера в дневные часы. Кроме того, наличие пика концентрации вниз по ветру за пределами города отражает воздействие города на дисперсию трассера за его пределами.
Благодарности
Работа была выполнена в рамках госзадания Института вычислительной математики и математической геофизики СО РАН № 0315-2016-0004 и при частичной финансовой поддержке РФФИ и Правительства Новосибирской области в рамках научного проекта 18-48-540005.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Fernando H., Lee S., Anderson J., Princevac M., Pardyjak E. and Grossman-Clarke S. Urban Fluid Mechanics: Air Circulation and Contaminant Dispersion in Cites // Environmental Fluid Mechanics. - 2001. - Vol. 1. - P. 107-164.
2. Martilli A. Numerical Study of Urban Impact on Boundary Layer Structure: Sensitivity to Wind Speed, Urban Morphology, and Rural Soil Moisture // J. Appl. Meteor. - 2002. - Vol. 41, No. 12. - P. 1247-1266.
3. Vu T., AshieY., Asaeda T. A к - s turbulence closure model for the atmospheric boundary layer including urban canopy // Bound.-Layer Meteor. - 2002. - Vol. 102. - P. 459-490.
4. Курбацкий А. Ф. Численное исследование воздействия поверхностного теплового пятна на структуру атмосферного пограничного слоя // Теплофизика и аэромеханика. -2005. - Т. 12. - № 1. - С. 41-60.
5. Курбацкий А. Ф. Улучшенная вычислительная модель для турбулентного атмосферного пограничного слоя с параметризацией городской шероховатости // Вычислительные технологии. - 2005. - Т. 10. - Часть 2. - С. 53-60.
6. Bornstein R., Johnson D. Urban-Rural Wind Velocity Differences // Atmos. Environ. -1977. - Vol. 11. - P. 597-604.
7. Курбацкий А. Ф., Курбацкая Л. И. Трехпараметрическая модель турбулентности для атмосферного пограничного слоя над урбанизированной поверхностью // Изв. РАН. Физика атмосферы и океана. - 2006. - Т.42. - №4. - С. 476-494.
8. Kurbatskii A. Computational modeling of the turbulent penetrative convection above the urban heat island in stably stratified environment. // J. Appl. Meteor. - 2001. - Vol. 40. - No. 10. -P.1748-1761.
9. Lu J., Arya P. et al. A Laboratory Study of the Urban Heat Island in a Calm and Stably Stratified Environment. Part I and II // J. Appl. Meteor. - 1997. - Vol. 36. - No. 10. -P.1377-1402.
10. Курбацкий А. Ф., Курбацкая Л. И. Численное исследование городского острова тепла: верификация эйлеровых моделей атмосферной диффузии // Оптика атмосферы и океана. - 2004. - Т. 17. - №5-6. - С. 470-477.
11. Snyder W. H., Lawson Jr. R.E., Shipman M.S., Lu J. Fluid modelling of atmospheric in the convective boundary layer // Bound. Layer Meteorol. - 2002. - Vol. 102. - No. 3. - P. 335-366.
© Л. И. Курбацкая, 2019