Исследование гидродинамического вибрационного фильтрования и разработка конструкции фильтра Текст научной статьи по специальности «Физика»

УДК 66.067 + 628.33 + 532.546 + 51-74

ИССЛЕДОВАНИЕ ГИДРОДИНАМИЧЕСКОГО ВИБРАЦИОННОГО ФИЛЬТРОВАНИЯ И РАЗРАБОТКА КОНСТРУКЦИИ ФИЛЬТРА

© 2012 В. А. Девисилов, И.А. Мягков, Е. Ю. Шарай

Московский государственный технический университет имени Н. Э. Баумана

Поступила в редакцию 06.10.2012.

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

Ключевые слова: гидродинамика, фильтрование, фильтр, неньютоновская жидкость, очистка жидкости.

1. ВВЕДЕНИЕ

В химической, нефтехимической, текстильной и пищевой промышленности остро стоит проблема очистки высоковязких жидкостей от механических примесей. Во многих случаях очищаемая жидкая дисперсная среда является неньютоновской жидкостью, эффективная вязкость которой зависит от интенсивности скоростей деформации. Особенностями фильтрования неньютоновской жидкости являются быстрое забивание пор твердыми частицами, зарастание пор вследствие облитерационных процессов, приводящее к снижению пропускной способности и ухудшению эксплуатационных свойств фильтрующего материала, что требует его частой замены. Высокая вязкость дисперсной среды приводит к необходимости создания большого перепада давления на фильтровальной перегородке (ФП). Для решения указанных проблем целесообразно использовать аппараты комбинированного действия, обеспечивающие высокую эффективность сепарационных процессов. Среди них можно выделить гидродинамический вибрационный фильтр (ГДВФ), позволяющий существенно сократить энергозатраты по сравнению с традиционным фильтрованием или аппаратами центробежного действия. ГДВФ отличается относительно простой конструкцией, надежностью и обеспечивает высокую степень очистки жидкости от твердых частиц при повышенном ресурсе работы [1]. Принцип гидродинамического фильтрования заключается в том,

Девисилов Владимир Аркадьевич, кандидат технических

наук, доцент. E-mail: devisil@mail.ru

Мягков Игорь Альбертович, ассистент.

E-mail: i.myagkov@bk.ru

Шарай Елена Юрьевна, аспирант.

E-mail: elena.sharai@yandex.ru

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

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

Таким образом, очистка жидкости в аппаратах на основе гидродинамического вибрационного фильтрования осуществляется за счёт сочетания следующих процессов:

• фильтрования через проницаемую структуру;

• центробежной сепарации;

• вибрационного разрушения осадка, образующегося на ФП;

• гидродинамического смыва осадка и удаления загрязнений из ГВФ.

Ко многим многофазным гетерогенным средам в широком диапазоне скоростей сдвига применимо реологическое уравнение состояния неньютоновской жидкости, записанное в виде степенной зависимости Оствальда-де Виля [3]. В работах [4, 5] с использованием компьютерного моделирования в программных комплексах Ansys CFX и Star-CD показано, что гидродинамическое фильтрование неньютоновской жидкости приводит к снижению нагрузки на фильтровальный материал, увеличению времени до регенерации фильтра, снижению перепада давления на ФП при её вибрации. Однако к настоящему времени отсутствует математическое описание гидродинамики процесса, происходя-

Рис. 1. Схема действующих сил при гидродинамическом вибрационном фильтровании

щего в гидродинамическом фильтре. Существует ограниченное количество работ, посвященных описанию отдельных механизмов, реализуемых в гидродинамическом фильтре. Большинство из них касается непосредственно фильтрования без учета влияния всех механизмов очистки или при других схемах движения потока жидкости [6, 7]. Следовательно, разработка аналитического аппарата, описывающего гидродинамические процессы в гидродинамическом фильтре для неньютоновской жидкости, подчиненной зависимости Оствальда-де Виля, является актуальной и практически важной задачей. Математическое описание гидродинамических процессов позволит подбирать основные режимные параметры фильтра для каждой жидкости при известных ее реологических характеристиках и параметрах загрязнений.

