Нелинейная динамика. 2017. Т. 13. № 2. С. 195-206. Полнотекстовая версия в свободном доступе http://nd.ics.org.ru Б01: 10.20537/па1702004
ОРИГИНАЛЬНЫЕ СТАТЬИ
УДК: 544.431.8+544.015.4
М8С 2010: 92Е99, 80А20, 74А50, 82С26, 70К50
Осциллирующие фазовые переходы раствор —газ и раствор — кристалл в каплях растворов с одним кристаллизующимся компонентом
В. Б. Федосеев
Экспериментально воспроизведены длительные апериодические несинхронные осциллирующие фазовые переходы газ-раствор и раствор - кристалл в ансамбле капель водного раствора ^С1. Наблюдения позволяют упростить ранее предложенную кинетическую модель осциллирующего режима фазовых превращений, исключив необходимость учета образования метастабильных фаз. Модель связывает скорость и направление потока растворителя в газовой фазе с давлением насыщенного пара у поверхности капли, температурой раствора, скоростями изменения температуры, объемов кристалла и капли. Анализ модели показывает, что присутствие в системе кристаллической фазы приводит к появлению по крайней мере двух особенностей (бифуркаций) зависимости химического потенциала летучего компонента от его количества в капле.
Ключевые слова: фазовые превращения, осцилляции, испарение, конденсация, кристаллизация, капли раствора
1. Введение
В работе [1] впервые описан эффект осциллирующих фазовых превращений в трехком-понентной трехфазной системе Н2О —^С1 —КС1, заключающийся в возникновении долговременных апериодических совместных фазовых превращений кристалл -раствор и рас-
Получено 07 января 2017 года После доработки 15 мая 2017 года
Работа поддержана грантом Российского научного фонда (проект №15-13-00137).
Федосеев Виктор Борисович [email protected]
Институт металлоорганической химии им. Г. А. Разуваева РАН 603137, Россия, г. Нижний Новгород, ул. Тропинина, д. 49
твор-пар. Эффект наблюдается в виде многократного самопроизвольного испарения летучего растворителя практически до полного исчезновения жидкой фазы вокруг кристалла с последующей конденсацией растворителя, растворением кристалла и восстановлением размера капель. В реальном времени эти осцилляции демонстрирует видеофайл, размещенный в дополнительных материалах к работе [1].
Мультистабильные состояния в столь простой по химическому составу системе наблюдаются впервые. Колебательный режим в гетерогенной трехфазной системе при стационарных внешних условиях представляет самостоятельный интерес. Возможность управлять процессом путем изменения внешних параметров (температура, состав газовой фазы) позволяет целенаправленно переводить систему между разными равновесными состояниями. Развитие кинетической модели этого процесса при простоте экспериментального воспроизведения результатов моделирования может быть актуально при создании материалов с памятью формы, создании искусственных нейронных сетей, моделировании устройств памяти на фазовых переходах.
Смеси, использованные в экспериментах, по составу близки к физиологическим растворам, поэтому результаты моделирования могут представлять интерес и при обсуждении взаимодействия коацерватов на стадии возникновения жизни.
Наиболее известным и изученным колебательным физико-химическим процессом является реакция Белоусова-Жаботинского [2]. Эта и подобные химические реакции наблюдаются в гомогенных системах. Основными периодически меняющимися во времени величинами в этих процессах являются концентрации компонентов реакционной смеси без изменения объема, геометрической формы и температуры системы.
Примером колебательных процессов с участием фазового превращения могут быть кольца Лизеганга [3], более века привлекающие внимание математиков [4]. В этом случае колебания концентрации компонентов сопровождают процесс перемещения межфазной границы и заканчиваются при затвердевании всего объема при охлаждении системы.
Осциллирующие фазовые превращения, при которых система многократно меняет объем и состав сосуществующих фаз, ранее не описаны. Теоретическое описание этого эффекта пока ограничено качественными предположениями в [1].
Одно из этих предположений основано на гипотезе о возникновении фазы твердого раствора ^С1-КС1, которая неустойчива в макроскопической системе при обычных условиях. Эта гипотеза опирается на термодинамические расчеты фазовых равновесий в системах малого объема [5, 6]. На основе наших оценок при малом объеме в смеси ^С1-КС1 —НО можно ожидать возникновения необычных в рассматриваемых условиях метастабильных фаз, в частности твердого раствора ^С1-КС1 [7]. Образованию этой фазы способствует большое пересыщение при условии рождения единичного центра кристаллизации. Возможность достижения необходимого пересыщения показана экспериментально в [7]. В [8] описаны некоторые термодинамические аспекты, связанные с энергией образования и размерами критического зародыша новой фазы при уменьшении объема системы.
2. Экспериментальная часть
Модель осциллирующего режима, основанная на гипотезе о роли метастабильного раствора, содержит большое число независимых переменных (размер и температура капли, размеры кристаллов каждого из нелетучих компонентов либо состав и размер метастабиль-ного кристалла) и термодинамических параметров (температура и давление пара летуче-
го компонента во внешней среде, мольные объемы и растворимости, энтальпии испарения и растворения компонентов).
Наблюдение осцилляций в двухкомпонентных системах могло бы существенно упростить и интерпретацию, и кинетическую модель явления. Для проверки этого предположения выполнены эксперименты с водными растворами, содержащими только один растворенный компонент.
Эксперименты с водными растворами ряда кристаллических соединений (NaCl, KCl, NaNÜ3, (NH2)2CO) показали, что осцилляции воспроизводятся и в более простых системах. Это позволило утверждать, что образование метастабильных фаз не является необходимым фактором возникновения осцилляций, и упростить модель процесса.
Экспериментальные исследования проводились на основе методики, описанной в [1]. Растворы наносились на стеклянную подложку с помощью распылителя. После испарения растворителя на подложке оставался ансамбль кристаллов разного размера. Средние размеры кристаллов, образующихся при высыхании капель, можно регулировать концентрацией распыляемого раствора. Подложка с ансамблем кристаллов помещалась на предметный столик микроскопа, затем на подложке создавался резервуар с чистым растворителем (дистиллированная вода). Описанные осцилляции можно наблюдать на отдельных кристаллах, находящихся в окрестности резервуара. В настоящей работе использована модифицированная методика, в которой резервуар с растворителем и ансамбль капель разделены (рис. 1). Ансамбль полученных распылением кристаллов находился на нижней поверхности предметного стекла на высоте от 2 до 5 мм от резервуара с растворителем. В такой схеме резервуар с растворителем может быть внесен, удален, перемещен или заменен на любой стадии опыта. Это делает возможным повторение опытов с одним и тем же ансамблем кристаллов или с растворителями разного состава. Изменение состава растворителя в резервуаре позволяет регулировать влажность среды, в которой находится наблюдаемый ансамбль капель.
1
Рис. 1. Схема установки для наблюдения осциллирующих фазовых превращений. 1 — объектив микроскопа, 2 — предметное стекло с нанесенными методом распыления каплями или кристаллами, 3 — боковые опоры, определяющие высоту расположения наблюдаемых объектов, 4 — растворитель, 5 — подвижная прозрачная подложка-резервуар для растворителя, 6 — прозрачное основание, перемещаемое объектоводителем микроскопа, 7 — нижняя подсветка.
Для устранения теплового влияния осветителя использовалась подсветка на основе светодиодного осветителя и зеркала. Для устранения влияния воздушных потоков в лабораторном помещении использовалась ширма. Следует отметить, что влияние внешних потоков нельзя считать критичным, так как в окрестности осциллирующих капель можно наблюдать капли и кристаллы разного размера и удаленности от резервуара, сохраняющие неизменным свое состояние на протяжении всего эксперимента.
Наблюдения состояния ансамбля капель и кристаллов сохранялись в виде видеофайлов либо серии кадров, сделанных с интервалами 2-5 секунд протяженностью до 20 минут.
При наблюдениях были зафиксированы осцилляции как отдельных, так и групп капель. На рисунке 2 показана выборка кадров, демонстрирующая состояние группы капель раствора ^С1 в отдельные моменты времени. В исходном состоянии до появления резервуара с водой на подложке присутствовали только кристаллы. После внесения резервуара с растворителем группа кристаллов, имеющих линейные размеры от 13 до 25 мкм, поглощает воду с образованием капель раствора. Для рисунка 2 отобраны моменты, в которые капли почти полностью теряют воду или восстанавливают размеры, достигая в диаметре 30-50 мкм. На полной кинограмме, приведенной на странице электронной версии статьи (видеофайл 1 — http://nd.ics.org.ru/nd1702004/#vidl), можно наблюдать и несинхронные изменения размеров отдельных капель с меньшей амплитудой изменения размера.
fF\ p—U'' IlMJ p 0 NaCl v ■ t 09 & »ц 4 _ 196 s NaCl P „ t * 9 9
Я 266« • ° NaCl Р * О© • _ 300» ° NaCl P v i, t * 9 * n 350S •J NaCl P " * 00 9 Jp 50ц
_ 370s 'О NaCl Р l * ф * j - л л 4005 О NaCl P t Ъ® * Я? »' J а «4» ' ° NaCl p * * Ф f J 50, J
Рис. 2. Состояние группы капель раствора NaCl на 62, 196, 266, 300, 350, 370, 400, 434 секунде от начала записи изображения. Верхний кадр слева демонстрирует исходное состояние (t = 0).
Осцилляции разной интенсивности и длительности по описанной методике удалось воспроизвести и с растворами других водорастворимых веществ (в частности, KCl, NaNÜ3, (NH2)2CO). В некоторых случаях затухающие апериодические колебания можно наблюдать в режиме перехода из состояния кристалла в состояние капли со стабильным размером. Рисунок 3 демонстрирует осцилляции, сопровождающиеся изменением формы кристалла. (Соответствующий рисунку 3 видеофайл приведен на странице электронной версии статьи (видеофайл 2 — http://nd.ics.org.ru/nd1702004/#vid2).) Это свидетельствует о том, что
Рис. 3. Изменение огранки кристалла УаС1 в процессе осциллирующих фазовых превращений. Состояния системы на 2, 67, 110, 116, 165, 179, 183, 187, 191 секунде от начала записи изображения.
процессы конденсации и испарения сопровождаются растворением и перекристаллизациеи частиц, на которых образуется капля. Кроме того, рисунок 3 демонстрирует, что осцилляции способствуют релаксации структуры кристалла к более устойчивой конфигурации.
3. Модель испарения и роста капель раствора, содержащих кристалл
Уменьшение числа химических компонентов и отсутствие метастабильных фаз существенно упрощает математическую модель.
Модель основана на рассмотрении потока растворителя в газовой фазе, вызванного разностью химических потенциалов летучего компонента в термостате (резервуар с растворителем или окружающая среда) и в окрестности капли
з = --> (ЗЛ)
где Т, Р — температура и давление летучего компонента в термостате, Т^, Р4 — температура и давление летучего компонента около поверхности капли, А1 — расстояние от резервуара с растворителем до поверхности осциллирующих капель, а — коэффициент [9]. В приближении идеального газа химический потенциал имеет вид у = Ц°(Т) + КТ 1п Р, здесь Я —
газовая постоянная, |°(Т) — стандартный химический потенциал растворителя в газовой фазе. Поток j, согласно (3.1), направлен в сторону понижения химического потенциала. В изотермических условиях уравнение (3.1) сводится к закону Фика, согласно которому поток направлен в сторону понижения давления пара.
Давление Pd летучего компонента в окрестности границы раствор-пар определено условием равенства химического потенциала летучего компонента в растворе и паре:
V0,8 = lo,s + RTd ln(xo) = Ho,gas = l°gas + RTd ln Pd.
Для систем с одним растворенным компонентом давление пара около поверхности капли зависит от радиуса r капли, концентрации (мольной доли) x растворенного компонента, размера L кристалла. Величины r, x, L связаны условиями сохранения (3.2)-(3.3), которые позволяют понизить число независимых переменных до трех, при этом Pd = P(Td,r,x).
Условия сохранения количества нелетучего компонента связывают количества молей растворителя no и растворенного компонента щ в растворе nis ив кристалле U\c с концентрацией раствора
х = —^—, п\ = nis + п\с. (3.2)
п0 + пи' v ;
Величины ni, Uis и no или x полностью определяют объемы Vd, Vc капли и кристалла соответственно. Геометрические параметры, в частности радиус кривизны поверхности капли Rd и размер кристалла L (длина ребра куба), определены, если задать форму капли (сферическая) и кристалла (куб):
1
Vd = noVo + niVi, VC = nicVi, L = Vc3, (3.3)
где Vo, Vi — мольные объемы растворителя и растворенного компонента соответственно. В приближении идеального раствора они не зависят от концентрации. Предполагается, что кристалл имеет кубическую форму и включен в объем капли.
Растворимость NaCl в рассматриваемых условиях практически не зависит от температуры [10], поэтому в качестве дополнительного приближения примем, что концентрация раствора в объеме капли в присутствии кристалла сопоставима с концентрацией xs насыщенного раствора и постоянна x = xs = const при L = 0. Это приближение пренебрегает массопереносом внутри капли и тем, что концентрация насыщенного раствора при малых объемах зависит от радиуса капли [8].
Кристалл полностью смачивается раствором, при этом содержащая кристалл капля
при Vd > 7ГL3 имеет сферическую форму. При условии Vd < тгЬ3 растущий кристалл
уже не вписывается в сферу и капля теряет сферическую форму. В этом случае форму капли можно приближенно представить в виде куба, грани которого являются основаниями сферических сегментов, заполненных раствором. Пренебрегая перекрывающимися и выступающими краевыми частями этих сегментов, объем покрывающего кристалл раствора можно оценить как
Vd-Vc = п0\о + пи\1 ~ бтг/?2 (Ra - = irh (h2 + |L2 ),
где h — высота сферического сегмента.
При уменьшении объема раствора (испарении) радиус кривизны поверхности опира-
Ь2 + 2Н2
ющегося на грань куба шарового сегмента Ка
4Н
возрастает (Ь растет, Н умень-
шается). Это несколько усложняет условия сохранения количества нелетучего компонента, связывающие объем капли Уа и размеры кристалла,
Уа = ±жВ% = п1У1(г1 + 1)-Ь3г1
л л/3 при г = Я(1> — Ь,
где п =
Уа = щУг (г? + 1) - Ь3г] и тг/г + | Ь2^ + Ь3 при /г < 1.
Уо 1 - х
V7! -г' '
(3.4)
В первом случае капля имеет сферическую форму, радиус капли совпадает с радиусом кривизны Ка. Во втором случае радиус кривизны Ка превышает радиус сферической капли объема Уа. Рисунок 4 схематически демонстрирует эти отличия.
Рис. 4. Несферическая капля, смачивающая кристалл, и сфера, описанная вокруг куба.
Предположение, что кристаллическая фаза представлена единственным включением кубической формы, допустимо для систем малых размеров. Для ^С1 (см. рис. 2 и 3) и некоторых других солей образование единичного кристалла характерно вплоть до на-нолитровых капель.
Давление насыщенного пара растворителя над сферической поверхностью идеального раствора можно описать, объединив уравнение Кельвина и закон Рауля
Р (Та, Ка, х) = Р (Та) (1 - х) ехр [т^щг ^о ),
(3.5)
где а — поверхностное натяжение поверхности капли, К — универсальная газовая постоянная.
Давление насыщенного пара над сферической поверхностью раствора, содержащего кристалл размера Ь, можно описать, объединяя (3.5) с уравнением Фрейндлиха, описывающим зависимость концентрации насыщенного раствора от размера кристаллов; получим
4а с.
Уо
2а
-Уо
Р (Та, Ка, Ь, х) = Р (Та) I 1 - х,еьвт* I ея*кт* ,
(3.6)
где х8 — концентрация насыщенного раствора при Ка ^ со, Ь ^ со, аС8 — поверхностная энергия межфазной границы раствор - кристалл [11].
Рассмотрим более детально поток растворителя в газовой фазе (3.1). Разность стандартных химических потенциалов растворителя ц° в небольшом интервале температур ограничим первыми членами полиномиального представления [12], при этом ц° (Та) — — (Т) = ЪАТ, где Ь — константа.
Используя (3.5) и уравнение Клапейрона — Клаузиуса 1п ^^ = > гДе ДьН —
Р (Т) И±а Т
энтальпия испарения растворителя, можно формально свести потоки массо- и теплоперено-са в газовой фазе к функциям градиентов концентрации и температуры [9]. С учетом (3.5) массоперенос в газовой фазе представим как
3=а[атУТ + ^1 ), (3.7)
где ат = Ъ + + Д 1п {Р (Г)), аР = + КГа 1п (1 - х) и УТ и
Т Я^ А1
Первое слагаемое определено состоянием термостата (резервуара) и свойствами растворителя. Ему соответствует направление потока в сторону понижения температуры. Второе слагаемое определено состоянием капли, оно определяет поток в сторону большего радиуса кривизны (для термостата радиус кривизны с) и высокой концентрации растворенного компонента (для термостата с чистым растворителем х = 0). При испарении сферической капли величина ар растет. Для несферических капель при испарении Ка — с и ар — ЯТа 1п(1 — х8) < 0. Коэффициенты, использованные в (3.1) и (3.7) для идеальных растворов, имеют простую связь с коэффициентами переноса (диффузии, термодиффузии, теплопроводности) [9].
Тепловые эффекты в испаряющейся капле раствора определяются энтальпией испарения растворителя АьН°, энтальпией растворения кристалла А3И° и теплопереносом ((х через поверхность капли
Яа = хг АьН° = зАлЬьН°, У0
V (3.8)
у1
(х » ХАаУТ,
где скорость изменения объема капли Уа = ]АаУо связана с массопереносом через ее по
3
верхность, Ad = 4nR2 — площадь поверхности капли, Л — коэффициент теплопроводности
в газовой фазе.
В приближении x = const скорости изменения объема кристалла Vc и капли Vd связаны уравнениями сохранения (3.2)-(3.3).
При конденсации растворителя происходит растворение кристалла, при испарении капли кристалл растет. Тепловые эффекты фазовых превращений имеют разный знак и отчасти компенсируют друг друга. Суммарный тепловой эффект Q = Qd + Qc + Qх определяет скорость изменения температуры капли
Td = Q Cp, (3.9)
где cp — теплоемкость капли с кристаллом. При грубой оценке cp = no (Rd, x) Cp,o (T) + + niCp,i (T), где Cp — молярные теплоемкости чистых компонентов.
Уравнения (3.2)—(3.9) связывают j и Т^ с изменениями объема кристалла Vc и капли Vd. Они составляют основу кинетической модели, способной воспроизвести осциллирующий режим совместных фазовых превращений кристалл - раствор и раствор -пар.
Ранее подобные кинетические модели были предложены для описания процесса испарения капель смесей летучих компонентов [13, 14]. В них методами численного моделирования обнаружена возможность немонотонного изменения радиуса капли во времени. При некоторых условиях такой эффект возникает при установлении стационарной концентрации раствора, после чего капля переходит в режим монотонного испарения [13]. В численных расчетах было получено немонотонное поведение температуры [14] и размера [15] капли испаряющейся многокомпонентной смеси летучих компонентов. Авторы этих работ связывают немонотонное поведение размера и температуры капли с переходными процессами. Подобный эффект наблюдался экспериментально [16].
Возможно, модели [13-16] с гомогенной конденсированной фазой могут воспроизвести и затухающие осцилляции, характерные для установления стационарного режима. В трехфазной системе осцилляции при релаксации системы к устойчивому состоянию можно наблюдать экспериментально. На рисунке 3 можно видеть, как в процессе осцилляций размера капли меняется огранка кристалла NaCl. При этом осцилляции продолжались и после перехода к более выгодной огранке кристалла.
Согласно [17], нет качественной связи между процессами массопереноса в газовой фазе и внутри капли. Это позволяет использовать более простую модель сферических капель, не рассматривающую влияние подложки и конвективных потоков в каплях. Однако и конвективные потоки могут быть причиной осцилляций капель. Экспериментально осциллирующий рост и движение капель в системе жидкость - жидкость наблюдались в процессе экстракции ацетона из водной среды каплями толуола [18]. В этом случае возникновение осцилляций, по мнению авторов, связано с конвекцией Марангони.
Простейшая модель, описывающая обратимый массоперенос между каплей раствора и газовой фазой (3.7), включает в себя условия локального равновесия на границе раствор-газ (3.5) и (3.6), условия баланса вещества (3.2)-(3.4), теплового баланса (3.8)-(3.9).
4. Особенности испарения и роста капель раствора, содержащих кристалл
Для трехфазной системы уравнения, описывающие процесс испарения капли, имеют две принципиальные особенности.
Одна из них связана с ростом кристалла и потерей каплей сферической формы. Согласно (3.4)-(3.5), при испарении капли с кристаллом (при х = х8) давление насыщенного пара Ра и химический потенциал растворителя уо = у0 даз + КТа 1п Ра возрастают до дости-
/3
жения каплей объема У^ = 7ГЬ3. При этом кубический кристалл вписан в сферу. При дальнейшем испарении капля теряет сферическую форму (рис. 4), после чего Ра и уо уменьшаются. Таким образом при У^ = 7ГЬ3 и конденсация, и испарение приводят к понижению давления пара Ра и химического потенциала уо (точка В на рис. 5). Это неустойчивое состояние при некоторых условиях можно рассматривать как одну из точек бифуркации, в которой меняется направление потока растворителя около капли.
Рис. 5. Зависимость давления насыщенного пара над каплями чистого растворителя (верхняя кривая); над каплями насыщенного раствора (xs) (верхняя штриховая линия); над каплями, содержащими кристалл размера L = 5, 10, 20 нм (штриховые линии); над испаряющимися каплями размера 100 нм, 1, 10, 100 мкм без кристалла (линии с маркером О); над испаряющимися каплями с кристаллом (линии с маркером ■); верхняя граница области существования сферических капель с кристаллом (жирная кривая). P(ж) — давление над плоской поверхностью чистого растворителя.
Другая особенность возникает при полном растворении кристалла. При L = 0, U\c = 0 концентрация x становится функцией объема капли (3.2). В соответствии с уравнениями (3.5)-(3.6) конденсация капли до растворения кристалла (x = xs) сопровождается уменьшением давления пара Pd и химического потенциала уо. После полного растворения конденсация понижает значение x, что сопровождается ростом Pd и уо.
В последнем случае нарушается симметрия процессов конденсации и испарения. При испарении зарождение кристаллической фазы происходит при достижении некоторого пересыщения (x > xs ). Степень пересыщения обычно имеет стохастический характер и зависит от размера капли [8].
Эволюцию испаряющихся (конденсирующихся) капель демонстрирует рисунок 5. Давление пара над каплями понижается при испарении летучего компонента в соответствии с (3.5). На примере капли с Rd = 10 мкм схематически показано изменение давления после возникновения зародыша кристаллической фазы.
При испарении капли раствора до момента зарождения кристалла (точка А на рисунке 5 соответствует образованию зародыша новой фазы размером L = 5 нм) концентрация нелетучего компонента растет. При этом в соответствии с (3.5) уменьшаются Pd и уо. Рост кристалла понижает пересыщение (процесс А ^ В), при этом давление повышается в соответствии с (3.5) до равновесного значения, соответствующего концентрации насыщенного раствора xs (или определяемого (3.6)). При дальнейшем испарении и росте кристалла капля
становится несферической (см. рис. 4). Этому состоянию соответствует точка В на рисунке 5. После этого при испарении радиус кривизны Rd возрастает до то при L ~ const, и на поверхности кристалла образуется пленка раствора. Если давление во внешней среде ниже, чем давление в точке B, испарение будет продолжаться с понижением давления пара над каплей (точка С на рисунке 5) путем образования островковых пленок и повышения концентрации в них за счет повышения растворимости, описанного в [8]. Состояние несферических капель (точка С на рисунке 5) показано схематически. Положение точки С определяется давлением летучего компонента во внешней среде, а радиус сопоставим с радиусом минимально возможной сферической капли, состоящей только из NaCl — это чуть меньше радиуса, при котором теряется сферическая форма.
Радиус капли при зарождении кристалла (точка A на рисунке 5) зависит от случайной величины (пересыщения). Неопределенность положения точек А и С может быть одной из причин нарушения периодичности наблюдаемых осцилляций.
Положение точек А, В и С (см. рис. 5) существенно зависит от температуры капли. В окрестности точек А и В давление Pd меняется при испарении или конденсации немонотонно. При некоторых условиях такое изменение Pd способно приводить к изменению направления потока j в окрестности капли.
5. Заключение
Упрощение модели [1] не отрицает возможности образования метастабильных фаз и их роли при возникновении осцилляций. Согласно наблюдениям, амплитуда изменения размера капли, воспроизводимость и устойчивость режима осцилляций в трехкомпонентной системе H2O-NaCl-KCl выше, чем в бинарных системах с единственным растворенным компонентом.
Описанная экспериментальная методика позволяет не только наблюдать эффект осцилляции фазовых превращений, но и управлять условиями, в которых эти осцилляции существуют. Эта возможность актуальна для создания и моделирования на основе описанного эффекта запоминающих устройств на фазовых переходах, сенсоров, элементов искусственных нейронных сетей.
Предлагаемая модель может представлять интерес для специалистов в области теории фазовых превращений и динамических систем.
Список литературы
[1] Федосеев В. Б., Максимов М.В. Осциллирующие фазовые превращения раствор - кристалл -раствор в системе состава KCl-NaCl-H2O // Письма в ЖЭТФ, 2015, т. 101, №6, с. 424-427.
[2] Ванаг В. К. Волны и динамические структуры в реакционно-диффузионных системах. Реакция Белоусова- Жаботинского в обращенной микроэмульсии // УФН, 2004, т. 174, №9, с. 991-1010.
[3] Liesegang R. E. Uber einige Eigenschaften von Gallerten // Naturwiss. Wochenschr., 1896, vol. 11, no. 30, pp. 353-362.
[4] Ильин А. М., Марков Б. А. Нелинейное уравнение диффузии и кольца Лизеганга // Докл. Акад. наук, 2011, т. 440, №2, с. 164-167.
[5] Tanaka T., Hara S. Thermodynamic evaluation of nano-particle binary alloy phase diagrams // Z. Metallkd., 2001, vol. 92, no. 11, pp. 1236-1241.
[6] Eichhammer Y., Heyns M., Moelans N. Calculation of phase equilibria for an alloy nanoparticle in contact with a solid nanowire // Calphad, 2011, vol. 35, no. 2, pp. 173-182.
[7] Федосеев В. Б., Шишулин А. В., Титаева Е. К., Федосеева Е.Н. О возможности образования кристалла твердого раствора NaCl-KCl из водного раствора при нормальной температуре в системах малого объема // ФТТ, 2016, т. 58, №10, с. 2020-2025.
[8] Федосеев В. Б., Федосеева Е.Н. Состояния пересыщенного раствора в системах ограниченного размера // Письма в ЖЭТФ, 2013, т. 97, №7, с. 473-478.
[9] Франк-Каменецкий Д. А. Диффузия и теплопередача в химической кинетике. Москва: Наука, 1987. 490 с.
[10] Справочник по растворимости: Т. 1. Бинарные системы: Кн. 2 / В. В. Кафаров (ред.). Москва-Ленинград: АН СССР, 1961. С. 108.
[11] Щукин Е. Л., Перцов А. В., Амелина Е. А. Коллоидная химия. 3-е изд. Москва: Высшая школа, 2004. 445 с.
[12] Рид Р., Праусниц Д., Шервуд Т. Свойства газов и жидкостей. Ленинград: Химия, 1982. 592 с.
[13] Кучма А. Е., Щёкин А. К., Куни Ф.М. Динамика изменения размера и состава закритической капли при бинарной конденсации // Коллоид. журн., 2011, т. 73, №2, с. 215-224.
[14] Mattila T., Kulmala M., Vesala T. On the condensational growth of a multicomponent droplet // J. Aerosol Sci., 1997, vol.28, no. 4, pp. 553-564.
[15] Мартюкова Д. С., Кучма А.Е., Щёкин А. К. Динамика изменения размера и состава бинарной капли в атмосфере двух конденсирующихся паров и пассивного газа при произвольных начальных условиях // Коллоид. журн., 2013, т. 75, №5, с. 625-632.
[16] Саверченко В. И., Фисенко С. П., Ходыко Ю.А. Кинетика испарения бинарной пиколитровой капли на подложке при пониженном давлении // Коллоид. журн., 2015, т. 77, № 1, с. 79-84.
[17] Тарасевич Ю.Ю., Исакова О. П., Кондухов В. В., Савицкая А. В. Влияние режима испарения на пространственное перераспределение компонентов в испаряющейся капле жидкости на твердой горизонтальной подложке // ЖТФ, 2010, т. 80, №5, с. 45-53.
[18] Engberg R.F., Wegener M., Kenig E.Y. The influence of Marangoni convection on fluid dynamics of oscillating single rising droplets // Chem. Eng. Sci., 2014, vol. 117, pp. 114-124.
Solution —gas and solution —crystal oscillatory phase transitions in the drops of aqueous solutions with one crystallizable component
Victor B. Fedoseev
Razuvaev Institute of Organometallic Chemistry of RAS
ul. Tropinina 49, GSP-445, Nizhny Novgorod, 603950, Russia
Long-term aperiodic oscillatory phase transitions of gas-solution and solution-crystal have been reproduced experimentally in the ensemble of aqueous droplets with a dissolved component. The new observations can simplify the dynamic model of phase transitions of the oscillating mode considered before. It reduces the number of independent variables and the need to consider the formation of metastable phases. The model establishes a relation between the speed and direction of solvent flow in the gas and the vapor pressure in the drop neighborhood, temperature of the drop, the rate of change in temperature, the rate of change in volumes of the drop and the crystal. The analysis has revealed that the presence of crystalline phases in the system causes at least two singularities (bifurcations) of the chemical potential of the volatile component with respect to its quantities in the drop.
MSC 2010: 92E99, 80A20, 74A50, 82C26, 70K50
Keywords: phase transitions, oscillation, evaporation, condensation, drop of solution, crystallization
Received January 07, 2017, accepted May 15, 2017
Citation: Rus. J. Nonlin. Dyn., 2017, vol. 13, no. 2, pp. 195-206 (Russian)