CC BY
96
УДК 66.021.1:532.62
Г. Н. Зиннатуллина, А. В. Давыдов, Р. Р. Байбеков, Н. Х. Зиннатуллин, В. В. Бронская
РАСЧЕТ ОСНОВНЫХ ТЕХНОЛОГИЧЕСКИХ ПАРАМЕТРОВ
ЦЕНТРОБЕЖНЫХ ПЛЕНОЧНЫХ АППАРАТОВ
Ключевые слова: реологические среды, центробежная пленка, толщина пленки, поле скоростей, давление, криволинейные
насадки.
В статье приведены результаты анализа гидродинамики центробежных жидких пленок, текущих по открытой поверхности насадок различной формы. Рассмотрены поведения вязких, аномально вязких, нелинейно упруго-вязких и вязко-пластичных жидких пленок.
Keywords: rheological liquidfilm thickness, centrifugal force, velocity field, pressure, curvilinear nozzle.
The results of the analysis of centrifugal hydrodynamics of liquid films flowing through the open surface of nozzles of different shapes in article. Viscous abnormally viscous, nonlinear viscoelastic and viscoplastic liquid film behavior was studied.
Для повешения удельной производительности технологических аппаратов давольно часто используется подвод внешней энергии в виде центробежной. Центробежные аппараты широко применяются в промышленности для проведения разнообразных гидромеханических, тепло- и массообмен-ных процессов [1-4]. В той работе рассмотрим случай, когда в аппарате реализуется тонкая центробежная ламинарная пленка со свободной поверхностью. Ротор этих аппаратов, состоящий из вала и рабочей насадки, изготавливается в виде плоского диска, конуса, сферы или иной криволинейной чаши или их комбинаций. Обрабатываемая жидкость, характеризующиеся различными реологическими свойствами, обычно подается в центр рабочего элемента - насадки или необходимый расход обеспечивается погружной частью насадок в жидкость. Рассмотрим законы сохранения импульса в центробежной жидкой пленке.
Вязкие (ньютоновские) жидкости
Изучению пленочного течения вязкой жидкости по открытой гладкой поверхности вращающегося плоского диска и конуса посвящен ряд теоретических и экспериментальных работ. В работе [5] в одномерной постановке определены среднее значение меридиональной скорости и толщины пленки жидкости. Аналогичная задача для плоского диска была решена в работе [6]. Влияние давления и поверхностного натяжения на процесс течения жидкости по конусу было проанализировано в работе [7]. При больших расходах маловязкой жидкости по поверхности ротора возникает значительное ее отставание. Этим вызвано появление работ в двумерной постановке. В работе [8] в приближенной постановке решена двумерная задача для плоского диска. Аналогичная задача решена в работе [9], предложенные в ней зависимости справедливы для больших радиусов. Наиболее полное изучение течения вязкой жидкости по поверхности плоского диска и конуса было выполнено К.В. Вачагиным [10]. Им учтены, кроме центробежной силы и силы вязкостного трения, также силы Кориолиса и давление. Получены выражения для меридиональной
скорости, скорости отстаивания, а также распределение давления по длине образующей. Они подтверждены экспериментальными данными автора. Однако необходимо подчеркнуть, что в вышеперечисленных работах рассматривался случай безволнового ламинарного течения вязкой жидкости.
Авторы работы [11] решают двумерную задачу ввода в профиль радиальной скорости безразмерной квазифункции тока. Задача сводится к определению некоторого симплекса, с помощью которого определяются все параметры течения. Авторы также выявили максимально возможный расход жидкости, подаваемой на диск, отвечающий условию «захлебывания» диска. Интересны работы, посвященные исследованиям особенностей гидродинамики перехода ламинарного течения в турбулентное. В работе [12] с помощью коалинового метода обнаружили на поверхности диска в области Яе=(1,9-3,0) 105 ряд треков в форме архимедовой спирали. Авторы предположили, что в центральной части диска существует ламинарный режим течения, а затем течение становится неустойчивым. Устойчивость пленочного течения вязкой жидкости на поверхности вращающегося диска была предметом исследования работ [12-13]. Были получены аналитические выражения, позволяющие определить границы устойчивости в зависимости от параметров течения. Изучение переходного режима течения на вращающемся диске было выполнено в работе [14]. Аналогичные исследования были выполнены авторами работы [15]. Авторы выделяют четыре зоны по радиусу вращающейся насадки: входная, первая ламинарно-волновая, турбулентная, ламинарно-волновая зона с мелкомасштабными возмущениями поверхности. Предложены аналитические зависимости для пределов зон. Влияние газового потока на характеристики жидкой центробежной пленки анализируется введением касательного напряжения на поверхности пленки в работе [16].
В центробежных пленочных аппаратах жидкость обычно в виде струи подается в центр рабочего элемента. В приосевой области складывается сложная трехмерная гидродинамическая обстановка. В приосевой области (ударная зона) конвективные
члены в уравнениях Навье-Стокса играют определяющую роль, а влияние центробежной силы на процесс постепенно нарастает от оси к переферии. Возникает проблема определения длины начального участка, когда гидродинамический пограничный слой полностью заполняет поперечное сечение потока. Вопросу определения длины начального участка Ьн посвящен ряд работ. В работе [17] длину начального участка Ьн предполагают определить в виде:
L,, =
Vy /V (v
1/2
(1)
где V - объемный расход жидкости, V - кинематическая вязкость жидкости, у - параметр, О - угловая скорость вращения ротора. Однако как определить параметр у неизвестно. Авторы
[18]определяют границу начального участка следующим образом:
LH = PV2 /(4л2цю)
10,25
(2)
Определению длины начального участка была посвящена работа [19]
dS
_ = i -i3s3-^L dl i
0
(3)
где о0 - толщина пленки.
Анализ этого уравнения позволил установить, что на безразмерном радиусе £ = 2,09 и далее все интегральные кривые полученные при различных начальных условиях, вырождаются в одну амплитуду, совпадающую с графиком Хинце-Мильборна. Это означает, что при £ = 2,09 заканчивается начальный участок. Численно определение Ьн по работам
[18] и [19] совпадают. Это тем более интересно, что выводы в работе [18] были получены в результате оценки значимости членов исходного уравнения, а в
[19] - путем прямого численного анализа. Длина начального участка Ьн была определена и для неподвижной насадки [20, 21]. Поскольку на начальном участке действие центробежной силы мало по сравнению с силой инерции, то результаты работ [20, 21] могут быть использованы для предварительной оценки гидродинамической ситуации. Подчеркнем, что работ, посвященных экспериментальному определению Ьн в известной нам литературе нет.
Течение тонких пленок вязкой жидкости по криволинейным вращающемся поверхностям было рассмотрено в работе [22]. Изучено стационарное течение жидкости при больших числах Экмана. Рассмотрена возможность получения пленки заданной толщины, ее устойчивость и возникновение на ее поверхности нелинейных волн. Показано, что течение жидкости на криволинейной поверхности более устойчиво, чем по плоскому диску. Неустоновив-шееся течение тонкого слоя вязкой жидкости по поверхности тела вращения, вращающейся с переменной угловой скоростью вокруг своей оси, при наличии внешней массовой силы, действующей вдоль оси тела рассмотрена в работе [23]. Задача решалась при начальных условиях, учитывающих
распределение толщины слоя и скоростей в начальный момент времени.
Исследование течения вязкой жидкости по поверхности вращающегося перфорированного конического диска выполнено в работе [24]. Система уравнений Навье-Стокса решается в двумерной постановке с учетом центробежных и Кориолисовых сил, однако без учета давления. Анализ работ в ламинарной постановке приведен и в работах [19, 25].
В экспериментальных работах, посвященных пленочному течению сплошной пленки по поверхности вращающегося диска, отмечается волнообразование [14,26,27] и переход ламинарного движения в турбулентное [15]. В работе [27] разработана номограмма определения средней толщины жидкой пленки в зависимости от определяющих параметров. Авторы работы [15] для характеристики перехода от ламинарно-волнового течения пленки к турбулентному предлагают использовать приведенное число Рейнольдса Ren, учитывающее влияние газового потока и центробежное число Вебера We„6. Изотермическое стабилизированное турбулентное течение вязкой жидкости при Ь>Ьн на поверхности вращающейся насадки в виде тонкой открытой сплошной пленки было изучено в работе [28]. Была рассмотрена криволинейная чаша в виде поверхности вращения второго порядка. Форма поверхности задается через sin а и cosa . Турбулентность течения жидкости была принята в виде полуэмпирической теории Прантля. Как показывают расчеты для стабилизированного участка (Ь>Ьн) скорость отставания жидкости значительно меньше меридиональной. Поэтому авторами работы было принято в виде
М- T
= P(0,45 о )2
dw„
dl
(4)
Двумерные уравнения движения Навье-Стокса свели к системе обыкновенных дифференциальных уравнений, которые решались численно методом Рунге-Кутта четвертого порядка. Определялись основные технологические параметры течения в безразмерных величинах. В работе обсуждаются различие между разными режимами течения. В некоторых центробежных пленочных аппаратах расход жидкости обеспечивается частично погруженными в жидкость коническими насадками, например аппараты Киршбаума-Штора, Р.Ш. Сафина и т.д. [29]. Ряд работ посвящен проблемам увлечения жидкой пленки вращающимися дисками [30-34]. В этих работах определены кинематические соотношения для определения контура свободной поверхности, предложены условия сшивки двух зон, рассчитаны профили скоростей как в зоне статического мениска, так и в динамической зоне. В данных работах считается, что количество увлекаемой телом жидкости полностью определяется процессами, протекающими выше уровня горизонтальной поверхности. Экспериментальные работы [35, 36] подтвердили зависимость расхода от глубины погружения усеченного конуса в жидкость. Первые попытки
теоретического анализа закономерностей работы аппарата с коническим, частично погруженным в жидкость, ротором имеются в работах [37-39]. Методом Кармана рассчитан пограничный слой на погруженной части конуса. Расчетные формулы получены в предположении постоянства толщины пограничного слоя. Экспериментально была выявлена зона перехода ламинарного режима в турбулентное. Однако решения полученные в работах [37-39], неоднозначно учитывают глубину погружения конуса и не объясняют факт отсутствия подъема жидкости при угловой скорости вращения, ниже некоторого критического значения. Расчет насосного эффекта вращающегося усеченного конуса, частично погруженного в вязкую жидкость в ламинарной и турбулентной постановке был выполнен в работе [40]. Были определены толщины пограничного слоя и
профили скоростей wx и wр и на их базе рассчитана производительность погружной части конуса:
• tcR— Ю253х 2L 4 1 2 V
V =-—-—s211--s +— s2 I sin a . (5)
4 v ^ 3 2 )
Предельная толщина слоя на свободной поверхности жидкости s, которая может подняться центробежной силой, определяется
3 . (6)
13 Г 1/2
-s3 -^s +1 = g/[co2KntgaJ .
Возможность насосного эффекта определяется соотношением:
(7)
ro2Rntga > g.
Турбулентный режим проанализирован на примере эмпирического соотношения «одной седьмой», которые с удовлетворительной точностью описывает поведение жидкости в окрестностях вращающегося конуса. Переход к турбулентному режиму течения сопровождается резким ростом расхода жидкости, понимаемый коническим ротором. Теоретические расчеты подтверждены экспериментальными данными [41-42].
Неньютоновские жидкости
Многие жидкости, обрабатываемые в центробежных пленочных аппаратах, являются реологическими средами. Для описания деформационного поведения вязких неньютоновских жидкостей широко применяется «степенной закон».
Ламинарное течение степенной жидкости по поверхности плоского диска в одномерной постановке было рассмотрено в работе [43]. Были получены зависимости для определения профиля меридиональной (радиальной) скорости и толщины пленки жидкости. В работах [44-47] была решена аналогичная задача для случая, когда рабочий элемент задан в виде криволинейной вращающейся
чаши типа z =
f (Vxr+y2
в специальной вра-
щающейся ортогональной системе координат в одномерной и двухмерной постановках. Были определены тангенциальные и меридиональные скорости,
давление и толщина пленки. Установлен безразмерный параметр по которому можно определить течение одномерное или двумерное.
Как известно, степенной закон не описывает кривую течения полностью, а только среднюю (вторую зону) часть. Для уравнения состояния Кросса, опиывающего все три зоны кривой течения, было проведено решение одномерной задачи численным методом [40]. Осесимметричное установившееся течение нелинейной упруго-вязкой жидкости, подчиняющейся реологическому закону:
1
CTij = -Р5ij + Н-эффВ1 + ^ P1B2
(8)
было рассмотрено в работах [49-52]. Здесь а у - тензор вязких напряжений, р - давление, 8 ¡у - единичный тензор, цэфф - эффективная вязкость, р1 - коэффициент нормальных напряжений, В1 ,В 2 - тензоры Уайта-Метцнера. Коэффициент нормальных напряжений был определен из данных вискозим-метрических исследований по методу Малкина. Была изучена одномерная и двумерная задачи, с учетом инерционных сил и без них. Использование метода интегральных соотношений позволило получить инженерные расчетные соотношения.
Разрыв центробежной пленки аномально-вязкой жидкости, подчиняющейся степенному закону был рассмотрен в работе [53]. Определены параметры процесса разрыва и минимальной плотности орошения вращающейся насадки на базе энергетической и моментной теорией. Определен параметр Я . Если Я < 1, следует пользоваться формулами, полученными из энергетической теории, если Я > 1 - формулами из моментной теории. Экспериментальные исследования по определению критического радиуса разрыва показали, что отклонение теории от эксперимента составляет в основном 25-30%. Разрыв пленки жидкости сопровождается образованием поверхностных струй. Распад этих струй на капли на краю центробежной насадки рассмотрен в работе [54].
В работе [55] разработана математическая модель движения концентрированной двухфазной средств по поверхности ротационного смесителя. Решение этой задачи позволило выбрать оптимальные размеры элементов ротора.
В работах [56-57] построена математическая модель течения степенной жидкости по внутренней поверхности вращающегося проницаемого конического ротора. Задача решена в трехмерной постановке. Автор установил закон изменения коэффициента проницаемости стенки конического ротора по длине образующей. Исследование процесса течения степенной жидкости по внешней поверхности криволинейной насадки выполнено в работе [58]. Получены безразмерные функции для определения меридиональной, осевой и тангенциальной компонент скорости, давления и толщины пленки. Полученные результаты были использованы при анализе процесса разделения тонкодисперсных суспензий. В работе [59-60] процесс течения степенной жидкости по внутренней проницаемой поверхности произволь-
ной формы. Разработан алгоритм и программа численного решения полных трехмерных уравнений движения. Установлены условия получения монодисперсных гранул. В работе [61] найдена радиальная координата пространственного пограничного слоя по поверхности пленки, что позволяет конструировать конические насадки, работающие без «захлебывания» при любом расходе степенной жидкости. Математическая модель течения жидкости по внутренней поверхности усеченного конуса записаны в приближении пограничного слоя. Течение двух несмешивающихся степенных жидкостей по внутренней поверхности конической насадки рассмотрено в работе [62-63].
Теоретическому и экспериментальному изучению течения вязко - пластических сред по оп-верхности вращающихся дисков посвящены ряд работ [64-67]. Вязко-пластичные среды обычно высоковязкие, поэтому задача решена в одномерном приближении. Получены параметры процесса течения, в частности, радиус формосохранения в виде:
r0 =-J— , (9)
Р
где х 0 - предельное напряжение сдвига.
Существование радиуса формосохранения обусловлено наличием относительного равновесия между х0 и напряжением, вызванным центробежной силой. Это обстоятельство позволило авторам работы [68] предложить центробежный метод определения реологических свойств вязкопластических сред.
В заключении отметим, что исследование центробежной жидкой пленки сводится к определению поля скоростей, давлений и толщины пленки, которые необходимы для анализа, гидромеханических, тепло- и массообменных процессов, происходящих в центробежных пленочных аппаратах.
Исследования центробежной жидкой пленки в изотермических условиях, приведенных в настоящей работе, являются отправной (началь-ной)точкой при анализе неизотермических задач.
Литература
1. Н.Х. Зиннатуллин, А.А. Булатов, С.Г. Николаева, Г.Н. Зиннатуллина. Вестник Казанского технологического университета, т.15,1, 125-127, 2012
2. Н.Х. Зиннатуллин, А.А. Булатов, В.Г.Кузнецов, И.М. Нафиков, Г.Н. Зиннатуллина. Вестник Казанского технологического университета, т.15,4, 107-109, 2012
3. А.В. Давыдов, В.В.Бронская,Н.Х. Зиннатуллин, Вестник Казанского технологического университета, т.15,23, 145-147, 2012
4. Н.Х. Зиннатуллин, А.А. Булатов, Р.Г.Галимуллин, А.И. Хайбуллина Вестник Казанского технологического университета, т.16,3, 66-68, 2013
5. J.O.Himze, H.Milbom. Journal of Applied Mechanics, 17,145.1950
6. A.G.Emslie, F.T.Bouner, L.G. Peck. Journal of Applied Physics, 29,858,1958
7. Р.Х.Мухутдинов. инженерно-физический журнал 4,8087, 1961
8. C.R.Adler, W.R.Marshall. Chem.Engng.Progr.,47, 515, 1951
9. А.В.Дятлов, С.Ф. Хохлов Труды Днепропетровского химико-технологического института, 10, 230-235, 1960
10. К.Ю. Вачагин, В.С. Николаев. Известия Вузов «Химия и химическая технология» т.3,6 1097-1102, 1960
11. Г.И.Лепехин, Г.В. Рябчук, Н.В. Тябин, Е.Р.Шульман Теоретические основы химической технологии, т.15,3,391-396, 1981
12. Ф.М.Гимранов, Н.Х.Зиннатуллин, Ф.А.Гарифуллин. Прикладная механика.т.12,7,85-90,1976
13. Г.М. Сысоев, В.Я. Шкадов, Инженерно-физический журнал, т 42, 6, 936-940, 1987
14. Б.И. Федоров, Г.З.Плавник, И.В. Прохоров, Л.Г.Жуховицкий. Инженерно-физический журнал. Т.31, 6, 1060-1067, 1976
15. А.И.Бутузов, И.И.Пуховой. Инженерно-физический журнал. Т.31, 12, 217-224, 1976
16. В.С.Николаев, К.Ю.Вачагин, Ю.Н.Барышев. Известия ВУЗов «Химия и химическая технология», т.10, 2, 237242,1967
17. N.Fraidenraich. Revista Mexicana de Fisika 25, 69-87, 1976
18. К.Гейзли, А.Чарваш Тепло- и массоперенос. Минск. Изв. АН БССР, Т.10, 401-419, 1968
19. О.Г. Николаева. Дисс. канд. техн. наук, КГТУ, Казань, 1998, 186 с.
20. А.А. Коптев, В.А. Юрченко. Тр. Тамбовского института хим.машиностроения, вып.2, 171-176, 1969
21. А. Ю. Глинкин, В.А. Рукавишников, Инженерно-физический журнал, т.39,1, 57-63, 1980
22. Е.И.Могилевский, В.Я. Шкадов. Механика жидкости и газа. 2, 18-32, 2009
23. О.Ф. Васильев, Н.С.Ханилов. ПМТФ, 3, 97-99, 1965
24. В.М.Готовцев, А.И.Зайцев, Н.Ф.Чупрынин, Ю.И. Макаров. Известия ВУЗов «Химия и химическая технология», т.24, 7, 920-923, 1981
25. А.А.Булатов Дисс. канд. техн. наук, КГТУ, Казань, 1987, 187 с.
26. О.А.Поваров, Е.Г.Васильченко, П.Г.Петров, Известия АН СССР Энергетика и транспорт.1,172-176, 1978
27. Г.А.Филлипов, О.А. Поваров, Е.Г.Васильченко, В.Н.Грушин, А.И.Николаевский В сб. Труды МЭИ.Москва, 1977, с.41-46
28. В.В.Бронская, А.В. Давыдов, Н.Х. Зиннатуллин, Вестник Казанского технологического университета, т.15,23, 145-147, 2012
29. Д.Т.Мусин Дисс. канд. техн. наук, КГТУ, Казань, 1987, 167 с.
30. K.Vijagraghvan. Ind. And Eng.Chem. Fundam. Т.21,4, 333-336, 1982
31. Н.И.Гельперин, Г.А. Носов, А.В. Мохоткин. Химическое и нефтеперерабатывающее машиностроение, 3, 1820
32. В.В.Пухначев ПМТФ,3,78-88,1977
33. В.И. Байков, З.П.Шульман, К.Энгельгард. ПМТФ,2 3843, 1984
34. В.И. Байков, ПМТФ, 3 70-78, 1978
35. Ю.И.Макаров В. Сб.: Труды МИХМа, т.19, 109-123
36. Ю.И.Макаров Газовая промышленность. 7, 28-31, 1961
37. В.А.Юрченко, А.А. Коптев, Г.С.Логосов. Химическое и нефтеперерабатывающее машиностроение. 12, 14-15, 1966
38. В.А.Юрченко, Г.С.Логосов, А.Г.Горст. Химическое и нефтеперерабатывающее машиностроение. М.: ЦИН-ТИХИМНЕФТЕМАШ. 1967, 58-61
39. В.А.Юрченко, А.А. Коптев, Г.С.Логосов. Химическое и нефтеперерабатывающее машиностроение. 4, 18-19, 1968
40. Н.Х. Зиннатуллин, А.А. Булатов, Ф.М. Гимранов, Д.Т. Мусин Теоретические основы химической технологии, т.32,6, 1-5, 1998
41. Д.Т.Мусин Дисс. канд. техн. наук, КГТУ, Казань, 1987, 167 с.
42. Н.Х. Зиннатуллин, Ф.М. Гимранов, А.А. Булатов, Изв. ВУЗов «Химия и химическая технология», т.41,3,65-68, 1998
43. A.Acrivos, M.J.Shan, E.E. Petersen. Journal of Applied Physics.v.31, 6, 963-968,1960
44. Н.Х. Зиннатуллин, К.Ю. Вачагин, Н.В.Тябин В.сб.: Труды КХТИ, вып. 35 146-153, 1965
45. Н.Х. Зиннатуллин, К.Ю. Вачагин, Н.В.Тябин Инженерно-физический журнал. Т.15, 2, 234-240, 1968
46. Н.Х. Зиннатуллин, И.М. Нафиков, А.А. Булатов, В.В. Антонов Инженерно-физический журнал. Т.69, 1, 112117, 1996
47. Н.Х. Зиннатуллин, К.Ю. Вачагин, Н.В.Тябин Изв. ВУЗов «Химия и химическая технология», т.12,3,1441-1445, 1969
48. В.П.Костромин, В.Г.Кузнецов, К.Ю.Вачагин Инженерно-физический журнал, т.30,1, 86-91, 1976
49. И.В.Флегентов, Н.Х. Зиннатуллин, В.сб.: Труды КХТИ, вып. 51,115-118, 1973
50. И.В.Флегентов, Н.Х. Зиннатуллин, В.сб.: Труды КХТИ, вып. 51,118-123, 1973
51. Ф.М. Гимранов, И.В.Флегентов, Н.Х. Зиннатуллин, В.сб.: Труды КХТИ, вып. 55,19-24, 1975
52. Н.Х. Зиннатуллин, Ф.М. Гимранов, И.В.Флегентов. Инженерно-физический журнал. Т.31, 2, 231-236, 1976
53. И.М. Нафиков Н.Х. Зиннатуллин. Инженерно-физический журнал. Т.39, 1, 47-50, 1980
54. И.М. Нафиков Н.Х. Зиннатуллин. Изв. ВУЗов «Химия и химическая технология», т.27,6,722-726, 1984
55. Ф.Г.Ахмадиев Дисс. канд. техн. наук, КХТИ, Казань, 1975, 147 с.
56. А.Э.Просвиров, Г.В.Рябчук. Известия РАН МЖГ, 6, 39-43, 1995
57. А.Э.Просвиров, Дисс. канд. техн. наук, ВГТУ ,Волгоград, 1996, 114 с.
58. А.С.Прокопенко, Г.В.Рябчук, Е.А. Селезнова Изв. ВУЗов «Химия и химическая технология», вып.5,121-122, 2002
59. Е.А.Смирнов А.С.Прокопенко, Г.В.Рябчук Изв. ВУЗов «Химия и химическая технология», вып.6,67-70, 2002
60. Е.А.Смирнов А.С.Прокопенко, Г.В.Рябчук Изв. ВУЗов «Химия и химическая технология», вып.2,162-163, 2003
61. В.А.Гордон, В.А.Осокин, Г.В. Рябчук, А.С. Чудин Изв. ВУЗов «Химия и химическая технология», т.48,9,112-115, 2005
62. А.С. Блинов, В.А.Гордон, А.Ю.Орешкин, В.А.Осокин, Г.А.Попович Г.В. Рябчук, Изв. ВУЗов «Химия и химическая технология», т.49,6,112-115, 2006
63. А.С. Блинов, В.А.Гордон, Ю.В.Грабельников, И.А.Никулин Г.В Рябчук, Изв. ВУЗов «Химия и химическая технология», т.50,7,92-95, 2007
64. Н.В.Тябин, Л.А.Шкляр, Е.П.Мосихин, Г.В. Виноградов В.сб.: Труды КХТИ, вып. 18, 123-130, 1953
65. Л.А.Шкляр Дисс. канд. техн. наук, КХТИ, Казань, 1954, 138 с.
66. .М.Гушкин Коллоидный журнал. 22, 573-580, 1960
67. .М.Гушкин Коллоидный журнал. 24, 263-268, 1962
68. Н.В. Тябин,Л.А.Шкляр, Е.П. Мосихин, Г.В.Виноградов.сб.: Труды КХТИ, вып. 16,47-51, 1951.
© Г. Н. Зиннатуллина - канд. ехн. наук, доц. каф. промышленной безопасности КНИТУ; А. В. Давыдов - асп. КГЭУ; Р. Р. Байбеков - ст. препод. каф. процессов и аппаратов химической технологии КНИТУ; Н. Х. Зиннатуллин - д-р техн. наук, проф. той же кафедры; В. В. Бронская - канд. техн. наук, доц. той же кафедры, dweronika@mail.ru.
