[论文解读] A sharp-interface model and its numerical approximation for solid-state dewetting with axisymmetric geometry

Quan Zhao|arXiv (Cornell University)|Nov 7, 2017

Fluid Dynamics and Thin Films参考文献 16被引用 4

一句话总结

本文提出了一种严格的轴对称尖锐界面模型，用于描述由表面扩散和接触线运动主导的固态脱湿过程，该模型通过能量变分原理推导得出。采用半隐式参数化有限元方法（PFEM），模拟了复杂的形态演化过程——包括边缘凸起形成、孔洞产生以及环形薄膜收缩——结果表明，空洞形成在很大程度上取决于各向异性参数σ和初始薄膜几何形状。

ABSTRACT

Based on the thermodynamic variation, we rigorously derive the sharp-interface model for solid-state dewetting on a flat substrate in the form of cylindrical symmetry. The governing equations for the model belong to fourth-order geometric curve evolution partial differential equations, with proper boundary conditions such that the total volume of the system is conserved and the total energy is dissipative during the time evolution. We propose a variational formulation for the sharp-interface model and then apply the parametric finite element method for solving it efficiently. Extensive numerical simulation results are presented lastly to demonstrate the morphological characteristics for solid-state dewetting.

研究动机与目标

• 通过热力学变分法推导轴对称几何下固态脱湿的数学上严谨的尖锐界面模型。
• 开发一种稳定且高效的数值格式，用于求解具有动态边界条件的四阶几何型偏微分方程。
• 研究三维固态脱湿中的复杂形态演化，如边缘回缩、断裂以及空洞形成。
• 分析表面能各向异性和初始薄膜几何形状对平衡态与动态形态的影响。
• 通过验证轴对称极限下的关键动力学行为，为未来全三维模拟奠定基础。

提出的方法

• 在轴对称约束下，通过对总表面能泛函进行第一变分，推导出轴向轮廓的四阶几何曲线演化偏微分方程。
• 边界条件确保体积守恒与能量耗散，接触线运动由松弛接触角条件控制。
• 引入变分格式以保证能量稳定性，并实现高效的数值处理。
• 在半隐式混合有限元框架中应用参数化有限元方法（PFEM），通过切向速度控制保持网格质量。
• 该方法隐式处理动态网格重分布，避免重新划分网格，相比显式标记粒子方法可采用更大的时间步长。
• 数值模拟在以弧长参数化的轴向曲线上进行，表面扩散与接触线运动实现完全耦合。

实验结果

研究问题

• RQ1 表面能的各向异性如何影响轴对称固态脱湿的形态演化？
• RQ2 在环形薄膜收缩过程中，何种条件会导致空洞形成？
• RQ3 初始薄膜几何形状与各向异性参数σ如何共同决定最终的平衡形态？
• RQ4 在轴对称固态脱湿中，主导不稳定性及转变机制（如边缘回缩与断裂）是什么？
• RQ5 PFEM能否以稳定且高效的方式准确捕捉如薄膜断裂与自融合等复杂拓扑变化？

主要发现

• 该模型正确地将薄膜/气相界面的动态演化描述为H⁻¹梯度流，接触线运动描述为L²梯度流，确保能量耗散与体积守恒。
• 数值模拟揭示了四种不同的形态演化类型：(1) 稳定单层薄膜，(2) 形成空洞的环形薄膜，(3) 不形成空洞的环形薄膜，(4) 分裂为中央薄膜与外环结构。
• 各向同性与各向异性表面能的相图表明，空洞形成强烈依赖于各向异性参数σ和初始薄膜高径比。
• PFEM在显著更大的时间步长下实现稳定模拟，并且网格质量优于传统标记粒子方法。
• 轴对称模型揭示了关键不稳定性，如轴向收缩、边缘凸起形成以及方位角方向的类似Rayleigh的断裂，尽管未完全捕捉全三维效应。
• 该模型预测仅当环形薄膜足够“肥厚”时，才能在方位角断裂前向中心迁移并合并，与实验观测一致。