Научная статья на тему 'О сдвиговой вязкости магнитных жидкостей'

О сдвиговой вязкости магнитных жидкостей Текст научной статьи по специальности «Физика»

CC BY
405
103
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
магнитная жидкость / КОЭФФИЦИЕНТ ТРЕНИЯ / сдвиговая вязкость / Magnetic liquid / Friction coefficient / shear viscosity

Аннотация научной статьи по физике, автор научной работы — Комилов К., Зарипов А. К.

Исследована сдвиговая вязкость магнитных жидкостей на основе керосина и додекана, когда коэффициент трения магнитных жидкостей является функцией термодинамических параметров состояния. Проведён численный расчёт коэффициента трения, времён релаксации и, а также коэффициента сдвиговой вязкости магнитных жидкостей на основе керосина и додекана в зависимости от температуры, плотности и концентрации. Полученные результаты находятся в удовлетворительном согласии с существующими экспериментальными данными.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Shear viscosity of magnetic liquids on the basis of kerosene and dodecane when the coefficient of the friction of magnetic liquids is functions of thermodynamic parameters of conditions is investigated. Is carried numerical calculation of coefficient of the friction, relaxation times, and also coefficient of shear viscosity of the magnetic liquids on the basis of kerosene and dodecane depending on temperature, density and concentration. The received results are in the satisfactory consent with existing experimental data.

Текст научной работы на тему «О сдвиговой вязкости магнитных жидкостей»

ДОКЛАДЫ АКАДЕМИИ НАУК РЕСПУБЛИКИ ТАДЖИКИСТАН _______________________________________2012, том 55, №5____________________________________

ФИЗИКА

УДК 532.7+537.84

К.Комилов, А.К.Зарипов

О СДВИГОВОЙ ВЯЗКОСТИ МАГНИТНЫХ ЖИДКОСТЕЙ

Таджикский национальный университет

(Представлено академиком АН Республики Таджикистан С.Одинаевым 07.02.2012 г.)

Исследована сдвиговая вязкость магнитных жидкостей на основе керосина и додекана, когда коэффициент трения магнитных жидкостей является функцией термодинамических параметров состояния. Проведён численный расчёт коэффициента трения Р, времён релаксации Тх и Т0, а

также коэффициента сдвиговой вязкости магнитных жидкостей на основе керосина и додекана в зависимости от температуры, плотности и концентрации. Полученные результаты находятся в удовлетворительном согласии с существующими экспериментальными данными.

Ключевые слова: магнитная жидкость - коэффициент трения - сдвиговая вязкость.

Благодаря таким уникальным свойствам, присущим только магнитным жидкостям, как текучесть, сжимаемость жидкой среды и значительная намагниченность, появляется возможность использовать их в различных областях техники и технологии. Например, в ультразвуковой дефектоскопии в качестве акустических контактных сред всё чаще используются магнитные жидкости, пропускающие не только продольные, но и поперечные ультразвуковые волны. Особенности прохождения и поглощения ультразвука в магнитных жидкостях определяются набором физических механизмов потерь и преобразования упругой энергии, к которым относятся вязкие потери, внутренний теплообмен, магнитогидродинамические процессы, дисперсия и другие нелинейные эффекты.

Коэффициент сдвиговой вязкости, который присутствует в любой жидкости, очень хорошо описывает механизм поглощения ультразвуковых волн и другие релаксационные процессы и его можно измерять экспериментально. Вязкостные свойства магнитных жидкостей исследуются различными методами, в том числе методом молекулярно-кинетической теории, который позволяет более существенно объяснить механизмы протекания ряда процессов в жидкости, а также их релаксации.

Поэтому можно заключить, что исследование явлений переноса, в частности определение коэффициента сдвиговой вязкости магнитных жидкостей методом молекулярно-кинетической теории, является актуальной задачей и оценка этих коэффициентов позволяет упростить расчёты технических устройств, а также предусмотреть перспективы их применения в различных областях технологии, медицины, сельском хозяйстве и т. д.

