Исследование движения воздуха и переноса примеси в уличном каньоне с использованием вихреразрешающей модели турбулентного течения

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2012 Математика и механика № 4

УДК 519.6

Е.А. Данилкин, Р.Б. Нутерман, А.А. Барт, Д.В. Деги, А.В. Старченко

ИССЛЕДОВАНИЕ ДВИЖЕНИЯ ВОЗДУХА И ПЕРЕНОСА ПРИМЕСИ В УЛИЧНОМ КАНЬОНЕ С ИСПОЛЬЗОВАНИЕМ ВИХРЕРАЗРЕШАЮЩЕЙ МОДЕЛИ ТУРБУЛЕНТНОГО ТЕЧЕНИЯ1

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

Ключевые слова: вихреразрешающее моделирование, турбулентность, уличные каньоны, параллельные вычисления.

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

На современном этапе развития теории турбулентности моделирование турбулентных течений в окружающей среде осуществляется, в основном, с использованием осредненных по Рейнольдсу уравнений Навье - Стокса и переноса (Reynolds Averaged Navier - Stokes - RANS-подход), для которых требуется решить проблему замыкания путем привлечения полуэмпирических моделей различного уровня сложности. Однако, несмотря на значительные успехи в развитии этих моделей, имеются определенные трудности при описании нестационарных турбулентных течений вблизи плохообтекаемых тел. В первую очередь это обусловлено определенными особенностями отрывных течений, а именно, наличием организованных когерентных структур, определяемых параметрами течения и геометрией области исследования [2, 11]. Метод крупных вихрей (Large Eddy Simulation - LES-подход) оказывается более предпочтительным при моделировании турбулентных

1 Работа выполнена при финансовой поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 (госконтракт № 14.B37.21.0667) и РФФИ в рамках научных проектов № 12-05-31341, № 12-01-00433а.

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

Необходимо отметить, что высокая вычислительная стоимость проведения расчетов на ЭВМ при реализации вихреразрешающего моделирования для описания турбулентности существенно ограничивает широкое применение ЬБ8-подхода для исследования течений в окружающей среде и является следствием, во-первых, сложной природы турбулентности, поскольку для явного разрешения энергосодержащих вихрей, характеризующих поведение потока, требуется использование очень подробной разностной сетки (порядка 106 - 108 ячеек в объеме, не превышающем размер наибольшего вихря в исследуемой системе). А во-вторых, в отличие от ЯЛ^-подхода, где решение исходной системы есть средние характеристики исследуемого течения, при использовании вихреразрешающего моделирования результатом решения является состояние течения в некоторый конкретный момент времени, и для получения средних характеристик потока необходимо провести осреднение по достаточно большому количеству временных шагов (величина порядка 104 - 106).

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

Физико-математическая постановка задачи

Рассматривается трехмерное нестационарное турбулентное движение несжимаемой среды с распределенными источниками примеси постоянной интенсивности.

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

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

f (x, t) = I G(x, x')f (x', t)dx', I G(x)dx = 1,

R3

R3

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

Полученная математическая модель трехмерного нестационарного турбулентного движения несжимаемой среды включает в себя отфильтрованные уравнения неразрывности и Навье - Стокса:

а также отфильтрованное дифференциальное уравнение, описывающее перенос пассивной газообразной примеси:

Здесь ui - фильтрованные компоненты мгновенного поля скорости; p - мгновенное значение давления; v - коэффициент кинематической вязкости; р -

плотность; C - концентрация примеси; F - функция, описывающая распределение источников; - тензор подсеточных напряжений, qj - подсеточный поток

массы; Sc - число Шмидта. По повторяющемуся индексу j проводится суммирование.

Система уравнений (2) - (6) остается незамкнутой, так как выражения (4) и (6) содержат исходные нефильтрованные компоненты скорости и концентрации.

