Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

NUMERICAL MODELING OF THE ATMOSPHERIC BOUNDARY LAYER WITH THE ACCOUNT OF ITS STRATIFICATION

Kuptsov A.I. 1 Akberov R.R. 1 Islamkhuzin D.Y. 1 Gimranov F.M. 1
1 FGBOU SEI HPE «KNITU» «Kazan National Research Technological University»
A computational model for a horizontally uniform atmospheric boundary layer accounting for various options of atmospheric stability is developed on the basis of Monin-Obukhov’s similarity theory. Numerical modeling of the atmospheric boundary layer is carried out using the Fluent package. For a turbulent closure, the standard k–ε turbulence model with modified turbulence closure constants and an additional source term in the equation for k are applied. Cases of neutral, stable and unstable stratification are considered. The computational model allows preserving horizontal uniformity of vertical profiles of velocity, pressure, density, temperature and atmospheric turbulent viscosity which are set as boundary conditions at the domain’s inlet. The obtained model can be applied to modeling of propagation of hazardous gas clouds in the atmosphere over the terrestrial surfaces including surfaces filled with buildings.
CFD modeling
Fluent 6.3
k–ε model
atmospheric boundary layer
atmospheric stability
1. Galeev A. D., Ponikarov S.I., Salin A.A. Modelirovanie posledstviy avariynogo proliva binarnogo rastvora s ispolzovaniem programmy Fluent [Modelling of the binary solution spill consequences using Fluent software]. Matematicheskoe modelirovanie, 2011, Vol. 23, no. 7, pp. 129–144.
2. Gorskiy V.G. Nauchno-metodicheskie aspekty analiza avariynogo riska. V.G. Gorskiy, G.A. Motkin, V.A. Petrunin, G.F. Tereschenko, A.A. Shatalov, T.N. Shvetsova-Shilovskaya. Moscow, Ekonomika i informatika publ., 2002, 260 p.
3. Monin A.S., Obuhov A.M. Osnovnye zakonomernosti turbulentnogo peremeshivaniya v prizemnom sloe atmosfery [Basic laws of turbulent mixing in the surface layer of the atmosphere]. Trudy Geofizicheskogo instituta AN USSR, 1954, no. 24. (151), pp. 163–187.
4. Shvyryaev A.A. Otsenka riska vozdeystviya zagryazneniya atmosfery v issleduemom regione: uchebnoe posobie. A.A. Shvyryaev, V.V. Menshikov. Moscow, MGU publ., 2004, 124 p.
5. Alinot C., Masson C. Aerodynamic simulations of wind turbines operating in atmospheric boundary layer with various thermal stratifications // ASME Conference Proceedings. 2002. Vol. 2002, no. 7476X. рр. 206–215.
6. Alinot C. Masson C. K-ε model for the atmospheric boundary layer under various thermal stratifications // Journal of Solar Energy Engineering. 2005. Vol. 127, no. 4. рр. 438–443.
7. Balogh M., Parente A., Benocci C. RANS simulation of ABL flow over complex terrains applying an enhanced k-ε model and wall function formulation: implementation and comparison for Fluent and OpenFOAM // Journal of Wind Engineering and Industrial Aerodynamics – 2012. Vol. 104, no. . 106. рр. 360–368.
8. Blocken B., Stathopoulos T., Carmeliet J. CFD simulation of the atmospheric boundary layer: wall function problems // Atmospheric Environment. 2007. Vol. 41, no.2. рр. 238–252.
9. Cheng A., Cheng D.T. Heritage and early history of the boundary element method // Engineering Analysis with Boundary Elements. 2005. no. 29. рр. 268–302.
10. Fluent Inc. Fluent 6.3. User’s Guide, Lebanon, 2006.
11. Parente A., Benocci C. On the RANS simulation of neutral ABL flows // In: Proceedings of the Fifth International Symposium on Computational Wind Engineering (Chapel Hill, North Carolina, USA, 23–27 May, 2010). рр. 1–8.
12. Pontiggia M., Derudi M., Busini V., Rota R. // Journal of Hazardous Materialsэ 2009. Vol. 171, no. 1. рр. 739–747.
13. Russell A. Computational fluid dynamics modeling of atmospheric flow applied to wind energy research // A thesis Master of Science in Mechanical Engineering (Boise State University, December 2009).
14. Versteeg H.K, Malalasekera W. An introduction to computational fluid dynamics. The finite volume method. Second edition. Harlow, 2007. рр. 517.