Целью данной работы являлась разработка математической модели магнитной жидкости, позволяющей определить коэффициент трения Р магнитных жидкостей, времена релаксации хх, т0, зависящие от параметров состояния, а также провести численные расчёты коэффициента сдвиговой

Адрес для корреспонденции: Зарипов Авзалшо Карамонович. 734025, Республика Таджикистан, г. Душанбе, пр. Рудаки, 17, Таджикский национальный университет. E-mail: afzal.z@mail.ru

вязкости % магнитных жидкостей в зависимости от плотности, концентрации и температуры. Полученные результаты демонстрируют повышение точности расчёта этих коэффициентов.

При определении коэффициента трения /3 встречаются некоторые трудности, поскольку он тесно взаимосвязан с параметрами состояния и межмолекулярными силами в жидкостях. Впервые Кирквуд [1] установил соотношение между коэффициентом трения и межмолекулярными силами в среде. Были предприняты попытки получения удовлетворительного соотношения для коэффициента трения 3 в различных приближениях. В [1] было предложено несколько моделей определения коэффициента трения 3.

Методом вычисления автокорреляционных функций силы было получено выражение для коэффициента трения 3 в виде [2]:

циал Леннард-Джонса с сильным отталкивательным членом как потенциал твёрдых сфер, т - масса частицы, п - числовая плотность.

Также в [2] было предложено выражение для коэффициента трения 3 на основе флуктуаци-онно-диссипативной теоремы, связывающей диссипативный коэффициент 3 с интегралом от автокорреляционной функции сил, действующих на молекулу со стороны окружающих молекул

Для проведения численного расчета 3 с помощью выражения (1) необходимо выбирать модельный потенциал взаимодействия частиц Ф(г) и радиальной функции распределения g(г) . В качестве потенциалов межчастичных взаимодействий были выбраны потенциалы, предложные в [3] и [4] соответственно:

(1)

радиальная часть оператора Лапласа, Ф(г) - «модифицированный» потен-

где С(Яс) = (а/Яс )12 - (а/Яс) и Яс = 21/ба.

В качестве радиальной функции распределения g (г) было использовано выражение, приведённое в [5]:

ё(г ) = у(р*)ехр[-], (4)

к!

где у(р ) = (2 - р ) / 2(1 - р )3 - контактная функция Карнахана-Старлинга,

р = жО п /6 = тгрЫ0О / 6М - приведённая плотность, Nо - число Авогадро, М - молярная масса.

Учитывая равенства (2)-(4) в выражении (1), приводим его к удобному виду для проведения численных расчётов:

