Поверхностно-ориентированные двумерные волновые структуры в ограниченных бинарных смесях

Язык труда и переводы:
УДК:
538.91
Дата публикации:
08 декабря 2022, 22:51
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Краснюк Игорь Борисович
Донецкий физико-технический институт имени А.А. Галкина
Аннотация:
Рассмотрен процесс формирования двумерных структур в бинарных смесях, к примеру, полимерных смесях или бинарных сплавах, который описывается уравнением Кана — Хиллиарда с динамическими граничными условиями с обратной связью. Изучаются асимптотические периодические импульсные структуры релаксационного, предтурбулентного и турбулентного типа с конечным, исчисляемым или неисчисляемым числом разрывов на периоде. Асимптотические структуры являются элементами аттрактора в начально-краевой задаче. Рассмотрено приложение модели к случаю бинарного полимера квадратной формы. Выполнено компьютерное моделирование, описывающее формирование кристаллитов в расплаве.
Ключевые слова:
уравнение Кана — Хиллиарда, нелинейные динамические граничные условия, ограниченные распределения, релаксация, предтурбулентный тип, турбулентный тип
Основной текст труда

Введение

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

Формулировка задачи

Мы рассматриваем краевую задачу, которая описывает процесс формирование двумерных пространственно-временных структур в симметричных диблок-сополимерах, ограниченных единичным кубом. Эволюция одной компоненты бинарной смеси моделируется двумерным конвективным уравнением Кана — Хиллиарда [1] с динамическими граничными условиями с обратной связью. Эта модель описывает распределение концентрации одной из компонент в бинарном расплаве после охлаждения поверхности. Мы полагаем, что скорости образования зародышей кристаллитов по сторонам квадрата являются нелинейными функциями концентрации и температуры, что соответствует имеющимся экспериментальным данным [2]. Частота зародышеобразования может быть аппроксимирована параболической функцией с отрицательной второй производной по температуре кристаллизации [2, Fig. 60].

Мы рассматриваем граничные условия в виде

u_{x}\left({t,0,y}\right)=F_{1}\left[{u_{x}\left({t,0,y}\right)}\right],\,\,\,\,u_{x}\left({t,l,y}\right)=F_{2}\left[{u_{x}\left({t,l,y}\right)}\right];                        (1)

u_{y}\left({t,x,0}\right)=G_{1}\left[{u_{y}\left({t,x,0}\right)}\right],\,\,\,\,u_{y}\left({t,x,l}\right)=G_{2}\left[{u_{y}\left({t,x,l}\right)}\right],                       (2)

где F_{1},\,\,F_{2},\,\,G_{1},\,G_{2}  — нелинейные функции.

Уравнение Кана — Хиллиарда может быть линеаризовано  в окрестности неупорядоченной фазы u=0 , и решение задачи может быть найдено в виде

u\left({t,x,y}\right)=f\left({t+a_{1}x}\right)+g\left({t+a_{2}y}\right),                                  (3)

гдеe f и g — неизвестные функции, a_{1} и a_{2} — параметры задачи. Она может быть сведена к разностному уравнению

f\left({t+{\frac {l}{a_{1}}}}\right)=\Phi _{1}\left[{f\left(t\right)}\right]                                         (4)

для функции f , и такому же уравнению для функции g , но на другом множестве \Phi _{2} . Решениями данных уравнений являются асимптотически периодические кусочно-постоянные импульсы [3]. Асимптотическое решение u(t,x,y)   является суммой этих импульсов.

Рассмотрим бинарную смесь, расположенную внутри квадрата стороной l . Положим, что на границе смесь переохлаждена до температуры T<T_{g} , где T_{g} — температура поверхностного распада. Мы рассматриваем объемную смесь полимеров, где в неупорядоченном состоянии вначале u=0  и u(t,x,y) — одна из компонент бинарной смеси.

Конвективное уравнение Кана — Хиллиарда имеет вид [1]

u_{t}-a_{1}u_{x}-a_{2}u_{y}=\left[{k_{1}u\left({1-u^{2}}\right)+u_{xx}}\right]_{xx}+\left[{k_{2}u\left({1-u^{2}}\right)+u_{yy}}\right]_{yy}.               (5)

Начальное условие имеет вид

u\left({x,0}\right)=u_{0}\left(x\right).                                            (6)

В настоящей работе рассматривается обобщенное граничное условие

\tau _{u}u_{x}=\Phi \left[{\partial _{n}u,\,\,\partial _{n}^{3}u,\,\,u}\right],                                           (7)

где \Phi — нелинейная функция с обратной связью.

