CC BY
13ISSN 1606-867Х (Print) ISSN 2687-0711 (Online)
■¿яеяеЕязэн^^И ИВнМм!11 иTüifiiJ
Конденсированные
https://journals.vsu.ru/kcmf/
среды и межфазные границы
Оригинальные статьи
Научная статья УДК 538.956
https://doi.org/10.17308/kcmf.2021.23/3683
Экспериментальное исследование и математическое моделирование автоколебательных процессов на межфазной границе электрод - магнитная жидкость в электрическом поле
В. С. Чекановш, Е. В. Кириллова2, А. В. Коваленко 3, Е. Н. Дискаева4
1Северо-Кавказский федеральный университет,
ул. Пушкина, 1, Ставрополь 355017, Российская Федерация
2Университет Прикладных Наук Рейн-Майн Kurt-Schumacher-Ring, 18, Висбаден 65197, Германия
3Кубанский государственный университет,
ул. Ставропольская, 149, Краснодар 350040, Российская Федерация
4МИРЭА-Российский технологический университет, пр. Кулакова, 8, Ставрополь 355000, Российская Федерация
В статье описана математическая модель автоколебаний в виде краевой задачи для нелинейной системы уравнений в частных производных, получено численное решение и показана адекватность модели путем сравнения с экспериментом. В модели использовалась классическая система дифференциальных уравнений материального баланса, Нернста-Планка и Пуассона без упрощений и подгоночных параметров. Целью статьи было исследование параметров автоколебаний концентрации в слое частиц дисперсной фазы магнитной жидкости на границе с электродом в электрическом поле на основе разработанной математической модели и ее обоснование с точки зрения физического механизма.
В результате численных экспериментов было определено критическое значение скачка потенциала, при котором начинаются автоколебания, найдено время установления периода колебаний и другие характеристики процесса. Разработан программный продукт «AutoWaveOl» с интуитивно понятным пользовательским интерфейсом и расширенным функционалом, предназначенный для изучения автоколебаний в тонком слое магнитного коллоида. Ключевые слова: магнитная жидкость, межфазная граница, приэлектродный слой, электрическое поле, автоколебания, математическая модель
Благодарности: исследование выполнено в рамках работы Международной лаборатории компьютерного и математического моделирования нелинейных процессов «Computer and mathematical modeling of nonlinear processes».
Для цитирования: Чеканов В. С., Кирилова Е. В., Коваленко А. В., Дискаева Е. Н. Экспериментальное исследование и математическое и моделирование автоколебательных процессов на межфазной границе электрод - магнитная жидкость в электрическом поле. Конденсированные среды и межфазные границы. 2021;23(4): 626-636. https://doi. org/10.17308/kcmf.2021.23/3683
For citation: Chekanov V. S., Kirillova E. V., Kovalenko A. V., Diskaeva E. N. Experimental study and mathematical modelling of self-oscillation at the electrode-magnetic fluid interface in an electric field. Kondensirovannye sredy i mezhfaznye granitsy = CondensedMatter and Interphases. 2021; 23 (4): 626-636. https://doi.org/10.17308/kcmf.2021.23/3683
И Чеканов Владимир Сергеевич, e-mail: [email protected] © Чеканов В. С., Кирилова Е. В., Коваленко А. В., Дискаева Е. Н., 2021
Контент доступен под лицензией Creative Commons Attribution 4.0 License.
Аннотация
1. Введение
Коллоидные растворы магнитных композиционных материалов, которые впоследствии получили название «магнитные жидкости (magnetic fluids)» были синтезированы в 60-е годы XX века и до сих пор привлекают внимание как экспериментаторов, так и теоретиков [1, 2].
Магнитные жидкости (МЖ) представляют собой ультрадисперсные устойчивые коллоиды ферро- или ферримагнитных однодоменных частиц, диспергированных в различных жидкостях и совершающих интенсивное броуновское движение. В качестве дисперсной фазы используют малые частицы таких металлов, как железо, кобальт, никель, гадолиний, их разнообразные ферриты, ферромагнитные окислы. Для предотвращения коагуляции коллоидного раствора, которая была бы неизбежной вследствие магнитного диполь-дипольного и ван-дер-ва-альсовского взаимодействий, и последующего укрупнения частиц в качестве стабилизаторов применяют поверхностно-активные вещества (ПАВ). Адсорбируясь на поверхности микрокристаллических дисперсных частиц, ПАВ образуют защитную оболочку, представляющую своеобразный структурно-механический барьер. Средний размер дисперсных частиц d ~ 10 нм, вследствие этого МЖ не расслаиваются и сохраняют свою однородность практически неограниченное время.
МЖ на основе диэлектрических жидкостей обладают уникальными электрическими и магнитооптическими свойствами, связанными, в том числе, и с приповерхностными явлениями. Одним из таких свойств МЖ является способность в электрическом поле образовывать на границе с электродом плотные слои, состоящие из частиц дисперсной фазы с защитной оболочкой. Эти слои существенно влияют на макроскопические свойства коллоидной системы и являются основой для получения материалов с заданными свойствами. Данное явление нашло применение в области отображения информации [3-5].
В приэлектродном слое наблюдалась и исследовалась самоорганизация - автоколебания и их распределенный аналог - автоволны [6].
Уникальность магнитной жидкости как экспериментальной среды состоит в том, что автоколебательный процесс в ней можно не только неограниченно долго наблюдать и исследовать, но и управлять им электрическим полем.
Математическое моделирование автоколебательного процесса позволит установить его фундаментальные закономерности и описать состояние процесса в любой точке пространства в произвольный момент времени.
В настоящей работе предложена математическая модель автоколебаний концентрации частиц МЖ, основанная на динамическом подходе, в котором процессы, происходящие в системе, определяются уравнениями в частных производных, описывающими законы сохранения материи и заряда. Модель не содержит подгоночных параметров. Химические реакции учитываются в граничных условиях.
2. Экспериментальная часть
При проведении эксперимента использовалась магнитная жидкость типа «магнетит в керосине». Объемная концентрация твердой фазы была рассчитана по плотности магнитной жидкости. Для МЖ плотностью р = 1,14 кг/м3 концентрация твёрдой фазы составила 3.4 об. %. Объемную концентрацию магнетита находили по величинам плотностей магнитной жидкости, дисперсионной среды (керосина) и дисперсной фазы (магнетита) по формуле:
р.....-р
Фт =-
Рт -Ро
где рМЖ, рО и рТ - плотности магнитной жидкости, керосина и магнетита. Плотность МЖ определялась весовым методом по выталкивающей силе, действующей на мерный цилиндр из латуни, погруженный в жидкость, и рассчитывалась по формуле:
р - Р2 Р =—, ¿V
где Р1 - вес цилиндра в воздухе; Р2 - вес цилиндра в жидкости; V - объем цилиндра; g - ускорение свободного падения.
Взвешивание цилиндра проводилось на аналитических весах марки ВЛР-200. Погрешность в определении р составляла ~ 0.1 %.
Относительная диэлектрическая проницаемость жидкости е = 2, удельная проводимость а = 3,8 • 10-7 Ом-м-1. Для определения диэлектрической проницаемости и проводимости прибором «Измеритель иммитанса LCR-821» (на частоте 1000 Гц) измерялись сопротивление и емкость ячейки с МЖ, затем вычислялись проводимость а и диэлектрическая проницаемость е. Средний размер частиц магнетита ~ 10 нм, тол-
В. С. Чеканов и др. Экспериментальное исследование и математическое моделирование...
щина защитной оболочки, состоящей из молекул олеиновой кислоты, ~1.5 нм [7, 8]
Для определения изменения отражательной способности (коэффициента отражения по терминологии [9]) от границы проводящий электрод - магнитная жидкость в электрическом поле была собрана экспериментальная установка, показанная на рис. 1. Установка представляет собой плоскопараллельную ячейку, состоящую из двух электродов. Один из электродов -стекло с проводящим прозрачным покрытием (1ТО). Для устранения засветки от поверхности ячейки на неё устанавливается стеклянная призма. Между электродами помещается магнитная жидкость. На грань призмы, приклеенной к стеклу иммерсионной жидкостью, падает луч диаметром 2 мм от лазерного источника света (1 = 650 нм) так, что угол падения на поверхность стекла составляет 45 градусов. Лучи, отраженные от границы призма-стекло и от слоя проводящее покрытие - МЖ, пространственно разделены (5 мм). Для нивелирования паразитной засветки фотодиодов перед поляроидом 7 была помещена диафрагма с отверстием ~ 1 мм. Лазер и поляроид могут вращаться относительно луча. Плоскость поляризации поляроида 10 сов -падает с плоскостью поляризации луча лазера. Поляроид служит для уменьшения деполяризованной составляющей света лазера. Вращением лазера с поляроидом добиваются того, чтобы s-компонента (ТЕ-волна) была перпендикулярна плоскости падения. Вращением поляро-
Рис. 1. Схема экспериментальной установки для измерения коэффициента отражения на границе электрод - магнитная жидкость: 1 - магнитная жидкость; 2 - пластина из фольгированного гети-накса; 3,11 - изолирующие прокладки из полистирола; 4 - стекло с проводящим прозрачным покрытием; 5 - прямоугольная равнобедренная призма; 6 - диафрагма; 7, 10 - поляроиды; 8 - фотодиод; 9 - лазерный источник света; 12 - электроды
ида 7 ослабляли свет, чтобы фотодиод работал в линейном режиме. Площадь поверхности электродов составляла S = 36х30 мм2, толщина слоя магнитной жидкости - 250 мкм. В экспериментах использовались образцы стекол с проводящим покрытием 1ТО (изготовитель ООО «Политех», г. Санкт-Петербург).
Методика эксперимента заключалась в следующем: на электроды подается напряжение и, и вследствие электрофореза на них образуются плотные слои толщиной 100-200 нм, состоящие из частиц магнетита. Стекло с проводящим покрытием 1ТО освещается монохроматическим источником света. Отраженный от поверхностей стекло - проводящее покрытие и проводящее покрытие - МЖ луч интерферирует, попадает через диафрагму 6 и поляроид 7 на фотодиод ФД-256. Фотодиод соединен с входом двухлучевого осциллографа GDS-71022. Так как фотодиод работает в линейном режиме, то напряжение на осциллографе прямо пропорционально интенсивности отраженных лучей (коэффициенту отражения).
Когда скачок потенциала между электродами равен фкр, наблюдается периодическое изменение интенсивности отраженного луча (рис. 2, луч 2), что связано с изменением толщины при-электродного плотного слоя.
При освещении поверхности ячейки светом с широким спектром длин волн изменение толщины слоя сопровождается периодическим изменением цвета поверхности электрода от зеленого до малинового. Зеленый цвет поверхно-
Рис. 2. Автоколебания тока (1) и оптического сигнала (2) в ячейке с магнитной жидкостью. Период колебаний Т * 0.5 с, скачок потенциала Дф = 5 В
сти соответствует спектру отражения в отсутствии слоя частиц, а малиновый - максимальной толщине приповерхностного слоя частиц. Периодическое изменение цвета поверхности или периодическое изменение интенсивности отраженного света (рис. 2, луч 2) интерпретируется нами как автоколебания концентрации частиц. Электрооптический эффект, заключающийся в изменении спектра отраженного света (при освещении белым светом) или интенсивности (при освещении монохроматическим светом) от изменяющегося в электрическом поле тонкого приповерхностного слоя магнитной жидкости, называется электроуправляемой интерференцией [10].
2.1. Физическая модель автоколебаний частиц в ячейке с магнитной жидкостью
Общая характеристика электропроводности магнитной жидкости. Электропроводность магнитной жидкости обусловлена несколькими механизмами. Во-первых, она связана с примесными ионами, концентрация которых зависит от качества отмывки дисперсного магнетита после соосаждения в процессе приготовления МЖ:
FeCl2• + 4Н20 + 2FeCl3•6H2O + 8Ш4ОН ^ ^ Fe3O4j + 8Ш4С1 + 20Н20.
Вследствие неполной диссоциации соли в жидкости также содержатся ионные комплексы - ионные пары, тройники [11, 12].
Второй тип электропроводности обусловлен контактными процессами на поверхностях электродов. Эти процессы называют инжекци-онными [11-14]. В средних высоковольтных полях порядка 106 В/м развивается электрохимическая инжекция, которая обуславливает электрогидродинамические эффекты как в жидких диэлектриках, так и в магнитной жидкости [14].
Кроме ионной проводимости перенос зарядов в магнитной жидкости может осуществляться магнитными частицами (МЧ) [16].
Образование объемного заряда в магнитной жидкости в электрическом поле. В электрическом поле в магнитной жидкости у электродов начинает формироваться объемный заряд, с которым связано много специфических эффектов [17-19], в том числе и автоволновой процесс. Механизм формирования в жидком диэлектрике объемного заряда под действием электрического поля в окрестности проводящих поверхностей подробно описан. В [13, 19] приводится выражение для определения толщины приэлектродного слоя объемного заряда:
о E FD£0eE0(b + ¡b2)
о ~ E-т„ ~-,
0 RT M s '
(1)
где тм - время максвелл-вагнеровской релаксации, D± - коэффициент диффузии ионов, Е0 -напряженность поля на границе зоны неравновесной реакции диссоциация-рекомбинация, Ь1, Ь2 - подвижности положительных и отрицательных ионов, F - постоянная Фарадея, Т -температура, R - универсальная газовая постоянная, а - омическая проводимость, е - диэлектрическая проницаемость, е0 - электрическая постоянная.
Соотношение (1) показывает, что толщину слоя объемного заряда можно оценить как расстояние, которое ион проходит за время релаксации. Именно поэтому эффект электризации тем больше, чем ниже электропроводность жидкости и, соответственно, больше время максвелл-вагнеровской релаксации. В работе [13] показано, что вследствие того, что у электродов одноименно заряженные ионы отталкиваются, происходит нарушение равновесия реакции диссоциация-рекомбинация примесных ионов, в результате чего и образуется объемный заряд.
Рассмотрим схему ячейки с магнитной жидкостью (рис. 3). Электрическое поле транспортирует анионы (катионы) из прикатодной (приа-нодной) области. Так как катод (анод) анионы (катионы) не поставляет, то прикатодная (приа-нодная) область обедняется анионами (анионами), что приводит к образованию прикатодной (прианодной) области с положительным (отрицательным) зарядом.
Рис. 3. Схема ячейки с магнитной жидкостью: 1, 2 - плотные приэлектродные слои магнитных частиц; 3,4 - области объемного заряда; А+ и В -примесные ионы; С+ и С- - заряженные частицы магнитной жидкости; Х+, X- инжектируемые ионы
Оценка величины области локализации объемного заряда для жидкости типа магнетит в керосине проведена с использованием следующих параметров: концентрация магнетита 3.4 об.%, удельная проводимость а» 10-7 Ом-1м-1. Расчет показал, что при скачке потенциала между электродами Дф = 5 В, расстоянии между электродами d = 5 • 10-5 м, подвижности примесных ионов Ь1 » Ь2 ~ 10-8 м2В-1с-1 [13] ширина области объемного заряда равна 8 ~ 4 • 10-7 м. Значение напряженности электрического поля в приэлек-тродной области изменяется согласно уравнению Пуассона из-за наличия объемного заряда плотностью р [21].
По данным работы [22] характерное время образования объемного заряда для жидкости средней концентрации магнетита (3.4 об.%) для образцов на основе керосина составляет величину порядка 0.05 с. Оценка времени образования объемного заряда в электрическом поле требует учета миграционной поляризации частиц магнетита, которая в рассматриваемом нами случае может быть описана в рамках теории Максвелла-Вагнера. Время масвелл-вагнеровской поляризации (2) много меньше времени образования объемного заряда:
e0 e 2 • 8.85 • 10-
10-
1.7 • 10-4 c.
(2)
Объемный заряд возникает также у поверхности дисперсной частицы. Область окружающего частицу объемного заряда, занимает ту часть жидкости, где ион, появившийся в результате диссоциации примеси у поверхности частицы, продолжает свое движение, не успев прореком-бинировать.
Концентрационные автоколебания магнитной жидкости. В отсутствие электрического поля частицы магнитной жидкости не заряжены [22]. В поле Е частицы заряжаются и приобретают индуцированный дипольный момент р. В результате действия поляризационной ^ = (р—)Е и кулоновской f = цБ сил частицы начинают движение к соответствующему электроду, в результате чего образуется концентрированный слой частиц. При этом визуально наблюдается изменение цвета поверхности ячейки с МЖ. После образования плотного приэлектродного слоя при контакте магнитных частиц с поверхностью электрода происходит электрохимическое взаимодействие молекул олеиновой кислоты, окружающих частицу, с электродом. Например, на катоде в результате реакции высво-
бождаются отрицательные ионы X , имеющие в плотном приэлектродном слое почти нулевую подвижность, что приводит к его зарядке. Когда кулоновская энергия отталкивания слоя и превысит энергию притяжения и , в слое возник-
г- Г-
нет неустойчивость, частицы перезарядятся и начнут движение к противоположному электроду - возникнет концентрационная волна. Дойдя до области объемного заряда II (рис. 3), частицы начнут движение в обратную сторону. Автоколебания ансамблей магнитных частиц удалось визуализировать с помощью метода электроу-правляемой интерференции [10] по периодическому изменению цвета поверхности ячейки с магнитной жидкостью и по колебаниям оптического сигнала (рис. 2). На противоположном электроде (аноде) также образуется плотный приэлектродный слой частиц. В результате при-электродных реакций образуются ионы Х+, что вызывает движение второй концентрационной волны навстречу первой.
Допустим, что все частицы в магнитной жидкости имеют одинаковый размер, электрический заряд и равную электрофоретическую подвижность Ь. Тогда в стационарном электрическом поле Е скорость частицы равна: V = Ь(Е, п), где п - направление движения.
Подвижность частицы Ь = 5 • 10-9 м2В-1с-1 была определена по формуле, представленной в [21]. Тогда скорость частицы в объеме ячейки равна у = 5 • 10-4 м-с-1.
С другой стороны, из равенства электростатической силы и силы вязкого трения:
qE = ^ • и, (3)
где d = 10-8 м- диаметр частицы, ^ = 1.5 • 10-3 Па-с - динамическая вязкость дисперсионной среды, Е - напряженность электрического поля, заряд частицы q равен заряду электрона [22]. По полученным нами экспериментальным данным период автоколебаний ~ 1 с, то есть время, за которое частица «оторвется» от электрода и вернется обратно.
За время t = 0.5 с частица пройдет путь ~ 540-5 м, то есть расстояние до противоположного электрода. Таким образом, при скачке потенциала Дф = 5 В наблюдается периодическое движение частиц от одного электрода к другому.
Поясним, почему наблюдаемые колебания концентрации являются автоколебаниями. Автоколебания — это незатухающие колебания в диссипативной динамической системе с нелинейной обратной связью, поддерживающиеся
t
M
7
за счет энергии постоянного, то есть непериодического внешнего воздействия [23]. В реализуемых лабораторных экспериментах источником постоянного (непериодического) воздействия является скачок потенциала на электродах. Роль нелинейного регулятора выполняет емкость приэлектродного слоя С, регулирующая накопление заряда в слое. В работе [24] было показано, что вследствие образования приэлектродного слоя ячейка с МЖ представляет собой конденсатор переменной емкости. Заряд в при-электродном слое накапливается за время t1, а затем после достижения критического значения напряжения на электродах и происходит перезарядка слоя за время Таким образом, постоянное воздействие преобразуется в переменное. Колеблющимся элементом системы является концентрация магнитных частиц. Колебания концентрации осуществляются через обратную связь, в качестве которой выступает сопротивление слоя, которое при и резко уменьшается. Диссипация энергии в автоколебательной системе возмещается за счет поступления в нее энергии из источника постоянного напряжения, благодаря чему автоколебания не затухают.
2.2 Одномерная математическая модель автоколебательного процесса
Нами разработана упрощенная математическая модель автоколебаний, которая позволяет определить зависимость характеристик процесса от исходных данных задачи, например, от скачка потенциала между электродами Лф, геометрических размеров и т. д. В предложенной модели учитывается перенос только положительно и отрицательно заряженных магнитных частиц с концентрациями С1 , С2 и потоками j1, ;'2, скачком потенциала Лф без учета химических реакций. Влияние примесных и инжектируемых ионов учитывается косвенно в граничных условиях. Основная идея построения модели заключается в исследовании вопроса о возможности возникновения колебательного движения из-за перезарядки частиц МЖ возле электродов или в областях локализации объемного заряда.
Обозначим х = 0 и х = Н - условные точки, в которых происходит процесс зарядки частиц МЖ, в частном случае это могут быть анод и катод. Начальное условие - СД0, х) = С10( х) и С2(0, х) = С20( х). Первоначальное распределение нейтральных частиц МЖ с известной концентрацией определяется диффузией [25]. В электри-
ческом поле в начальный момент времени t = 0 происходит зарядка частиц МЖ. В окрестности х = Н частицы МЖ заряжаются отрицательно, а в окрестности х = 0 заряжаются положительно. В электрическом поле отрицательно заряженные магнитные частицы начинают двигаться в сторону анода, а положительно заряженные в сторону катода. В результате у электродов образуются плотные слои.
Использована одномерная нелинейная математическая модель с классической системой дифференциальных уравнений материального баланса, Нернста-Планка и Пуассона [26-29] без упрощений и подгоночных параметров:
^ -J± i - 12
dt ~ dx'' =
h --
d j
—zp.qj-DiC, i =1,2,
RT dx ' dx
F
£
dx 2 - e (Z1C1 + Z2Cl) ,
(4)
(5)
(6)
где С, ji - концентрации и потоки заряженных частиц МЖ, zi - зарядовые числа, ф - потенциал, ег - диэлектрическая проницаемость МЖ, F -постоянная Фарадея, D¡ - коэффициенты диффузии заряженных частиц МЖ, Т0 - температура, R - универсальная газовая постоянная.
Поскольку в представленной Ш модели автоколебаний частиц в магнитной жидкости рассматривается два типа заряженных магнитных частиц, то процесс перезарядки и соответствующие колебания происходят от одной до другой области объемного заряда.
Предположим, что процесс перезарядки волны заряженных частиц МЖ в точках х = 0 и х = Н происходит постепенно по мере подхода заряженных частиц к ним. Например, отрицательно заряженные частицы в точке х = 0 получают положительный заряд и становятся положительно заряженными, а положительно заряженные частицы МЖ (если они есть) доходят до точки х = Н и заряжаются отрицательно (рис. 4б). Отметим, что процесс перезарядки описывается краевыми условиями при х = 0 (7) и при х = Н (8), скачок потенциала Лф задается в виде (9).
RT
DC j-Di 2 dx
dC2 dx
- hlA ,
-^DC j- D2 d-C-
RT dx 2 dx
(7)
- h
x-0
x-0
В. С. Чеканов и др. Экспериментальное исследование и математическое моделирование...
U 0.2 0.4 0.6 0.8 1 о 0.3 0.4 0.6 0.8 I
В Г
Рис. 4. Распределение концентраций положительно (С1) и отрицательно (С2) заряженных частиц МЖ в моменты времени: а) t = 0.01 с, б) t = 0.02 с, в) t = 0.03 с, г) t = 0.04 с
Функция (11) соответствует условиям (10) и закону Ома в момент времени t = 0.
3. Обсуждение результатов численного исследования
Численные исследования проведены при следующих параметрах задачи: расстояние между электродами Н = 5-10-5 м, расчет проведен для t е[0,10] секунд с шагом запоминания 0.01 с. Параметры сетки: п = 20 000 элементов. Скачок потенциала Дф исследовался в диапазоне от 1 до 20 В. В расчетах принято некоторое начальное распределение концентрации магнитных частиц, соответствующее экспериментальным данным:
С1(0,х) = С2(0,х)=3.4 об. %
Для численного решения краевой задачи математической модели использовался метод конечных элементов в среде моделирования в Comsol Multiphysics.
Объемная концентрация частиц магнетита, выраженная в моль/м3, равна:
С = 0.73593 • 10-3 моль/м3,
С (0, х) = С10 (х) = 0.0074 • е-х(0/01Н) моль/м3,
С2(0, х) = С20( х) = 0.0074 • е-(н - х/(0/01Н).
v RTo
^ „ ^ dC., -ДС^- - Д 2
dj
IX
\RTo
1 dx
-ДС^- - Д^ dx
dj
'2VJ1 dX
= JlK ,
j2K,
(8)
j(t,0) = Dj, j(i, H) = 0 ,
(9)
где ]1А, ]2А, ]1К, ]2К - потоки перезаряженных частиц у анода и катода.
Начальные условия (t = 0) определяют исходное состояние системы, а ее дальнейшая эволюция определяется уравнениями и краевыми условиями. Так как тонкие приповерхностные слои образуются у обоих электродов, будем считать, что в начальный момент времени частицы МЖ сосредоточены у каждого из электродов (10):
С1(0, х) = С10( х), С2(0, х) = С20( х) (10)
Как будет показано ниже, такое начальное состояние приводит к появлению двух движущихся навстречу друг к другу концентрационных волн. При этом волны взаимодействуют друг с другом. Для изменения потенциала начальное условие определим линейной функцией (11):
H - Х
j(0, x) = j0( x) = —— Dj.
(11)
Согласно этим начальным условиям частицы МЖ первоначально распределены достаточно плотно у электродов, то есть приповерхностный слой уже образовался.
3.1. Результаты численного анализа для двух концентрационных волн
На рис. 4а представлено начальное распределение концентраций частиц СД£,х) и С2(Ь,х). В начальный момент времени они заряжаются при координате х = 0 - положительно, при х = Н -отрицательно. Положительно заряженные частицы МЖ начинают двигаться по направлению к х = Н, а отрицательно заряженные магнитные частицы к х = 0 (рис. 4б). Пунктирной линией (С2) показаны отрицательно заряженные магнитные частицы, которые идут к положитель-
но заряженному электроду (аноду), а сплошной линией (С1) -положительно заряженные магнитные частицы, которые идут к отрицательно заряженному электроду (катоду). Первоначально две разные по амплитуде концентрационные волны двигаются навстречу друг другу. Затем их амплитуды постепенно становятся равными (рис. 4г). После того, как концентрационные волны достигают соответствующих электродов, частицы перезаряжаются. Далее цикл повторяется.
В момент времени t = 0.05 с (рис. 5а) образовался фронт положительно заряженных магнитных частиц. Один цикл полностью завершился (рис. 5б), начался следующий цикл (рис. 5в) и т. д. В моменты времени t = 0.14 с, t = 0.17. (рис. 5д, е) количество положительно и отрицательно заряженных частиц становится почти постоянным.
Рис. 5. Распределение концентраций положительно (С1) и отрицательно (С2) заряженных частиц МЖ в момент времени: а) t = 0.05с., б) t = 0.06 с, в) t = 0.07 с, г) t = 0.1 с, д) t = 0.14 с, е) t = 0.17 с
В. С. Чеканов и др.
Экспериментальное исследование и математическое моделирование...
В момент времени 0.29 с (рис. 6а) видно, как волна положительно заряженных магнитных частиц пошла к катоду, а волна отрицательно заряженных частиц движется навстречу, к аноду. При t = 0.33 с (рис. 6б) у анода волна максимально выросла и дальше начнет падать (рис. 6б).
Таким образом, частицы совершают движение от одного электрода к другому и обратно. Перезарядка происходит на электродах.
3.2. Сравнительный анализ численного и лабораторного экспериментов
Анализ выявил, что период колебаний в численном эксперименте варьируется от 0.5 до 0.6 с, что согласуется с экспериментальными данными. Численное исследование проведено при различных значениях скачка потенциала Лф в диапазоне от 1 до 20 В. Получено, что чем больше скачок потенциала, тем меньше период колебаний, что также подтверждается экспериментально. Кроме того, необходимо отметить, что чем меньше скачок потенциала, тем более процесс диффузии преобладает над электромиграцией, что подавляет автоколебательный процесс.
Показано, что при скачке потенциала меньше 3 В концентрационные волны не образуются при компьютерном моделировании и не наблюдаются экспериментально. Были рассмотрены различные значения зарядовых чисел частиц от ±1 до ±4 и выявлено, что чем больше величина заряда, тем более хаотичным становится движение и форма волны.
3.3. Программная реализация
Создан программный продукт «AutoWave01» с интуитивно понятным пользовательским интерфейсом и расширенным функционалом, предназначенный для изучения автоволновых
процессов в тонком слое магнитного коллоида. Использование приложения «AutoWave01» существенно сокращает затраты на дорогостоящие эксперименты и позволяет провести предварительный анализ с целью нахождения оптимальных параметров, давая возможность сохранения всех результатов исследования в удобном формате.
4. Заключение
В работе проведено математическое и компьютерное Ш моделирование автоколебательного процесса частиц в МЖ, которое на качественном уровне полностью отвечает лабораторному эксперименту. При повышении напряжения на электродах до определенного порогового значения и возникают колебания частиц дисперсной фазы между электродами. Это явление связано с зарядкой и перезарядкой частиц. Предсказательной ценностью данной модели являются возможность определения критического значения скачка потенциала, нахождения времени установления периода колебаний и другие характеристик процесса.
Заявленный вклад авторов
Все авторы сделали эквивалентный вклад в подготовку публикации.
Конфликт интересов
Авторы заявляют, что у них нет известных финансовых конфликтов интересов или личных отношений, которые могли бы повлиять на работу, представленную в этой статье.
Список литературы
1. Мурадова А. Г., Шарапаев А. И., Зайцева М. П., Кузнецова С. А., Юртов Е. В. Наноструктурирован-
Рис. 6. Распределение концентраций положительно (С1) и отрицательно (С2) заряженных частиц МЖ в момент времени а) t = 0.29 с, б) t = 0.33 с (в увеличенном масштабе)
В. С. Чеканов и др. Экспериментальное исследование и математическое моделирование...
ные оксиды железа. Успехи в химии и химической технологии. 2019;1(211): 77-78. Режим доступа: https://www.elibrary.ru/item.asp?id=38098384
2. Zakinyan A. Electrical conductivity of inverse magnetic fluid emulsion. Magnetohydrodynamics. 2018;54(1-2):163-166. https://doi.org/10.22364/ mhd.54.1-2.29
3. Чеканов В. В., Кандаурова Н. В., Рахмани-на Ю. А., Чеканов В. С. Индикатор ультразвука 2. Патент РФ № 2446384; опубл. 27.03.2012; бюл. № 9.
4. Чеканов В. В., Кандаурова Н. В., Рахмани-на Ю. А., Чеканов В. С. Индикатор теплового излучения. Патент РФ № 2446422; опубл. 27.03.2012; бюл. № 9.
5. Чеканов В. В., Кандаурова Н. В., Рахмани-на Ю. А., Чеканов В. С. Индикатор разности потенциалов. Патент РФ № 2449382; опубл. 27.04.2012; бюл. № 12.
6. Chekanov V. V., Kandaurova N. V. Chekanov V. S. Observation of the autowave process in the near-electrode layer of the magnetic fluid. Spiral waves formation mechanism. Journal of Molecular Liquids. 2018;272: 828-833. https://doi.org/10.1016/j.molliq.2018.10.073
7. Шлиомис М. И. Магнитные жидкости. Успехи физических наук, 1974;112(3): 427. https://doi. org/10.3367/UFNr.0112.197403b.0427
8. Rosenzweig R. E. Ferrohydrodynamics. Cambridge University Press; 1985. 344 p.
9. Борн М., Вольф Э. Основы оптики. М.: Наука; 1973. 720 с.
10. Чеканов В. В., Кандаурова Н. В., Чеканов В. С., Романцев В. В. Использование электроу-правляемой интерференции в электроперестраи-ваемом гаетофильтре и для наблюдения автоволнового процесса в приэлектродном слое магнитной жидкости. Оптический журнал. 2019;86(1): 21-26. https://doi.org/10.17586/1023-5086-2019-86-01-21-26
11. Жакин А. И. Ионная электропроводность и комплексообразование в жидких диэлектриках. Успехи физических наук. 2003;173(1): 51-68. https:// doi.org/10.3367/ufnr.0173.200301c.0051
12. Жакин А. И. Кинетика агрегирования в неполярных жидких диэлектриках. Электронная обработка материалов. 2015;51(4): 49-60. Режим доступа: https://eom.ifa.md/ru/journal/ shortview/1206
13. Жакин А. И. Приэлектродные и переходные процессы в жидких диэлектриках. Успехи физических наук. 2006;176(3): 279-310. https://doi. org/10.3367/UFNr.0176.200603d.0289
14. Жакин А. И. Электрогидродинамика. Успехи физических наук. 2012;182(5): 495-520. https://doi. org/10.3367/UFNr.0182.201205b.0495
15. Стишков Ю. К., Богданов Д. В. Влияние неравновесных приэлектродных слоев на структуру
ЭГД течений в треионной модели диэлектрической жидкости. Журнал технической физики. 2017;87(8): 1144. https://doi.org/10.21883/JTF.2017.08.44719.1882
16. Кожевников В. М., Ларионов Ю. А., Демин М. С. Перенос и накопление заряда в слое магнитодиэлектрического коллоида с наноразмер-ными частицами Вестник Северо-Кавказского государственного технического университета. 2008;3: 46 50. Режим доступа: https://www.elibrary.ru/item. asp?id=11674845
17. Kolesnikova A., Zakinyan A. Rotating magnetic field induced structure formation in a magnetic fluid emulsion. Magnetohydrodynamics. 2018;54(1-2): 4548. https://doi.org/ 10.22364/mhd.54.1-2.8
18. Kuz'ko A. E., Chekanov V. S. Anisotropy of magnetic fluid conductivity in constant magnetic fields. Surface Engineering and Applied Electrochemistry. 2020;56(6): 727-733. https://doi.org/10.3103/ S1068375520060095
19. Kozhevnikov V. M., Larionov Y. A., Chuenko-va I. Y. Features of electrical properties in a structured thin magnetic fluid layer. Magnetohydrodynamics. 2018;54(1-2): 55-59. https://doi.org/10.22364/ mhd.54.1-2.10
20. Духин С. С., Эстрела Льопис В. Р., Жолковский Э. К. Электроповерхностные явления и электрофильтрование. Киев: Наукова думка; 1985. 288 с.
21. Erin K. V. Study of the kinetics of space charge formation in colloidal magnetic nanoparticles in liquid dielectrics. Surface Engineering and Applied Electrochemistry. 2017;53(4): 327-332. https://doi. org/10.3103/S 1068375517040044
22. Ерин К. В. Электрооптические исследования электрофореза коллоидных частиц магнетита в керосине в поле приэлектродного объемного заряда. Коллоидный журнал. 2015;77(1): 24-29. https:// doi.org/10.7868/S0023291215010073
23. Вавилин В. А. Автоколебания в жидкофаз-ных химических системах. Природа. 2000;5Ж 19-24. Режим доступа: https://www.elibrary.ru/item. asp?id=35128874
24. Chekanov V. V. Kandaurova N. V. Chekanov V. S. Effect of a nearsurface nanolayer formation on the magnetic fluid electrical properties. Acta Tech-nica. 2018;63(4): 555-562.
25. Чеканов В. С., Коваленко А. В., Дискае-ва Е. Н., Кириллова Е. В. Математическое моделирование автоволнового процесса в тонком слое магнитного коллоида. Экологический вестник научных центров Черноморского экономического сотрудничества. 2020;17(4): 57-67. https://doi. org/10.31429/vestnik-17-4-57-67
26. Urtenov M. K., Uzdenova A. M., Kovalenko A. V. Nikonenko V. V., Pismenskaya N. D., Vasil'eva V. I., Pourcelly G. Basic mathematical model of overlimiting
В. С. Чеканов и др. Экспериментальное исследование и математическое моделирование...
transfer enhanced by electroconvection in flow-through electrodialysis membrane cells. Journal of Membrane Science. 2013; 447: 190-202. https://doi. org/10.1016/j.memsci.2013.07.033
27. Nikonenko V. V., Vasil'eva V. I., Akberova E. M., Uzdenova A. M., Urtenov M. K., Kovalenko A. V., Pis-menskaya N. D., Mareev S. A., Pourcelly G. Competition between diffusion and electroconvection at an ion-selective surface in intensive current regimes. Advances in Colloid and Interface Science. 2016;235: 233-246. https://doi.org/10.10Wj.cis.2016.06.014
28. Gudza V. A., Pismenskiy A. V., Urtenov M. K., Shkorkina I. V., Chubyr N. O., Коваленко А. V. The influence of water dissociation/recombination on transport of binary salt in diffusion layer near ion exchange membrane. Journal of Advanced Research in Dynamical and Control Systems. 2020;12(4): 923-935. https://doi.org/10.5373/JARDCS/V12SP4/20201563
29. Kovalenko A. V., Yzdenova A. M., Sukhinov A. I., Chubyr N. O., Urtenov M. Kh. Simulation of galvanic dynamic mode in membrane hydrocleaning systems taking into account space charge. AIP Conference Proceedings. 2019;2188(1): 050021. https://doi. org/10.1063/1.5138448
Информация об авторах
Чеканов Владимир Сергеевич, к. т. н., доцент кафедры информационных систем и технологий, Северо-Кавказский Федеральный Университет, Ставрополь, Российская Федерация; e-mail: [email protected]. ORCID iD: https://orcid.org/0000-0002-2680-2883.
Кириллова Евгения Вадимовна, к. ф.-м. н., профессор Университета Прикладных Наук Рейн-Майн города Висбаден, Германия; e-mail: kirillova@web. de. ORCID iD: https://orcid.org/0000-0002-6797-0920.
Коваленко Анна Владимировна, д. ф.-м. н., доцент, заведующая кафедрой анализа данных и искусственного интеллекта, Кубанский государственный университет, Краснодар, Российская Федерация; e-mail: [email protected]. ORCID iD: https://orcid.org/0000-0002-3991-3953.
Дискаева Елена Николаевна, к. ф.-м. н., доцент кафедры промышленных технологий, Филиал федерального государственного бюджетного образовательного учреждения высшего образования «МИРЭА — Российский технологический университет» в г. Ставрополе, Ставрополь, Российская Федерация; e-mail: [email protected]. ORCID iD: https://orcid.org/0000-0002-5185-6023.
Поступила в редакцию 01.04.2021; одобрена после рецензирования 02.09.2021; принята к публикации 15.11.2021; опубликована онлайн 25.12.2021.
