CC BY
22
DOI 10.23859/1994-0637-2017-1-76-6 УДК 542.63:621.929
Кожевников Сергей Олегович
Кандидат технических наук, доцент, Ивановский государственный политехнический университет (Иваново, Россия) E-mail: [email protected]
Калинин Евгений Николаевич
Доктор технических наук, профессор, Ивановский государственный политехнический университет (Иваново, Россия) E-mail: [email protected]
К ПОСТАНОВКЕ ЗАДАЧИ МОДЕЛИРОВАНИЯ ДВИЖЕНИЯ ЖИДКОЙ СРЕДЫ В ПЕРЕМЕШИВАЮЩЕМ АППАРАТЕ
Аннотация. В статье представлены результаты моделирования жидкой среды в аппарате, оснащенном вращающейся мешалкой с образованием воронки. Разработана математическая модель движения идеальной жидкости в аппарате, оснащенном мешалкой, равномерно вращающейся относительно вертикальной оси рабочей камеры.
В результате модельного эксперимента средствами вычислительного комплекса MathCaD в форме численной модели получены зависимости составляющих скоростей движения жидкости в аппарате и получены профили радиальной, окружной и осевой составляющей абсолютной скорости потока жидкости. Получены функциональные зависимости, дающие возможность оптимизировать скоростной режим движения жидкости в аппарате с вращающейся мешалкой и решить задачи по геометрическому и динамическому синтезу аппарата.
Ключевые слова. Компьютерное моделирование, вращающаяся мешалка, образование воронки.
© Кожевников С.О., Калинин Е.Н., 2017
Kozhevnikov Sergey Olegovich
PhD (Technology), Associate Professor, Ivanovo state Polytechnic University (Ivanovo, Russia) E-mail: [email protected]
Kalinin Evgeniy Nikolayevich
Doctor of technical Sciences, Professor, Ivanovo state Polytechnic University (Ivanovo, Russia) E-mail: [email protected]
ABOUT THE TASK FOR SIMULATION OF LIQUID MOTION IN THE KNEADING MACHINE
Abstract. The article presents the results of the fluid motion simulation in the apparatus equipped with a rotary stirrer with formation of funnel. A mathematical model of the motion of an ideal fluid in an apparatus equipped with a stirrer, uniformly rotating about a vertical axis of the working chamber is developed.
As a result of the model experiment by means of MathCaD computing system, correspondence between components of liquid velocities in the machine and received pro-profiles of the radial, axial and circumferential component of the absolute velocity of the liquid flow were obtained in the form of numerical-term model. The functional relationship provided in the paper enable to optimize the speed of fluid movement in the apparatus equipped with a rotary stirrer agitator and solve problems on the geometric and dynamic synthesis apparatus.
Key words. Computer modeling, rotary stirrer, funnelling.
Введение
Выполненный анализ конструкций перемешивающих устройств, предназначенных для получения материалов высокой дисперсности, показал, что они требуют
дальнейшего совершенствования [2], [3]. Применяемые для этой цели устройства имеют высокие показатели потребляемой мощности [2].
Наши исследования [2] показали, что для получения высокодисперсных суспензий можно использовать аппарат, оснащенный рабочим органом, состоящим из ротора с трубками переменного сечения. При прохождении жидкости через трубки переменного сечения, которые в свою очередь совершают вращательно движение, скорость среды на выходе многократно возрастает. Моделирование процесса движения жидкой среды позволит определить распределение скоростей движения жидкости в аппарате, выполнить его конструктивную, режимную оптимизацию и представляет особый интерес.
Основная часть
Перспективным направлением развития химико-технологических процессов является получение высокодисперсных материалов. Применяемые в промышленности конструкции перемешивающих устройств для получения материалов высокой дисперсности характеризуются небольшим разнообразием рабочих органов. Анализ способов получения материалов высокой дисперсности и конструктивного оформления процесса показал, что существующие перемешивающие устройства имеют высокую потребляемую мощность и требуют дальнейшего совершенствования рабочих органов смесителей [2], [3].
Собственные исследования [2], а также исследования других авторов [1] показали, что для получения высокодисперсных материалов наилучшим является применение вращающихся мешалок с дополнительным прохождением жидкости через трубки переменного сечения. Особый интерес представляет моделирование распределения скорости движения жидкости в аппарате с мешалкой. Моделирование процесса движения жидкой среды в аппарате позволит определить оптимальный скоростной режим движения жидкости, предотвратить образование застойных зон и решить задачи его геометрического и динамического синтеза.
Представим движение жидкости в аппарате с глубиной погружения Н мешалки (рис. 1). При этом движение жидкости симметрично относительно оси вращения мешалки.
Рис. 1. Расчетная схема: 1 - аппарат; 2 - мешалка
Для упрощения анализа влияние вязкости среды не учитывается. Свободную поверхность будем полагать плоской горизонтальной, пренебрегая ее понижением у оси.
В общем случае, когда на оси потока окружные скорости достигают бесконечно больших значений и в реальном потоке образуется воронка, задача приводится к неоднородному дифференциальному уравнению [4]:
+ к2у = -кс, (1)
д2 дг ( г дг )
где у - функция тока жидкости, м; г - текущий радиус, м; 2 - координата по высоте аппарата, м; к - постоянная, равная удвоенному отношению модуля угловой скорости вращения частицы к модулю скорости ее перемещения; С - постоянная.
Начало цилиндрических координат совместим с вертикальной осью мешалки и с горизонтальной осью центров лопаток. Граничные условия приняты в соответствии с рис. 1:
у(0, г) = у „, у(Н, г) = 0,
у( 2,0) = 0,
где постоянное значение функции тока у0 связано с насосным эффектом, создаваемым мешалкой, зависимостью
V *
Р
У0 = 2% • 2%
где V* - насосный эффект мешалки, м3/с, определяемый по формуле [4]:
VР = J ■ п ■ ёъм,
где J - коэффициент, зависящий от типа мешалки [4]; п - частота вращения мешалки, с-1; - диаметр мешалки, м.
Решим уравнение (1) при помощи рядов Фурье. Для этого введем вместо функции тока жидкости функцию
и = У-У0 (1 , (2)
С учетом (2) уравнение (1) примет вид
д2и д (1 ди\ 72гт ,/г к2у0
ТГ + г-дГ\-~д~ 1 + к и = -к(ку0 + С)+ ^2 , (3)
д2 дг (г дг ) Н
с граничными условиями и(0,г) = 0; и(Н,г) = 0; и(2,0) = -у0(1 -(2/Н}).
Решение уравнения (3), удовлетворяющее указанным граничным условиям, можно представить в виде ряда Фурье. Функция тока винтового потока
1 V-X1
Н ) к П=1 п
1+
с
1 - (-1)п—
ПК| -1
кН,
1+
1+-
с
1 - (-1)п— к^о
кН
-1
О,
Н
„ /(д. л
где £ 2п - постоянная, при-> 1:
п ^ кН
= ^т I -1 ■к ■г ■ К ■
кН
\
-1 ■ к ■ г
пк
при-= 1:
кН
^п = 1;
пк
при-< 1:
кН
О =-К 1 .к .г-7, ■
п 2\ I кН 1 1
п-, пк|кг
Л
где 71 (х) - функция Бесселя второго рода; К1 (х) - функция Бесселя мнимого аргумента.
При практических расчетах значения аргумента бесселевых функций х оказываются обычно во втором и последующих членах ряда столь большими, что можно воспользоваться асимптотическими представлениями этой функции при х^-да [4]
7 (х > "¿Х; к(х (Ч2Г
Окружная скорость определяется выражением
к у + С
уф="
где С - коэффициент, зависящий от скорости движения жидкости в аппарате, м /с. Коэффициент С определяется из уравнения
С = Ур/ г.
Радиальная скорость винтового потока
г
V
г Нг
2Е
1 +
С
1 - (-1)п — ку0
кН
-1
1 +
1 +-
С
1 - (-1)п — ку0
О.
кН
-1
П%2 ,
еоэ--1
Н
Осевая скорость потока
2к2у0 1
V2 =--^
% п=1 П
1 +
1 +-
С
1 - (-1)" —
кУ0
П"%\ -1
кН,
/ . "%2
II 81П-
" Н
"% 1 гл/
где при-> 1 постоянная Ип:
кН
О" =
п%
при-= 1:
кН
/ ч2 1
( п% \ -1 ( ■ К ■
оп =(
кН
Л
-1 ■ к ■ г
п%
при-< 1:
кН
П/ =-%1 -(^ ■ у ■ [ 1 -Ы' ■ к ■ г .
" 2 \ кН) 1 V (кН)
( )
где п - нечетное число, п =1,3,5....
Расчет параметров винтового потока жидкости в аппарате с перемешивающим устройством без учета физических свойств движущейся среды реализован в среде МаШСЛО в соответствии с приведенными уравнениями при следующих входных параметрах: ём = 0,05 м, J = 10, С = 0,2 м2/с, Н = 0,11 м, г = 0,14 м.
Распределение скоростей в аппарате с мешалкой показано на рис. 2, а профиль образующейся воронки - на рис. 3.
Рис. 2. Профили скоростей жидкости в аппарате с мешалкой
Рис. 3. 3-D модель профиля воронки жидкости в аппарате с мешалкой: 1 - аппарат; 2 - воронка жидкости; 3 - мешалка
Выводы
1. Разработана математическая модель движения идеальной жидкости в аппарате, оснащенном мешалкой, равномерно вращающейся относительно вертикальной оси рабочей камеры, реализованная средствами вычислительного комплекса MathCAD в форме численной модели с возможностью анализа кинематических параметров жидкой среды.
2. В результате модельного эксперимента получены зависимости составляющих скоростей движения жидкости в аппарате при входном значении скорости вращения n = 10 c-1 и получены профили радиальной, окружной и осевой составляющей абсолютной скорости потока жидкости.
3. Полученные функциональные зависимости дают возможность определить оптимальный скоростной режим движения жидкости в аппарате с вращающейся мешалкой и, как следствие, решить задачи по геометрическому и динамическому синтезу аппарата.
Литература
1. Васильев О.Ф. Основы механики винтовых и циркуляционных потоков. М. ; Л., 1958. 144 с.
2. Кожевников С.О., Гуюмджян П.П. Методы диспергирования жидких сред // Информационная среда вуза: Материалы XIII Междунар. науч.-техн. конф. Иваново, 2008. С. 865-868.
3. Кожевников С.О. Обзор основных типов перемешивающих устройств для получения технологических жидкостей // Вестник научно-промышленного общества. Вып. 6. М., 2003. С. 67-69.
4. Кожевников С.О. Определение конструктивных параметров перемешивающего устройства непрерывного действия // Вестник научно промышленного общества. Вып. 19. М., 2012. С. 8-12.
References
1. Vasil'ev O.F. Osnovy mekhaniki vintovykh i tsirkuliatsionnykh potokov [Fundamentals of mechanics of screw and flow patterns]. Mostow ; Leningrad, 1958. 144 p.
2. Kozhevnikov S.O., Guiumdzhian P.P. Metody dispergirovaniia zhidkikh sred [Methods of dispersing liquid media]. Informatsionnaia sreda vuza: Materialy KhIII Mezhdunar. nauch.-tekhn. konf. [University Infomedia: Proceedings of VIII International conference. scientific and engineering. Conf.]. Ivanovo, 2008, pp. 865-868.
3. Kozhevnikov S.O. Obzor osnovnykh tipov peremeshivaiushchikh ustroistv dlia polucheniia tekhnologicheskikh zhidkostei [Review of the main types of mixing devices for processing liquids].
Vestnik nauchno-promyshlennogo obshchestva. Vyp. 6. [Bulletin of the scientific and industrial community], 2003, Iss. 6, pp. 67-69.
4. Kozhevnikov S.O. Opredelenie konstruktivnykh parametrov peremeshivaiushchego ust-roistva nepreryvnogo deistviia [Determination of structural parameters continuous mixing apparatus]. Vestnik nauchno promyshlennogo obshchestva. Vyp. 19 [Journal of Industrial Research Associations], 2012, Iss. 19, pp. 8-12.
Кожевников С.О., Калинин Е.Н. К постановке задачи моделирования движения жидкой среды в перемешивающем аппарате // Вестник Череповецкого государственного университета. 2017. №1. С. 40-47.
For citation: Kozhevnikov S.O., Kalinin E.N. About the task for simulation of liquid motion in the kneading machine. Bulletin of the Cherepovets State University, 2017, no. 1, pp. 40-47.