Данная работа посвящена созданию математической модели течения неньютоновской жидкости в гидродинамическом фильтре при вращающейся ФП с использованием уравнений На-вье-Стокса.

2. ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ

Рассматриваемая схема течения жидкости в гидродинамическом фильтре показана на рис. 2. Аппарат состоит из конического корпуса 2, в котором расположена цилиндрическая ФП 1. В кольцевой зазор сверху подается суспензия. Под действием перепада давления очищаемая жидкость фильтруется через перегородку 1, обладая окружной V, радиальной юг и продольной составляющими скорости. Часть жидкости перепускается вдоль перегородки и выходит через

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

Математическая модель процесса течения неньютоновской жидкости в конфузорной области фильтра представляет собой совокупность уравнений Навье-Стокса в компонентах тензора напряжений, уравнения неразрывности, реологического уравнения среды, условия равномерности расхода жидкости через ФП, граничных условий, записанных в цилиндрической системе координат.

Полагаем, что: распределение скоростей не зависит от угловой координаты, режим течения

Рис. 2. Схема течения в гидродинамическом фильтре: фильтровальная перегородка; 2 - корпус фильтра

1

ламинарный, жидкость несжимаема, ФП представлена эквивалентным гидравлическим сопротивлением.

К многофазным высоковязким средам применимо уравнение состояния неньютоновской жидкости, которое записывают в виде степенной зависимости Оствальда-де Виля:

где

dvr vr dvz

r- + — +—- = 0

r„

dr

kAn~l dv

dr

dz

(9)

t = k\n

или

% = k Hn_1r. (1)

Эффективная вязкость неньютоновской

жидкости определяется как

in—1

M = k\r\n .

(2)

В случае сложного деформирования степенной закон Оствальда-де Виля (1) представляется в виде:

Ъ = кЛ"-1уг]. (3)

В цилиндрической системе координат интенсивность тензора скоростей деформаций А записывается в виде

A=

t j+2Cvr j+#J+(dp+p #r+(!

(4)

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

^ = У<Р(г'7) ,

\г = \г {г, г), Уг = (г, г) .

Тогда формула (3) примет следующий вид:

% =k

'М+ДТ+/М +|4+vvl Jdv,

dr ) +1r J +\ dz J +1 dr + r I +1 dz I +! dr + '

(5)

Уравнения Навье-Стокса и неразрывности в компонентах тензора деформации для рассматриваемого случая имеют вид:

vdvL + v dvL — vl =—\_ dp + 1 f d%rr +d%rz +%rr —%ррЛ

r dr z dz r pdr p

dr dz r

(6)

—= — 1 dp +1 f % +dTzL +

dv dv vm ,

z rz zz rz

vr — + vz . .

dr dz r p dz p\ dr dz r

(7)

dvp vrvp "'p vr —— +-— + v

dr r

dvp 1 f d%rp T %rp

dz p

dr dz

%pp = kA

n—1

1 dvp

r + _ p

r dp

% = kAn— dv

dz

(10)

kA -

1 1 f dvr dv

+ ■ dz dr

rp

kA -

kAn—1 -

1 1 f dvV

1 dvr

dr

■ + —

r dp

2

dv

z + dvp Л

dp dz

-- у

Подставив в систему уравнений (6)-(9) значения компонентов тензора напряжений ъ, определенных согласно соотношениям (10), получим систему уравнений гидродинамики неньютоновской псевдопластической жидкости. Уравнения в конечном виде не приводятся в силу их громоздкой записи.

Система уравнений (6)-(9) имеет следующие граничные условия:

Уг = 0 при г = г у, Г = Я - ^(а) ;

уг = 0 при Г = Я - ztg(a) ;

Уг = Уу при г = г у;

^ = 0 при г = Я - ztg(a) ; (11)

= гу со при г = гу ; Р = РЯ (г) при 7 = 0 ; Р = Рь (г) пРи г = Ь.

