Исследование процесса движения волокнистой суспензии в разгонном узле дезинтегратора мокрого помола

УДК 621.926.7

В. С. Севостьянов, А. Е. Качаев, С. А. Михайличенко, Т. Л. Сиваченко, А. А. Фарафонов

ИССЛЕДОВАНИЕ ПРОЦЕССА ДВИЖЕНИЯ ВОЛОКНИСТОЙ СУСПЕНЗИИ В РАЗГОННОМ УЗЛЕ ДЕЗИНТЕГРАТОРА МОКРОГО ПОМОЛА

UDC 621.926.7

V. S. Sevostyanov, А. Y. Kachayev, S. А. Mikhailichenko, Т. L. Sivachenko, А. А. Farafonov

STUDY OF MOVEMENT OF FIBROUS SUSPENSION IN THE DISPERSAL NODE OF A DISINTEGRATOR FOR WET GRINDING

Аннотация

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

Ключевые слова:

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

Abstract

Materials reinforced by fibers are used in the construction and the formation of their homogeneous structure is an important task. Uniform mixing of long fibers in the liquid and a stable movement of such suspension in the nodes having a complex shape requires the development of analytical models of these processes. The study of the movement of fibrous suspension in the dispersal node of a disintegrator is important and urgent.

Key words:

fibrous suspension, disintegrator, two-phase flow, absolute speed, dispersal node, wet milling.

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

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

Основная часть

Известна конструкция дезинтегратора [2, 3], в котором, помимо «сухого» помола, возможен процесс диспергирования волокнистых суспензий. На рис. 1 представлена камера помола дезинтегратора, оснащенная разгонным узлом.

© Севостьянов В. С., Качаев А. Е., Михайличенко С. А., Сиваченко Т. Л., Фарафонов А. А., 2016

Именно этот узел, трехмерная модель которого приведена на рис. 2, определяет производительность агрегата, его

эффективность на стадии разгона суспензии перед её комплексной гидродинамической обработкой.

Рис. 1. Конструкция камеры помола дезинтегратора, оснащенная разгонным узлом: 1 - наружный

ротор; 2 - внутренний ротор; 3 - разгонный узел

Рис. 2. Трехмерная модель разгонного узла с тороидально-эллиптическими патрубками, установленного на валу внутреннего ротора дезинтегратора: 1 - вращающийся кольцевой канал; 2 - тороидально-

эллиптический канал

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

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

В решении задачи рассмотрим стационарное осесимметричное закрученное течение вязкой несжимаемой жидкости, содержащей незначительное по массе количество примеси в виде волокон. Движение суспензии осуществляется в соответствии с рис. 2: на входе в разгонный узел - во вращающемся кольцевом канале постоянного сечения - участок 1; на выходе из разгонного узла - в сужающемся тороидально-эллиптическом канале - участок 2.

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

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

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

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

йи йу 1 йр о

и— + у— =---— + у-V -и;

йг р

и--+ V---

й2 йг Г

2

1 йр =---—+ V

р йг

с у ^

V2,

г У

йм йм V - м

и--+ V--I--

йг Г

2

= и

с ЛХ, ^

V м- — г 2

V

а уравнение неразрывности потока представим как

йи 1 й(V - г) Л

— +----- = 0,

г йг

где г - цилиндрические координаты; и, V, м - осевая, радиальная и окружная скорости потока в выбранной системе координат, м/с; р, р - плотность и давление волокнистой суспензии соответст-

венно, кг/м ; Па; V - плоский оператор Лапласа,

д2 д

1 д

V2 + + -Т".

дz дг г дr

Решим выражения (1)...(4) последовательно, используя метод подстановки и простейшее автомодельное решение задачи (все распределения которого зависят от одной переменной координаты) о закрученном течении вязкой жидкости в кольцевом канале, для чего воспользуемся равенствами: и(г) = U(г), г) = V(г), м>(г) = W(ог).

Тогда получим, что давление потока

р

где о - угловая скорость вращения разгонного узла, рад/с; р^г) - линейная функция давления от г; г1 - внутренний радиус кольцевого пространства, равный диаметру вала внутреннего ротора гв; г2 - внешний радиус кольцевого пространства.

Сила внутреннего трения потока, действующего на г2 в направлении течения, зависит от осевой скорости потока и равна:

2п • ИТ • ц

с1и (г)4 йг

г + йг

йи (г)4

йг

= 2п- Ит  ц— йг