Моделирование пограничного слоя атмосферы при помощи CFD пакетов позволяет достаточно точно прогнозировать распространение облаков опасных газов, возникающих в случае аварий на химических опасных объектах. При этом крайне важным является адекватное описание поведения распространения газа в зависимости от устойчивости атмосферы. Следовательно, обязательным условием является корректное задание граничных метеоусловий (профили скорости, температуры, турбулентные характеристики и т.д.), которые потом не будут искажены и изменены по длине расчетной области в процессе численного моделирования. Еще одной сложной задачей является выбор граничных условий для турбулентных характеристик, от которых непосредственно зависит коэффициент турбулентной диффузии, и, следовательно, реальное распределение газа в любой точке атмосферы.

При моделировании распространения облаков опасных газов в атмосфере для учета атмосферной турбулентности на основе вертикального градиента температуры чаще всего применяется классификация по Паскуиллу [2, 4]. При использовании этой классификации профиль температуры задается простейшим линейным уравнением, в котором отсутствует влияние теплового потока от поверхности земли с воздухом.

В качестве других граничных условий турбулентных характеристик, некоторые исследователи использовали значения масштаба и интенсивности турбулентности, которые считали заранее известными [1]. Профиль скорости и параметры турбулентности, используемые в расчетах, определялись методом последовательных приближений, до получения одинаковых профилей по всей длине расчетной области. Такой подход позволяет проводить оценочные расчеты, но не учитывает варианты устойчивости атмосферы в конкретных условиях.

В целом при использовании CFD пакетов необходимо учитывать, что решаемые уравнения не вполне адекватно описывают физику турбулентности в атмосферном пограничном слое. Так, например, наиболее часто применяемые уравнения k–εpsilon модели для замыкания турбулентности (для нахождения турбулентной вязкости) не содержат членов, учитывающих генерацию турбулентности силами плавучести. Различные варианты устойчивости атмосферы (нейтральная, устойчивая, неустойчивая) должны учитываться с помощью различных граничных условий. При моделировании используют масштаб длины Монина ‒ Обухова, абсолютное значение которого определяет высоту, до которой влияние сил трения на турбулентность превышает влияние сил плавучести. При положительном значении этой длины атмосфера считается устойчивой, другими словами, данная ситуация подавляет турбулентность, что должно учитываться в модели. Соответственно, при отрицательном масштабе длины Монина ‒ Обухова атмосфера неустойчива, и данная ситуация должна усиливать турбулентность. Нейтральная атмосфера – это адиабатическая атмосфера без притока тепла извне, то есть с нулевым потоком тепла от поверхности земли. При нейтральной атмосфере масштаб длины Монина – Обухова равен бесконечности, что означает, что влиянием сил плавучести на турбулентность можно пренебречь. Моделирование всех этих трех сценариев в CFD пакетах представляет большую трудность, особенно, когда ставится цель – создание горизонтально-однородного пограничного слоя атмосферы, при котором вертикальные профили всех параметров потока должны сохраняться на протяжении длины порядка нескольких километров, на которой практически происходит рассеивание облаков опасных газов до безопасных для человека концентраций. В литературе предлагались различные подходы для численного моделирования пограничного слоя атмосферы, включая изменение констант турбулентности, использование различных граничных условий для верхней и выходной границы, добавление источниковых членов и т.д. [5–7, 11–13]. Большая часть этих работ была посвящена вопросам моделирования нейтральной и устойчивой атмосферы. Предлагались уравнения как для сжимаемой, так и для несжимаемой среды. При этом важную роль играет выбор того или иного CFD пакета, поскольку решение уравнений, граничные условия, параметры турбулентности (высота шероховатости и т.д.) реализованы в различных CFD пакетах по-разному. Наиболее часто в работах используются пакеты Fluent, CFX и Star CD (или разновидность Star CCM+).

