Твёрдотельное смачивание осесимметричной тонкой плёнки на осесимметричных криволинейных подложках: моделирование и симуляция

Основная информация

• ID статьи: 2501.00783
• Название: Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: modeling and simulation
• Авторы: Zhenghua Duan, Meng Li*, Chunjie Zhou (Школа математики и статистики, Чжэнчжоуский университет)
• Классификация: math.NA cs.NA
• Дата публикации: 3 января 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2501.00783

Аннотация

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

Научный контекст и мотивация

Значимость проблемы

Твёрдотельное смачивание (Solid-State Dewetting, SSD) — это явление нестабильности твёрдых тонких плёнок при температурах значительно ниже точки плавления, при котором плёнка подвергается смачиванию или агломерации, образуя сложные морфологические изменения и формируя малые частицы на подложке. Это явление находит широкое применение в оптических и магнитных устройствах, тонкоплёночных системах, датчиках и формировании катализаторов.

Ограничения существующих исследований

1. Ограничения по размерности: Существующие теоретические исследования сосредоточены главным образом на плоских подложках с ограниченным вниманием к топологически структурированным подложкам
2. Двумерные ограничения: Имеющиеся исследования сосредоточены на двумерных сценариях с ограниченными исследованиями SSD на трёхмерных криволинейных подложках
3. Недостаток численных методов: Отсутствуют эффективные численные методы, одновременно сохраняющие энергетическую устойчивость и сохранение объёма

Научная мотивация

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

Основные вклады

1. Теоретическое моделирование: Посредством термодинамической вариационной формулировки выведена модель с острым интерфейсом для SSD трёхмерной осесимметричной плёнки на осесимметричной криволинейной подложке, охватывающая случаи слабой и сильной анизотропии
2. Инновации в численных методах:
   • Введение симметризованной матрицы поверхностной энергии и вывод новой симметризованной вариационной формулировки
   • Проектирование безусловно энергетически устойчивой параметрической аппроксимации методом конечных элементов посредством тщательной дискретизации граничных членов
   • Разработка структурно-сохраняющего метода, одновременно сохраняющего сохранение объёма и энергетическую устойчивость
3. Верификация алгоритма: Посредством обширных численных экспериментов верифицирована сходимость и структурно-сохраняющие свойства предложенного метода
4. Исследование физических явлений: Изучены различные интересные явления, включая миграцию малых частиц на подложках с положительной/отрицательной кривизной, события защемления и сжатие краёв

Подробное описание методов

Определение задачи

Исследование процесса эволюции твёрдотельного смачивания осесимметричной плёнки на осесимметричной криволинейной подложке, включая:

• Входные данные: начальная морфология плёнки, геометрия подложки, параметры материала
• Выходные данные: траектория эволюции плёнки во времени
• Ограничения: сохранение объёма, убывание энергии, граничные условия контактной линии

Математическая модель

Полная свободная энергия

Полная свободная энергия системы выражается как: $W = \int\int_S \gamma_{FV}(N) dS + (\gamma_{FS} - \gamma_{VS}) A(\Gamma_o/\Gamma_i)$

где $S$ — интерфейс плёнка/пар, $\Gamma_i$ и $\Gamma_o$ — внутренняя и внешняя контактные линии.

Осесимметричное упрощение

При осесимметричном предположении поверхность параметризуется как: $S(s,\phi) := (r(s)\cos\phi, r(s)\sin\phi, z(s))$

Полная энергия упрощается до: $W = 2\pi\int_0^L \gamma(\theta)r ds + 2\pi(\gamma_{FS} - \gamma_{VS})\int_{c_l}^{c_r} \hat{x}(c)dc$

Модель с острым интерфейсом

Посредством термодинамической вариационной формулировки получено управляющее уравнение: $X_t \cdot n = \frac{1}{r}(r\mu_s)_s, \quad 0 < s < L(t), t > 0$

