CC BY
66
© H.H. Арефьев, C.M. Штин, 2007
УДК 553.543
H.H. Арефьев, С.М. Штин
МЕТОД ОПРЕДЕЛЕНИЯ РЕОЛОГИЧЕСКИХ ХАРАКТЕРИСТИК САПРОПЕЛЯ
Л ля расчета течения сапропеля по каналам гидроаппаратов (насосов) и по трубопроводам необходимо знать его реологические характеристики. Известно, что сапро-пели обладают свойствами неньютоновских жидкостей, а течение их с достаточной степенью приближения может быть описано уравнением Бингама - Шведова [1, 2]:
6Ч (1)
Р = Ро +Пп^~^ , (1)
6И
где р - напряжение сдвига; р0 - предельное напряжение сдвига; г)т - коэффициент вязкости; 64 / 6И - градиент скорости сдвига.
Характер реологического поведения той или иной жидкости или дисперсии и ее реологические параметры можно с достаточной точностью опре делить, лишь изучив процесс ее течения. В практике реологических исследований широко применяют капиллярный и ротационный вискозиметры [3, 4, 5, 6, 7, 8, 9, 10]. Наиболее широкое распространение получили ротационные вискозиметры с коаксиальными цилиндрами. Исследования реологических параметров вязкопластичных жидкостей чаше всего проводят на ротационных вискозиметрах системы М.П. Воларовича РВ-4, РВ-7, РВ-8, которые по величине зазора относятся к группе интегральных приборов. Стандартные методи-
ки вискозиметрических измерении позволяют определить значения эффективной вязкости при фиксированных скоростях вращения ротора вискозиметра, а также определять предельные напряжения сдвига при течении однородных консистентных смазок, коллоидных растворов, высокомолекулярных соединении или тонкодисперсных систем [11]. Однако сапропели содержат грубодисперсные частицы и структурные агрегаты, превышающие по размерам зазоры вискозиметров, в результате при измерениях наблюдается заклинивание роторов вискозиметров. Необоснованное увеличение зазоров приводит к нарушению однородности полеИ напряжении и скоростеи деформации, что приводит к искажению ожидаемых результатов. Кривые течения дисперсных систем, полученные в результате опытов по стандартным методикам измерении на вискозиметрах системы М.П. Воларовича, имеют расхождения: при заданноИ величине касательного напряжения большим радиусам соответствуют меньшие значения градиентов скоростеИ.
По мнению П.А. Ребиндера, Н.Н. Кулакова, М.П. Воларовича, М. РеИ-нера, Дж. ОлдроИда основноИ причиной расхождения получаемых значениИ реологических параметров и зависимости этих параметров от геометрических размеров вискозиметров
является пристенное скольжение. На необходимость учета пристенного скольжения при расчетах течения вязкопластичных жидкостей указывается в работах М. Муни, Грина, Скотт-Блера, А.П. Юфина, А.Х. Кима, Д.М. Толстого, В.В. Синицыны, Г.В. Виноградова. Эффект пристенного скольжения можно представить себе как результат повышения текучести в пристенном слое по сравнению с текучестью в объеме. Первая экспериментальная работа, в которой было обнаружено пристенное скольжение дисперсных систем, была выполнена в 1920 г. Грином [12]. Воларович М.П. и Кулаков Н.Н., проводившие наблюдения для торфомассы, получили значение пристенного слоя воды порядка 0,15 мкм [13, 14]. Причем в стальных трубах эффект оказался значительно большим, чем в стеклянных трубах. Однако коэффициент пристенного скольжения значим только для течения дисперсных систем. Гипотеза прилипания к стенке для низкомолекулярных систем завоевала сейчас всеобщее признание. Влияние проскальзывания на расход воды оказывается едва уловимым (< 5 %) в самом благоприятном случае в опытах Шнел-ла [15]. Поэтому во всех исследованиях эффекта пристенного скольжения [13, 14, 16, 17] предполагается, что дисперсионная среда прилипает к стенке. Пристенное скольжение дисперсных систем по мнению Толстого Д.М. [13] связано с существованием пристенных слоев дисперсионной среды, средняя толщина которых может достигать порядка микрометра, а в грубодисперсных системах и значительно большей величины. На практике для устранения проскальзывания применяют роторы вискозиметров с рифлеными поверхностями рабочих цилин-
дров. Однако применение рифленых измерительных поверхностей не является гарантией полного устранения пристенного эффекта. Это следует уже хотя бы из того, что сушествен-ную роль играет профиль рифлений, особенно профиль внешних частей выступов. Кроме того, опыты на рифленых поверхностях допустимо проводить только при невысоких скоростях деформаций, так как после разрушения структуры в быстро движушемся потоке испытуемой жидкости рифления могут вызывать сильные возмуше-ния [6]. Поэтому такой путь устранения скольжения не является строгим и вполне надежным. Пристенное скольжение можно оценить по методу М. Муни для вискозиметров с коаксиальными цилиндрами [18]. В этом случае необходимо проводить опыты на приборах с набором не менее четырех пар коаксиальных цилиндров различных радиусов, что очень усложняет и удорожает проведение опытов.
Наиболее широко исследователи применяют метод анализа потерь напора в трубах, основанный на теории капиллярной вискозиметрии. При этом зависимость гидравлического уклона от скорости движения исследуемой жидкости в трубопроводе аппроксимируют уравнением Букингэма, принимают линейное распределение касательных напряжений сил трения по живому сечению трубопровода [19, 20]. Большинство исследователей отмечают условный характер получаемых реологических характеристик, их зависимость от диаметров трубопроводов, а также от конструкционного материала трубопровода [10, 21, 22]. Неинвариантность реологических характеристик можно объяснить тем, что метод анализа потерь напора не учитывает эффект пристенного скольжения. Букингэм впервые ввел по-
правку на пристенное скольжение в капиллярном вискозиметре. Она выражается добавочным объемом AQ истечения в секунду, обусловленным скольжением вязкопластичной («бин-гамовской») жидкости по стенке [14]. На режимах установившегося течения пристенное скольжение может быть количественно оценено по методу М. Муни, который в случае капиллярной вискозиметрии дает блестящие результаты [6, 18]. Такой подход к оценке реологических характеристик кормовых смесей использовали B.C. Мкртумян и B.B. Переведенцев [21] и рекомендовали его к практическому применению при исследовании вязкопластичных жидкостей. Сущность метода состоит в следующем. Необходимо экспериментально изучить напорно-расходные характеристики при течении исследуемого материала в горизонтальных трубах минимум четырех - пяти диаметров не менее
0,05 м. По экспериментальным данным следует построить зависимости эффективного градиента скорости
D Q й
De =----з от касательных напряжений
nR
сдвига tr на стенках трубопровода, то есть интегральные реологические кривые. Их расхождение для различных диаметров трубопроводов свидетельствует о наличии пристенного скольжения.
По интегральным кривым находят зависимости De(1/R) при tr = const и по их виду судят о характере зависимости скорости скольжения от напряжений сдвига. Экстраполированием функции De(1/R) при tr = const до 1/ R = 0 находят градиент скорости в объемных слоях исследуемого материала. Путем многократных повторений последней операции при разных tr находят полную реологическую кривую течения исследуемого
материала [21]. В результате удается определить инвариантные реологические характеристики, независяшие от размеров трубопроводов, а также оценить скорости скольжения исследуемого материала относительно стенок трубопроводов.
Очевидно, что такой метод обладает рядом достоинств: получение инвариантных реологических характеристик, возможность оценки скорости скольжения исследуемого материала по стенке трубопровода, а следовательно, возможность учета этой скорости в расчетных методиках трубопроводного гидротранспорта вязкопластичных жидкостей. Но наряду с вышеперечисленными достоинствами очевидны и недостатки этого метода: громоздкость оборудования,
высокая трудоемкость исследовательских работ.
Основываясь на теоретических исследованиях течения вязкопластичных жидкостей с гидросмазкой под действием перепада давления в круглом канале, выполненных в [23], а также на экспериментальных исследованиях течения вязкопластичных жидкостей [21], достаточно точно описываемых реологическим уравнением Бингама -Шведова, разработаем упрошенный метод определения реологических характеристик.
Можно предположить, что проскальзывание при течении по трубам вязкопластичных жидкостей обусловлено образованием тонкой пленки из жидкости, менее вязкой, чем исследуемая. Эта пленка образуется на внутренней поверхности трубы и толшина ее составляет 8 . Также можно предположить, что эта пленка образована дисперсионной средой и является ньютоновской жидкостью с коэффициентом динамической вязкости Л . Например, для сапропеля или торфомассы дисперсионной средой
является вода. На основании этих предположений и результатов исследований [23] можем написать выражение для определения расхода вязкопластичной жидкости с гидросмазкой при течении по трубе:
Q = пИ ^ V;
пИ3
4Ппп
1 -
4 Ро_ 3 Рш
Ро
(2)
где Рщ - напряжение сдвига на стенке трубы; - скорость скольжения.
Если обе части выражения (2) раз, то получим
делить на пИ3
Рщ
П3 И 4Пп
1 -4 р. +1
3 Рш 31 Р,
Ро
Vs =■
м
Р =|-дР|И
щ 1 дг) 2 ’
(4)
(5)
где -
др
дг
оси канала.
Если (4) с учетом (5) ввести в первый член уравнения (2), то получим известное уравнение Букингэма [14] добавочного расхода с учетом пристенного скольжения:
ДQ =
пИ 38 ДР
Это свидетельствует о верности теоретических основ для разработки нового метода реологических исследований.
На основании экспериментальных данных, полученных при исследовании реологических характеристик смесей [21], которые можно отнести к жидкостям Бингама, было получено, что скорость скольжения можно определить по формуле:
V; = И? Р - Ь), (6)
где а и Ь - величины, зависящие от свойств вязкопластичной жидкости.
Приравнивая правые части выражений (4) и (6), получим
(3)
Левая часть выражения (3) представляет собой эффективный градиент скорости О, вычисленный по среднеобъемной скорости потока, первый член правой части - градиент скорости в пристенном слое, второй -величину, пропорциональную градиенту скорости в объемных слоях.
На основании исследований [23] можно написать выражение для определения и Рщ:
Р„д
8 = м 11 - А
И2 I Рщ
(7)
Таким образом, по формуле (7) можно определить толщину смазывающего пристенного слоя 8 в зависимости от свойств вязкопластичной жидкости, радиуса трубы и напряжения сдвига на стенке трубы. Из уравнения (3) с учетом (6) получим:
О = ;|г (Рщ - Ь) +
Р.
4п
IП
1 - 4 Р. +
3 Рщ
1 Г Ро 4 V Рщ у
(8)
градиент давления вдоль
Уравнение (8) справедливо для ламинарного установившегося изотер-миического течения. На основании уравнения (8) разработан новый метод экспериментальных исследований, которые проводятся при движении вязкопластичной жидкости по двум трубам с радиусами Их и И2. При этом замеряются эффективные градиенты скорости (ОИ)! и напряжения сдвига на стенках Рщ/. Тогда на основании уравнения (8) можно написать для одинаковых напряжений Рщ1 для обоих диаметров труб:
),= Д3 ( - Ь) +
РШ,
4п
і п
Р
•* ш.
4п
ІП
щ
1 _ 4 Рь +1 (£и
3 РШ1 3 1 РШ1
Рщ _ Ь
(°щ1 )р«,„ (°щ2 )р«,„ Ршк Ь
(11)
а =
(Д )рщ (0Я2 )рщ
щ3 _ Щ Ї (ш] _ Ь
(12)
(Ря2 )рШк щ3 (РШк Ь) Ршк _ 3 Ро + 33 ррг
Рш
Ш]
1 4 Ро + 1 ' Ро Ї4 Пп/
_ 3 Рш, 3 V Рш, у (9) 4
1 _3 Р-+1
( V Ро
3 Р
Ш]
V РШ] у
(14)
Из системы уравнений (9) получим для } ф к
(ОИ1 2 ). = а Г -1 - -1 )( - Ь ),
(0и1 Р.к ~(0и2 )Р„к = а [й3 “ И3 - Ь )•
(10)
Если поделить правые и левые части уравнений системы (1о) соответственно, то получим
( 0и1 к“ ( 0и 2 Рр
Из уравнения (11) определяется значение Ь. Из одного из уравнений (1о) определяется значение а:
Из системы уравнений (9) можем написать для ] ф к :
(Вщ1 к _ (Рш _ Ьр Рш] _ 3 Ро + 3
(13)
Из уравнения (13) определяется значение ро- Значение пПё можно определить из выражения:
Следовательно, вновь разработанный метод определения реологических характеристик вязкопластичной жидкости позволяет сократить объем исследовательских работ по сравнению с известным методом М. Муни. При этом необходимо отметить, что в прежнем методе М. Муни применяется графо-аналитический метод обработки экспериментальных данных для получения реологических характеристик [21], в то время как во вновь разработанном - чисто аналитический, что свидетельствует о большей достоверности получаемых результатов по новому методу.
Для проверки разработанного метода были использованы результаты экспериментальных исследований течения кормовых смесей [21], которые относятся к жидкостям Бингама. В работе [21] на лабораторном стенде были исследованы реологические характеристики кормовых смесей по методу М. Муни: каждую кормовую смесь транспортировали по пяти трубам с радиусами 25, 40, 52, 65 и 75 мм. При этом на разных режимах выполнялись следующие замеры: потери давления на отрезке трубы и объемный расход жидкости. По известным формулам авторами [21] были вычислены напряжения сдвига на стенке трубы Рш и эффективные градиенты скорости В Результаты замеров представлены в виде зависимостей В = ¡(Рш) для всех типоразмеров труб и смесей.
Для проведения проверки разработанного метода исследований были взяты экспериментальные данные из [21] для двух смесей (№3 и №6) при
а) Смесь А/3
Ъ, '/с
90 30 9а Р„, Па
Рис. 1
1 - труба диаметром 80 мм
2 - труба диаметром 130 мм
течении их по двум трубам радиусами 25 и 52 мм. На основании опытных данных с помощью вновь разработанного метода исследований были рас-
считаны реологические характеристики смесей №3 и №6. По рассчитанным реологическим характеристикам смесей №3 и №6 были вычислены эффективные градиенты скоростей В по (8) для труб радиусами 4о и 65 мм. Теоретически полученные параметры В = {(.Рш) сравнивались с имеющимися зависимостями, определенными экспериментально на лабораторном стенде [21]. На рисунке показаны теоретически полученные кривые зависимости В = {(.Рш) для труб радиусами 4о и 65 мм и экспериментальные значения [21] (рис.
1,а - смесь №3, рис. 1б -смесь №6). Сравнение показывает, что максимальная относительная погрешность между теоретически и экспериментально полученными зависимостями не превышает о,5 %, откуда можно сделать вывод, что разработанный метод определения реологических характеристик вязкопластичной жидкости по результатам исследований на двух трубах позволяет получать достоверные данные, обеспечивающие расчет режимов течения вязкопластичной жидкости по трубам различных диаметров с высокой степенью точности.
------------ СПИСОК ЛИТЕРАТУРЫ
1. Методические указания по расчету гидравлического транспорта сапропелей. Министерство мелиорации и водного хозяйства СССР. Главнечерноземводстрой. ВНИИГиМ им. А.Н. Костякова. (Утвержде-
ны НТС Главнечерноземводстроя при Мин-водхозе СССР, протокол № 216 от
20.11.80 г.). - М., 1981. - 53 с.
2. Лопотко М.З., Лецко А.П., Дубинин С.К. Рекомендации по технологии промыш-
ленной добычи сапропелей из открытых водоемов для удобрений. (Академия наук Белорусской ССР. Институт торфа). - Минск: Наука и техника, 1981. - 78 с.
3. Воларович М.П. Исследование реологических свойств дисперсных систем. // Коллоидный журнал, т. XVI, №3, 1954, с. 227 - 240.
4. Лишгван И.И. Микро- и макрореология дисперсных систем. - Минск: 1975. - 39 с.
5. Косаревич И. В. Структурообразова-ние в дисперсиях сапропелей. - Минск: Наука и техника, 1990. - 248 с.
6. Белкин И.М. и др. Измерение вязкости и физико-механических характеристик материалов. - М.: Машиностроение, 1967. -272 с.
7. Уилкинсон У. Л. Неньютоновские жидкости. Гидромеханика, перемешивание и теплообмен. - М.: Мир, 1964. - 216 с.
8. Воларович М.П. Применение методов исследования вязкости и пластичности в прикладной минералогии. // Тр. ин-та прикл. минералогии. Вып. 66, 1934, с. 46 - 52.
9. Шишенко Р. И. Гидравлика глинистых растворов. - Баку: Госнаучтехиздат нефтяной и горнотопливной промышленности, 1951. - 133 с.
10. Полишук В.Г. Разработка метода определения реологических характеристик текучих сельскохозяйственных материалов для проектирования гидромеханизации на животноводческих фермах и комплексах. // Реф. дис. на соиск. учен. степ. канд. техн. наук. - Саратов, 1978. - 18 с.
11. Воларович М.П., Лазовская Н.В. Ротационные вискозиметры для исследования реологических свойств дисперсных систем высокомолекулярных соединений. // Коллоидный журнал, т. 28, вып. 2, 1966, с. 28
- 47.
12. Greer H. Industrial Rheology and Rheologieal struktures. New-Jork, 1949. - p. 56.
13. Толсгой Л.М. Скольжение жидкостей и дисперсных систем по твердым поверхностям. Дис. на соиск. учен. степ.
верхностям. Дис. на соиск. учен. степ. докт. физ.-матем. наук. - М., 1953. - 349 с.
14. Толстой Л.М. Об эффекте пристенного скольжения дисперсных систем. Происхождение, размеры и значение эффекта. // Коллоидный журнал, 1947, т. 9, №6, с. 450 - 461.
15. Шульман З.П., Берковский Б.М. Пограничный слой неньютоновских жидкостей.
- Минск: Наука и техника, 1966. - 238 с.
16. Олдройд Лж. Г. Неньютоновское течение жидкостей и твердых тел. //В кн. Реология: Теория и приложения. Под редакцией Ф. Эйриха. - М.: Изд-во ИД, 1962, с. 757 - 793.
17. Metzner A.B. et al. Jnhomogeneous flows of non-newtonian fluids: generation of spatial concentration gradients. // Journal of non-newtonian fluids Mechanies, 1979, V.5, p. 449-462.
18. Mooney M. Explisit formulas of slip and fluidity. // Journal Reology, V.2, №2, 1931, p. 210.
19. Смолдырев A.E., Сафонов Ю.К. Трубопроводный транспорт концентрированных гидросмесей. - 2 изд. -
М.: Машиностроение, 1989. - 256 с.
20. Юфин А.П. Гидромеханизация. - М.: Стройиздат, 1974. - 222 с.
21. Мкртумян В.С., Переведенцев В.В. Реологические исследования полужидких кормовых смесей. // Механизация и электрофикация социалистического сельского хозяйства, 1971, №1, с. 39 - 42.
22. Чижов A.E. Исследование перемещения вязкопластичных кормосмесей и воды по пластмассовым трубам и резинотканевым рукавам. // Реф. дис. на соиск. учен. степ. канд. техн. наук. - Курск, 1974.
- 19 с.
23. Арефьев H.H. Течение вязкопластичной жидкости в круглом канале с гидросмазкой. // Деп. в ЦБНТИ концерна Росречфлот 14.03.91, № 291-рф. (Опубл. в НТИО вып. 5, 1991). - М., 1991. - 12 с.
— Коротко об авторах-----------------------------------------------------------
Арефьев Николай Николаевич - кандидат технических наук, Октябрьский ССРЗ, Н. Новгород,
Штин Сергей Михайлович - кандидат технических наук, Московский государственный горный университет.
