ФИЗИКА
УДК 541. 182. 022: 532. 135
DOI: 10.18384/2310-7251-2022-4-6-16
НЕНЬЮТОНОВСКОЕ ТЕЧЕНИЕ НАНОЖИДКОСТИ НА ОСНОВЕ ОКСИДА ТИТАНА
Вековищев М. П., Кирсанов Е. А.
Государственный социально-гуманитарный университет
140411, Московская обл., г. Коломна, ул. Зелёная, д. 30, Российская Федерация
Аннотация
Цель: рассмотреть реологическое поведение наножидкости, полученной на основе нано-частиц оксида титана в воде и этиленгликоле.
Процедура и методы. Проведена аппроксимация экспериментальных данных уравнениями структурной реологической модели на отдельных интервалах скорости сдвига. Результаты. Показана связь коэффициентов реологических уравнений с характером изменения структуры наножидкости, а именно формирования и разрушения агрегатов на-ночастиц.
Теоретическая и/или практическая значимость. Предложены уравнения, которые способны аппроксимировать экспериментальные данные на отдельных интервалах скорости сдвига, соответствующих определённому структурному состоянию наножидкости. Ключевые слова: наножидкость, наночастицы оксида титана, структурная реологическая модель, реологические кривые
NON-NEWTONIAN FLOW OF A NANOFLUID BASED ON TITANIUM OXIDE
M. Vekovishchev, E. Kirsanov.
State University of Humanities and Social Studies
ul. Zelenaya 30, Kolomna 140411, Moscow Region, Russian Federation
Abstract
Aim. We consider the rheological behavior of a nanofluid obtained on the basis of titanium oxide nanoparticles in water and ethylene glycol.
Methodology. The experimental data are approximated using equations of the structural rheological model on separate intervals of the shear rate.
© CC BY Вековищев М. П., Кирсанов Е. А., 2022.
W
Results. A relationship is obtained between the coefficients of rheological equations and the nature of changes in the structure of the nanofluid, namely, the formation and destruction of nanoparticle aggregates.
Research implications. Equations are derived that make it possible to approximate experimental data at individual shear rate intervals corresponding to a certain structural state of the nanofluid. Keywords: nanofluid, titanium oxide nanoparticles, structural rheological model, rheological curves
Введение
Наножидкости представляют собой ультрадисперсные системы с жидкой дисперсионной средой [1-4]. В случае твёрдых наночастиц систему часто называют золем или коллоидным раствором наночастиц. Некоторые исследователи [2] считают такие наножидкости особым объектом, отличным от традиционных систем коллоидной химии. Одним из важных применений наножидкости является увеличение теплопередачи в жидких средах.
При получении наножидкостей обычно используют высокодисперсные порошки стабильных металлов, их оксидов с характерными размерами частиц от 1 до 100 нм. В качестве жидкой фазы чаще всего выступают водные растворы, содержащие различные органические добавки и соли, препятствующие слипанию частиц. После смешивания нанопорошка с несущей жидкостью образец нано-жидкости подвергают длительному воздействию ультразвука [5-8].
Наибольшее внимание исследователей направлено на зависимость вязкости наножидкости от концентрации и размера наночастиц, а также от температуры. Однако при достаточно высокой концентрации наночастиц наножидкости демонстрируют неньютоновское поведение [5; 7-10]. Поэтому имеет смысл рассмотреть реологическое поведение наножидкостей с точки зрения структурной реологической модели [11; 12].
Наножидкость на основе оксида титана в воде
В работе [8] исследовано стационарное течение наножидкости на основе наночастиц оксида титана в воде (рис. 1, а). Наночастицы имеют размеры от 7 до 20 нм и диспергированы в деионизированной воде.
а б
Рис. 1 / Fig. 1. Реологические кривые наножидкости на основе наночастиц оксида титана в воде при температуре 25° С и объёмных концентрациях Ф= 0,12 (1), 0,10 (2), 0,08 (3), 0,05 (4): а - кривые вязкости в двойных логарифмических координатах; б - кривые течения в корневых координатах / Curves of nanofluids based on titanium oxide nanoparticles in water at a temperature of 25° C and volumetric concentrations Ф= 0.12 (1), 0.10 (2), 0.08 (3), and 0.05 (4): (а) viscosity curves in double logarithmic coordinates; (b) flow curves in root coordinates Источник: [8].
По существующим представлениям [5] присутствие частиц в несущей «ньютоновской» жидкости увеличивает вязкость системы в результате возмущения сдвигового течения вращающимися частицами, а в особенности агрегатами частиц. В теории наножидкостей принято наименование - нанокластеры. Поэтому вязкость наножидкостей увеличивается с ростом концентрации дисперсной фазы. Значительное уменьшение сдвиговой вязкости при увеличении скорости сдвига (сдвиговое разжижение) объясняют разрушением нанокластеров при высоких напряжениях сдвига. Обычно полагают, что при достаточно низких концентрациях зависимость вязкости от скорости сдвига описывается степенным законом (уравнением Оствальда).
В рамках структурной реологической модели [11] вязкость увеличивается пропорционально количеству агрегированных частиц, т. е. возрастает с концентрацией. Разрушение агрегатов частиц объясняется действием разрывающих гидродинамических сил, сходных по природе с силами Стокса. Судя по виду кривых вязкости (рис. 1, а), степенной закон не выполняется достаточно строго при всех рассмотренных концентрациях, кроме того, имеются «ньютоновские» участки при низких или при очень высоких скоростях сдвига.
W
Аппроксимация на выделенных интервалах скоростей сдвига проводилась нами с помощью обобщённого уравнения течения [11].
Т 1/2у1/2
т1/2 +V2 у1/2. (1)
<у!/2 + X
Физический смысл коэффициента компактности коэффициента агрегации Тс1/2 и коэффициента вязкости Кэссона П1/2с описан в работе [11]. Результаты аппроксимации удобно представить в корневых координатах (рис. 1, б), поскольку рассчитанная кривая переходит в прямолинейную зависимость при высоких скоростях сдвига. Коэффициенты реологического уравнения представлены в табл. 1.
Суммарная вязкость неагрегированной системы описывается коэффициентом П, и, в целом, уменьшается с уменьшением концентрации, т. е. количества частиц в единице объёма. Структурная часть вязкости ТС/2 / X обусловлена количеством и величиной агрегатов и также уменьшается с уменьшением концентрации. Ньютоновское поведение соответствует динамическому равновесию структуры, при котором количество агрегированных частиц не зависит от скорости сдвига.
Таблица 1 / Table 1
Коэффициенты обобщённого уравнения течения, величина корня предельной нулевой вязкости, корень из структурной вязкости, рассчитанные для наножидкости на основе наночастиц оксида титана в воде при температуре 25° С и различных объёмных концентрациях / Coefficients of the generalized flow equation, the value of the root of the limiting zero viscosity, and the root of the structural viscosity calculated for an aqueous nanofluid based on titanium oxide nanoparticules at a temperature of 25° C and various volume fractions
концентрация 0,12 0,10 0,08 0,05 0,018
Т1/2, Па1/2 10,41 5,63 4,24 2,70 0,125
пС/2,(Пас)1/2 0,088 0,107 0,075 0,047 0,030
X , с-1/2 0,265 0,371 0,567 0,451 0
т1/2/ X 39,35 15,18 7,48 5,98 -
П1/2(0), (Па с)1/2 39,44 15,29 7,55 6,03 -
Источник: по данным авторов.
При очень низких концентрациях наночастиц ньютоновское поведение при высоких скоростях более выражено (рис. 2, а). В работе [9] использовались на-ночастицы ТЮ2 размером около 20 нм, диспергированные в воде при объёмной концентрации 1,18 об.%.
о -
-2 -
lgп, Пас
уравнение 1
-2
-3
lgп, Пас
о о
II скачок скорости 1 уменьшением вязкости 0_ lg т, Па
чхоэоррсооосоооо.
-1
о
а
б
о
0
Рис. 2 / Fig. 2. Реологические кривые наножидкости на основе наночастиц оксида титана в воде при объёмной концентрациях Ф= 0,018: а - зависимость вязкости от скорости сдвига; б - зависимость вязкости от напряжения сдвига / Rheological curves of a nanofluid based on titanium oxide nanoparticles in water at
volume concentrations Ф= 0.018: (а) dependence of viscosity on shear rate; (b) dependence of viscosity on shear stress Источник: [9].
При такой низкой концентрации уравнение (1) можно использовать только на малом интервале низких скоростей сдвига. В дальнейшем, при крайне малом увеличении напряжения сдвига т происходит резкое увеличение скорости сдвига Y (и соответствующее резкое уменьшение вязкости), которое связывают со сдвиговым расслоением вещества. Этот эффект хорошо виден на кривой п(т) в двойных логарифмических координатах (рис. 2, б) и на кривой т( Y) в корневых координатах (рис. 3). Скачкообразный переход от одного режима течения к другому показан на рис. 2, а пунктирной линией (п ~1/ Y). Участок высоких скоростей сдвига можно считать «ньютоновским», поскольку зависимости п(т) и П( Y) примерно параллельны оси абсцисс. Экспериментально полученные значения вязкости П превышают рассчитанную предельную вязкость Пс максимально разрушенной структуры (штрихпунктирная прямая на рис. 2, а).
1
0,8 0,6 0,4 0,2 0
Ц1/2,Па1/2
' уравнение 1 скачок скорости
y = 0,0455x - 0,1469 R2 = 0,9994 уравнение 2 ^
y = 0,0387x - 0,0171 R2 = 0,9997
T уравнение 2
•,1/2 -1/2 Y ,с
0
5
10
15
20
25
Рис. 3 / Fig. 3. Кривая течения в корневых координатах для наножидкости на основе наночастиц оксида титана в воде при объёмной концентрациях Ф= 0,018 / Flow curve in root coordinates for a nanofluid based on titanium oxide nanoparticles in water at volume
concentrations Ф= 0.018
Источник: [9].
Однако, на рис. 3 можно видеть, что в корневых координатах с достаточно хорошей точностью выполняется линейная зависимость вида
ТШ =л1/2уШ -Т1/2 . (2)
В рамках структурной реологической модели [12] коэффициенты этого реологического уравнения имеют следующий вид:
= П~2 +BNk з/(кз +k0; TÍ = BN (k, - k2)/k + h). (3)
В модели предполагается, что агрегаты и индивидуальные частицы являются элементами структуры, которые обеспечивают вязкость при сдвиговом течении. Если в единице объёма содержится N частиц, то количество агрегированных частиц обозначим как N2, а количество индивидуальных частиц как Ni.
Выражение (3) содержит следующие коэффициенты: k 2 — константа скорости формирования агрегатов при случайных столкновениях частиц; ко - константа скорости спонтанного разрушения агрегатов, например, в результате теплового движения; k1 - константа скорости разрушения агрегата под действием растягивающих гидродинамических сил, к3 - константа скорости формирования агрегатов частиц из одиночных частиц под действием сдвига. Коэффициент B не зависит от скорости сдвига.
Коэффициент TÍ является положительным, если к0 > к2, и отрицательным, если к0 < к2. Поэтому на данном интервале скоростей возможны три вида реологического поведения. В первом случае (к0 > к2,) вязкость увеличивается с ростом скорости сдвига (явление сдвигового затвердевания). Во втором случае (к0 < к2) вязкость уменьшается с ростом скорости сдвига, причём коэффициент т^2 меньше или сравним по величине с коэффициентом ní. В третьем случае коэффициент TÍ близок к нулю при условии ко ~ к2, тогда реологическое поведение
похоже на «ньютоновское» течение с практически постоянным значением сдвиговой вязкости Цсу .
Наножидкость на основе оксида титана в этиленгликоле
Реологическое поведение наножидкости на основе оксида титана в этиленгликоле рассмотрено в работе [10]. Реологические кривые несколько отличаются по виду от тех, что получены для водной наножидкости. Однако, на графиках (рис. 4, 5) имеются аналогичные участки кривых: режим, близкий к «ньютоновскому»; участок снижение вязкости за счёт разрушения агрегатов (уравнение 1). Кроме того, на рис. 4 жирной линией выделено реологическое поведение, описанное уравнением (2). Штриховая горизонтальная линия изображает вязкость этиленгликоля.
lgП, Пас
уравнение 1
равнение
_1_I_I_I_I_1_
-1 -0,5 0 0,5 1 1,5 2 2,5 3 lg у, с
-1
0
Рис. 4 / Fig. 4. Кривые вязкости в двойных логарифмических координатах для наножидкости на основе наночастиц оксида титана в этиленгликоле при массовых концентрациях 25% (1), 15% (2), 10% (3), 5% (4) / Viscosity curves in double logarithmic coordinates for a nanofluid based on titanium oxide nanoparticles in ethylene glycol at mass concentrations of 25% (1), 15% (2), 10% (3), 5% (4)
Источник: [10].
Результаты аппроксимации более ясно выражены на графиках в корневых координатах (рис. 5). Рассчитанная по уравнению (1) зависимость асимптотически приближается к прямой при высоких скоростях сдвига (рис. 5, а). Структура наножидкости постепенно разрушается при увеличении скорости сдвига (ki>0, k3=0). Коэффициенты обобщённого уравнения течения приведены в табл. 2.
Прямая, соответствующая уравнению (2), пересекает ось ординат при положительных значениях (рис. 5, б). Вязкость уменьшается с ростом скорости
сдвига; коэффициент т^2 меньше коэффициента п"2. В этом случае одновременно происходят два процесса: формирование агрегатов под действием сдвига и разрушение агрегатов под действием сдвига, но доминирует процесс разрушения (ко < к2, к1 > 0, кз > 0).
т1/2,Па1/2
уравнение 1
10
15 20
а
25
■X> 1 0,25
^ 3 0,2
4
0,15
0,1
1/2 -1 / 2
У ,с 0,05
30 35 0
' уравнение2
Т1/2,Па1/2
уравнение 1
Р4
У = 0,1415x + 0,0619 R2 = 0,9946
•,1/2 -1/2 Y ,с
-I_I_I_L.
0 0,2 0,4 0,6 0,8 1 1,2 1,4 б
Рис. 5 / Fig. 5. Кривые течения в корневых координатах для наножидкости на основе наночастиц оксида титана в этиленгликоле при массовых концентрациях 25% (1),
15% (2), 10% (3), 5% (4): а - кривые течения на полном интервале скоростей сдвига; б - кривая течения для массовой концентрации 5% на интервале низких скоростей / Flow curves in root coordinates for a nanofluid based on titanium oxide nanoparticles in ethylene glycol at mass concentrations of 25% (1), 15% (2), 10% (3), 5% (4): (a) flow curves over the full range of shear rates; (b) flow curve for 5% mass concentration in
the low velocity range
Источник: [10].
0
5
Таблица 2 / Table 2
Коэффициенты обобщённого уравнения течения, величина корня предельной нулевой вязкости, корень из структурной вязкости, рассчитанные для наножидкости на основе наночастиц оксида титана в этиленгликоле при различных массовых концентрациях / Coefficients of the generalized flow equation, the value of the root of the limiting zero viscosity, and the root of the structural viscosity calculated for an aqueous nanofluid based on titanium oxide nanoparticules at various weight concentrations
концентрация 25% 15% 10% 5%
ТГ, Паш 1,430 0,603 0,371 0,162
Щ2 ,(Па с)1/2 0,140 0,128 0,125 0,123
X , с-1/2 0,055 0,239 0,945 0,917
концентрация 25% 15% 10% 5%
тС2/ X 26,10 2,52 0,393 0,176
П1/2(0), (Па с)1/2 26,24 2,65 0,518 0,299
Источник: по данным авторов.
Выводы
Рассмотрено реологическое поведение наножидкостей на основе наночастиц оксида титана в воде и в этиленгликоле. Аппроксимация реологических кривых в случае стационарного сдвигового течения проведена с помощью уравнений структурной реологической модели. Обнаружены характерные особенности реологического поведения, которые встречаются в других структурированных жидкостях: участок ньютоновского течения; сдвиговое разжижение; сдвиговое затвердевание, эффект резкого увеличения скорости сдвига при незначительном увеличении напряжения сдвига.
Статья поступила в редакцию 14.10.2022 г. ЛИТЕРАТУРА
1. Mahbubul I. M., Saidur R., Amalina M. A. Latest developments on the viscosity of nanoflu-ids // International Journal of Heat and Mass Transfer. 2012. Vol. 55. Iss. 4. P. 874-885. DOI: 10.1016/J.IJHEATMASSTRANSFER.2011.10.021.
2. Рудяк В. Я. Современное состояние исследований вязкости наножидкостей // Вестник Новосибирского государственного университета. Серия: Физика. 2015. Т. 10. № 1. С. 5-22.
3. Haisheng Chen, Yulong Ding, Chunqing Tan. Rheological behaviour of nanofluids // New Journal of Physics. 2007. Vol. 9. P. 367-391. doi:10.1088/1367-2630/9/10/367.
4. Investigation of thermal conductivity and rheological properties of nanofluids containing graphene nanoplatelets / Mehrali M., Sadeghinezhad E., Tahan Latibari S., Kazi S. N., Mehrali M., Zubir M. N., Metselaar H. S. C. // Nanoscale Research Letters. 2014. Vol. 9. P. 15-24. DOI: 10.1186/1556-276X-9-15.
5. Hasan S. M. Non-Newtonian rheological characteristics of oil-based metal oxide nanoflu-ids: a thesis for the degree master of science. Department of Mechanical Engineering, Northern Illinois University. De Kalb, Illinois, 2017. 55 p.
6. Saeedi A. H., Akbari M., Toghraie D. An experimental study on rheological behavior of a nanofluid containing oxide nanoparticle and proposing a new correlation // Physica E: Low-dimensional Systems and Nanostructures. 2018. Vol. 99. P. 285-293. DOI: 10.1016/j.physe.2018.02.018.
7. Rheological non-Newtonian behaviour of ethylene glycol-based Fe2O3 nanofluids / Pastoriza-Gallego M. J., Lugo L., Legido J. L., Pineiro M. M. // Nanoscale Research Letters. 2011. Vol. 6. P. 560-566. DOI: 10.1186/1556-276X-6-560.
8. Tseng W. J., Lin K. Rheology and colloidal structure of aqueous TiO2 nanoparticle suspensions // Materials Science and Engineering: A. 2003. Vol. 355. Iss. 1-2. P. 186-192. DOI: 10.1016/S0921-5093(03)00063-7.
9. Heat transfer and flow behavior of aqueous suspensions of TiO2 nanoparticles (nanofluids) flowing through a vertical pipe / He Y., Chen H., Ding Y., Cang D., Lu H. // International
Journal of Heat and Mass Transfer. 2007. Vol. 50. Iss. 11-12. P. 2272-2281. DOI: 10.1016/j.ijheatmasstransfer.2006.10.024.
10. Rheological and volumetric properties of TiO2-ethylene glycol nanofluids / Cabaleiro D., Pastoriza-Gallego M. J., Gracia-Fernandez C., Pineiro M. M., Lugo L. // Nanofluids. Na-noscale Research Letters. 2013. Vol. 8. Article number: 286. DOI: 10.1186/1556-276X-8-286.
11. Кирсанов Е. А., Матвеенко В. Н. Неньютоновское течение дисперсных, полимерных и жидкокристаллических систем. Структурный подход. М.: Техносфера, 2016. 384 с.
12. Кирсанов Е. А., Матвеенко В. Н. Вязкость и упругость структурированных жидкостей. М.: Техносфера, 2022. 284 с.
REFERENCES
1. Mahbubul I. M., Saidur R., Amalina M. A. Latest developments on the viscosity of nanoflu-
ids. In: International Journal of Heat and Mass Transfer, 2012, vol. 55, iss. 4, pp. 874-885. DOI: 10.1016/J.IJHEATMASSTRANSFER.2011.10.021.
2. Rudyak V. Ya. [State-of-the-art study of nanofluid viscosity]. In: Vestnik Novosibirskogo
gosudarstvennogo universiteta. Seriya: Fizika [Vestnik NSU. Series: Physics], 2015, vol. 10, no. 1, pp. 5-22.
3. Haisheng Chen, Yulong Ding, Chunqing Tan. Rheological behaviour of nanofluids. In: New
Journal of Physics, 2007, vol. 9, pp. 367-391. DOI: 10.1088/1367-2630/9/10/367.
4. Mehrali M., Sadeghinezhad E., Tahan Latibari S., Kazi S. N., Mehrali M., Zubir M. N., Met-
selaar H. S. C. Investigation of thermal conductivity and rheological properties of nanofluids containing graphene nanoplatelets. In: Nanoscale Research Letters, 2014, vol. 9, pp. 15-24. DOI: 10.1186/1556-276X-9-15.
5. Hasan S. M. Non-Newtonian rheological characteristics of oil-based metal oxide nanoflu-
ids: a thesis for the degree master of science. Department of Mechanical Engineering, Northern Illinois University. De Kalb, Illinois, 2017. 55 p.
6. Saeedi A. H., Akbari M., Toghraie D. An experimental study on rheological behavior of a
nanofluid containing oxide nanoparticle and proposing a new correlation. In: Physica E: Low-dimensional Systems and Nanostructures, 2018, vol. 99, pp. 285-293. DOI: 10.1016/j.physe.2018.02.018.
7. Pastoriza-Gallego M. J., Lugo L., Legido J. L., Piñeiro M. M. Rheological non-Newtonian
behaviour of ethylene glycol-based Fe2O3 nanofluids. In: Nanoscale Research Letters, 2011, vol. 6, pp. 560-566. DOI: 10.1186/1556-276X-6-560.
8. Tseng W. J., Lin K. Rheology and colloidal structure of aqueous TiO2 nanoparticle suspen-
sions. In: Materials Science and Engineering: A, 2003, vol. 355, iss. 1-2, pp. 186-192. DOI: 10.1016/S0921-5093(03)00063-7.
9. He Y., Chen H., Ding Y., Cang D., Lu H. Heat transfer and flow behavior of aqueous sus-
pensions of TiO2 nanoparticles (nanofluids) flowing through a vertical pipe. In: International Journal of Heat and Mass Transfer, 2007, vol. 50, iss. 11-12, pp. 2272-2281. DOI: 10.1016/j.ijheatmasstransfer.2006.10.024.
10. Cabaleiro D., Pastoriza-Gallego M. J., Gracia-Fernandez C., Pineiro M. M., Lugo L.
Rheological and volumetric properties of TiO2-ethylene glycol nanofluids. In: Nanofluids. Nanoscale Research Letters, 2013, vol. 8, article number: 286. DOI: 10.1186/1556-276X-8-286.
11. Kirsanov Ye. A., Matveenko V. N. Nen'yutonovskoe techenie dispersnykh, polimernykh i
zhidkokristallicheskikh sistem. Strukturnyi podkhod [Non-Newtonian flow of dispersed, polymeric and liquid crystal systems]. Moscow, Tekhnosfera Publ., 2016. 384 p.
12. Kirsanov Ye. A., Matveyenko V. N. Vyazkost' i uprugost' strukturirovannykh zhidkostey [Viscosity and elasticity of structured liquids]. Moscow, Tekhnosfera Publ., 284 p.
ИНФОРМАЦИЯ ОБ АВТОРАХ
Вековищев Михаил Петрович - кандидат физико-математических наук, доцент кафедры физики и химии Государственного социально-гуманитарного университета; e-mail: mpv.71@mail.ru;
Кирсанов Евгений Александрович - кандидат физико-математических наук, доцент кафедры физики и химии Государственного социально-гуманитарного университета; e-mail: Kirsanov47@mail.ru.
INFORMATION ABOUT THE AUTHORS
Mikhail P. Vekovishchev - Cand. Sci. (Phys.-Math.), Assoc. Prof., State University of Humanities and Social Studies; e-mail: mpv.71@mail.ru;
Evgeny A. Kirsanov - Cand. Sci. (Phys.-Math.), Assoc. Prof., State University of Humanities
and Social Studies;
e-mail: Kirsanov47@mail.ru.
ПРАВИЛЬНАЯ ССЫЛКА НА СТАТЬЮ
Вековищев М. П., Кирсанов Е. А. Неньютоновское течение наножидкости на основе оксида титана // Вестник Московского государственного областного университета. Серия: Физика-математика. 2022. № 4. С. 6-16. DOI: 10.18384/2310-7251-2022-4-6-16.
FOR CITATION
Vekovishchev M. P., Kirsanov E. A. Non-Newtonian flow of a nanofluid based on titanium oxide. In; Bulletin of the Moscow Region State University. Series: Physics and Mathematics, 2022, no. 4, pp. 6-16.
DOI: 10.18384/2310-7251-2022-4-6-16.