и qj характеризуют влияние мелкомасштабных вихрей на эволюцию крупномасштабных вихрей, и их необходимо моделировать, как в подходе Буссинеска через установление их связи со скоростями щ и концентрацией C. Последнее и составляет суть методики подсеточного моделирования [1, 6, 7, 11, 20].

В настоящее время существует большое количество подходов к моделированию подсеточных масштабов [1, 20], но наиболее широкое распространение получили модели, основанные на использовании турбулентной вязкости (Eddy Viscosity Models, EVM) для не разрешенных явно турбулентных пульсаций. В таких моделях тензор напряжений и турбулентный поток массы вычисляются по формулам

' діїі Л дт

dt 3x}- р dxi dx}- V dx}- j

uiu l -uiu l, i, j = 1,2,3 ,

1 J 1 J 7 7 J 7 7 7

dC dUjC д f v дС Л dqj

dt dxj dxj Sc dx j dxj

j j V j j j

qj = Cu  - CU:, j = 1,2,3 .

1 j j j J 77

+ F;

1 fe+duLN

2 dxj dxi

V j * j

Здесь - тензор скорости деформации, построенный по фильтрованному полю скорости щ , Кт = Кт (м, х, t) > 0 - коэффициент турбулентной вязкости, зависящий от решения, Бсх - турбулентное число Шмидта.

Выбор зависимости Кт = Кт (м, х, t) чрезвычайно разнообразен. Например, для модели Смагоринского и ее модификаций турбулентная вязкость определяется на основе локальных параметров потока, а для динамических и дифференциальных моделей (и их комбинаций) - из предыстории развития течения. Более подробную информацию по моделям турбулентной вязкости можно найти в работах [1, 3, 12, 20].

В данной работе использовалась модель Смагоринского [21], в которой влияние мелкомасштабных вихрей на эволюцию крупномасштабных аппроксимируется выражением тЦш°8 = -2С<?Д1|8|8у, где Д^ - шаг сетки модели, |8| = ф8~8~ -норма тензора скорости деформации, С3 - постоянная Смагоринского.

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

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

Пренебрегая отклонением профиля скорости от логарифмического закона, скорость в первом расчетном узле над шероховатой поверхностью можно определить, используя следующую зависимость:

- параметр шероховатости, к и 0,4 - постоянная Кармана, г - расстояние до поверхности. Предположение о логарифмическом распределении профиля скорости использовалось при оценке влияния напряжений трения у стенки.

Аппроксимация дифференциальной задачи и численный метод решения

Численное решение представленной выше системы дифференциальных уравнений в частных производных осуществляется на основе метода конечного объема с использованием разнесенной разностной сетки, когда значения компонент скорости определяются на гранях конечных объёмов, а скалярные характеристики

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

где мт = / р - динамическая скорость, %М1 - напряжения трения на стенке, г0

При вычислении интегралов применяется кусочно-полиномиальная интерполяция для зависимых от х1,x2,х3 величин. Аппроксимация конвективных членов уравнения переноса выполняется с использованием направленной схемы QUICK [16]. Аппроксимация диффузионных членов осуществляется с использованием центрально-разностной схемы второго порядка. Для решения уравнений переноса применяется явная схема по времени (Адамса - Бэшфорда). Результатом дискретизации является явная разностная схема второго порядка аппроксимации по времени и пространству, являющаяся условно устойчивой.

В гидродинамической части модели для согласования полей скорости и давления использовалась схема предиктор-корректор, в соответствии с которой явная схема Адамса - Бэшфорда для уравнений движения выполняла функцию предиктора, а коррекция поля скорости, удовлетворяющего уравнению неразрывности на новом временном слое, выполнялась на основе решения разностного уравнения для давления. Для решения системы линейных алгебраических уравнений для отыскания давления использовался метод сопряженных градиентов (CG) с использованием предобуславливания методом верхней релаксации с красно-черным упорядочиванием [5].

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

Параллельная реализация