3. РЕЗУЛЬТАТЫ РАСЧЕТА

Для нахождения гидродинамических характеристик потока необходимо численно решить систему дифференциальных уравнений в частных производных (6)-(9) при заданных граничных условиях (11) для соответствующей геометрии корпуса фильтра. Решение системы дифференциальных уравнений реализовано на ПК с помощью программного комплекса Maple как последовательная итерационная процедура. Численное решение выполнялось для неньютоновской псевдопластичной жидкости в интервале показателя степени n от 0.6 до 1.0 при угле конусности а= 11° и соотношениях геометричес-

r

2

v

p

2

r

2

2

r

ких размеров фильтра R/L = 0.45, r/L = 0.20 (рис. 1). В процессе расчета определялись давление p, окружная V, радиальная V и продольная г скорости потока в расчетной области.

На рис. 3а, 3б и 3в приведены зависимости безразмерных окружной , радиальной ут и продольной V составляющих скорости течения неньютоновской жидкости от безразмерной радиальной координаты г для различных значений показателя степени п при z/L = 0.75. Соотношения между размерными и безразмерными величинами следующие:

<Р -

ОТ,

!

Дп=1

И

1.0

0.0

1.0

1.1

1.2

1.3

1.4

1.5

Рис. 4. Зависимости безразмерной эффективной вязкости Д неньютоновской жидкости от радиальной координаты для различных показателей степени п (1 - 0.6; 2 - 0.8; 5 - 1.0) при z/L = 0.75: линии - расчет по уравнениям; маркеры - расчет в Ли8у8 СРХ

г/ уо Уо '/

Установлено, что окружная скорость V практически не зависит от индекса течения п (рис. 3а). Радиальная и продольная составляющие скорости уменьшаются с увеличением индекса течения п (рис. 3б, 3в).

На рис. 4 (кривые 1-3) приведены зависимости безразмерной эффективной вязкости Д неньютоновской жидкости от безразмерной радиальной координаты г для различных значений показателя степени при z/L = 0.75. Безразмерная эффективная вязкость определяется как:

0.0 а)

Уг

Из анализа полученных расчетных зависимостей рис. 4, следует, что при уменьшении показателя степени п неньютоновской псевдопластической жидкости эффективная вязкость уменьшается.

Для определения адекватности полученных результатов проведен расчет в программном

0.0 б)

в)

1.1 1.2

1.4

1.5 г

1.1 1.2

1.3

1.4

1.5 г

1.1

1.2

1.3

1.4

1.5

7

Рис. 3. Зависимости безразмерных составляющих скорости течения неньютоновской жидкости от радиальной координаты т для различных показателей степени п (1 - 0.6; 2 - 0.8; 3 - 1.0) при z/L = 0.75: а - окружная составляющая ; б - радиальная составляющая ; в - продольная составляющая \2

Г

V

г ^

vz =

комплексе Ли8у8 СРХ для рассматриваемой модели при аналогичных граничных условиях. На рис. 4 квадратными маркерами показаны значения эффективной вязкости при п = 0.8, определенные в модуле Ли8у8 СРХ. Как видно из рисунка, получено хорошее согласование результатов при решении двумя методами.

При решении системы уравнений (6)-(9)

V <р/ у?ншх 1.0

а) 1.0

1.4

б)

1.4

1.8 Г

0.0 в)

1.4

1.

Рис. 5. Зависимости безразмерных составляющих скорости течения неньютоновской жидкости от радиальной координаты для различных чисел Рейнольдса Иеп при z/L = 0.75 (1 - 2.8; 2 - 6.5; 3 - 10.5) и при z/L = 0.50 (4 - 2.8; 5 - 6.5; 6 - 10.5): а - окружная составляющая ; б - радиальная составляющая ; в - продольная составляющая

