2025-11-15T19:55:11.527544

Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: modeling and simulation

Duan, Li, Zhou
In this work, we consider the solid-state dewetting of an axisymmetric thin film on a curved-surface substrate, with the assumption that the substrate morphology is also axisymmetric. Under the assumptions of axisymmetry, the surface evolution problem on a curved-surface substrate can be reduced to a curve evolution problem on a static curved substrate. Based on the thermodynamic variation of the anisotropic surface energy, we thoroughly derive a sharp-interface model that is governed by anisotropic surface diffusion, along with appropriate boundary conditions. The continuum system satisfies the laws of energy decay and volume conservation, which motivates the design of a structure-preserving numerical algorithm for simulating the mathematical model. By introducing a symmetrized surface energy matrix, we derive a novel symmetrized variational formulation. Then, by carefully discretizing the boundary terms of the variational formulation, we establish an unconditionally energy-stable parametric finite element approximation of the axisymmetric system. By applying an ingenious correction method, we further develop another structure-preserving method that can preserve both the energy stability and volume conservation properties. Finally, we present extensive numerical examples to demonstrate the convergence and structure-preserving properties of our proposed numerical scheme. Additionally, several interesting phenomena are explored, including the migration of 'small' particles on a curved-surface substrate generated by curves with positive or negative curvature, pinch-off events, and edge retraction.
academic

Solid-state dewetting of axisymmetric thin film on axisymmetric curved-surface substrates: modeling and simulation

基本信息

  • 论文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
  • 发表时间: 2025年1月3日 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2501.00783

摘要

本文研究轴对称薄膜在轴对称弯曲表面基底上的固态去润湿现象。在轴对称假设下，弯曲表面基底上的表面演化问题可简化为静态弯曲基底上的曲线演化问题。基于各向异性表面能的热力学变分，作者严格推导了由各向异性表面扩散控制的尖锐界面模型及相应边界条件。连续系统满足能量衰减和体积守恒定律，这为设计结构保持数值算法提供了动机。通过引入对称化表面能矩阵，推导出新颖的对称化变分公式。通过精心离散化变分公式的边界项，建立了轴对称系统的无条件能量稳定参数有限元近似。应用巧妙的修正方法，进一步开发了能够同时保持能量稳定性和体积守恒性质的结构保持方法。

研究背景与动机

问题的重要性

固态去润湿(Solid-State Dewetting, SSD)是指固体薄膜在远低于熔点的温度下发生的不稳定现象，薄膜会发生去润湿或团聚，形成复杂的形态变化并在基底上形成小颗粒。这一现象在光学和磁性器件、薄膜、传感器和催化剂形成等领域具有广泛应用。

现有研究的局限性

  1. 维度限制：现有理论研究主要集中在平面基底上，对拓扑图案化基底关注较少
  2. 二维局限：已有研究主要关注二维场景，对三维弯曲表面基底上的SSD研究有限
  3. 数值方法不足：缺乏能够同时保持能量稳定性和体积守恒的高效数值方法

研究动机

  • 三维薄膜和弯曲表面基底通常具有旋转对称性，可将复杂的三维SSD简化为弯曲基底上基于曲线的系统
  • 需要发展能够处理强各向异性效应的数值方法
  • 弯曲基底的几何特性对薄膜演化行为有重要影响

核心贡献

  1. 理论建模：通过热力学变分推导了三维轴对称薄膜在轴对称弯曲表面基底上SSD的尖锐界面模型，涵盖弱和强各向异性情况
  2. 数值方法创新：
    • 引入对称化表面能矩阵，推导出新颖的对称化变分公式
    • 通过精心离散化边界项，设计了无条件能量稳定的参数有限元近似
    • 开发了同时保持体积守恒和能量稳定的结构保持方法
  3. 算法验证：通过广泛的数值实验验证了所提方法的收敛性和结构保持性质
  4. 物理现象探索：研究了多种有趣现象，包括小颗粒在正/负曲率基底上的迁移、夹断事件和边缘收缩

方法详解

任务定义

研究轴对称薄膜在轴对称弯曲表面基底上的固态去润湿演化过程，包括：

  • 输入：初始薄膜形态、基底几何、材料参数
  • 输出：薄膜随时间的演化轨迹
  • 约束：体积守恒、能量衰减、接触线边界条件

数学模型

总自由能