В качестве основного подхода распараллеливания выбрана двумерная геометрическая декомпозиция сеточной области, которая подразумевает выделение каждому процессорному элементу некоторой сеточной подобласти со всеми принадлежащими этой подобласти значениями неизвестных сеточных функции (рис. 1).

Ахз

х2№>

Х1

Пэ Т-

(0, Ц-

1.’

...;..и.п.Э . 3...L.

-1-1ПЭ-4

"МтШ"

Рис. 1. Иллюстрация организации пересылок граничных сеточных значений при 2Б-декомпозиции сеточной области для обеспечения однородности вычислительного процесса

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

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

Ыу+1

Ыу

- внутренние узлы;

- граничные узлы, необходимые для расчетов соседнему процессорному элементу;

1,2

ЫхЛх+1

Рис. 2. Иллюстрация принадлежности узлов вычислительной сетки для одной подобласти

Использование многопроцессорной вычислительной техники (кластер ТГУ СКИФ СуЪепа, 358 узлов/696 шестиядерных процессоров 1Пе1Хеоп 5670, 2,93ГГ ц) позволило существенно сократить время вычислений при решении описанного класса задач. В таблице представлено время работы (в часах) параллельной программы для различного количества процессорных элементов, используемых при решении уравнений Навье - Стокса для течения над шероховатой пластиной (5 000 шагов по времени). Заметим, что при проведении вычислений для турбулентных течений требуется около 25 000-250 000 шагов по времени.

Время счета для случая моделирования течения над шероховатой пластиной

на сетке 120x120x30

Число процессов 1 4 16 25 64 100

Время счета, ч 38,68 23,34 2,86 1,12 0,83 0,81

Результаты математического моделирования

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

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

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

- условия прилипания. Для задания начальных условий выбраны аналитические выражения, предложенные в работе [18] для трех компонент скорости (ReT= uTh/v = 180,h - полуширина канала):

u1(x1,x2,x3)/Uin = C(— (x3 /h)8 ) + e2nsin(nx3 /h)cos( /h)sin(x2 /h);

U2

Uj(x1,x2,x3)/Uin = —e(1 + cos(nx3 /h))sin(x1 /h)sin(x2 /h);

где C = 7,764ReT 7 • uT / Uín и e = 0,1C. Эти начальные условия удовлетворяют граничным условиям и по своей сути являются возмущением, наложенным на среднее течение для задания начального уровня турбулентности в потоке.

Условия проведения расчетов, геометрия области исследования и параметры течения выбраны в соответствии с работой [13]. Размеры области исследования составляли 4nh, 2nh, 2h , расчеты проводились на равномерной сетке 96 х 80 х 64 . Результаты моделирования течения жидкости показали хорошее согласование с экспериментальными данными [19] и данными прямого численного моделирования [15] (рис. 3, 4).

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

Рис. 3. Профиль осредненной по времени продольной компоненты скорости

3 -

о

л

§

о

5

5

С

х

л

и

2 -

.V

я

і

3

а'

і

_ь

>

•О '

о и'- БШ О и'- Ехрі

- дитек

>

>.о '

«••г.

о 4-

о

-1—

0,2

0,4

—і—

0,6

—і—

0,8

Х3

Рис. 4. Флуктуации продольной компоненты скорости

Эксперимент с обтеканием цилиндра [17] также использовался для апробации построенной математической модели турбулентного течения. Моделируется турбулентный поток вокруг цилиндра квадратного сечения, расположенного в канале (рис. 5). Рассматриваемый случай имеет довольно простую геометрию, но представляет большой интерес с точки зрения моделирования процесса образования

Рис. 5. Расчетная область

1

крупных вихреи с их последующим распадом и диссипации в возникающем неустойчивом следе за цилиндром. Размеры геометрии области исследования представлены на рис. 5, в расчетах использовалась область Nxx х Nx2 х Nx3 = 132 х 122 х 22 . Выбор области течения и значения числа Рейнольдса Re = UinD / v соответствуют эксперименту Lyn et al. [17].

Движение жидкости поддерживалось за счет задания постоянной скорости на входной границе, направленной вдоль оси Ox1. В качестве начальных данных использовались постоянные значения скорости U1 = Uin = 0,55 м/с , U2 = U3 = 0 . Расчеты (и осреднение величин) проводились в течение двенадцати периодов отрыва вихря после времени инициализации, равного 50D / Uin .

Эксперимент с обтеканием цилиндра квадратного сечения хорошо описан, имеются данные измерений, представленные в виде таблиц. Этот эксперимент часто используется для апробации различных подходов моделирования турбулентности. Набегающий ламинарный поток рассеивается на наветренной стороне цилиндра, а с боковых поверхностей отрываются вихри. Экспериментально установленная частота отрыва вихрей f соответствует значению безразмерного числа Струхаля Sh = fD / Uin = 0,13 . Такое же значение частоты отрыва вихрей получено численно в результате расчетов.

Вычислительный эксперимент показал наличие хорошего уровня согласования расчетов с экспериментальными данными [17]. Профиль компоненты скорости (Uj), осредненной по времени, хорошо накладывается на данные измерений. На рис. 6 представлены результаты расчетов изменения осредненной по времени продольной скорости потока вдоль оси, ориентированной по потоку и проходящей через середину поперечной плоскости расчетной области.

<“ і)