варьировался расход жидкости через фильтр, при этом обобщенное число Рейнольдса Иеп, определенное во входном сечении при z = 0, изменялось от 2.8 до 10.5. На рис. 5 представлены зависимости окружной V радиальной V и продольной V составляющих скорости течения неньютоновской жидкости от радиальной координаты г в безразмерных осях при z/L = 0.75 (кривые 1-3) и z/L = 0.50 (кривые 4-6) для Иеп = 2.8 (кривые 1, 4), Иеп = 6.5 (кривые 2, 5), Иеп = 10.5 (кривые 3, 6).

Установлено, что максимум радиальной (рис. 5б) и продольной (рис. 5в) составляющей скорости смещается в сторону ФП, что объясняется отсосом жидкости через проницаемую перегородку. С увеличением критерия Иеп кривая окружной составляющей скорости (рис. 5а) приобретает более пологий вид.

На рис. 6 показаны графики безразмерной эффективной вязкости для критерия Рейнольдса Иеп = 2.8 (кривые 1, 4), Иеп = 6.5 (кривые 2, 5), Иеп = 10.5 (кривые 3, 6) при z/L = 0.75 (кривые 13) и z/L = 0.50 (кривые 4-6). Из анализа полученных кривых следует, что с увеличением обобщенного числа Рейнольдса эффективная вязкость увеличивается в кольцевой области фильтра, но при этом вблизи ФП для всех кривых ее значение одинаково, что объясняется равенством радиальной скорости для всех рассматриваемых случаев при г = 1.

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

0.7 п

0.6

0.5

0.4

1.0

1.4

1.?

Рис. 6. Зависимости безразмерной эффективной вязкости неньютоновской жидкости от радиальной координаты для различных чисел Рейнольдса Иеп при z/L = 0.75 (1 - 8.5; 2 - 22.4; 3 - 30.6) и при z/L = 0.50 (4 - 8.5; 5 - 22.4; 6 - 30.6)

м 1.2 п

0.3

0.0 4-Т-1-1-1

1.0 1.2 1.4 1.6 1.8

Г

Рис. 7. Зависимости безразмерной эффективной вязкости неньютоновской жидкости от радиальной координаты для различных угловых скоростей вращения щ при z/L = 0.50 (1 - 12 рад/с; 2 - 24 рад/с; 3 - вращение отсутствует) при п = 0.8 и z/L = 0.75

о ФП 12 рад/с (рис. 6, кривая 1), 24 рад/с (рис. 7, кривая 2) и при неподвижной ФП (рис. 7, кривая 3) для индекса течения п = 0.8 и z/L = 0.50.

Кривая 3 рис. 7 имеет экстремум (максимум), соответствующий максимальному значению скорости в данном сечении кольцевой области, за счет увеличения эффективной вязкости в зонах с малыми скоростями сдвига.

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

На основе проведенных расчетов окружной, радиальной и продольной составляющих скоростей жидкости по уравнениям (6)-(9) можно получить линии тока в кольцевом зазоре между корпусом фильтра и фильтровальной перегородкой. На рис. 8 показаны линии тока при обобщенном числе Рейнольдса Иеп = 2.8, определенном во входном сечении фильтра при п = 0.8 и о= 12 рад/с, для описанной выше геометрии. В программном комплексе Ли8у8 СРХ получены траектории твердых частиц при аналогичных граничных условиях. На рис. 7 изображены линия тока (поз. 3), построенные на основе расчета скоростей по уравнениям (6)-(9) и промежуточные положения частиц (поз. 4) диаметром 8

мкм при р/р = 0.4, полученные в программном комплексе Ли8у8 СРХ, которые хорошо накладываются на линии тока, полученные из численного решения уравнений (6)-(9). Анализ многочисленных расчетов показал возможность считать с приемлемой достоверностью, что траектории движения частиц размером до 20 мкм совпадают с линией тока жидкой фазы.

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

4. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ

Для проведения экспериментальных исследований создан стенд на базе модельного ГДВФ. Основные конструктивные узлы ЭС, методика измерений, результаты экспериментов детально представлены в [1].