系统的总自由能表示为： W=∫∫SγFV(N)dS+(γFS−γVS)A(Γo/Γi)W = \int\int_S \gamma_{FV}(N) dS + (\gamma_{FS} - \gamma_{VS}) A(\Gamma_o/\Gamma_i)

其中SS为薄膜/蒸汽界面，Γi\Gamma_i和Γo\Gamma_o为内外接触线。

轴对称简化

在轴对称假设下，表面可参数化为： S(s,ϕ):=(r(s)cos⁡ϕ,r(s)sin⁡ϕ,z(s))S(s,\phi) := (r(s)\cos\phi, r(s)\sin\phi, z(s))

总能量简化为： W=2π∫0Lγ(θ)rds+2π(γFS−γVS)∫clcrx^(c)dcW = 2\pi\int_0^L \gamma(\theta)r ds + 2\pi(\gamma_{FS} - \gamma_{VS})\int_{c_l}^{c_r} \hat{x}(c)dc

尖锐界面模型

通过热力学变分推导得到控制方程： Xt⋅n=1r(rμs)s,0<s<L(t),t>0X_t \cdot n = \frac{1}{r}(r\mu_s)_s, \quad 0 < s < L(t), t > 0

化学势定义为： μ=(γ(θ)+γ′′(θ))κ−γ(θ)zs+γ′(θ)rsr\mu = (\gamma(\theta) + \gamma''(\theta))\kappa - \frac{\gamma(\theta)z_s + \gamma'(\theta)r_s}{r}

边界条件包括：

  1. 接触线条件：确保接触线保持在弯曲基底上
  2. 松弛接触角条件：描述接触角的动力学演化
  3. 零质量通量条件：保证体积守恒

对称化变分公式

引入对称矩阵： B(θ)=(γ(θ)−γ′(θ)γ′(θ)γ(θ))(cos⁡2θsin⁡2θsin⁡2θ−cos⁡2θ)+S(θ)[12I−12(cos⁡2θsin⁡2θsin⁡2θ−cos⁡2θ)]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μn=∂s[rB(θ)∂sX]−γ(θ)e1r\mu n = \partial_s[rB(\theta)\partial_s X] - \gamma(\theta)e_1

数值方法

参数有限元方法

采用参数有限元方法(PFEM)求解变分公式：

  • 时间离散：半隐式格式
  • 空间离散：线性有限元
  • 边界处理：精心设计的边界项离散化

能量稳定方法

设计的数值格式满足离散能量稳定性： W(Xm+1)−W(Xm)≤0W(X^{m+1}) - W(X^m) \leq 0

结构保持方法

通过修正项δfm+12\delta f^{m+\frac{1}{2}}实现精确体积守恒： V(Xm+1)−V(Xm)=0V(X^{m+1}) - V(X^m) = 0

实验设置

测试案例

论文考虑三种典型初始配置：

  1. 球形薄膜：半径1.5的球形薄膜在半球基底上
  2. 环形薄膜：厚度0.5的轴对称环形薄膜在半球基底上
  3. 球形薄膜：半径1.8的球形薄膜在球形基底上

材料参数

  • 接触线迁移率：η=100\eta = 100
  • 表面能参数：σ=−32\sigma = -\frac{\sqrt{3}}{2}
  • 各向异性表面能：γ(θ)=1+βcos⁡(4θ)\gamma(\theta) = 1 + \beta\cos(4\theta)
    • β=0\beta = 0：各向同性
    • 0≤β≤1150 \leq \beta \leq \frac{1}{15}：弱各向异性
    • β>115\beta > \frac{1}{15}：强各向异性

数值参数

  • 网格尺寸：h=2−7h = 2^{-7} 到 2−82^{-8}
  • 时间步长：Δt=2−9\Delta t = 2^{-9} 到 2−102^{-10}
  • 收敛容差：10−810^{-8}

实验结果

收敛性验证

通过流形距离度量数值误差： eh,Δt(t):=∣(Ωh,Δt(t)\Ωr(t))∪(Ωr(t)\Ωh,Δt(t))∣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−310^{-3}量级
  • 结构保持方法：实现机器精度的体积守恒(误差∼10−15\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理论建模、数值方法和实验研究的主要工作，文献调研全面且权威。

总体评价：这是一篇高质量的数值分析论文，在理论建模、数值方法和应用验证方面都有重要贡献。方法具有良好的数学基础和数值性质，为三维弯曲基底上的固态去润湿研究提供了有力工具。