Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

MODELING OF PROCESSES OF VENTELATION

Datsyuk T.A. 1 Sauts A.V. 1 Yurmanov B.N. 1 Taurit V.R. 1
1 Saint-Petersburg State University of Architecture and Construction
The work deals with the modeling of ventilation and dispersion of its emissions in the atmosphere. Simulation is based on solving a system of Navier-Stokes equations, conservation laws of mass, momentum, heat. Cover various aspects of the numerical solution of the given equations. Proposed system of equations, which allows to calculate the value of the background of the coefficient of turbulence. For hypo sonic approach proposed solution together with those given in article equations of fluid dynamics equations standing of the ideal of a real gas and steam. This equation is a modification of the equation of Van-der-Waals forces and more accurately take into account the size of the molecules of a gas or vapor and their interaction. On the basis of the terms of the thermodynamic stability of the obtained ratio, which allows to exclude physically impossible roots when solving the equation of the volume. Is performed the analysis of the known settlement models, and computational packages fluid dynamics.
modeling
ventilation
turbulence
equations of heat and mass transfer
the equation of state
the real gas
steam
dissipation
Введение

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

Цель статьи - использование аналитических зависимостей, с помощью которых выполняются расчеты количества вредных выделений, определить размеры каналов, воздуховодов, шахт и выбор способа обработки воздуха и т.д. При этом целесообразно использовать программный продукт «Поток» с модулем «VSV». Для подготовки исходных данных необходимо наличие схем проектируемых вентиляционных систем с указанием длин участков и расходов воздуха на концевых участках. Входными данными для расчета являются описание систем вентиляции и требования, предъявляемые к ней. Используя математическое моделирование, решаются следующие вопросы:

  • выбор оптимальных вариантов подачи и удаления воздуха;
  • распределение параметров микроклимата по объему помещений;
  • оценка аэродинамического режима застройки;
  • выбор мест для воздухозабора и удаления воздуха.

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

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

Для моделирования турбулентных потоков необходимо решать систему уравнений сохранения массы и Рейнольдса (сохранения импульса):

(1)

(2)

где t - время, X=Xi,j,k - пространственные координаты, u=ui,j,k - компоненты вектора скорости, р - пьезометрическое давление, ρ - плотность, τij - компоненты тензора напряжений, sm - источник массы, si - компоненты источника импульса.

Тензор напряжений выражается в виде:

(3)

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

Для получения информации о полях температуры Т и концентрации с вредных веществ система дополняется следующими уравнениями:

уравнение сохранения количества тепла

(4)

уравнение сохранения пассивной примеси с

(5)

где Cр - коэффициент теплоёмкости, λ - коэффициент теплопроводности, k=ki,j,k - коэффициент турбулентности.

Базовый коэффициент турбулентности kбаз определяется с помощью системы уравнений:

 (6)

где kф - фоновый коэффициент турбулентности, kф =1-15 м2/с; ε = 0,1-04;

Коэффициенты турбулентности определяются с помощью уравнений:

 (7)

На открытой территории при малой диссипации значение kz определяется по уравнению:

kk = k0 z /z0; (8)

где k0 - значение kk на высоте z0 (k0 = 0,1 м2/с при z0 = 2 м).

На открытой площадке профиль скоростей ветра не деформирован, т.е.

При неизвестной стратификации атмосферной на открытой площадке профиль скоростей ветра можно определить:

 ; (9)

 (10)

где z0 - заданная высота (высота флюгера); u0 - скорость ветра на высоте z0; B = 0,15.

При соблюдении условии (10) локальный критерий Ричардсона Ri определяется как:

 (11)

Продифференцируем уравнение (9), приравняем уравнения (7) и (8), оттуда выразим kбаз

 (12)

Приравняем уравнение (12) с уравнениями системы (6). В полученное равенство подставим (11) и (9), в окончательном виде получим систему уравнений:

 (13)

Пульсационный член , следуя идеям Буссинеска, представляется в виде:

(14)

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

 (15)

 (16)

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

Для турбулентных потоков, изучаемых в вентиляционной практике, целесообразно использовать или гипотезу Буссинеска о малости изменений плотности, или так называемое «гипозвуковое» приближение. Напряжения Рейнольдса считаются пропорциональными осредненным по времени скоростям деформаций. Вводится коэффициент турбулентной вязкости , данная концепция выражается как:

. (17)