Мы ограничимся изучением конвективного уравнения Кана-Хиллиарда в ограниченной области \Omega \subset R^{3} . Для простоты положим \Omega — куб со стороной l . На каждой поверхности в прострастве x,y имеют место динамические граничные условия с дополнительными классическими  стационарными граничными условиями

\partial _{n}^{3}u=\Psi \left[{\partial _{n}u,\,\,\partial _{n}^{3}u,\,\,u}\right]                                            (8)

на \partial \Omega ,  где  \Psi \left[\cdot \right] — данный функционал.

Приведение задачи к разностному уравнению непрерывного типа

Определим функцию параметра порядка u=\left(N_{A}-N_{B}\right)/\left(N_{A}+N_{B}\right) , где цепь субэлементов N_{A}  ковалентно связана с цепью N_{B}  и N_{A}+N_{B}=1 ; a^{-1} — безразмерная скорость конвекции.   k:=k-k_{c} — параметр Флори — Хиггинса, k_{c} — критическое значение фазового распада. При k>k_{c}  в окрестности неупорядоченной фазы u=0  наблюдаются синусоидальные флуктуации. При  флуктуации параметра порядка монотонны.

Вблизи среднего значения u=0   (при N_{A}=N_{B} ), линеаризованное уравнение Кана — Хиллиарда

u_{t}+a_{1}u_{x}+a_{2}u_{y}=\left[{k_{1}u+u_{xx}}\right]_{xx}+\left[{k_{2}u+u_{yy}}\right]_{yy}                                       (9)

имеет решение

u\left({t,\,x,\,y}\right)=f\left(t-x/a_{1}\right)+g\left(t-y/a_{2}\right)                                       (10)

при a_{1},a_{2}>0 . Тогда мы получаем два обычных разностных уравнения:

\lambda _{1}f_{\zeta \zeta }\left(\zeta \right)+f_{\zeta \zeta \zeta \zeta }\left(\zeta \right)=0;                                                            (11)

\lambda _{2}g_{\eta \eta }\left(\eta \right)+g_{\eta \eta \eta \eta }\left(\eta \right)=0                                                              (12)

с \zeta =t+x/a_{1} , \eta =t-{y/{a_{2}}}   и   \lambda _{1}=k_{1}a_{1}^{2},\,\,\,\lambda _{2}=k_{2}a_{2}^{2} .

Положим \lambda _{1}<0 . По определению, \lambda _{1}=\chi -\chi _{c} , где \chi — параметр Флори — Хиггинса, который характеризует взаимодействие между атомами в полимерных смесях. \chi <\chi _{c} , где \chi _{c} — критический параметр фазового расслоения неупорядоченной фазы полимерной смеси на две упорядоченные Если \chi -\chi _{c} , имеются монотонные распределения концентраций компонент бинарной смеси. В этом случае доминируют поверхностные возмущения и можно говорить о поверхностно-индуцированном упорядочении в объеме расплава при , где T_{c} — критическая температура фазового расслоения на две упорядоченные фазы.

Для нахождения точного решения положим a_{1},a_{2}>0 . Тогда функция f(x(t),t)  является константой вдоль характеристики {\frac {dx(t)}{dt}}=a_{1} , функция g(y(t),t) — вдоль {\frac {dy(t)}{dt}}=a_{2} . Интегрируя уравнения (11) и (12) от точек \zeta _{0}=t,\,\,\,\eta _{0}=t  до \zeta _{0}=t+{\frac {1}{a_{1}}},\eta _{0}=t+{\frac {1}{a_{2}}} вдоль характеристик {\frac {dx(t)}{dt}}=a_{1} и {\frac {dy(t)}{dt}}=a_{2} , получим

f'\left({l,t}\right)-f'\left({0,t}\right)=C_{1}\exp \left({\alpha _{1}t}\right)+C_{2}\exp \left({\alpha _{2}t}\right),                             (13)

где \alpha _{1}={\sqrt {\lambda _{1}}},\,\,\,\,\,\alpha _{2}={\sqrt {\lambda _{2}}}  и C_{1},\,\,C_{2}\in R .  Ограничимся изучением случая C_{1}=0 . Тогда f'\left({l,t}\right)-f'\left({0,t}\right)=C_{2}\exp \left({\alpha _{2}t}\right)  где a_{2}<0 , и f'\left({l,t}\right)-f'\left({0,t}\right)\Rightarrow 0  при t\to +\infty   для всех точек t\in R^{+} . Далее, интегрируя (11) от \zeta =t  до \zeta =t-l/a_{1} , получаем