й ( йи (г) ^

йг

йг,

где ИТ - высота кольцевого канала, м; и(г) - осевая скорость потока суспензии в канале, м/с; йг - изменение радиуса, м; Л - динамическая вязкость волокнистой суспензии, Па-с.

При стационарном течении двухфазного потока сумма обеих сил обра-

щается в нуль. Поэтому на основании предположений Гогена и закономерностей Пуазеля [4] можно записать следующее уравнение:

йг

йи (г) ~ йг

РЬР г,

Ит-ц '

где Р2 - давление суспензии в нижней (конечной) части кольцевого пространства, МПа; Р1 - давление суспензии в верхней (начальной) части кольцевого пространства, МПа.

Решение (8) обращается в нуль при г = г1 и при г = г2, следовательно,

и

Р - Р 1 2 1 1

4ИТц

[

+

+

2 2 г - г

С

Мп

1п

г

V '1 У

( г ^

V г2 У

Радиальная составляющая скорости не зависит от координаты г и прямо пропорциональна координате г. Тогда из (1) получим

V(г) ^^ - V -V2 •V(г) = 0. (10) йг

Решение уравнения (10) в плоскости ХУ приводится к виду (при усло-

д2

вии, что — = 0 )

дг

V(г) = 2v -V2 •г .

Из уравнения (4) определяем значение окружной скорости потока на радиусе г:

W (г) = и (г) -о- г

И

.

Из выражения (12) очевидно, что при увеличении высоты кольцевого канала разгонного узла окружная скорость недиспергированной волокнистой сус-

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

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

Подставляя в (12) уравнение (9) для определения осевой скорости потока, получим

ж ( г ) =

1п

с л г

V г2 J

к

ю-г

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

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

Рис. 3. Эпюры скоростей волокнистой суспензии в разгонном узле дезинтегратора

В потоке суспензии на волокно действует сила сопротивления которая пропорциональна первой степени скорости поперечного смещения:

р = 21-т-В8-г,

й„

где т - масса волокна, кг; I - длина волокна, м; йц - диаметр волокна, м; г - поперечная координата волокна в

. йг

точке центра тяжести массы, г = —, м;

йг

Бц - коэффициент сопротивления боковой поверхности волокна во вращающемся потоке [4],

Б, =

4п • ц

I • 1п(7,4 Яе) где Яе - число Рейнольдса, = {V)ds Р

Яе

ц

где V - абсолютная скорость потока

суспензии, м/с, которую можно найти как

V = у1 и2

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

т •

й2г Н2 п йг

= т-^ - 2Б— (18) йг г йг

где О^ - центробежная скорость волок-

на, м/с, определяемая согласно выражению

О я =о-г,

где тв - угловая скорость вращения ориентированного волокна в кольцевом пространстве, рад/с; г - произвольный радиус вращения волокна в кольцевом канале, м.

Используя в (18) выражение (19), при этом сокращая (18) на массу волокна, получим

й 2г йг2

= оО •г - 2 Б

йг „ ч

>и

Аналитическое решение (20) ищем в виде

г = еп .

Подставляя (21) в (20), получим

у2 •е* + 2у-еп-о2 •еп = 0.

Алгебраическими корнями уравнения (22) будут действительные числа

= -Б ±^1Б2 + о2 .

Общее решение уравнения (22) будет представлено следующим распределением:

г

где С1,2 - константы интегрирования, определяемые из начальных условий:

йг Л

г = 0; г = г,; — = 0. 0 йг

Значит, из (24) при наличии начальных условий следует, что

-Ух+С2-У2 = 0; , С1 + С2 = г0.

Подставляя (24) в (25), определяем постоянные константы интегрирования:

с = го -

2 2 5 + Ю,

л/В

2 2 5 + Ю,

С = ^ В5 С2 2 -

2 +

В

2 . ,,2 5 + Ю,

На основании аналитических исследований (18)...(26) получим уравнение для определения радиальной скоро-

сти движения волокна во вращающемся потоке суспензии на любом радиусе кольцевого канала

V = Л - С - е»' + Г2 - с2 - е^.

После подстановки в (26) всех компонентов уравнения окончательно имеем

V = го-v 2

2

в2 1 г -

2 2 5 +Ю

В, - Ув,2 + ю,2

ТВ, е

2 2 2 + ю

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

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

О п-р-(р2 -р) х у 8/- кТ

2 2 / 4 4 \ г2 - 1 (г2-

2

1п

V г1 У

Осредненная скорость суспензии в трубе кольцевого сечения

_ 1

V =— V,

0'

где v0 - максимальная скорость в середине сечения кольцевого канала, м/с.

Тогда, с учетом (30), выражение (29) примет вид:

о=п-р-2

2 2 г2 - г1

г г ^

_2

V г1 У

1п

.

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

о, > иэ -Оэ,

где о, - расход суспензии на верти-

кальном кольцевом участке, кг/с; пэ -количество тороидально-эллиптических патрубков, шт.; Qэ - расход суспензии через тороидально-эллиптический патрубок, кг/с.

С помощью вероятного метода [5] для определения скорости V подберем уравнение и коэффициенты в этом решении так, чтобы удовлетворялось граничное условие на стенке трубы V = 0. Направим координатные оси У и 2 вдоль главных нормальных осей поперечного сечения тороидально-эллиптической трубы (см. рис. 3) и получим решение в виде

V = А-У2 + Б-22 + V,

0'

где А и В - геометрические коэффициенты.

Выражению (33) соответствует уравнение

Р2 - Р

2 А + 2 Б =—2-1

Ит-Ц

На внутренней поверхности тороидально-эллиптической трубы V = 0, значит,

А-У2 + Б-22 + v0 = 0.

Приведем (35) к уравнению эллиптического сечения трубы

жидкости на оси трубы эллиптического сечения

У0 =

Р2 - Р1 а2 - Ь2

2ИТ - ц а + Ь

2

Значит, выражение (39) на основании (36) примет вид:

V = V,,

т2 \

- у! - 2_

ч а2 Ь2,

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

Теперь определяем расход суспензии через трубу эллиптического сечения. Поверхности, на которых скорость V постоянна, - эллиптические цилиндры:

У1 21 = 1

а'2 + Ь'2 " '

полуоси которых определяются соотношениями:

а'2 = а2 ^; Ь'2 = Ь2 ^,

у2 гу1 -+-1 -0

а Ь

Тогда коэффициенты можно будет представить и вычислить как

А = -

а

Б = -2.

Ь2

Решая последовательно уравнения (36).(38), получаем выражение для определения скорости течения вязкой

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

Qэ =р1 v•йS = л - а-Ь

-Р-

IV У-йу.

После интегрирования (43) окончательно имеем

„ п-р-а-Ь ^ =---= V0.

Выражение (44) является частным случаем для расчета расхода тороидально-эллиптического патрубка и получено методом интегрирования уравнений движения суспензии внутри разгонного узла дезинтегратора.

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

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

Работа выполнена в рамках гранта РФФИ и Правительства Белгородской области, проект № 14-41-08054 р_офи_м.

СПИСОК ЛИТЕРАТУРЫ

1. Отходы горнодобывающих предприятий как сырье для производства мелкозернистого бетона, армированного фибрами / А. В. Клюев [и др.] // Вестн. БГТУ им. В. Г. Шухова. - 2010. - № 4. - С. 81-84.

2. Пат. 2377070 РФ, МПК7 В 02 С 13/22. Дезинтегратор / А. М. Гридчин, В. С. Севостьянов, Е. С. Глаголев, А. Е. Качаев ; заявитель и патентообладатель Белгород. гос. технолог. ун-т им. В. Г. Шухова. -№ 2008118079/03 ; заявл. 06.05.08 ; опубл. 27.12.09, Бюл. № 36. - 9 с.

3. Технологические модули и агрегаты для комплексной переработки природных и техногенных материалов / В. С. Севостьянов, Е. В. Синица, М. В. Севостьянов, А. В. Уральский, А. Е. Качаев, Д. Н. Перелыгин // Актуальные проблемы в строительстве и архитектуре. Образование. Наука. Практика : сб. докл. 66-й Всероссийской науч.-техн. конф. - Самара : Самар. гос. архитектурно-строит. ун-т, 2009. -С. 210-213.

4. Мазо, А. Б. Гидродинамика : учеб. пособие / А. Б. Мазо, К. А. Поташев. - Казань : КГУ, 2008. -

126 с.

5. Ерош, И. Л. Дискретная математика : учеб. пособие / И. Л. Ерош, М. Б. Сергеев, Н. В. Соловьев. -СПб. : СПбГУАП, 2005. - 144 с. : ил.

Статья сдана в редакцию 18 января 2016 года

Владимир Семенович Севостьянов, д-р техн. наук, проф., Белгородский государственный технологический университет им. В. Г. Шухова. E-mail: sevostynov_00@mail.ru.

Александр Евгеньевич Качаев, канд. техн. наук, гл. конструктор, ООО «Воскресенский завод «Машиностроитель».

Сергей Анатольевич Михайличенко, канд. техн. наук, проф., Белгородский государственный технологический университет им. В. Г. Шухова. E-mail: prorector@intbel.ru.

Татьяна Леонидовна Сиваченко, соискатель, Белгородский государственный технологический университет им. В. Г. Шухова. E-mail: tatsianamail.86@gmail.com.

Александр Александрович Фарафонов, аспирант, Белгородский государственный технологический университет им. В. Г. Шухова. Тел.: +7-4722-36-50-30.

Vladimir Semenovich Sevostyanov, DSc (Engineering), Prof., Belgorod State Technological University named after V. G. Shukhov. E-mail: sevostynov_00@mail.ru.

Aleksandr Yevgenyevich Kachayev, PhD (Engineering), Chief Designer of OOO Vosktesensky «Mashinostroitel» plant.

Sergey Anatolyevich Mikhailichenko, PhD (Engineering), Prof., Belgorod State Technological University named after V. G. Shukhov. E-mail: prorector@intbel.ru.

Tatiana Leonidovna Sivachenko, соискатель, Belgorod State Technological University named after V. G. Shukhov. E-mail: tatsianamail.86@gmail.com.

Aleksandr Aleksandrovich Farafonov, PhD student, Belgorod State Technological University named after V. G. Shukhov. Phone: +7-4722-36-50-30.