А2 = 11 р * |г ~6 (22г -6 - 5) gl (г1 , (5)

1 2

ад

На основе формул (5) и (6) проведён численный расчёт. В табл. 1 представлены результаты численного расчёта коэффициента трения А магнитной жидкости на основе керосина и додекана с частицами магнетита Fe3O4 при T=297 К и различных концентрациях без учёта вклада внешнего магнитного поля. Значения для температуры, плотности, концентрации и намагниченности насыщения заимствованы из [6] и [7].

Таблица 1

Зависимость коэффициента трения магнитных жидкостей на основе керосина и додекана от термодинамических параметров состояния

Т, К р, кг/м3 ф, % М^10-3, А/м Рг1012, кг/с Т11-1013, с Т01-109, с Р2-1012, кг/с І12-1013, с Т02-109, с

Магнитная жидкость на основе керосина

[6] формула (5) формула (6)

297 800 0 0 0.703 2.670 0.050 0.811 2.314 0.057

860 1.37 6.3 0.744 2.530 0.052 0.858 2.192 0.061

970 3.85 14.6 0.820 2.301 0.058 0.946 1.995 0.067

1090 6.35 23 0.907 2.087 0.064 1.046 1.809 0.074

1140 7.94 27.9 0.943 2.007 0.067 1.089 1.740 0.077

1230 9.75 36.9 1.013 1.872 0.072 1.169 1.622 0.083

Магнитная жидкость на основе додекана

[7] формула (5) формула (6)

273 1137.3 13.3 35 4.741 0.345 0.700 5.358 0.305 0.792

293 1118.2 12.6 34 4.517 0.361 0.622 5.011 0.325 0.690

303 1099.0 11.8 32.5 4.286 0.378 0.570 4.716 0.344 0.628

323 1079.9 10.9 28.8 4.120 0.392 0.514 4.468 0.362 0.558

Из табл. 1 видно, что значения коэффициента трения А с увеличением концентрации в отсутствии внешнего магнитного поля возрастает. Из-за отсутствия экспериментальных данных по опре-

делению коэффициента трения ( мы не имеем возможности провести сравнение результатов численного расчёта с экспериментальными результатами.

Для проведения численного расчёта зависимости сдвиговой вязкости % магнитной жидкости

от параметров состояния воспользуемся ранее полученным выражением для коэффициента сдвиговой вязкости [8]:

, ч пкТт 2жп а г 3 дФ г , ч 5^0 / (дИ)

гА®) =--------+-----------I ёгг----I О(г,г ,®) 1 + -!-2-1МVII---

^ ( 1 1 + ®Т )2 15 { дг 1 Г )_ 2(( \дЗ)рт _

Г/Я- , (7)

дг

Т ( 2 ^ ^

где О (г, Г1,«) = Т1 --- [шф-СОБф) ехр (-ф)-(8тф2 - ) ехр (-ф )], ф (г, Г, а>) =

2 ®о ) '

= (®т0 / 2)12 (г + г ), Т = (а1 / 2кТ, т1 = ш / 2( , ^0 - магнитная проницаемость вакуума.

В отсутствии внешнего магнитного поля, выражение (8) принимает следующий вид:

, ч пкТт 2жп2аъ Г з дФ \ ( ,дgo „

г, (®) =------Цт +---------I Лт — I О, (г, г,,®)—0 . (8)

/Л ® 1 + (®т)2 1 5 I дг I ^ , 1 , Удг 1 1 ()

Согласно выражению (8), сдвиговая вязкость магнитной жидкости описывается как трансляционной, так и структурной релаксацией. Времена структурной релаксации характеризуется коэффициентом трения магнитной жидкости ( .

На основе выражения (8) с учётом (5) проведён численный расчёт. В табл. 2 представлены результаты численного расчёта зависимости сдвиговой вязкости г г магнитной жидкости на основе керосина с частицами магнетита Бе304 при Т=297 К от плотности и концентрации.

Таблица 2

Зависимость коэффициента сдвиговой вязкости магнитной жидкости на основе керосина от температуры, плотности и концентрации.

Т, К р , кг/м3 ф, % г,, Па ■ с г,, Па ■ с

[6] формула (8)

297 800 0 0.0013 0.0008

860 1.37 0.0014 0.0010

970 3.85 0.0017 0.0013

1090 6.35 0.0023 0.0018

1140 7.94 0.0028 0.0020

1230 9.75 0.0037 0.0025

В табл. 3 представлены результаты численного расчёта зависимости сдвиговой вязкости г магнитной жидкости на основе додекана от температуры, плотности и концентрации.

Таблица 3

Зависимость коэффициента сдвиговой вязкости магнитной жидкости на основе додекана от температуры, плотности и концентрации.

T, К p , кг/м3 9, % Г, Па ■ с rs, Па ■ с

[7] формула (8)

273 1137.3 13.3 2.264 1.744

293 1118.2 12.6 1.492 1.308

303 1099.0 11.8 1.064 1.040

323 1079.9 10.9 0.803 0.813

Сравнение результатов, приведённых в табл. 3, проиллюстрировано на рисунке, где приведена графическая зависимость коэффициента сдвиговой вязкости магнитной жидкости на основе доде-кана от концентрации.

Рис. Зависимость сдвиговой вязкости магнитной жидкости на основе додекана от концентрации:

I - расчёт по формуле (8), II - данные работы [7].

Из сравнения кривых видно, что значения коэффициента сдвиговой вязкости магнитной жидкости на основе додекана с увеличением плотности и концентрации возрастают. С возрастанием температуры значение сдвиговой вязкости уменьшается. Согласно приведённым данным, полученные теоретические результаты на основе выражения (8) хорошо согласуются с результатами работы [7]. При малых концентрациях эти результаты с результатами работы [7] почти совпадают, а при высоких концентрациях наблюдается их отклонение.

Таким образом, можно заключить, что теоретическое исследование магнитных жидкостей методом молекулярно-кинетической теории, который учитывает все необратимые процессы внутри жидкости, позволяет повысить точность численных расчётов, которые можно успешно использовать как для современной теории жидкостей, так и для многочисленных технических и технологических применений.

Поступило 07.02.2012 г.

ЛИТЕРАТУРА

1. Rice S.A., Gray P. Statistical mechanics of simple liquids. - New York, 1965.

2. Грэй П. Физика простых жидкостей, т.1 /Под ред. Г.Темперли и др. - М.: Мир, 1971, с. 136-192.

3. Zuowei Wang, Christian Holm, Hanns Walter Müller. - Phys. Rev., 2002, E 66, pp. 021405-1-02140513.

4. Роулинсон Дж. - Физика простых жидкостей, т. 1 / Под ред. Г.Темперли и др. - М.: Мир, 1971, с. 63-80.

5. Юхновский И.Р., Головко М.Ф. Статистическая теория классических равновесных систем. - Киев: Наукова думка, 1980, 372 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

6. Полунин В. М. Акустические эффекты в магнитных жидкостях. - М.: Физматлит, 2008, 207 с.

7. Виноградов А.Н.- Вестн. Моск. университета. Сер. 2. Химия, 1999, т. 40, № 2, с. 90-93.

8. Одинаев С., Комилов К., Зарифов А. - ЖФХ, 2006, т.80, № 5, с. 864-871.

^.Комилов, А.К.Зарипов

ОИД БА ЧАСПАКИИ ЛАГЗИШИ МОЕЪ^ОИ МАГНИТЙ

Донишго^и миллии Тоцикистон

Дар макола часпакии лагзиши моеъх,ои магнитй дар асоси керосин ва додекан, мавриди функсияи параметрх,ои термодинамикии х,олат будани зариби соиши моеъх,ои магнитй, тахдик карда шудааст. Х,исобх,ои ададии вобастагии зариби соиш /3, вактх,ои релаксатсия хх ва х0,

инчунин зариби часпакии лагзиши моеъхри магнитй дар асоси киросин ва додекан аз хдрорат, зиччй ва консентратсия гузаронида шудааст. Натичах,ои хрсилшуда бо маълумотх,ои тачрибавии мавчуда ба таври каноатбахш мувофикат мекунанд.

Калима^ои калиди: моеъи магнитй - зариби соиш - часпакии лагзиш.

K.Komilov, A.K.Zaripov ABOUT SHEAR VISCOSITY OF MAGNETIC LIQUIDS

Tajik National University Shear viscosity of magnetic liquids on the basis of kerosene and dodecane when the coefficient of the friction of magnetic liquids is functions of thermodynamic parameters of conditions is investigated. Is carried numerical calculation of coefficient of the friction 3, relaxation times xx, x0 and also coefficient of

shear viscosity of the magnetic liquids on the basis of kerosene and dodecane depending on temperature, density and concentration. The received results are in the satisfactory consent with existing experimental data.

Key words: magnetic liquid - friction coefficient - shear viscosity.

i Надоели баннеры? Вы всегда можете отключить рекламу.