Химический потенциал определяется как: $\mu = (\gamma(\theta) + \gamma''(\theta))\kappa - \frac{\gamma(\theta)z_s + \gamma'(\theta)r_s}{r}$

Граничные условия включают:

1. Условие контактной линии: обеспечение остаётся контактная линия на криволинейной подложке
2. Условие релаксации контактного угла: описание кинетики эволюции контактного угла
3. Условие нулевого потока массы: гарантирование сохранения объёма

Симметризованная вариационная формулировка

Введена симметричная матрица: $B(\theta) = \begin{pmatrix} \gamma(\theta) & -\gamma'(\theta) \\ \gamma'(\theta) & \gamma(\theta) \end{pmatrix} \begin{pmatrix} \cos 2\theta & \sin 2\theta \\ \sin 2\theta & -\cos 2\theta \end{pmatrix} + S(\theta)\left[\frac{1}{2}I - \frac{1}{2}\begin{pmatrix} \cos 2\theta & \sin 2\theta \\ \sin 2\theta & -\cos 2\theta \end{pmatrix}\right]$

Посредством этой матрицы уравнение химического потенциала может быть записано как: $r\mu n = \partial_s[rB(\theta)\partial_s X] - \gamma(\theta)e_1$

Численный метод

Параметрический метод конечных элементов

Применён параметрический метод конечных элементов (PFEM) для решения вариационной формулировки:

• Дискретизация по времени: полунеявная схема
• Пространственная дискретизация: линейные конечные элементы
• Обработка границ: тщательно спроектированная дискретизация граничных членов

Метод энергетической устойчивости

Спроектированная численная схема удовлетворяет дискретной энергетической устойчивости: $W(X^{m+1}) - W(X^m) \leq 0$

Структурно-сохраняющий метод

Посредством корректирующего члена $\delta f^{m+\frac{1}{2}}$ достигается точное сохранение объёма: $V(X^{m+1}) - V(X^m) = 0$

Экспериментальная установка

Тестовые случаи

В статье рассмотрены три типичные начальные конфигурации:

1. Сферическая плёнка: сферическая плёнка радиусом 1,5 на полусферической подложке
2. Кольцевая плёнка: осесимметричная кольцевая плёнка толщиной 0,5 на полусферической подложке
3. Сферическая плёнка: сферическая плёнка радиусом 1,8 на сферической подложке

Параметры материала

• Подвижность контактной линии: $\eta = 100$
• Параметр поверхностной энергии: $\sigma = -\frac{\sqrt{3}}{2}$
• Анизотропная поверхностная энергия: $\gamma(\theta) = 1 + \beta\cos(4\theta)$
  • $\beta = 0$: изотропный случай
  • $0 \leq \beta \leq \frac{1}{15}$: слабая анизотропия
  • $\beta > \frac{1}{15}$: сильная анизотропия

Численные параметры

• Размер сетки: $h = 2^{-7}$ до $2^{-8}$
• Временной шаг: $\Delta t = 2^{-9}$ до $2^{-10}$
• Допуск сходимости: $10^{-8}$

Результаты экспериментов

Верификация сходимости

Посредством метрики расстояния многообразия численная ошибка измеряется как: $e^{h,\Delta t}(t) := |(\Omega^{h,\Delta t}(t)\backslash\Omega^r(t)) \cup (\Omega^r(t)\backslash\Omega^{h,\Delta t}(t))|$

Результаты показывают:

• Порядок пространственной сходимости: второй порядок, соответствующий теоретическим предсказаниям
• Порядок временной сходимости: первый порядок
• Сохранение хорошей сходимости в случаях изотропии, слабой и сильной анизотропии

Верификация структурно-сохраняющих свойств

Энергетическая устойчивость

Численные эксперименты подтверждают монотонное убывание дискретной энергии:

• Кривые энергии при различных временных шагах показывают убывающий тренд
• Интенсивность анизотропии не влияет на энергетическую устойчивость

Сохранение объёма

• Метод энергетической устойчивости: относительная ошибка объёма порядка $10^{-3}$
• Структурно-сохраняющий метод: достижение сохранения объёма с машинной точностью (ошибка $\sim 10^{-15}$)

Исследование физических явлений

Явление миграции частиц

На криволинейной подложке наблюдаются:

• Подложка с положительной кривизной: малые частицы стремятся мигрировать в область с меньшей кривизной
• Подложка с отрицательной кривизной: поведение частиц зависит от интенсивности анизотропии
• Эффект анизотропии: сильная анизотропия значительно влияет на скорость миграции

Явление защемления

Длинная плёнка подвергается защемлению в процессе эволюции:

• После защемления переинициализируется как две независимые части
• В конечном итоге образуются более мелкие островковые и кольцевые плёнки

Сжатие краёв

Сжатие краёв полубесконечной ступенчатой плёнки:

• Плёнка постепенно переползает через угол подложки
• Чем больше интенсивность анизотропии, тем медленнее скорость сжатия

Связанные работы

Теоретическое моделирование

• Srolovitz & Safran (1986): первое предложение модели с острым интерфейсом
• Jiang и др. (2012): метод фазового поля для обработки топологических изменений
• Zhao и др. (2019): модель редукции размерности для осесимметричной геометрии

Численные методы

• Barrett, Garcke & Nürnberg: параметрический метод конечных элементов (метод BGN)
• Bao & Zhao (2021): структурно-сохраняющий PFEM
• Li и др. (2023): симметризованный метод

Исследования криволинейных подложек

• Jiang и др. (2018): SSD на двумерной криволинейной подложке
• Bao и др. (2024): структурно-сохраняющий метод для двумерной криволинейной подложки

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

Заключение и обсуждение

Основные выводы

1. Успешно установлена полная математическая модель SSD на осесимметричной криволинейной подложке
2. Разработан эффективный численный метод, одновременно сохраняющий энергетическую устойчивость и сохранение объёма
3. Выявлено важное влияние геометрии криволинейной подложки на поведение эволюции плёнки

Ограничения

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

Направления будущих исследований

1. Расширение на неосесимметричную геометрию
2. Рассмотрение эффектов многофизического взаимодействия
3. Разработка методов адаптивной сетки
4. Применение к практическим промышленным задачам

Глубокая оценка

Преимущества

1. Теоретическая строгость: вывод из первых принципов термодинамики, полный и строгий
2. Инновационность методов: симметризованная вариационная формулировка и структурно-сохраняющая схема обладают инновационностью
3. Численная устойчивость: безусловная энергетическая устойчивость избегает ограничений на временной шаг
4. Физическая значимость: численные результаты раскрывают богатые физические явления
5. Достаточность экспериментов: включает верификацию сходимости, устойчивости и применимости

Недостатки

1. Область применения: осесимметричное предположение ограничивает универсальность метода
2. Вычислительная эффективность: вычислительные затраты для крупномасштабных трёхмерных задач могут быть высокими
3. Чувствительность параметров: выбор некоторых параметров материала требует более подробного обсуждения
4. Экспериментальная верификация: отсутствует сравнение с данными реальных экспериментов

Влияние

1. Научный вклад: предоставляет важные теоретические и численные инструменты для исследования SSD на криволинейных подложках
2. Практическая ценность: может быть применён к проектированию наноустройств с осесимметричными характеристиками
3. Методологическое значение: структурно-сохраняющие численные методы могут быть обобщены на другие задачи геометрической эволюции

Сценарии применения

• Проектирование тонкоплёночных устройств с осесимметричными характеристиками
• Подготовка наноструктур на криволинейных поверхностях
• Контроль морфологии частиц катализаторов
• Поверхностная инженерия оптических и магнитных устройств

Библиография

Статья содержит 55 высококачественных ссылок, охватывающих основные работы в теоретическом моделировании SSD, численных методах и экспериментальных исследованиях, с полным и авторитетным обзором литературы.

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