Рис. 8. Картина течения в гидродинамическом фильтре: 1 - фильтровальная перегородка; 2 - корпус фильтра; 3 - линии тока неньютоновской жидкости при п = 0.8, полученные при численном решении системы уравнений (13)-(15); 4 - траектории частиц, полученные при расчете в Лшуэ СРХ

Пневмогидравлическая схема ЭС представлена на рис. 9. Согласно схеме очищаемая жидкость из исходной емкости 22 шестеренчатым насосом 11 подается в систему ГВФ 8, где очищается от механических примесей. Вращение ФП ГВФ обеспечивается электродвигателем 21 через клиноременную передачу; вибрация ФП - электродвигателем 6 с кулисным механизмом 7. Все электродвигатели снабжены датчиками числа оборотов.

Работа ЭС возможна как по замкнутому контуру, так и по разомкнутому При работе по замкнутому контуру очищенная в ГВФ 8 жидкость из приемной емкости 28 шестеренчатым насосом 23 подается в исходную емкость 22 для дальнейшей рециркуляции. В результате многократной рециркуляции через ГВФ 8 очищаемой жидкости механические примеси откладываются в емкости 29 с концентратом загрязнений, а в исходной емкости 22 образуется очищенная жидкость. Количество проходов жидкости при рециркуляции через ГВФ определяется требуемой концентрацией механических примесей в очищенной жидкости. Необходимо также отметить, что в системе ЭС установлен волокнистый фильтр 26, который позволяет при необходимости уже сразу после первого прохода уда-

лить механические примеси, которые остались после ГВФ 8. Задание одного из двух вышесказанных способов очистки достигается переключением трехходового клапана 27. При работе по разомкнутому контуру очищенная в ГВФ 8 жидкость выводится из системы через трехходовой клапан 24.

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

Результаты экспериментального исследования гидродинамического вибрационного фильтра с хорошей степенью совпадения (15-20 %) подтвердили спрогнозированные теоретическими расчетами и машинным моделированием се-парационные и гидравлические параметры

Рис. 9. Пневмогидравлическая схема экспериментального стенда: 1 - манометр образцовый; 2 - клапан трехходовой; 3 - расходомер индукционный; 4 - клапан обратный; 5 -насос гидроструйный; 6 - генератор вибрации ГДВФ (электродвигатель); 7 - механизм кулисный; 8 - ГДВФ; 9 - манометр образцовый; 10 - электродвигатель шестеренчатого насоса; 11 - насос шестеренчатый; 12 -расходомер турбинный; 13 - клапан балансировочный; 14 - расходомер переменного перепада давления; 15-клапан балансировочный; 16 - расходомер переменного перепада давления; 17 - электродвигатель шестеренчатого насоса; 18 - клапан балансировочный; 19 - компрессор; 20 - емкость с раствором загрязнений; 21 - привод вращения ГДВФ (электродвигатель); 22 - емкость исходная; 23 - насос шестеренчатый; 24 - клапан трехходовой; 25 - клапан обратный; 26 - фильтр волокнистый; 27 - клапан трехходовой; 28 - емкость приемная; 29 - емкость с концентратом загрязнений

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

Были проведены экспериментальные исследования по определению концентраций частиц на выходе из кольцевой области на модельных неньютоновских средах, представляющих собой водные растворы сополимера акриловой кислоты и алли-ловых эфиров пентаэритрита (САКАПа). Отличие расчетных и экспериментальных значений по количеству частиц с медианным диаметром 15 мкм и среднеквадратическим отклонением 0.6 при р/ рч = 0.4, присутствующих на выходе из фильтра в перепускаемом потоке, составляло не более 15 %. Эти исследования подтверждают основные выводы, сделанные в работе.

Некоторые более детальные результаты представлены в работе [8], в расширенном варианте опытные данные будут опубликованы позднее.

Отвод

очищенной жидкости

5. РАЗРАБОТКА ОПЫТНОГО ОБРАЗЦА ФИЛЬТРА