В настоящей работе разработан подход с использованием уравнений для сжимаемой и несжимаемой среды через решение уравнения состояния идеального газа с использованием пакета Fluent 6.3. В модели на верхней границе используется условие Дирихле [9] с известными расчетными значениями всех параметров потока, на входной границе задаются вертикальные профили через применение пользовательских функций (UDF), на нижней границе задается значение температуры для уравнения энергии, а для компонент скорости и турбулентных характеристик используются стандартные пристеночные функции, полученные из логарифмического закона стенки для скорости. На выходной границе используется условие Pressure Outlet [10]. При этом особое внимание уделяется профилю давления, т.к. задание однородного профиля давления на выходе приводит к проблемам со сходимостью и к искажению действительного поля течения под действием силы тяжести.

Физико-математическая модель

Система дифференциальных уравнений, определяющих трехмерный нестационарный поток газовоздушной смеси, перенос энергии, газа, турбулентной кинетической энергии и скорости диссипации турбулентной кинетической энергии, имеет следующий вид [14]:

Уравнение неразрывности:

kuptsov001.wmf (1)

где ρ – плотность воздуха или газовоздушной смеси, кг/м3; xi – координаты x, y, z в метрах (x – горизонтальная координата в направлении ветра; y – горизонтальная координата в направлении поперек ветра; z – вертикальная координата); ui – компоненты осредненной скорости ветра u, v, w, м/с.

Уравнение для переноса импульса:

kuptsov002.wmf (2)

где p – давление, Па; kuptsov003.wmf – пульсационные скорости, м/с.

Уравнение для переноса энергии:

kuptsov004.wmf (3)

а энтальпия, в единицах кДж/кг, равна:

h = cpT, (4)

где cp – удельная теплоемкость при постоянном давлении, кДж/(кг·К); T – температура воздуха, К; λ – коэффициент теплопроводности, Вт/(м·К); Prt – турбулентное число Прандтля; Prt = 0,85.

Уравнение для переноса газа:

kuptsov005.wmf (5)

где Ys – массовая доля компонента s газовоздушной смеси; Sct – турбулентное число Шмидта; D – коэффициент молекулярной диффузии, зависящий от состава смеси, м2/с; Sc – источниковый член, определяющий генерацию газовой примеси, кг/(м3·с). В случае включения этого уравнения в модель через соответствующую опцию в Fluent («Species Transport on»), во всех уравнениях переноса значения коэффициентов переноса и всех переменных потока будут значениями для смеси, определяемые используемыми по умолчанию правилами для смешения.

Уравнение для переноса k (турбулентной кинетической энергии):

kuptsov006.wmf (6)

где Sk – источниковый член.

Генерация (подавление) турбулентности, обусловленная температурной (плотностной) стратификацией, определяется из выражения [10]:

kuptsov007.wmf (7)

где gi – компонент вектора гравитации в i-м направлении, м/с2; β – коэффициент термического расширения, 1/К.

Для идеальных газов уравнение выше приводится к виду

kuptsov008.wmf (8)

Уравнение для переноса ε (скорости диссипации турбулентной кинетической энергии):

kuptsov009.wmf (9)

где Sε – источниковый член, кг/(м·с4), который здесь принимается равным нулю; тензор Eij определяется стандартными выражениями k–ε модели.

В уравнении (9) имеется связанный с силами плавучести коэффициент С3ε, рассчитанный из следующего соотношения:

kuptsov010.wmf (10)

где v и u – компоненты вектора скорости газовоздушного потока, м/с, параллельный и перпендикулярный вектору силы тяжести соответственно.

Уравнение для турбулентной вязкости, кг/(м·с):

kuptsov011.wmf (11)

Уравнение для коэффициента турбулентной диффузии:

kuptsov012.wmf (12)

Уравнение состояния идеального газа:

ρ = (pраб + p)/(R·T/Mω), (13)

где pраб – рабочее давление, Па; p – избыточное давление, рассчитанное пакетом Fluent, Па; R – универсальная газовая постоянная, Дж/(моль·К); T – температура газа, К; Mω – молекулярный вес газа, г/моль.

