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

إزالة الرطوبة بالحالة الصلبة للأغشية الرقيقة محورية التماثل على الركائز ذات الأسطح المنحنية محورية التماثل: النمذجة والمحاكاة

المعلومات الأساسية

معرّف الورقة: 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) تشير إلى ظاهرة عدم الاستقرار التي تحدث في الأغشية الصلبة في درجات حرارة أقل بكثير من نقطة الانصهار، حيث تخضع الأغشية لإزالة الرطوبة أو التكتل، مما يشكل تغييرات شكلية معقدة وتشكيل جزيئات صغيرة على الركيزة. تتمتع هذه الظاهرة بتطبيقات واسعة في الأجهزة البصرية والمغناطيسية والأغشية والمستشعرات وتشكيل المحفزات.

قيود البحث الحالي

القيود البعدية: يركز البحث النظري الحالي بشكل أساسي على الركائز المستوية، مع اهتمام أقل بالركائز ذات الأنماط الطوبولوجية
القيود ثنائية الأبعاد: يركز البحث الموجود بشكل أساسي على السيناريوهات ثنائية الأبعاد، مع بحث محدود عن SSD على الركائز ذات الأسطح المنحنية ثلاثية الأبعاد
نقص الطرق الرقمية: نقص الطرق الفعالة التي يمكنها الحفاظ على الاستقرار في الطاقة وحفظ الحجم في نفس الوقت

دافع البحث

عادة ما تتمتع الأغشية ثلاثية الأبعاد والركائز ذات الأسطح المنحنية بتماثل دوراني، مما يسمح بتبسيط SSD ثلاثي الأبعاد المعقد إلى نظام قائم على المنحنيات على ركيزة منحنية
الحاجة إلى تطوير طرق رقمية قادرة على التعامل مع التأثيرات الموجهة القوية
للخصائص الهندسية للركيزة المنحنية تأثير مهم على سلوك تطور الأغشية

المساهمات الأساسية

النمذجة النظرية: اشتقاق نموذج الواجهة الحادة لـ SSD للأغشية الرقيقة محورية التماثل ثلاثية الأبعاد على الركائز ذات الأسطح المنحنية محورية التماثل من خلال التباين الديناميكي الحراري، يغطي الحالات الموجهة الضعيفة والقوية
ابتكار الطرق الرقمية:
- إدخال مصفوفة الطاقة السطحية المتماثلة واشتقاق صيغة تباينية متماثلة جديدة
- من خلال التقطيع الحذر لحدود الصيغة التباينية، تم تصميم تقريب عنصر محدود معاملي مستقر في الطاقة بشكل غير مشروط
- تطوير طريقة محافظة على البنية قادرة على الحفاظ على حفظ الحجم والاستقرار في الطاقة في نفس الوقت
التحقق من الخوارزمية: التحقق من التقارب والخصائص المحافظة على البنية للطريقة المقترحة من خلال تجارب رقمية واسعة
استكشاف الظواهر الفيزيائية: دراسة ظواهر مثيرة للاهتمام متعددة، بما في ذلك هجرة الجزيئات الصغيرة على الركائز ذات الانحناء الموجب/السالب وأحداث القطع والانكماش الحافة

شرح الطريقة

تعريف المهمة

دراسة عملية تطور إزالة الرطوبة بالحالة الصلبة للأغشية الرقيقة محورية التماثل على الركائز ذات الأسطح المنحنية محورية التماثل، بما في ذلك:

المدخلات: الشكل الأولي للغشاء، الهندسة الركيزة، معاملات المواد
المخرجات: مسار تطور الغشاء بمرور الوقت
القيود: حفظ الحجم، تناقص الطاقة، شروط الحدود لخط التلامس

النموذج الرياضي

الطاقة الحرة الكلية

يتم التعبير عن الطاقة الحرة الكلية للنظام كما يلي: $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}$

تتضمن الشروط الحدية:

شرط خط التلامس: ضمان بقاء خط التلامس على الركيزة المنحنية
شرط زاوية التلامس المرتخية: وصف الديناميكا الحركية لتطور زاوية التلامس
شرط تدفق الكتلة الصفري: ضمان حفظ الحجم

الصيغة التباينية المتماثلة