Результаты теоретических исследований, аналитических расчетов и экспериментов позволили выбрать оптимальную конструкцию фильтра. Информация о методике расчета конструктивных и режимных параметров носит характер "ноу-хау" и может быть предоставлена заинтересованным организациям на договорных основах. На рис. 10, 11 показана конструкция опытного образца ГДВФ с системами привода, подачи жидкости и узлом генерирования вибрации.

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

Узел генерирования колебаний в ГДВФ

Подвод загрязненной жидкости

16_ 17

Рис. 10. Опытный образец ГДВФ: 1 - корпус сварной; 2 - перегородка фильтровальная из ПСМ; 3 - штуцер ФП зажимной; 4 - патрубок для отвода концентрата; 5 - штуцер зажимной (толкатель); 6 - кольцо разбрызгивающее с лопатками; 7 - вал торсионный; 8 - стакан корпусной; 9 - крышка корпусная; 10 - винт установочный; 11 - патрубок выходной; 12 - гайка шлицевая регулировочная с лапчатой шайбой; 13 - кольцо вибрационное; 14 - втулка скольжения биметаллическая; 15 - патрубок входной; 16 - втулка скольжения биметаллическая; 17 - крышка корпусная; 18 - винт установочный; 19 - узел генерирования колебаний в ГДВФ; 20 - шпилька регулировочная; 21 -манжета армированная; 22 - гайка корончатая регулировочная со шплинтом; 23 - кольцо опорное; 24 - подшипник шариковый упорный; 25 - гайка шлицевая регулировочная с лапчатой шайбой; 26 - пружина сжатия; 27 - гайка корончатая регулировочная со шплинтом; 28 - шпилька направляющая; 29 - выступ (зуб) толкателя; 30 - выступ дорожки вибрационного кольца

Рис. 11. Схема компоновки опытного образца ГДВФ: 1 - струбцина; 2 - столешница; 3 - электродвигатель ГДВФ; 4 - муфта ГДВФ с крестовым элементом; 5 -рама ГДВФ сварная; 6 - генератор колебаний ГДВФ; 7 - ГДВФ; 8 - патрубок ГДВФ для отвода концентрата загрязнений; 9 - патрубок ГДВФ выходной; 10 - трубопровод соединительный; 11 - патрубок насоса входной; 12 - насос; 13 - патрубок ГДВФ входной; 14 - патрубок насоса выходной; 15 - муфта насоса с крестовым элементом; 16 - рама насоса сварная; 17 - электродвигатель насоса; 18 - столешница; 19 - струбцина; 20 -опоры вибрационные для горизонтального монтажа ГДВФ; 21 - опоры вибрационные для вертикального монтажа ГДВФ; 22 - винт установочный для стыкового крепления блоков

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

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

моторных и трансформаторных масел. Нами разработаны конструкции всех аппаратов установки (ГДВФ, аппарат для удаления легких фракций, дополнительной очистки и, при необходимости, введения присадок). На рис. 12 представлена компоновка установки на автомобиле.

Рис. 12. Компоновка установки по регенерации отработанных масел в кузове автомобиля

6. ЗАКЛЮЧЕНИЕ

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

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

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

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

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

ность их использования для практических расчетов гидродинамических и сепарационных процессов в гидродинамическом фильтре.

6. Разработана конструкция опытного образца фильтра, готовая к испытаниям в процессе эксплуатации.

Работа выполнена в рамках ФЦП "Научные и научно-педагогические кадры инновационной России"на 2009 - 2013 годы.

Обозначения

интенсивность тензора скоростей деформаций, с-1;

