employee from 01.01.2007 until now
Belgorod, Russian Federation
Gubkin, Belgorod, Russian Federation
VAC 05.17.00 Химическая технология
VAC 05.23.00 Строительство и архитектура
GRNTI 20.53 Технические средства обеспечения информационных процессов
Recently in the Russian Federation observed the development of systems of alternative energy - biogas technology. At the same time to intensify the process of obtaining biogas installation equipped with heating systems and mixing. A computer model of a bioreactor with mechanical stirring of the substrate, the most widely used in biogas plants. Produced simulation of mixing biomass in the multifunctional complex ANSYS, in which determine the direction and flow rate of components under mechanical mixing method. It was established that in the top of the reactor biomass hardly moves, while at the bottom of the reactor there is a twofold overspeed.
computer modeling, mechanical stirring, bioreactor
Введение. В последнее время на территории РФ в сельском хозяйстве начали внедряться биогазовые станции и установки, основной задачей которых является утилизация органических отходов [1, 2]. Утилизация происходит за счет активной жизнедеятельности микроорганизмов, которые переваривают белки, жиры, углеводы. Что бы обеспечить комфортную среду для размножения и жизнедеятельности микроорганизмам, необходимо поддерживать определенный, равномерный уровень температуры и распределение бактерий по всему объему биомассы. Для решения данных проблем применяют различные системы перемешивания, при этом в биогазовых установках наибольшее применение получил механический способ [2].
Известно, что для комфортного размножения бактерий, накладываются определенные ограничения, а именно: скорость движения потоков биомассы не должна превышать 0,6 м/с, а направление потоков рекомендуется создавать за счет тангенциальных и аксиальных скоростей [3]. Это обосновывается градиентом плотностей в вертикальной плоскости реактора, а также возможностью создания плавающей корки, препятствующей выходу газа. Поэтому актуальным является исследование скорости и направления движения биомассы при механическом перемешивании. Эффективным методом исследования гидродинамики двухфазных потоков, получившим в последнее время широкое применение, является компьютерное моделирование [4–6].
Методология. В работе использовались методы математического и компьютерного моделирования. Для создания модели биореактора использовался программный комплекс SolidWorks 2015 – САПР. Моделирование процесса перемешивания производилось методом Sliding-mesh MRF вращения в программном пакете ANSYS FLUENT 15.0.7.
Основная часть. Разработан биореактор, представляющий собой резервуар цилиндрической формы с расположенной внутри мешалкой (рис. 1).
Биореактор заполнен до высотной отметки 2,2 метра. В качестве перемешиваемой среды была выбрана жидкость с плотностью ρ = 1016,04 м3/кг и динамической вязкостью µ = 0,45 Па·с, что соответствует навозу крупного рогатого скота влажностью 92 % при температуре 37 °С.
Рис. 1. Модель биореактора с мешалкой
Математическая постановка.
Движение вязкой, несжимаемой жидкости описывается уравнением Навье-Стокса и уравнением неразрывности:
(1)
(2)
где i,j = 1,2,3, предполагается суммирование по одинаковым индексам, x1, x2, x3 – пространственные координаты, t – время, fi – массовые силы.
В этой системе уравнений независимыми искомыми параметрами являются 3 компоненты скорости и давление p.
В качестве граничного условия задается условие прилипания на всех твердых стенках: поверхность емкости и мешалки:
(3)
Аналитическое решение уравнений (1, 2) возможно только для локальных задач (течение Пуазейля и плоское течение Куэтта), поэтому более сложные задачи выполняют численным решением.
Из множества программ, для решения задач CFD, был выбран коммерческий пакет ANSYS FLUENT 15.0.7.
Из физической постановки следует, что рассматриваемая задача перемешивания нестационарная, т.к. имеется подвижный элемент-мешалка, которая вращается с постоянной угловой скоростью в неподвижной системе отсчета.
Наличие пакета FLUENT в ANSYS R15 позволяет смоделировать движущуюся систему отсчета и тем самым решить задачи, в которых имеются движущиеся части. Это делается путем установки в выбранных зонах ячеек движущейся системы отсчета. Когда используется движущаяся система отсчета, уравнения движения преобразуются путем введения дополнительных членов ускорения, которые появляются из-за преобразования неподвижной системы отсчета в движущуюся систему отсчета. Решая эти уравнения установившемся способом, может быть смоделирован поток вблизи движущихся частей.
В нашей задаче использовалось MRF моделирование (Multiple Reference Frame), позволяющее смоделировать вращение движущихся частей с одновременным выходом из них газа.
Рассмотрим систему координат, которая вращается с постоянной угловой скоростью относительно неподвижной системы отсчета как показано на рисунке 2. Вектор
показывает начало вращающейся системы координат.
Ось вращения задается единичным вектором направления
так, что бы:
(4)
В CFD задачах вычислительная область позиционируется относительно вращающейся системы координат так, чтобы любая точка в области вращения определялась вектором
из начала вращающейся системы.
|
Y |
|
Z |
|
X |
|
Стационарная система |
|
координат |
|
z |
|
x |
|
y |
|
Подвижная система координат |
|
CFD |
|
Область |
|
Ось вращения |
|
r |
|
0 |
|
r |
|
R |
|
w |
Рис. 2. Вращающаяся система отсчета относительно неподвижной
Переместить скорость движения жидкости из неподвижной системы в подвижную можно за счет применения следующего отношения:
(5)
(6)
где
– относительная скорость, рассматриваемая относительно вращающейся системы координат,
– абсолютная скорость относительно стационарной системы,
– скорость «вихря», возникающего из-за движущейся системы отсчета.
Число Рейнольдса Re при течении в трубе определяется соотношением:
(7)
где u – средняя скорость движения жидкости в трубе; d0 – внутренний диаметр трубы; ρ – плотность.
Для процессов, осуществляемых в аппаратах для перемешивания жидкостей, трудно определить критерии, в состав которых входит скорость u. Распределение скоростей в аппаратах с мешалками носит очень сложный характер, и в этом случае применять среднюю скорость u как при течении жидкости в трубе невозможно. Поэтому используют условную скорость, равную окружной скорости конца лопастей мешалки, опуская множитель π = 3,14.
Для течения жидкости в аппаратах с мешалками принимают условно 𝝎=πndm, и d0 = dm (диаметр мешалки) [19]. После исключения множителя π = 3,14 как постоянной величины получаем:
(8)
Для создания модели использовался программный комплекс SolidWorks 2015 – САПР, предназначенная для разработок изделий любой степени сложности и назначения.
Расчетная модель состоит из 3-х элементов:
– перемешивающее устройство;
– емкость биореактора;
– область вращения, обусловленная принципом работы sliding-mesh вращения в FLUENT.
Емкость биореактора разбита на две части. В нижней части расположена область вращения. Это позволяет детализировать расчетную сетку в зоне возможной турбулентности, возникающей вне области вращения тем самым уточнив расчет.
Наиболее важным этапом при численном расчете любой задачи является качество построенной сетки. В нижней части емкость биореактора и в зоне вращения использовались ячейки в форме tetraider, это обеспечивает лучшую сходимость при коэффициенте турбулизации выше 1000. Минимальный размер ячеек 2,3е-2.
В верхней части емкости биореактора использовались ячейки в форме hexaeder, обеспечивающие лучшую сходимость при низких значениях турбулизации. Максимальное значение ячеек – 5e-2.
Число ячеек составило более 2 млн. единиц. При построении сетки использовалась размерная функция, которая позволила наложить более детальную сетку вблизи мешалки, так как в этих областях предполагаются наибольшие градиенты скорости и давления.
Был произведен расчет на сходимость заданных параметров, одним из показателей качества моделирования является плавное понижение экспоненты, без резких скачков и провалов (см. рис. 3).
Рис. 3. График сходимости невязок методом конечных объемов
Расчет производился для Re = 4077,12, в результате которого было получено значение угловой скорости ω = 3 рад/с. В качестве режима движения жидкости использовалось уравнение Эйлера. Сходимость достигнута на 430 итерации.
Результаты моделирования распределения скоростей представлены на рисунке 4.
Рис. 4. Скорости движения биомассы
при механическом перемешивании
Как видно из рисунков в верхней части биореактора биомасса практически не движется, вследствие чего возникают застойные зоны, уменьшающие эффективность работы биореактора и препятствующие выходу газа. Кроме того в нижней части биореактора наблюдается двукратное превышение допустимых скоростей (0,8-0,9 м/с), что негативно сказывается на жизнедеятельности метанообразующих бактерий.
Выводы. Разработана компьютерная модель биореактора получения биогаза с механическим перемешиванием. Произведено моделирование процесса перемешивания биомассы в биореакторе в многофункциональном комплексе ANSYS. Получены составляющие скорости движения потока биомассы при механическом способе перемешивания.
1. Evstyunichev M.A., Il'ina T.N. Osobennosti syr'evoy bazy Belgorodskoy oblasti dlya proizvodstva biogaza // Vestnik Belgorodskogo gosudarstvennogo tehnologicheskogo universiteta im. V.G. Shuhova. 2013. № 5. S. 170-173.
2. Suslov D.Yu. Poluchenie biogaza v bioreaktore s barbotazhnym peremeshivaniem: dissertaciya ... kandidata tehnicheskih nauk: 05.17.08 // Ivanovskiy gosudarstvennyy himiko-tehnologicheskiy universitet. Belgorod, 2013.
3. Eder B., Shul'c H. Biogazovye ustanovki // perevod s nem. kompanii «Zorg Biogas»: 2008. 268 s.
4. Feoktistov A. Yu., Staroverov S. V., Gol’tsov A. B., Kireev V. M. Design of aspiration shelters for constricted conditions // Chemical and Petroleum Engineering. 2013. T. 49. №. 3-4. pp. 261-264.
5. Shaptala V.G., Shaptala V.V., Suslov D.Yu. Voprosy modelirovaniya i rascheta barbotazhnyh reaktorov // Vestnik Belgorodskogo gosudarstvennogo tehnologicheskogo universiteta im. V.G. Shuhova. 2013. №5. S. 189-192.
6. Kafarov V.V., Vinarov A.Yu., Gordeev L.S. Modelirovanie biohimicheskih reaktorov. M., Lesnaya prom-st', 1979. 344 s.