Коэффициент эффективной вязкости вычисляется как сумма молекулярного и турбулентного коэффициентов:

 (18)

«Гипозвуковое» приближение предполагает решение совместно с приведенными выше уравнениями уравнения стояния идеального газа:

ρ = p/(RT) (19)

где p - давление в окружающей среде; R - газовая постоянная.

Для более точных расчётов плотность примеси можно определить, используя модифицированное уравнение Ван-дер-Ваальса для реальных газов и паров

 (20)

где константы N и M - учитывают ассоциацию/диссоциацию молекул газа или пара; а - учитывает прочее взаимодействие; b' - учитывающая размеры молекул газа; υ=1/ρ.

Выделяя из уравнения (12) давление р и дифференцируя его по объёму (учёт термодинамической устойчивости) получится следующее соотношение:

 . (21)

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

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

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

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

  • прямое моделирование «DNS» (решение нестационарных уравнений Навье - Стокса);
  • решение осредненных уравнений Рейнольдса «RANS», система которых, однако, незамкнута и нуждается в дополнительных замыкающих соотношениях;
  • метод крупных вихрей «LES», который основан на решении нестационарных уравнений Навье - Стокса с параметризацией вихрей подсеточного масштаба;
  • метод «DES», который является комбинацией двух методов: в зоне отрывных течений - «LES», а в области «гладкого» потока - «RANS».

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

В настоящее время успешно применяются для моделирования задач вентиляции сертифицированные пакеты, такие как «STAR-CD», «FLUENT» или «ANSYS/FLOTRAN». При правильно сформулированной задаче и рациональном алгоритме решения получаемый объем информации позволяет на стадии проектирования выбрать оптимальный вариант, но выполнение расчетов с использованием данных программ требует соответствующей подготовки, и некорректное их использование может привести к ошибочным результатам.

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

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

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

Наиболее простыми являются однопараметрические модели, использующие концепцию турбулентной вязкости «nt», а турбулентность предполагается изотропной. Модифицированный вариант модели «nt-92» рекомендуется при моделировании струйных и отрывных течений. Хорошее совпадение с результатами эксперимента дает также однопараметрическая модель «S-A» (Спаларта - Альмараса), которая содержит уравнение переноса для величины .

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

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

Хорошо известна модель, включающая уравнения для переноса энергии турбулентности k и скорости диссипации турбулентной энергии ε. Модели типа «k- могут использоваться как для пристеночных течений, так и для более сложных отрывных течений.

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

В настоящее время к числу наиболее перспективных относят модели «SSG» и «Gibson-Launder», где используется нелинейная связь тензора турбулентных напряжений Рейнольдса и тензора осредненных скоростей деформаций. Они разрабатывались для улучшения прогнозирования отрывных течений. Поскольку в них рассчитываются все компоненты тензоров, они требуют больших компьютерных ресурсов по сравнению с двухпараметрическими моделями.

Для сложных отрывных течений некоторые преимущества выявило применение однопараметрических моделей «nt-92», «S-А» по точности предсказания параметров течения и по скорости счета по сравнению с двухпараметрическими моделями.

Например, в программе «STAR-CD» предусмотрено использование моделей типа «k-e», Спаларта - Альмараса, «SSG», «Gibson-Launder», а также метод крупных вихрей «LES», и метод «DES». Два последних метода лучше подходят для расчета движения воздуха в условиях сложной геометрии, где будут возникать многочисленные отрывные вихревые области, но они требуют больших вычислительных ресурсов.

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

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

Задание с достаточной степенью точности геометрических особенностей исследуемого объекта. Для этих целей можно рекомендовать для построения трехмерных моделей такие пакеты, как «SolidWorks», «Pro/Engeneer», «NX Nastran». При построении расчетной сетки количество ячеек выбирается так, чтобы получить достоверное решение при минимальном времени расчета. Выбрать следует одну из полуэмпирических моделей турбулентности, являющейся наиболее эффективной для рассматриваемого течения.

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

Рецензенты:

  • Воликов Анатолий Николаевич, доктор технических наук, профессор кафедры теплогазоснабжения и охраны воздушного бассейна, ФГБОУ ВПОУ «СПБГАСУ», г. Санкт-Петербург.
  • Полушкин Виталий Иванович, доктор технических наук, профессор, профессор кафедры отопления, вентиляции и кондиционирования воздуха, ФГБОУ ВПО «СПбГАСУ», г. Санкт-Петербург.