U2

0,4.

0 -

-0,4

-4

h

р

*

Expr.[17]  QUICK

x1/D

(“22> ]

1 - - -

Л* « *

0,6 - » ' п п *

1° v  1 ч

0,4 -  /

£ '

і

0,2 - і і 1

0 1 1 і " * " Г " 1 1 1

=Expr.[17] • QUICK

0

x1/D

Рис. 6. Сравнение профиля средней скорости иь полученного на экспериментальном стенде, с результатами моделирования, полученными в данной работе

Анализ результатов показывает, что использование направленной схемы РШК позволяет достичь довольно точного согласования результатов численного эксперимента с измерениями [17] для средних характеристик потока при использовании замыкания Смагоринского с постоянным коэффициентом. На рис. 6 также

представлен профиль осредненных по времени нормальных турбулентных напряжений вниз по потоку за цилиндром. Из рисунков видно, что получено достаточно хорошее соответствие с данными эксперимента [17] для относительных пульсаций компоненты скорости (й2 ) .

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

200

150

100

50

0

-50

-100

-150

-200

-250

-300

-350

-400

-450

-500

-550

-600

¡200 І 150 І 100 50 0

-50

-100

-150

-200

-250

-300

-350

-400

-450

-500

-550

1-600

200

150

100

50

0

-50

-100

-150

-200

-250

-300

-350

-400

-450

-500

-550

-600

08 -0,06 -0,04 -0,02 0 0,02 0,04 0,06 0,08 0,1 0,12 0,14

Рис. 7. Векторное поле скорости и карта давления в моменты времени 11 < ¿2 < Хз, где Х2 = ¿1 +0,2Т, Хз = ¿1 +0,4Т, Т - период отрыва вихря

На основе построенной математической модели турбулентного течения несжимаемой среды проведен ряд расчетов для трехмерной модели уличного каньона. Геометрия исследуемой области представлена на рис. 8. В данной работе исследуется характер распространения примеси в зависимости от соотношения ширины и высоты уличного каньона, расположения источника примеси и скорости потока.

Рис. 8. Расчетная область