إدخال مصفوفة متماثلة: $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.5 على ركيزة نصف كروية
غشاء حلقي: غشاء حلقي محوري التماثل بسمك 0.5 على ركيزة نصف كروية
غشاء كروي: غشاء كروي بنصف قطر 1.8 على ركيزة كروية

معاملات المواد

معامل هجرة خط التلامس: $\eta = 100$
معامل الطاقة السطحية: $\sigma = -\frac{\sqrt{3}}{2}$
الطاقة السطحية الموجهة: $\gamma(\theta) = 1 + \beta\cos(4\theta)$ $γ (θ) = 1 + β cos (4 θ)$
- $\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): طريقة محافظة على البنية للركائز المنحنية ثنائية الأبعاد

توسع هذه الورقة هذه الطرق إلى الركائز المنحنية محورية التماثل ثلاثية الأبعاد، وهو تقدم مهم في هذا المجال.

الاستنتاجات والمناقشة

الاستنتاجات الرئيسية

تم إنشاء نموذج رياضي كامل لـ SSD على الركائز المنحنية محورية التماثل بنجاح
تم تطوير طريقة رقمية فعالة قادرة على الحفاظ على الاستقرار في الطاقة وحفظ الحجم
كشفت عن التأثير المهم للهندسة المنحنية للركيزة على سلوك تطور الأغشية

القيود

الافتراض محوري التماثل: يحد من الأشكال الهندسية التي يمكن معالجتها
الانتشار السطحي السائد: يتجاهل التأثيرات المرنة والتفاعلات الكيميائية
التعقيد الحسابي: لا تزال تكلفة الحساب للمشاكل ثلاثية الأبعاد مرتفعة نسبياً

الاتجاهات المستقبلية

التوسع إلى الهندسة غير محورية التماثل
الأخذ في الاعتبار تأثيرات الفيزياء المتعددة
تطوير طرق الشبكة التكيفية
التطبيق على المشاكل الصناعية الفعلية

التقييم المتعمق

المزايا

الصرامة النظرية: الاشتقاق من المبادئ الأولى للديناميكا الحرارية، كامل وصارم
ابتكار الطريقة: الصيغة التباينية المتماثلة والشكل المحافظ على البنية مبتكران
الاستقرار الرقمي: استقرار الطاقة غير المشروط، يتجنب قيود خطوة زمنية
المعنى الفيزيائي: تكشف النتائج الرقمية عن ظواهر فيزيائية غنية
كفاية التجارب: تتضمن التحقق من التقارب والاستقرار والتطبيقات

أوجه القصور

نطاق التطبيق: الافتراض محوري التماثل يحد من عمومية الطريقة
الكفاءة الحسابية: قد تكون تكلفة الحساب مرتفعة للمشاكل ثلاثية الأبعاد الكبيرة
حساسية المعاملات: اختيار بعض معاملات المواد يفتقر إلى مناقشة مفصلة
التحقق التجريبي: نقص المقارنة مع بيانات التجارب الفعلية

التأثير

المساهمة الأكاديمية: توفير أدوات نظرية ورقمية مهمة لبحث SSD على الركائز المنحنية
القيمة العملية: يمكن تطبيقها على تصميم أجهزة النانو والهندسة الرقيقة
المعنى المنهجي: يمكن تعميم طريقة المحافظة على البنية الرقمية على مشاكل تطور الهندسة الأخرى

السيناريوهات المطبقة

تصميم أجهزة الأغشية ذات الخصائص محورية التماثل
تحضير الهياكل النانوية على الأسطح المنحنية
التحكم في شكل جزيئات المحفزات
الهندسة السطحية للأجهزة البصرية والمغناطيسية

المراجع

تتضمن الورقة 55 مرجعاً عالي الجودة، تغطي الأعمال الرئيسية في النمذجة النظرية وطرق SSD والبحث التجريبي، مع مراجعة أدبية شاملة وموثوقة.

التقييم الإجمالي: هذه ورقة عالية الجودة في التحليل الرقمي، مع مساهمات مهمة في النمذجة النظرية والطرق الرقمية والتحقق التطبيقي. تتمتع الطريقة بأساس رياضي جيد وخصائص رقمية، وتوفر أداة قوية لبحث إزالة الرطوبة بالحالة الصلبة على الركائز المنحنية ثلاثية الأبعاد.