Константы турбулентности согласно работам [5–6] допускается принимать

Cμ = 0,0333; σk = 1,00; σε = 1,3;

C = 1,176; C = 1,92.

Данные константы турбулентности в случае моделирования пограничного слоя атмосферы отличаются от констант турбулентности, используемых в Fluent по умолчанию. Результаты расчетов показали, что коэффициенты турбулентной диффузии, полученные при использовании k-ε модели с указанными константами хорошо согласуются с коэффициентами турбулентной диффузии, полученными по данным Монина и Обухова. При этом если принять турбулентное число Шмидта за единицу, то коэффициент турбулентной диффузии, полученный в результате численного моделирования, будет полностью совпадать с коэффициентом Kz, найденным по теории подобия Монина – Обухова [3]:

Kz = (K·u*·z)/(1 + β*·z/L), (14)

где K – константа Кармана (в расчетах использовалось значение 0,41); u* – динамическая скорость, м/с; z – вертикальная координата, м; β* – некоторая универсальная постоянная, подлежащая определению на основе эмпирических данных; L – масштаб длины Монина – Обухова, м.

Уравнения для измененной «стандартной» k–ε модели для несжимаемой среды остаются теми же, за исключением уравнения состояния идеального газа (13). В этом случае оно выглядит уже как

ρ = (pраб)/(R·T/Mω). (15)

В данной работе не рассматривалась модель Буссинеска для плотности, из-за того, что у этой модели есть существенный недостаток: ее нельзя применять в случае больших перепадов температуры по высоте, так как в этом случае линейная аппроксимация не будет удовлетворять приемлемой точности [10].

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

Моделирование же несжимаемой среды возможно как с фиксацией профилей турбулентных характеристик, так и без фиксации. При этом сохранение профиля атмосферной турбулентной вязкости во всей расчетной области происходит за счет использования в уравнении для переноса кинетической энергии турбулентности k (6) источникового члена [5–6]:

kuptsov014.wmf (16)

Таким образом, применение модели несжимаемой среды без фиксации профилей позволяет учесть изменения турбулентности в случае препятствий. Главное отличие от работ [5–6] состоит в том, что источниковый член (16) используется в случае моделирования несжимаемой среды и при этом источниковый член в уравнении для переноса скорости диссипации кинетической энергии турбулентности ε (9) равен нулю. Благодаря этому затраты машинного времени на выполнение расчетов существенно сокращаются.

Особое внимание при моделировании следует обратить на создание поля давления под действием гравитационных сил земли.

Граничные условия

Моделирование горизонтально-однородного пограничного слоя атмосферы с учетом его стратификации производится до начала выбросов. При этом рассматривался стационарный процесс. Следовательно, в вышеприведенных уравнениях (1–3, 5–6, 9) производные по времени равны нулю; однако в Fluent они используются для получения стационарного решения методом установления по времени.

Для входной (левой) и верхней границ расчетной области использовалось граничное условие «Velocity Inlet», требующее задания значений или профилей всех компонент скорости, температуры, k и ε в соответствии с моделируемым состоянием атмосферной стратификации [3, 5–7, 11–13]:

kuptsov015.wmf kuptsov016.wmf kuptsov017.wmf C = 0;

kuptsov018.wmf kuptsov019.wmf

Вход и верхняя граница расчетной области для различных состояний атмосферной стратификации задавались как граничное условие – Velocity Inlet (скорость на входе [3, 5–7, 11–13]):

1. Моделирование неустойчивого состояния атмосферы, L < 0:

kuptsov020.wmf (17)

kuptsov021.wmf (18)

kuptsov022.wmf (19)

kuptsov023.wmf (20)

kuptsov024.wmf (21)

kuptsov025.wmf (22)

где u(z) – скорость ветра на высоте z, м/с; z0 – высота шероховатости, м; kuptsov026.wmf kuptsov027.wmf – универсальные безразмерные функции; T(z) – температура на высоте z, К; T0 – температура поверхности (земли), К; T* – масштаб температуры, К; g – ускорение свободного падения, м/с2; ε – скорость диссипации турбулентной кинетической энергии, м2/с3; k – турбулентная кинетическая энергия, м2/с2;