f'''\left(t\right)-f'''\left(t-{\frac {l}{a_{1}}}\right)+\lambda _{1}\left(f'\left(t\right)-f'\left(t-{\frac {l}{a_{1}}}\right)\right)=0.              (14)

Мы используем двойные граничные условия Неймана. Из разностного уравнения следует функциональное

G_{1}\left[f\left(t\right)\right]-F_{2}\left[f\left(t-l/{a_{1}}\right)\right]+\lambda _{1}\left(G_{1}\left[f\left(t\right)\right]-F_{2}\left[f\left(t-{l}/{a_{1}}\right)\right]\right)=0.                (15)

Его решение

f(t)=\Phi _{1}\left[f\left(t-{l}/{a_{1}}\right)\right],                                        (16)

где -l/a_{1}\geq t<0 . Множество \Phi _{1}\in C^{2}\left({I,I}\right) ,   I — открытый интервал. Рассматриваем только класс унимодальных отображений [3], с одной особой точкой.

Аналогично получаем разностное уравнение для g(\eta )

g(t)=\Phi _{2}\left[g\left(t-{l}/{a_{2}}\right)\right],                                         (17)

где a_{2}>0 . Решение уравнения  (17) может быть найдено пошаговой итерацией начальной функции h_{1}\left(t\right) на -l/a_{1}\leq t<0 , \Phi _{1}:I\to I . Обычно предельное решение p_{1}(t)  является кусочно-постоянной апериодической функцией с ограниченным, счетным или несчетным числом точек разрыва \Gamma  на периоде. Если \Gamma конечно, имеем дело с решением релаксационного типа. Если \Gamma  счетно, имеем предельное распределение пре-турбулентного типа. Если \Gamma несчетно — турбулентного (рис. 1).

Рис. 1. Предельное распределение релаксационного типа в двумерном случае

Здесь p_{1}\left(t\right)\in P^{+}  почти для всех точек t\in R^{+}  исключая набор точек нулевой меры Лебега, где P^{+}   набор фиксированных точек отображения \Phi _{1}  [3]. Подвижность атомов k_{1}  и скорость конвекции a_{1}  являются параметрами бифуркации. Подвижность k_{1}=k_{1}\left(\theta \right)  может быть монотонной функцией безразмерной температуры \theta . Особый случай наблюдается при k_{1}:\sim k_{1}\theta ^{-1} . Если смесь охлаждается на плоских стенах, на поверхности появляются осциллирующие структуры с немонотонными амплитудами. Эти волноподобные структуры распространяются в объем смеси с экспоненциально затухающими амплитудами. При t\to +\infty , пространственно-темпоральные структуры стремятся к

p\left({\zeta ,\eta }\right)=p_{1}\left(\zeta \right)+p_{2}\left(\eta \right),                                               (18)

где  \zeta =t-x/a_{1} и   \eta =t-{y/{a_{2}}} . Здесь p_{1}\left(\zeta \right)\in P_{1}^{+}  и p_{1}\left(\zeta \right)\in P_{2}^{+} , где P_{1}^{+},\,\,P_{2}^{+}   — притягивающие фиксированные точки отображений \Phi _{1}  и \Phi _{2} . На рис. 2 представлены результаты компьютерного моделирования предельных распределений концентраций.

Рис. 2. Предельные распределения пре-турбулентного типа со счетным числом особых точек на периоде

Предельные распределения пре-турбулентного типа имеют утолщенные линии. Для предельных распределений турбулентного типа имеются Канторовы множества линий разрыва волновых фронтов распространения концентрации по объему. Результатом является направленное упорядочение, например в направлениях Ox  и Oy  мы имеем решения релаксационного типа, или один из следующих случаев: релаксация \oplus  пре-турбулентный тип, релаксация \oplus турбулентный тип, претурбулентный \oplus турбулентный тип и т. д.

Выводы

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

Литература
  1. Eden A., Kalantarov V.K. The convective Cahn-Hilliard equation. Applied Mathematics Letters, 2007, vol. 20 (4), pp. 455—461.
  2. Скрипов В.П., Коверда В.П. Спонтанная кристаллизация переохлажденных жидкостей. Москва, Наука, 1984, 232 с.
  3. Sharkovsky, A.N., Maistrenko Yu.L., Romanenko E.Yu. Difference equations and their applications. Ser. Mathematics and Its Applications. Dordrecht, Kluwer Academic, 1993, 250 p.
  4. Krasnyuk I.B. Spatial-temporal oscillations of the order parameter in confined diblock copolymer mixtures. International Journal of Computational Materials Science and Engineering, 2013, vol. 2 (1), pp. 1350006-1–1350006-11.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.