Задача решалась в следующей постановке. В поперечном основному потоку направлении использовались условия периодичности. В этом случае необходимо располагать боковые границы расчетной области так, чтобы длина области в поперечном направлении была больше, чем размер самого большого вихря в каньоне. Для рассматриваемого случая крупнейший вихрь ограничен шириной каньона Ж = 20 м, тогда поперечную длину каньона можно взять равной 30 м (Ь / Ж = 1,5) [14]. Периодические граничные условия задавались в продольном направлении с целью имитировать бесконечную серию каньонов. Расчеты проводились на сетке 182 х 54 х180. Источник поступления примеси постоянной интенсивности располагался вблизи поверхности на высоте Иг = 0,125 м над центром основания рассматриваемой расчетной области.

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

На основе численных экспериментов показано, что увеличение скорости основного потока способствует более интенсивному выносу примеси из уличного каньона. Также показано, как геометрические параметры уличного каньона влияют на вид течения и уровень загрязнения. Так, в случае уменьшения высоты Н уличного каньона или увеличения ширины Ж до соотношения Н / Ж = 0,5 центр основного вихря смещается к наветренному заданию и вихрь растягивается во всю длину канала (рис. 9). Это приводит к уменьшению скорости вращательного движения воздушных масс у подветренной стороны каньона и, как следствие,

Рис. 9. Изолинии концентрации примеси, осредненной вдоль оси Ох2 , векторное поле скорости, и1п = 5 м/с ; источник примеси

z/H

0,8

0,4

C/Cm

Рис. 10. Нормированная и осреднен-ная концентрация примеси на наветренной стороне уличного каньона. Точками обозначены экспериментальные данные [14]

примесь менее интенсивно выносится из каньона и возрастают локальные значения концентрации примеси. Дальнейшее увеличение расстояния между зданиями Н / Ж = 0,125 приводит к образованию двух рециркуляционных зон: большой вихрь у подветренной стороны и малый - у наветренной. В этом случае примесь, поступающая от источника, расположенного в центре уличного каньона, уносится в сторону подветренного здания, где она циркулирует в турбулентном вихре.

Помимо проведения параметрических расчетов для случая Ж = 20 м, Н = 24 м, ЦУ1П = = 5 м/с (источник примеси располагается в центре каньона в точке х = 30 м, у = 15 м, г = 0,125 м) выполнено сравнение с экспериментальными данными из работы [14]. Резуль-

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

Заключение

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

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

ЛИТЕРАТУРА

1. Волков К.Н., Емельянов В.Н. Моделирование крупных вихрей в расчетах турбулентных течений. М.: Физматлит, 2008. 368 с.

2. Гарбарук А.В., Стрелец М.Х., Шур М.Л. Моделирование турбулентности в расчетах сложных течений: учебное пособие. СПб: Изд-во Политехи. ун-та, 2012. 88 с.

3. Глазунов А.В. Вихревое моделирование турбулентности с использованием смешенного динамического локального замыкания // Известия РАН. Физика атмосферы и океана. 2009. Т. 45. № 1. С. 7-28.

4. Данилкин Е.А., Старченко А.В. К выбору способа декомпозиции при численном решении систем связанных дифференциальных уравнений на многопроцессорной технике с распределенной памятью // Третья Сибирская школа-семинар по параллельным вычислениям. Томск : Изд-во Том. ун-та, 2006. С. 95-101.

5. Данилкин Е.А., Старченко А.В. Параллельная реализация численного метода решения системы уравнений Навье - Стокса при моделировании крупных вихрей турбулентных течений // Вестник Новосибирского государственного университета. Сер. Информационные технологии. 2009. Т. 7. № 2. С. 49-61.

6. Колмогоров А.Н. Локальная структура турбулентности в несжимаемой вязкой жидкости при очень больших числах Рейнольдса // ДАН СССР. 1941. Т. 30. С. 299-303.

7. КурбацкийА.Ф. Лекции по турбулентности: в 2 ч. Введение в турбулентность. Новосибирск, 2000. 118 с.; Моделирование турбулентных течений. Новосибирск, 2001. 136 с.