2. Моделирование нейтрального состояния атмосферы, L = ∞:

kuptsov028.wmf (23)

kuptsov029.wmf (24)

kuptsov030.wmf (25)

kuptsov031.wmf (26)

3. Моделирование устойчивого состояния атмосферы, L > 0:

kuptsov032.wmf (27)

kuptsov033.wmf (28)

kuptsov034.wmf (29)

kuptsov035.wmf (30)

kuptsov036.wmf (31)

kuptsov037.wmf (32)

При включении уравнения для переноса газа на входной границе необходимо задание массовых долей всех компонентов газо-воздушной смеси.

На нижней границе области задавалось условие прилипания («wall»), подразумевающее, что на этой границе u = v = w = 0. Однако при использовании этой опции в сочетании со стандартной k–ε моделью граничное условие задается не на самой стенке, а в центре первой пристеночной ячейки расчетной сетки через использование логарифмического закона стенки для скорости, температуры и концентрации. Для k и ε задание граничных условий на нижней границе не требуется (для стандартной k–ε модели Fluent использует для k условие ∂k/∂z = 0, а значение ε в центре первой пристеночной ячейки расчетной сетки рассчитывается по стандартной формуле).

Эквивалентная шероховатость:

kuptsov038.wmf (33)

где Ks – эквивалентная высота элементов шероховатости на подстилающей поверхности, используемая в Fluent вместо z0 для логарифмического закона стенки, м; Cs – константа шероховатости, Cs = 0,5 [8].

Значение температуры на нижней границе области задавалось постоянным, в соответствии с вертикальным профилем температуры, выбранным для рассматриваемого условия атмосферной стратификации.

На выходной (правой) границе области использовалось граничное условие «Pressure Outlet», требующее задания вертикального профиля давления. При этом градиенты всех остальных переменных потока в нормальном направлении к выходу принимаются равными нулю. В случае возникновения втекания воздуха через эту границу на некоторых итерациях при решении уравнений Fluent рассматривает те участки границы, через которые произошло втекание воздуха, в качестве входной границы. Следовательно, пользователь должен заранее предусмотреть эту ситуацию и задать корректные вертикальные профили для температуры, k и ε на выходной границе для ускорения сходимости решения. Для этой цели в данной работе выходные профили этих переменных потока задавались по тем же зависимостям, что и вертикальные профили на входе.

Отличие от предыдущих работ [5–7, 11–13] заключается в том, что в этих работах на выходной границе использовалось условие либо «Outflow», которое не требует задания вертикального профиля для давления, либо «Pressure Outlet» с заданием однородного вертикального профиля по высоте или с использованием профиля на основе уравнения гидростатики. В данной работе используется условие «Pressure Outlet», но давление задается неоднородным по высоте, когда профиль давления на выходной границе устанавливается неоднородным, через использование опции «Boundary Profile» (граничный профиль). Этот профиль получается итерационным путем. Сначала задается однородный профиль, и расчет ведется до сходимости. Затем записывается профиль из середины области, где профиль давления более-менее установился, и используется в качестве граничного условия на выходе. Затем снова записывается граничное условие из середины области и используется в качестве граничного условия на выходе. Данный итерационный процесс повторяется до тех пор, пока разница между итерациями будет неощутимой.

Расчетная сетка

Описанная выше физико-математическая модель описывает трехмерный процесс. Однако для сокращения времени вычислений в данной работе использовано двухмерное приближение к реальной задаче. Для этого достаточно использования двухмерных конечных элементов для дискретизации расчетной области задачи.

Таким образом, двухмерная расчетная сетка была построена из элементов четырехугольной формы и сгущалась в сторону всех четырех границ области. Размеры области двухмерной расчетной сетки менялись в зависимости от моделирования определённой устойчивости атмосферы в диапазоне от 3000 до 5000 м по длине. Особое внимание было уделено размеру первого ряда ячеек сетки вблизи нижней границы области, т.к. для этого ряда ячеек должно выполняться сразу несколько условий:

