CC BY
48
Исследование режимов виброкипящего гранулированного слоя с использованием пакета ОрепРОАМ
''2'3 Н.С. Орлова <norlova.umi.vnc@gmail.com> 2 Я.Н. Качалкина <yankevich-2502@mail.ru>
1 Южный математический институт Владикавказского научного центра
РАН и Правительства РСО-А, Россия, Владикавказ, ул. Маркуса, 22;
2 Северо-Кавказский горно-металлургический институт (государственный
технологический университет), Россия, Владикавказ, ул. Николаева, 44; 3 Финансовый университет при Правительстве РФ, Россия, Владикавказ, ул.
Молодежная, 7.
Аннотация. Виброкипение широко используется в различных технологических процессах. В связи с этим исследование режимов виброкипения является актуальным. В работе проведено исследование режимов виброкипящего слоя с использованием пакета ОрепРОАМ. Представлены результаты моделирования динамики виброкипящего слоя частиц стекла при различных значениях амплитуды и частоты колебаний, толщины слоя засыпки. Выявлены режимы, при которых слой теряет устойчивость, и образуются всплески гранулированного материала и фонтанирующие каналы.
Ключевые слова: виброкипящий слой, гранулированный материал, математическое моделирование, ОрепРОАМ, 1\уоР1ш8еЕи1егРоат, режимы виброкипения
1. Введение
В настоящее время имеющиеся модели виброкипения описывают поведение движения двухфазных сред только при определенных условиях [1-3], т.е. имеют различные области применения, которые зависят от входных параметров (толщина слоя засыпки материала, размер твердых частиц, амплитуда и частота колебаний полки, на которой располагается материал, и т.д.). Разработка универсальной трехмерной модели виброкипящего слоя имеет практическую значимость, так как за счет увеличения площади поверхности контакта фаз виброкипение широко используется в различных технологических процессах: очистка газов, химические технологии, сушка гранулированного материала. Несмотря на огромное практическое значение и
143
разнообразие работ в этом направлении, общепринятой универсальной модели пока нет.
В данной работе представлены результаты моделирования динамики виброкипящего слоя с использованием свободно распространяемого пакета для решения прикладных задач гидро и аэромеханики ОрепРОАМ при поддержке программы "Университетский кластер" с удаленным доступом к консоли на управляющем узле вычислительного кластера ВЬ2х220 [4]. Для описания процесса виброкипения был доработан решатель 1\уоРЬа5еЕи1егРоат, который использовался для моделирования динамики кипящего гранулированного слоя [5]. В решателе 1\уоРЬа5еЕи1егРоат реализована двухжидкостная модель кипящего (ожиженного) слоя на основе континуального подхода (подхода Эйлера), при котором движение слоя рассматривается как движение двух взаимодействующих континуумов, связанных с газом и частицами. Основные уравнения двухжидкостной модели - уравнения неразрывности и уравнения количества движения для обеих фаз [5-8].
2. Постановка задачи
Для моделирования динамики виброкипящего слоя был доработан решатель 1\уоРЬа5еЕи1егРоат, в котором вместо использования абсолютной системы отсчета для пространственных координат вводилась относительная система отсчета, движущаяся вместе с контейнером (полкой, на которой располагается слой материала). Предполагалось, что стенки контейнера не деформируются и перемещаются как твердое тело. В этой движущейся системе отсчета нет движения стенок. Таким образом, вычислительная сетка является статической в движущейся неинерциальной системе отсчета, что значительно упрощает процедуру численного решения уравнений. Концептуально, результаты должны быть эквивалентны результатам, полученным из решения основных уравнений в инерциальной системе координат с использованием сетки, которая движется со стенками. Поскольку рассматриваются вертикальные вибрации, считается, что будет меняться только вертикальная координата у. В связи с этим уравнения количества движения для обеих фаз примут следующий вид [9]:
где С1 - ускорение колебаний полки (стенки), на которой располагается слой материала. Т.к. рассматриваются только вертикальные колебания (вибрации), то ускорение колебаний определяется вертикальной составляющей
д
(a Ut) + V • (а ифиф) + V • (а ) = VP + а (g - а) +МФ;
рф
(1)
dt
а = ) = (о,а ,о) , где С1у = —А в} 5ш(й/). О) = 2т$ (А- амплитуда
колебаний, / - частота колебаний).
Кроме того, решалось уравнение неразрывности (2) для обеих фаз.
(2)
Индекс «^означает принадлежность к фазе (твердой «а» или газовой «/>»): а(р
- объемная доля соответствующей фазы; р^ - плотность фазы; - вектор
скорости фазы; - тензор эффективных напряжений; Р - давление
газовой фазы; g - ускорение свободного падения; М ,Р - член, моделирующий обмен импульсом между фазами. Выражения для нахождения коэффициентов и членов, входящих в уравнения (1-2), подробно описаны в литературе [6-8].
В решателе 1\уоРЬа5еЕи1егРоат обе фазы считаются несжимаемыми. Для газовой фазы реализована полуэмпирическая двухпараметрическая модель турбулентности к - е [5,6]. Член, моделирующий обмен импульсом между фазами, определяется силами трения, возникающими между двумя фазами:
где К - функция сопротивления, которая зависит от объемной доли газовой фазы. При С(/1 <0.8 функция сопротивления определяется по формуле (3), а
при аь >0.8 - по формуле (4) [6-8].
Мер =-К! drag ,
1
Mdrag = аааъКи г,
(3)
(4)
где 11 г — и а и ь - относительная скорость фаз, с1а - диаметр частиц, [1Ъ -динамическая вязкость газа, Св - коэффициент сопротивления.
Для учета эффектов, обусловленных взаимодействием частиц друг с другом, используется кинетическая теория (по аналогии с кинетической теорией газа), с помощью которой можно выразить эффективные напряжения, возникающие в дисперсной фазе за счет движения частиц и за счет столкновений частиц друг с другом. По аналогии с термодинамической температурой, вводится гранулярная температура д, как средняя энергия флуктуаций скорости частиц. В решателе для расчета гранулярной температуры используется алгебраическое уравнение (5) [7].
В уравнении (5) /)л - тензор скоростей деформации твердой фазы, коэффициенты — К^ определяются следующим образом:
в
(
(5)
2 а К
\
Ds =I|w/a+(wa)T ,
K,=2(\ + e)pag0,
^3 =
daPa
ЗО-е)
\ + -(\ + е\Ъе-\ )aag0
Л
/
\2{l-ef
Pag 0
dn4n
gо =-
1-
ос„
1
V
ос
V атах у
где g0 - радиальная функция контакта, ССатах - максимальное значение объемной доли частиц, е - коэффициент восстановления в случае столкновений частица-частица. Более подробное описание коэффициентов представлено в работах [7,8].
3. Начальные и граничные условия
Задача решалась в трехмерном приближении. В связи с тем, что введена относительная система отсчета, используются следующие начальные условия:
U a = 0; Ub = 0 ; Р = 0; Яа=0,6; 0 = 0.
Далее представлены граничные условия на левой, правой, передней и задней стенках, а также на нижней стенке (полке):
Üa=0;Üb=0;VP = 0;Vaa=0;Ve = 0.
Граничные условия на верхней свободной поверхности:
Ша = 0; VUb =0-УР = 0усса=0ув = 0.
Размеры вычислительной области: высота - 0,4 м, ширина - 0.4 м, толщина -0.4 м. Процесс виброкипения рассчитывался за 2 с. При этом использовался шаг по времени, равный 1 х 10"4 с. Шаг по координате х равен 0.005 м, шаг по координате у - 0.005 м, шаг по координате z равен 0.005 м. Проводилось распараллеливание вычислений на 12 ядрах. В таблице 1. представлены
значения входных параметров задачи. В качестве твердой фазы были выбраны сферические частицы стекла.
Табл. 1. - Значения входных параметров задачи.
Описание Значение
1 плотность твердой фазы (частицы стекла), ра [кг/м3] 2600
2 плотность газовой фазы (воздуха), рь [кг/м3] 1,2
3 вязкость газовой фазы (воздуха), ///, [Пахс] 1,5х10"5
4 коэффициент восстановления частица-частица, С [-] 0,9
5 начальное значение объемной доли частиц в слое, ам [-] 0,6
6 начальная высота слоя, //,, [м] 0,05-0,1
7 ширина слоя, /,|м| 0,4
8 толщина слоя, Щм] 0,4
9 диаметр частиц, с1а[м] 0,0003
4. Результаты расчетов
С использованием разработанной трехмерной модели виброкипящего слоя было проведено исследование режимов виброкипения. Расчеты проводились при значениях амплитуды колебаний в диапазоне А = 1,5-9 мм, частоты колебаний в диапазоне /= 20 - 60 Гц при значениях начальной толщины слоя частиц 50 - 100 мм. В результате серии трехмерных вычислительных экспериментов с распараллеливанием на супер-ЭВМ наблюдается волнообразная поверхность гранулированного материала. С увеличением амплитуды и частоты колебаний слой частиц теряет устойчивость и образуются всплески гранулированного материала.
На рис. 1 представлены результаты моделирования динамики виброкипящего слоя в момент времени 1,5 с при различных значениях амплитуды колебаний А и при частоте колебаний/= 20 Гц для слоя с толщиной засыпки 50 мм.
Рис. 1. Структура виброкипящего слоя при разных значениях амплитуды колебаний.
Из рис. 1 видно, что с увеличением амплитуды колебаний А > Змм начинают появляться всплески гранулированного материала, а при амплитуде А = 9 мм всплески больше похожи на фонтанирующие каналы.
На рис. 2 представлены аналогичные результаты, полученные при частоте колебаний/= 40 Гц.
Рис. 2. Структура виброкипящего слоя при разных значениях амплитуды колебаний.
В отличие от предыдущего случая (рис. 1) всплески гранулированного материала появляются уже при амплитуде колебаний . I = 3 мм. А при амплитуде А = 9 мм наблюдаются явно выраженные фонтанирующие каналы. Очевидно, что образование всплесков в большей степени зависит от амплитуды колебаний, чем от частоты. Они начинают появляться при значениях амплитуды А > 2 мм. При значениях амплитуды колебаний. I > 6 мм в слое наблюдаются локализованные фонтанирующие каналы. С дальнейшим увеличением частоты колебаний (до 60 Гц) при высоких значениях амплитуды (> 6 мм) количество фонтанирующих каналов
возрастает, зазор между нижней частью слоя и полкой (на которой располагается материал) становится больше.
Вычислительные эксперименты проводились и для более толстых слоев (с толщиной засыпки до 100 мм). На рис. 3 представлены результаты сравнения степени расширения виброкипящего слоя с толщиной засыпки 50 мм (кривая 1) и 100 мм (кривая 2) при частоте колебаний 20 Гц и при амплитуде колебаний 1,5-9 мм.
Нтах/Н
Рис. 3. Степень расширения виброкипящего слоя
Качественно результаты расчетов, полученные при различных значениях амплитуды и частоты колебаний для слоя с толщиной засыпки 100 мм мало отличаются от результатов, полученных для слоя с толщиной 50 мм. С увеличением амплитуды и частоты колебаний в слое наблюдаются всплески (при А > 2 мм) и фонтанирующие каналы (при А > 6 мм). Но степень расширения более толстого слоя (100 мм) меньше, чем степень расширения относительно тонкого (50 мм), как видно из рис. 3.
5. Заключение
В результате серии вычислительных экспериментов наблюдается волнообразная поверхность гранулированного материала. С увеличением амплитуды и частоты колебаний слой частиц теряет устойчивость и образуются всплески гранулированного материала. Образование всплесков в большей степени зависит от амплитуды колебаний, чем от частоты. Они начинают появляться при значениях амплитуды А > 2 мм. При значениях 150
амплитуды колебаний А > 6 мм в слое наблюдаются локализованные фонтанирующие каналы.
Таким образом, в результате моделирования виброкипящего слоя с использованием пакета ОрепРОАМ исследованы режимы виброкипения. Выявлены режимы, при которых слой теряет устойчивость и образуются всплески гранулированного материала и фонтанирующие каналы.
Список литературы
[1]. Н.С. Орлова. Тестирование двух моделей виброожиженного слоя. Известия вузов. Северо-Кавказский регион. Технические науки. №2, 2012 г. стр.42-45.
[2]. Н.С. Орлова. Сравнение расчетов по двухжидкостной модели виброожиженного слоя с экспериментальными данными. Инженерно - физический журнал, том 85, №6, 2012 г. стр.1202-1207.
[3]. Н.С. Орлова. Сравнение результатов экспериментального исследования виброкипящего слоя с расчетами по гидродинамической модели гранулярного газа. Инженерно- физический журнал. Т. 87, № 2, 2014. стр. 429-435.
[4]. Консоль на управляющий узел кластера (2012), "BL2x220 Cluster Console" // URL: https://unihub.ru/resources/bl2x220cc (дата обращения: 10.06.2014).
[5]. Е.С. Каменецкий, Н.С. Орлова, М.В. Волик, Д.Г. Минасян. Исследование динамики кипящего гранулированного слоя с использованием пакета OpenFoam. Известия вузов. Северо-Кавказский регион. Естественные науки. №5, 2014 г. стр. 37-42.
[6]. Н. Rusche. Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions. Thesis submitted for the degree of Doctor of Philosophy of the University of London and Diploma of Imperial College. 2002. 343 p.
[7]. B. van Wachem. Derivation, Implementation, and Validation of Computer Simulation Models for Gas-Solid Fluidized Beds - Dissertation at Delft University of Technology. 2000. 222 p.
[8]. D. Gidaspow. Multiphase flow and fluidization: Continuum and kinetic theory descriptions. Boston: Academic Press Inc., 1994. 211 p.
[9]. N.S. Orlova, M.V. Volik. Modelling of vibrofluidized bed dynamics using OpenFoam. Waves and vortices in complex media: 5-th International Scientific School of Young Scientists, November 25-28, 2014, Moscow. M.: MAKS Press, 2014. pp. 72-74.
Investigation of vibrofluidized granular layer modes using OpenFOAM
''2'3N.S. Orlova <norlova.umi.vnc@gmail.com> 2 Ya.N. Kachalkina <yankevich-2502@mail.ru> 1 South Mathematical Institute of the Vladikavkaz Scientific Center of the Russian
Academy of Sciences and the Government of Republic of North Ossetia-Alania, Russia, Vladikavkaz, Marcus
street, 22
2 North-Caucasian Mining and Metallurgical Institute (State Technological University), Russia, Vladikavkaz, Nikolaev street, 44 3 Financial University under the Government of the Russian Federation, Russia, Vladikavkaz, Molodezhnaya street, 7
Abstract. Vibrofluidization is widely used in various industrial processes. So the study of vibrofluidization modes is current research. In this paper we investigated the modes of vibrofluidized layer using the package OpenFOAM. The results of modelling of the vibrofluidized layer dynamics at different values of the amplitude and frequency of oscillations, the thickness of the filling, were presented. The modes in which layer becomes unstable were identified. Wavelets formed from granular material and the jetting channels appear in these modes.
Key words: vibrofluidized layer, granular material, mathematical modeling, OpenFOAM, twoPhaseEulerFoam, vibrofluidization modes
References
[1]. N.S. Orlova. Testirovanie dvuh modeley vibroozhizhennogo sloya. [Testing of two models of vibrofluidized layer], Izvestiya vuzov. Severo-Kavkazskiy region. Tehnicheskie nauki. [Proceedings of the universities. North Caucasus region, technical sciences] 2012, no. 2, pp. 42-45 (in Russian)
[2]. N. S. Orlova. Comparison of calculations by the two-field vibrofluidized-bed model with experimental data. Journal of Engineering Physics and Thermophysics vol. 85, no. 6,pp 1305-1310. doi 10.1007/sl0891-012-0775-x
[3]. N. S. Orlova. Comparison of the Results of Experimental Investigations of a Vibrofluidized Bed with Calculations by a Granular Gas Hydrodynamic Model. Journal of Engineering Physics and Thermophysics vol. 87, no. 2, pp 443-449. doi 10.1007/sl0891-014-1030-4
[4]. Konsol na upravlyayuschiy uzel klastera (2012), "BL2x220 Cluster Console" // URL: https://unihub.ru/resources/bl2x220cc (data obrascheniya: 10.06.2014). [Console to the
control node of the cluster (2012), "BL2x220 Cluster Console" // URL: https://unihub.ru/resources/bl2x220cc (date of treatment: 10.06.2014).]
[5]. E.S. Kamenetskiy, N.S. Orlova, M.V. Volik, D.G. Minasyan. Issledovanie dinamiki kipyaschego granulirovannogo sloya s ispolzovaniem paketa OpenFoam. [The study of the dynamics of fluidized granular layer using the package OpenFoam.] Izvestiya vuzov. Severo-Kavkazskiy region. Estestvennyie nauki. [Proceedings of the universities. North Caucasus region, natural sciences] 2014, no. 5, pp. 37-42 (in Russian)
[6]. H. Rusche. Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions. Thesis submitted for the degree of Doctor of Philosophy of the University of London and Diploma of Imperial College. 2002. 343 p.
[7]. B. van Wachem. Derivation, Implementation, and Validation of Computer Simulation Models for Gas-Solid Fluidized Beds - Dissertation at Delft University of Technology. 2000. 222 p.
[8]. D. Gidaspow. Multiphase flow and fluidization: Continuum and kinetic theory descriptions. Boston: Academic Press Inc., 1994. 211 p.
[9]. N.S. Orlova, M.V. Volik. Modelling of vibrofluidized bed dynamics using OpenFoam. Waves and vortices in complex media: 5-th International Scientific School of Young Scientists, November 25-28, 2014, Moscow. M.: MAKS Press, 2014. pp. 72-74.