£ - индекс консистентности, Па сп; Ь - длина корпуса фильтра, м; п - степень нелинейности кривой течения (показатель неньютоновского поведения, п < 1;

рк(г) - функция распределения давления в радиальном сечении при 2 = 0, Па;

рь(г) - функция распределения давления в радиальном сечении при 2 = Ь, Па; р - давление, Па;

Я - радиус конической части фильтра при 2 = 0, м; Яеч - эквивалентный гидравлический радиус входного сечения фильтра, м;

гу - радиус фильтровальной перегородки, м; г - радиальная цилиндрическая координата, м; V у - средняя скорость фильтрования, м с-1;

средняя скорость во входном сечении

фильтра, м с-1;

0

V ^ V г, Vг - окружная, радиальная и продольная составляющие скорости в цилиндрической системе координат, м с-1;

V рпах - максимальная окружная составляющая скорости, соответствующая Иеп = 10.5, определенном в сечении при г/Ь = 0.75;

V пып - минимальная (в силу отрицательного знака) радиальная составляющая скорости, соответствующая Иеп = 10.5, определенном в сечении при г/Ь = 0.75;

V япах - максимальная продольная составляющая скорости, соответствующая Иеп = 10.5, определенном в сечении при г/Ь = 0.75;

Г - безразмерная радиальная координата;

Vф - безразмерная радиальная составляющая скорости;

Уг - безразмерная радиальная составляющая скорости;

У2 - безразмерная радиальная составляющая скорости;

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

а- угол конусности корпуса фильтра, град;

Уу - компоненты тензора скоростей деформации, с-1;

ц-эффективная вязкость неньютоновской жидкости, Па с;

Цп=1 - вязкость жидкости при п = 1, Па с;

Ц - безразмерная эффективная вязкость неньютоновской жидкости;

р- плотность жидкости, кг м-3;

рч - плотность частиц, кг м-3;

ТТ - компоненты тензора напряжений, Па;

со- угловая скорость вращения, рад с-1;

Иеп = рЯе^ V 02-п/к - обобщенное число Рейнольдса.

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

1. Девисилов В.А., Мягков И.А. Мобильная установка для восстановления отработанных масел / / Безопасность в техносфере. 2007. № 5. С. 36.

2. Мочалин Е.В., Мочалина И.Г. Эффективность сепарации взвешенных примесей вращающимся фильтроэ-лементом // Вестник НТУ "ХПИ". 2011. № 10. С. 3.

3. Романков П.Г. Гидродинамические процессы химической технологии // Теорет. основы хим. технологии. 1972. Т. 6. № 6. С. 855.

4. Девисилов В.А., Шарай Е.Ю. Моделирование работы гидродинамического фильтра с использованием программного комплекса // Безопасность в техносфере. 2009. № 5. С. 21.

5. Девисилов В.А., Шарай Е.Ю. Моделирование течения неньютоновской жидкости вблизи вибрирующей фильтровальной перегородки гидродинамического фильтра // Безопасность в техносфере. 2010. № 5. С. 23.

6. Холпанов Л.П., Ибятов Р.И. Математическое моделирование гидродинамики на проницаемых поверхностях // Теорет. основы хим. технологии. 2003. Т. 37. № 3. С. 227

7. Холпанов Л.П., Ибятов Р.И. Математическое моделирование динамики дисперсной фазы // Теорет. основы хим. технологии. 2005. Т. 39. № 2. С. 190.

8. Девисилов В.А., Мягков И.А. Экспериментальное исследование гидродинамического вибрационного фильтра // Химическое и нефтегазовое машиностроение. 2011. № 8. С. 13 -17.

INVESTIGATION OF HYDRODYNAMIC VIBRATIONAL FILTERING AND DEVELOPMENT OF CONSTRUCTION OF A FILTER

© 2012 V.A. Devisilov, I.A. Myagkov, E.Yu. Sharay

Moscow State Technical University named after N.E. Bauman

Mathematical model of current of non-newton liquid in hydrodynamic filter during rotation of filtering partition have been developed. Using numerical methods system of equations of reodynamics have been decided. The sample of filter have been created.

Key words: hydrodynamics, filtering, filter, non-newton liquid, purification of liquid

Vladimir Devisilov, Candidate of Technics, Associate Professor. E-mail: devisil@mail.ru Igor Myagkov, Assistant. E-mail: i.myagkov@bk.ru Elena Sharay, Postgraduate. E-mail: elena.sharai@yandex.ru