УДК 66.063.8
Ю. В. Караева, И. А. Трахунова, Г. Р. Халитова,
С. И. Исламова
ВЛИЯНИЕ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ МЕТАНТЕНКА НА ЭФФЕКТИВНОСТЬ ПРОЦЕССА МЕТАНОВОГО БРОЖЕНИЯ
Ключевые слова: метантенк, гидравлическое перемешивание.
Исследовано влияние геометрических параметров метантенка на эффективность процесса метанового брожения биоотходов при использовании двух систем гидравлического перемешивания. Проведены численные исследования для аппаратов с соотношением высоты к радиусу H/R от 1 до 2 и объемом от 3 до 100 м3.
Keywords: anaerobic digester, hydraulic mixing.
We studied the influence of geometrical parameters on the efficiency of the anaerobic methane fermentation digester for two hydraulic mixing systems. Numerical investigations were performed for tanks with height-to-radius ratio H/R from 1 to 2 and capacities from 3 to 100 m3.
Введение
В процессе метанового брожения биоотходов важную роль играет перемешивание. Перемешивание способствует равномерному распределению популяции бактерий, рН и температуры; повышает контакт между микроорганизмами и субстратом; предотвращает формирование осадка на дне и образование корки на свободной поверхности; позволяет предотвратить формирование пустот и скоплений, уменьшающих рабочую площадь метантенка; помогает освободить пузырьки биогаза [1, 2, 3]. Согласно ГОСТ Р 537902010, оптимальное перемешивание сырья в метантенке увеличивает выход биогаза на 50%.
Особенностью системы гидравлического перемешивания содержимого метантенка является то, что она позволяет создать оптимальные условия для метаногенного консорциума, представляющего собой уникальную группу бактерий, чувствительную к изменениям температуры и различного рода механическим воздействиям.
Способы перемешивания
Учитывая особенности процесса метанового брожения, разработана установка анаэробной переработки органических отходов с системой гидравлического перемешивания (патент на полезную модель № 115350). Отличительной
особенностью системы перемешивания метантенка является наличие модифицированного контура, в котором перемешивание осуществляется рециркуляцией субстрата через два подводящих и один отводящий патрубок, расположенный в верхней части метантенка. Биомасса под давлением подается через нижний подводящий патрубок в придонную область метантенка, препятствуя образованию осадка, и через верхний подводящий патрубок на уровень зеркала жидкости, препятствуя формированию корки на свободной поверхности в метантенке.
На основе описанного в патенте принципа работы биогазовой установки смоделирован метантенк с системой гидравлического
перемешивания, представленный на рис.1 (способ
7111^1/у/^ггпг>п идпд ллґ>ниі/?пгліа 1\
Рис. 1 - Способ гидравлического перемешивания 1
Разработанный способ гидравлического перемешивания, представляет собой комбинацию трех внешних контуров циркуляции органического субстрата, описанных в [4, 5], и одного
модифицированного контура, описанного в патенте [6].
Предложенный способ сравнивается со способом гидравлического перемешивания,
представленным в [4, 5], как широко
распространенным на практике (способ
гидравлического перемешивания 2).
Рециркуляция органического субстрата осуществляется с помощью внешнего циркуляционного контура, представляющего собой трубу с центробежным насосом, при этом подача органического субстрата осуществляется через боковую стенку в придонную область аппарата, а отвод - из верхней части метантенка. В работе рассмотрен метантенк с системой гидравлического перемешивания, состоящей из четырех циркуляционных контуров (рис.2).
Таким образом, в работе рассмотрено два способа гидравлического перемешивания: модифицированный способ 1 и стандартный способ 2.
Целью численных исследований является оценка эффективности модифицированной системы гидравлического перемешивания при различных геометрических параметрах метантенка.
Рис. 2 - Способ гидравлического перемешивания 2
Математическая модель гидравлического перемешивания
Для оценки качества перемешивания рассмотрим процесс гидравлического
перемешивания в реакторе метанового брожения. Математическое моделирование происходящих при этом процессов проведено при следующих допущениях.
a) Субстрат является гомогенной жидкостью, вязкость и плотность которой зависят от концентрации.
b) Процессы диффузии пренебрежимо малы по сравнению с конвективным переносом массы.
c) Температура постоянна.
dV ... -
p-----= divT + p ■ c ■ g
dT
Tdlv(v )=0 = 0
dT
где V - вектор скорости, м/с; т - время, с;
(1) (2) (3)
бт
субстанциональная (полная) производная по
времени; Т - тензор напряжений, Т = -р! + 2м{° )о; р - давление, Па; I - единичный тензор; й - тензор скоростей деформации; /и - вязкость субстрата, Пах; д - вектор ускорения свободного падения.
При этом плотность смеси определяется плотностью фаз, а эффективная вязкость связана с объемной концентрацией дисперсной среды формулой Эйнштейна [7].
Для получения решения системы дифференциальных уравнений (1) - (3) необходимо задать краевые условия. Начальное распределение объемной концентрации при т = 0 задается для условия наличия осадка на дне метантенка.
Начальные условия: V = о , с0 = f(х,у,г). Граничные условия для скорости:
- на твердых границах резервуара задаются условия прилипания жидкости;
- на входных отверстиях задаются профили скорости, соответствующие сформировавшемуся профилю скорости ньютоновской жидкости в
круглой трубе, концентрация на входе равна концентрации на выходе;
- полагая, что выход из резервуара представляет трубу такой длины, на выходе из которой можно считать движение жидкости установившимся, принимаем в качестве граничных условий на выходе
- условие стабилизации скорости;
- на открытой границе: (рд - р)п + 2рйП = 0 , где
рд - давление газа, п - нормаль к поверхности субстрата
В качестве критерия перемешивания предлагается безразмерная величина [8]:
K T)= Vc (г)/Vc0
(4)
где Vc (г) - средний по объему метантенка выход биогаза, зависящий от распределения поля концентрации органического субстрата; Vc0 -
предельный выход биогаза.
Результаты численных исследований
Для решения поставленной задачи
использовался пакет моделирования COMSOL Multiphysics 3.5, позволяющий решать системы нелинейных дифференциальных уравнений в частных производных методом конечных элементов. На основе возможностей СО.MSOL Multiphysics 3.5 была создана новая пользовательская модель,
адаптированная к условиям рассматриваемой задачи.
Достоверность полученных результатов обеспечивается тем, что в основе математической модели лежат классические уравнения механики сплошных сред. Кроме того, для верификации метода решения были проведены численные расчеты с использованием пакета прикладных программ COMSOL Multiphysics 3.5 задачи
истечения круглой струи, имеющей аналитическое решение [9].
В данной работе исследованы метантенки с объемом рабочей области от 3 до 100 м3 и соотношением 1 < H / R < 2 [10, 11],
R0/R = 0.025 , где Н - высота метантенка, R -радиус метантенка, R0 - радиус входных и
выходных патрубков.
Соотношение высоты циркуляционных патрубков и высоты метантенка постоянно и составляет:
hjH = 0.95 ; h2/H = 0.05 ; h3/H = 0.2,
где hi - высота верхних циркуляционных патрубков; h2 - высота нижних циркуляционных патрубков; h3 - высота бокового циркуляционного патрубка в модифицированном контуре.
Мощность на прокачку органического субстрата в рассматриваемых системах гидравлического перемешивания одинакова.
В результате расчетов получены данные о качестве перемешивания для метантенков с различными геометрическими размерами.
На рис. 3 представлены значения
отношения К1/К2, где К - критерий оценки качества перемешивания для метантенков со способом перемешивания 1, К2 - критерий оценки качества перемешивания для аппаратов со способом перемешивания 2.
К,/К2
Рис. 3 - Эффективность перемешивания в зависимости от соотношения И/Я и объема метантенка
Значения отношения К1/К2 больше 1 указывает на то, что применение модифицированного способа перемешивания эффективно. Как видно из рисунка, с увеличением объема аппарата качество перемешивания в метантенке с модифицированным способом повышается относительно стандартного способа перемешивания для всех аппаратов с рассмотренными соотношениями Н/Я.
Максимальный эффект использования
модифицированной системы перемешивания наблюдается для метантенков с соотношением Н/Я = 1. Использование модифицированной системы гидравлического перемешивания в случае Н/Я = 1.5 целесообразно для метантенков больших объемов. Для метантенков с соотношением Н/Я=2 значения отношения К^К2 близки к 1, что объясняется тем, что верхняя подающая струя не достигает нижних слоев и не влияет на перемешивание осадка в центральной части придонной области.
Предложенный критерий оценки качества перемешивания отражает интенсивность выхода биогаза - целевого продукта процесса метанового брожения. На рис. 4 представлено количество биогаза, получаемое при оптимальном режиме перемешивания, для метантенков различных объемов с соотношением Н/Я = 1 и Н/Я = 1.5 .
Количество биогаза, получаемое в аппаратах с соотношением Н/Я = 2 не представлено, т.к. разница в объеме получаемого биогаза при рассмотренных способах гидравлического перемешивания незначительна.
Из диаграммы видно, что при модифицированном способе гидравлического перемешивания возможен больший выход биогаза для всех рассмотренных объемов. Следует отметить, что при одинаковом объеме метантенков,
соотношение H/R = 1 обеспечивает больший выход биогаза.
W, м3
1600п
Рис. 4 - Выход биогаза в зависимости от способа перемешивания и соотношения H/R
Заключение
Исследовано влияние геометрических пропорций аппарата на эффективность процесса метанового брожения биоотходов при использовании двух систем гидравлического перемешивания. Проведенные численные исследования показали, что использование модифицированной системы гидравлического перемешивания эффективно при соотношении H/R < 1.5. Максимальный эффект от использования предложенного способа
перемешивания будет достигаться при соотношении H/R < 1.
Работа выполнена в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России на 2009-2013 годы» (Соглашение №14.B37.21.0299 между
Министерством образования и науки Российской Федерации и ФГБОУ ВПО «КНИТУ»).
Литература
1. M. Elnekave, N. Tufekci, S. Kimchie, G. Shelef Tracing the mixing efficiency of a primary mesophilic anaerobic digester in a municipal wastewater treatment plant. Fresen. Environ. Bull., 15, 9B, 1098-1105 (2006).
2. A.P. Borole, K.T. Klasson, W. Ridenour, J. Holland, K. Karim, M.H. Al-Dahhan. Methane production in a 100-L upflow bioreactor by anaerobic digestion of farm waste. Appl. Biochem. Biotechnol, 131, (1-3), 887-896 (2006).
3. А.З. Миндубаев, С.Т. Минзанова, Е.В. Скворцов, В.Ф. Миронов, В.В. Зобов, Ф.Ю. Ахмадуллина, Л.Г. Миронова, Д.Е. Белостоцкий, А.И. Коновалов, Оптимизация параметров выработки биогаза в лабораторном масштабе, Вестник Казан. технол. ун-та, 4, 233 - 240 (2009).
4. Веденев А.Г. Биогазовые технологии в Кыргызской Республике / А.Г. Веденев, Т.А. Веденеева - Бишкек: Типография «Евро», 2006. - 90 с.
5. Веденев А.Г. Строительство биогазовых установок. Краткое руководство / А.Г. Веденев, А.Н. Маслов -Бишкек: Типография «Евро», 2006. - 28 с.
6. Пат. 115350 Российская Федерация. МПК U 1,
C02F3/28 (2006.01) Установка анаэробной переработки органических отходов / Караева Ю.В., Халитова Г.Р., Трахунова И.А.; заявитель и патентообладатель:
Учреждение Российской академии наук Казанский научный центр РАН. -№ 2011124487/05, заявл.
16.06.2011; опубл. 27.04.2012.
7. Р.И. Нигматуллин, Динамика многофазных сред. Т. 1, Наука, Москва, 1987. 464 с.
8. Караева Ю.В., Трахунова И.А., Халитова Г.Р., Даминов А.З. Оценка качества перемешивания в реакторах метанового брожения. Вестник Казан. технол. ун-та 10, 226-230 (2012).
9. Г. Шлихтинг. Теория пограничного слоя. Наука, Москва, 1974. 713 с.
10. Яковлев С.В. Канализация / С.В. Яковлев, Я.А. Карелин, А.И. Жуков, С.К. Колобанов. - М.: Стройиздат, 1975. - 632 с.
11. Пааль Л.Л. Справочник по очистке природных и сточных вод // Л.Л.Пааль, Я.Я.Кару, Х.А. Мельдер, Б.Н. Репин. - М.: Высш.шк., 1994. - 336 с.
© Ю. В. Караева - к.т.н., ст. науч. сотр. лаб. ЭТПИЭ Исследовательского центра проблем энергетики КазНЦ РАН; вед. науч. сотр. каф. химической технологии древесины КНИТУ, [email protected]; И. А. Трахунова - мл. науч. сотр. лаб. ЭТПИЭ Исследовательского центра проблем энергетики КазНЦ РАН, мл. науч. сотр. каф. химической технологии древесины КНИТУ, [email protected]; Г. Р. Халитова - к.т.н., ст. науч. сотр. лаб. ЭТПИЭ Исследовательского центра проблем энергетики КазНЦ РАН; вед. науч. сотр. каф. химической технологии древесины КНИТУ, [email protected]; С. И. Исламова - науч. сотр. лаб. ЭТПИЭ Исследовательского центра проблем энергетики КазНЦ РАН.