8. Нутерман Р.Б., Бакланов А.А., Старченко А.В. Моделирование аэродинамики и распространения выбросов от автотранспорта в городском подслое // Математическое моделирование. 2010. Т. 22. № 4. С. 3-22.

9. Ортега Дж. Введение в параллельные и векторные методы решения линейных систем. М.: Мир, 1991. 364 с.

10. Турбулентные сдвиговые течения: пер. с англ. / ред. А.С. Гиневский. М.: Машиностроение, 1982. 432 с.

11. Хлопков Ю.И. Лекции по теоретическим методам исследования турбулентности. М.: МФТИ, 2005. 178 с.

12. Germano M., Piomelli U., Moin P., Cabot W.H. A dynamic subgrid-scale eddy viscosity model // Phys. Fluids. A. 1991. V. 3. P. 1760-1765.

13. Gokarn A., Battaglia1F., Fox R.O. Large eddy simulations of incompressible turbulent flows using parallel computing techniques // Int. J. Numer. Meth. Fluids. 2008. V. 56. No. 10. P. 1819-1843.

14. Hoydysh W.G., Dabberdt W.F. Kenematics and dispersion characteristics of flows in asymmetric street canyons // Atmospheric Environment. 1988. V. 22. P. 2677- 2689.

15. Kim J, Moin P., Moser R. Turbulence statistics in fully developed channel flow at low Reynolds number // J. Fluid Mechanics. 1987. V. 177. P. 133-166.

16. Leonard B. A Stable and Accurate Convective Modeling Procedure Based on Quadratic Upstream Interpolation // Computer Methods in Applied Mechanics and Engineering. 1979. V. 19. P. 59-98.

17. Lyn D., Einav S., Rodi W., et al. A laser-Doppler velocimetry study of ensemble averaged characteristics of the turbulent near wake of a square cylinder // J. Fluid Mech. 1995. V. 304. P. 285-319.

18. Moin P., Kim J. On the numerical solution of time dependent viscous incompressible fluid flows involving solid boundaries // J. Computational Physics. 1980. V. 35. P. 381-392.

19. Niederschulte M.A., Adrian R.J., Hanratty T.J. Measurements of turbulent flow in a channel at low Reynolds numbers // Experiments in Fluids. 1990. V. 9. P. 222-230.

20. Sagaut P. Large eddy simulation for Incompressible Flow. 3rd ed. An Series: Scientific Computation, 2006. 556 p.

21. Smagorinsky J. General circulation experiments with the primitive equations. I: The basic experiment // Monthly Weather Review. 1963. V. 91. No. 3. P. 99-165.

22. Walton A., Cheng A.Y.S. Large-eddy simulation of pollution dispersion in an urban street canyon. Part 2: idealised canyon imulation // Atmospheric Environment. 2002. V. 36. P. 3615-3627.

Статья поступила 22.10.2012 г.

Danilkin E.A., Nuterman R.B., Bart A.A., Degi D.V., Starchenko A.V. STUDY OF AIRFLOW AND POLLUTANT TRANSPORT IN AN URBAN STREET CANYON USING LARGE EDDY SIMULATION OF THE TURBULENT FLOW. In this paper, the large eddy simulation model is presented, as well as a modification of the algortihm for numerical solving the system of Navier-Stokes equations. The modification is oriented to the supercomputer technique. After approbation of the constructed model on test problems, the proposed mathematical model was applied to studying the airflow and pollutant transport in an urban street canyon.

Keywords: large eddy simulation, turbulence, urban street canyons, parallel computing

DANILKIN Evgeniy Alexandrovich

E-mail: ugin@math.tsu.ru

NUTERNANRoman Borisovich

E-mail: nutrik@math.tsu.ru

BART Anrdey Andreevich

E-mail: baza@math.tsu.ru

DEGI Dmitrii Vladimirovich

E-mail: dimadegi@sibmail.com

STARCHENKO Alexander Vasil’evich

E-mail: starch@math.tsu.ru