УДК 541. 182. 022: 532. 135 Ю. Н. Тимошин, Е. А. Кирсанов
АНАЛИЗ КОЭФФИЦИЕНТОВ ОБОБЩЁННОГО УРАВНЕНИЯ ТЕЧЕНИЯ ANALYSIS OF COEFFICIENTS OF THE GENERALIZED FLOW EQUATION
Московский государственный областной социально-гуманитарный институт, Коломна, Московская область. E-mail: Kirsanov47@mail.ru
Поведение коэффициентов реологического уравнения при изменении физикохимических характеристик системы обсуждается в рамках обобщённой модели течения структурированных систем. Объясняется зависимость величины коэффициентов от температуры, вязкости дисперсионной среды, Z- потенциала.
Ключевые слова: обобщённая модель течения, реологическое уравнение, характеристики дисперсной системы.
Behavior of rheological equation coefficients is discussed in the framework of the generalized flow model of structured systems under change of physical and chemical characteristics of the system. The dependence of the coefficients on temperature, viscosity of dispersion medium and Z - potential is explained.
Key words: generalized flow model, rheological equation, characteristics of the dispersed system.
Введение
Для проверки справедливости новой реологической модели совершенно недостаточно оценить точность аппроксимации экспериментальных данных. Как уже отмечалось [1], большое число эмпирических и полуэмпирических уравнений способно удовлетворительно описать пластичное или псевдопластичное поведение дисперсных систем в силу сходства математической формы уравнений. Увеличение числа коэффициентов (реологических параметров) значительно улучшает качество «подгонки», не облегчая понимания физико-химических процессов течения.
Наличие микрореологической модели, возможность получить уравнения, связывающие коэффициенты со структурными и физико-химическими характеристиками системы, придают реологической модели необходимую достоверность. Обобщённая модель течения [2] позволяет предсказать поведение коэффициентов реологического уравнения при изменении свойств дисперсной фазы и дисперсионной среды. Обобщённое уравнение течения связывает напряжение сдвига со скоростью сдвига:
12
1/2 _ гс 1/2 .1/2.
1 +-*-.1/2 V Y
~Л+^с 'У'
© Тимошин Ю. Н., Кирсанов Е. А., 2012
Коэффициент вязкости Кэссона лс соответствует предельной эффективной вязкости при у ^ да, в этом предельном случае агрегаты полностью разрушены. Если коэффициент х равен нулю, то размеры агрегата неограниченно возрастают при у ^ 0 . В этом случае наблюдается нелинейное пластичное течение, коэффициент тприобретает смысл предельного динамического напряжения и имеется тенденция к образованию сплошной сетки (каркаса) при нулевой скорости сдвига. Если х больше нуля, то размеры агрегата остаются конечными при бесконечно малой скорости сдвига, при этом наблюдается псевдопластичное течение. Коэффициент т^ характеризует степень агрегации
системы, он уменьшается при уменьшении объёмной концентрации Ф и равен нулю, если агрегаты отсутствуют. В последнем случае наблюдается ньютоновское течение.
Получены [2, 3] уравнения для коэффициентов:
1/2 1/2 „,1/2/1 7 — А
Лс = Л» =Ло (1 — кФ)
т 1/2 = т 1/2
‘'с 0с
1
\А2
V1 — к 2Ф/
(2)
(3)
-12 =
2 ( F
Аа2
.1/2
где
а 2а — 1
а2 =
2 А2 +1 а
4
а = —. 3
Коэффициенты к и А описывают гидродинамическое взаимодействие отдельных частиц в сдвиговом течении, коэффициенты к2 и ^2, в общем, связаны с гидродинамическим взаимодействием агрегатов При этом коэффициенты к и А не обусловлены какой-либо реальной упаковкой частиц в агрегатах.
Таким образом, коэффициенты связаны с вязкостью дисперсионной среды л0 и с
силой сцепления между частицами . В уравнение (3) входит величина FA , т. е. сила, препятствующая разрыву агрегата и нормированная на единицу площади сечения агрегата ( FA = F 8 / п Г ). В общем виде:
л!2 = лГ МФ), т/2 = РТ, Г,_ (Ф).
Функциональные зависимости от объёмной концентрации Ф, показанные в уравнениях (2) и (3), подтверждены в работах [3, 4].
Можно предположить, что характер гидродинамического взаимодействия (т. е. значения коэффициентов к и А) не зависит существенно от вязкости жидкости, в которой диспергированы частицы, и от размеров частиц, тогда величина (лс / л0 ) не должна
зависеть от вязкости ло и от размеров частиц (среднего радиуса г). Для суспензий мо-нодисперсных латексов, например, показано [5], что величина лУ2 практически не зависит от концентрации электролита в водной дисперсионной среде или от размеров частиц при постоянной объёмной концентрации Ф, а величина тУ2 увеличивается с уменьшением размеров частиц или с уменьшением концентрации электролита.
1
3
3
Рассмотрим более подробно характер изменения коэффициентов реологического уравнения при изменении физико-химических характеристик системы.
Водная суспензия синтетического латекса GR-S X-667
Обобщенная модель течения [2, 3] использована для описания течения водной суспензии синтетического латекса, исследованного Мароном и Фоком [6]. Суспензия GR-S Х-667 представляет собой сферические частицы сополимера бутадиен-стирола. При объемной концентрации Ф > 0,25 суспензия является неньютоновской жидкостью. Средний диаметр частиц Ds = 0,14 мкм.
На рис. 1 показаны экспериментальные данные и теоретические кривые течения, построенные с помощью обобщенного уравнения течения (1) по табличным данным Марона и Фока [6]. Значения реологических параметров приведены в табл. 1.
При низкой объемной концентрации Ф система обладает предельным динамическим напряжением сдвига т с и может быть отнесена к нелинейным пластичным жидкостям (х = 0). При более высокой концентрации суспензия обнаруживает псевдопла-стичное поведение (х ф 0).
7
6
5
4
3
2
1
0
.1/2
50
100
150 у1'
а
б
Рис. 1. Реологическое поведение водной суспензии GR-S Х-667 в корневых координатах
при различных температурах: а - при минимальной концентрации Ф = 0,2990, б - при наибольшей концентрации Ф = 0,6017. Сплошные кривые построены по уравнению (1). Пояснения в табл. 1
0
Таким образом, не существует принципиального отличия между пластичным и псевдопластичным поведением в водной суспензии синтетического латекса. При увеличении концентрации дисперсной фазы возникают некоторые ограничения на размеры агрегатов, которые препятствуют их неограниченному росту и образованию сплошной сетки при бесконечно низкой скорости сдвига. Возможно, эти ограничения связаны с увеличением числа столкновений агрегатов частиц.
Поведение коэффициентов реологического уравнения (табл. 1) позволяет сделать выводы об увеличении значений г°./2и г]12 с ростом объёмной концентрации Ф,
уменьшении г]12 с ростом температуры. Коэффициент х характеризует рыхлость или компактность упаковки частиц внутри агрегатов, его нулевое значение указывает на образование сетки при нулевой скорости сдвига. Большие положительные значения х свидетельствуют о возможности спонтанного отрыва частиц от агрегатов и о конечном размере агрегатов при у ^ 0 .
Таблица 1
Коэффициенты обобщённого уравнения течения водной суспензии синтетического латекса GR-S X-667 при различных объёмных концентрациях дисперсной фазы и при различных температурах образца
Ф 20,13°С (ромб) 29,94°С (треугольник)
ГЇ'1 1/2 Лс ж гГ _,1/2 Лс Ж
0,2990 0,0856 0,0549 0 0,0926 0,0491 0
0,3970 0,227 0,0720 0 0,213 0,0648 0
0,4398 0,338 0,0839 0 0,370 0,0747 0
0,4876 0,556 0,100 0 0,534 0,091 0
0,5390 1,861 0,121 8,33 1,88 0,108 8,73
0,5603 1,82 0,148 3,82 1,72 0,137 3,48
0,5866 2,15 0,192 1,41 2,05 0,181 1,34
0,6017 2,51 0,232 0,95 2,65 0,208 1,30
0 - 0,0316 - - 0,0282 -
Ф 39,74°С (квадрат) 50,00°С (крест)
ГЇ1 _,1/2 Лс X гГ _,1/2 Лс X
0,2990 0,0933 0,0446 0 0,0755 0,0411 0
0,3970 0,247 0,0583 0 0,237 0,0537 0
0,4398 0,355 0,0683 0 0,364 0,0628 0
0,4876 0,521 0,0835 0 0,603 0,0747 0
0,5390 2,13 0,0942 10,6 1,71 0,0943 7,29
0,5603 2,14 0,118 6,19 2,11 0,109 5,97
0,5866 2,27 0,159 1,95 2,52 0,142 2,65
0,6017 2,73 0,189 1,37 2,79 0,175 1,59
0 - 0,0256 - - 0,0234 -
В табл. 1 также приведены значения корня вязкости воды л°12 при соответствующих температурах.
На рис. 2 представлены результаты проверки уравнения (2). Приемлемая аппроксимация получается при исключении из расчёта значений коэффициентов при объёмной концентрации Ф = 0,539. Эти точки представлены на графике, но выделены «тенью». По нашему мнению [7] некоторая неоднозначность рассчитанных значений коэффициентов связана с некорректностью обратной задачи и необходимо повышение
точности измерений, чтобы получить «гладкие» зависимости коэффициентов от физико-химических характеристик системы.
0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 Ф
Рис. 2. Зависимость величины (лУ2 ) от объемной концентрации Ф для водной суспензии латекса GR-S Х-667 при температурах: 20,13 °С (ромб); 29,94 °С (треугольник);
39,74 С (квадрат); 50,00 С (крест). Аппроксимация с помощью уравнения (2). Значения (Л,°/2) при концентрации 0,5390 явно занижены и исключены из аппроксимации. Коэффициенты уравнения (2) представлены в табл. 2
Рис. 3. Зависимость величины (г°./2) от объемной концентрации Ф для GR-S Х-667
при температурах: 20,13 °С (ромб); 29,94 °С (треугольник); 39,74 °С (квадрат); 50,00 °С (крест). Аппроксимация выполнена с помощью уравнения (3), коэффициенты указаны в табл. 2
Таблица 2
Коэффициенты уравнений (2) и (3) для водной суспензии синтетического латекса GR-S X-667 при различных температурах образца и рассчитанное значение силы сцепления между частицами FS
t , °С k А k2 А2 т1'2 ‘'0 Fs 1016 Символы
20,13 1,49 0,87 0,051 222,9 0,0026 4,2 Ромб
29,94 1,47 0,92 0,054 222,9 0,0017 1,8 Треугольник
39,74 1,49 0,88 0,051 222,9 0,0029 5,2 Квадрат
50,00 1,50 0,86 0,051 222,9 0,0029 5,2 Крест
Коэффициенты к и А уравнения (2) почти не зависят от температуры системы, причём увеличение к всегда связано с уменьшением А, что явно обусловлено особенностями аппроксимации методом суммы квадратов разностей. Поэтому кривые на рис. 2 расположены параллельно друг другу. Изменение величины лУ2 с температурой, ско-
рее всего, связано с изменением вязкости воды (лУ2).
1/2 , 1/2
Л С /Л0
о-
-0
і і і і
15 20 25 30 35 40 45 50
3,1
3,2
3,3
3,4
3,5
3
а
б
Рис. 4. Температурные зависимости коэффициента вязкости Кэссона водной суспензии GR-S Х-667: отношение лУ 2 / лУ 2 при разных температурах при увеличении концентрации от 0,299 до 0,6017 (а); ^ цс как функция обратной температуры (1 / Т) (б)
На рис. 4, а показана вариация отношения лУ2 /лУ2 при изменении температуры. С учётом неоднозначности коэффициентов реологического уравнения, можно считать, что это отношение остаётся постоянным, что вполне согласуется с постоянством коэффициентов к и А. В соответствии с законом Аррениуса получаются прямые в системе координат (^ щс - 1/Т), причём коэффициент наклона прямых, полученных при увеличении концентрации дисперсной фазы, практически совпадает с наклоном нижней прямой, соответствующей вязкости дистиллированной воды. На графике (рис. 4, б) указан интервал концентраций. Таким образом, изменение коэффициента вязкости Кэс-сона лс с температурой целиком определяется температурной зависимостью вязкости дисперсионной среды л0 . Это представляется логичным, поскольку коэффициент Лс соответствует вязкости системы с полностью разрушенными агрегатами.
На рис. 2, б представлена зависимость коэффициента агрегации от объёмной концентрации и выполнена аппроксимация по уравнению 3 отдельно для каждой температуры. Судя по разбросу точек на графике, регулярной зависимости т1°2 от температуры не наблюдается, что отражается на величине коэффициентов т°./2, к2 и А2 уравнения 3 в табл. 2. Можно считать, что А2 = 222,94; к2 = 0,051 и оценить величину силы сцепления, т. е. силы, необходимой для разрыва связи между соседними частицами FS = 5-10 - Н. Легко допустить, что сила сцепления практически не изменяется в столь малом интервале температур (от 20,13 до 50 °С).
Водная суспензия оксида алюминия
Рассмотрим поведение коэффициентов реологического уравнения на примере водной суспензии порошка оксида алюминия а-А12 03, стационарное течение которой исследовано в [8]. Частицы имеют почти сферическую форму, средний диаметр 0,43 мкм, вязкость воды при температуре измерения 15 °С равна 1,14 мПа с. Концентрация электролита №С1 равна 10"2 М. £ - потенциал частиц был измерен в разбавленной суспензии в соответствии с величиной pH раствора (данные в табл. 3). Экспериментальные данные [8] и кривые, рассчитанные по уравнению (1), показаны на рис. 5.
Таблица 3
Коэффициенты обобщенного уравнения течения для суспензии оксида алюминия при объёмной концентрации Ф = 0,175 и различном значении £ - потенциала
£ - потенциал, мВ 18 21,5 23,5 25,5 29,5
pH 7,98 7,75 7,58 7,41 6,97
гУ2, Па1/2 1,98 0,97 0,73 0,58 1,38
л1/2,( Пас)1'2 0,065 0,069 0,067 0,064 0,044
X , с^ 0,209 0 0 0,042 18,8
Символы кривой ромб треуг. откр. треуг. квадрат круг
а б
Рис. 5. Реологическое поведение водной суспензии оксида алюминия: а - зависимость вязкости от скорости сдвига в двойных логарифмических координатах; б - кривые течения в корневых координатах. Пояснения в табл. 3
Как указывалось авторами [8], увеличение величины £ - потенциала приводит к увеличению электростатического отталкивания частиц и к уменьшению силы взаимодействия между частицами. Действительно, увеличение электрокинетического потенциала от 18 до 25,5 мВ приводит к значительному уменьшению коэффициента агрегации тс, который прямо пропорционален силе сцепления между частицами. Величина
коэффициента вязкости Кэссона г)с практически постоянна, поскольку его значение не
зависит от степени агрегации частиц в системе, а связано с потерями энергии при обтекании отдельных частиц (табл. 3).
Исключением является поведение суспензии при £ - потенциале частиц, равном 29,5 мВ, кривая течения которого резко отличается от кривых, расположенных на графике выше.
Неводные суспензии поливинилхлорида (пластизоли)
ПВХ-пластизоли представляют собой неводные суспензии поливинилхлорида, полученные эмульсионной полимеризацией. Пластизоли демонстрируют неньютоновское поведение, которое существенно зависит от природы дисперсионной среды [9]. При нагревании частицы смолы могут растворяться в пластификаторе, кроме того, происходит набухание смолы, что приводит к увеличению объёма дисперсной фазы. Коллоидная стабильность ПВХ-пластизолей обусловлена их стерической стабилизацией. Роль стабилизирующей добавки играет растворённая часть полимерных цепей ПВХ, которая необратимо закреплена на ПВХ-частице.
Полагают, что частицы ПВХ в минеральном масле или глицерине способны к образованию агрегатов (значительно флокулированная суспензия). На рис. 6 показаны кривые течения в корневых координатах для частиц ПВХ большого размера. Реологические характеристики системы приведены в табл. 4.
а б
Рис. 6. Кривая течения пластизоля РУС в корневых координатах при объёмной концентрации 0,2 и среднем диаметре частиц 1154 нм в минеральном масле (а) и в глицерине (б)
Пластификатор ди-2-этилгексил 0-фталат (ООР) имеет ньютоновскую вязкость
0,059 Па с и частично растворяет частицу поливинилхлорида. Такая дисперсионная среда (ООР) обеспечивает флокуляционную стабильность пластизоля. Кривые течения в корневых координатах сходны с теми, что показаны на рис. 6, и так же описываются обобщённым уравнением течения. Структурно-реологические параметры различных суспензий пластизолей ПВХ представлены в табл. 4.
Таблица 4
Реологические характеристики ДВ’Л'-пластизолей
Система Ф 1172 с т!/2 п X 1 /і 'с 'о
PVC-C (1154 нм) в минеральном масле (по=0,03 Пах) 0,2 0,255 2,00 0,187 2,16
PVC-C (1154нм) в глицерине (0,103 Па-с) 0,2 0,458 3,73 0,221 2,04
PVC-A (293нм) в DOP (0,059 Па-с) 0,2 0,354 0,291 0,03 2,12
РVC - В(619нм)в DOP (0,059 Па-с) 0,2 0,377 0,144 0,03 2,40
PVC - В (619нм) в DOP (0,059 Па-с) 0,54 4,04 6,62 0 276,8
Для частиц разного размера в жидкости с разной вязкостью, но при одинаковой объёмной концентрации Ф, величина отношения гг I го является практически постоянной. Это означает, что коэффициенты к и А в уравнении (2) мало изменяются при переходе от одной системы к другой, вне зависимости от вязкости дисперсионной среды или размера частицы.
Коэффициент ^существенно выше в тех жидкостях, которые являются плохими растворителями, это вполне понятно, поскольку плохая растворимость способствует флокуляции ПВХ-частиц в этих жидкостях.
Величина т\п увеличивается с уменьшением размера частиц при фиксированной объёмной концентрации Ф, поскольку растёт число частиц и тем самым увеличивается возможность агрегирования. Это происходит в соответствии с уравнением [10]
т1п = BN~2I ~ = B (Ф /4/3nr 3)(~2/ ~), (5)
в котором N - полное число частиц в единице объема, к2 - константа скорости форми-
~ ■ 1/ 2
рования агрегатов, к1 у - константа скорости разрушения агрегата под действием растягивающих гидродинамических сил.
Увеличение объёмной концентрации Ф приводит к резкому увеличению обоих
коэффициентов Кэссона (т1/2 и г112).
п П
Таким образом, характер изменения коэффициентов обобщённого уравнения течения согласуется с предсказаниями обобщённой модели течения.
Список использованной литературы
1. Кирсанов Е. А., Тимошин Ю. Н. II Жидкие кристаллы и их практическое использование.
2010. Вып. 4 (34). С. 38 - 47.
2. Кирсанов Е. А. II Жидкие кристаллы и их практическое использование. 2006. Вып. 1/2. С. 21 - 30.
3. Кирсанов Е. А., Ремизов С. В., Новоселова Н. В., Матвеенко В. Н. II Вестн. Моск. ун-та. Сер. 2 : Химия. 2007. Т. 48. № 1. С. 22 - 26.
4. Матвеенко В. Н., Кирсанов Е. А., Ремизов С. В. II Вестн. Моск. ун-та. Сер. 2 : Химия. 2006. Т. 47. № 6. С. 393 - 397.
5. Кирсанов Е. А., Тимошин Ю. Н., Новоселова Н. В., Матвеенко В. Н. II Вестн. Моск. ун-та. Сер. 2 : Химия. 2006. Т. 47. № 6. С. 387 - 392.
6. Maron S. N., FokS. M. II J. Colloid Sci. 1955. Vol. 10. P. 482 - 493.
7. Кирсанов Е. А., Тимошин Ю. Н. II Жидкие кристаллы и их практическое использование.
2011. Вып. 1 (35). С. 48 - 56.
8. Ishihara T., Katsuki H., Kuno H. II Rheologica Acta. 1987. Vol. 26. P. 172 - 181.
9. Willey S. J., Macosko C. W. II J. Rheology. 1978. Vol. 22. № 5. P. 525 - 545.
10. Тимошин Ю. Н., Кирсанов Е. А. II Вестник МГОСГИ. Математические и естественные науки. 2010. № 1 (9). С. 97 - 102.
Поступила в редакцию 29.11.2011 г.