1) 30 < y+ < 500;

2) yp > y0.

В расчетах y принималась как координата высоты расчетной области.

Первое условие означает, что центр первой ячейки сетки должен находиться в развитой турбулентной (логарифмической) части турбулентного пограничного слоя атмосферы, имеющего масштаб, приблизительно равный высоте элементов шероховатости на подстилающей поверхности. Второе условие означает, что центр первой ячейки сетки yp должен быть выше элементов шероховатости y0 на подстилающей поверхности [8].

Результаты расчетов

Результаты расчетов приведены на рисунке.

На рисунке представлены результаты численного моделирования в двухмерной расчетной области. В качестве сохраняемых параметров вдоль расчетной области были выбраны: атмосферная турбулентная вязкость, скорость и температура. Вертикальные профили на входной границе x = 0 м строились с учетом теории подобия Монина – Обухова и сравнивались с профилями в сечении x = 500 м при различных стратификациях атмосферы в случаях сжимаемой и несжимаемой сред с учетом фиксации или с учетом использования источникового члена.

pic_35.tif

Численное моделирование атмосферного пограничного слоя: 1 – Граничное условие на входе; 2 – Профиль через 500 м, сжимаемая среда, фиксация; 3 – Профиль через 500 м, сжимаемая среда, фиксация, стандартные константы; 4 – Профиль через 500 м, несжимаемая среда, фиксация; 5 – Профиль через 500 м, несжимаемая среда с источником

Согласно результатам численного моделирования лучше всего профили атмосферной турбулентной вязкости сохраняются в случае сжимаемой среды и при фиксации турбулентных характеристик. При этом профили в расчетах с использованием несжимаемой среды отличаются не намного: при использовании фиксации в среднем на 2,7 % от граничных условий на входе, а в случае использования источникового члена на 7,1 %. Лучшее сохранение профилей турбулентной атмосферной вязкости в случае фиксации турбулентных характеристик объясняется малым изменением плотности атмосферного воздуха в ходе моделирования по длине расчетной области, причем в случае сжимаемой среды речь идет о практически полном сохранении входного профиля. При использовании источникового члена профили отличаются не только из-за изменения плотности и влияния, оказываемого самим членом, но еще из-за заданного граничного условия на верхней границе. Для того чтобы снизить влияние граничного условия на верхней границе, необходимо увеличивать расчетную сетку по высоте. В таком случае возможно приемлемое моделирование до такой определенной высоты, когда профиль атмосферной турбулентной вязкости будет не сильно отличаться от входного профиля. Наглядно видно, что использование стандартных констант турбулентности существенно завышает значение атмосферной турбулентной вязкости по сравнению с граничными условиями, полученными по теории подобия Монина – Обухова. Это объясняется тем, что в уравнении (11) для атмосферной турбулентной вязкости стандартная константа турбулентности Cμ принимается равной 0,09.

Профили скорости и профили температуры при разных устойчивостях от граничных заданных условий в среднем отличаются не более чем на 0,1 %.

Тем самым расчеты показали что, несмотря на то, что плотность воздуха в случаях сжимаемой и несжимаемой среды задается различными уравнениями, они не сильно влияют на соответствующие профили скорости и температуры.

Исходя из результатов расчетов, было замечено, что с дальнейшим увеличением продольной координаты x профили начинают отличаться от входного профиля. Для сохранения профилей на еще более дальние расстояния необходимо увеличивать расчетную область в этом направлении. Увеличение сетки в направлении x положительно влияет на сохранение профилей в определенной точке, но увеличивает время расчета.

Такие малые расхождения в профилях атмосферной турбулентной вязкости, скорости и температуры позволяют адекватно учесть распространение опасных газов с учетом стратификации атмосферы во всей расчетной области.

Заключение

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

Рецензенты:

Зиннатуллин Н.Х., д.т.н., профессор кафедры ПАХТ Казанского национального исследовательского технологического университета, г. Казань;

Строганов В.Ф., д.х.н., профессор, зав. кафедрой ХиИЭС Казанского государственного архитектурно-строительного университета, г. Казань.

Работа поступила в редакцию 05.08.2014.