A phase field model with arbitrary misorientation dependence of grain boundary energy

1. بيان المشكلة

تُستخدم نماذج المجال الطوري (phase-field models)، وتحديداً نماذج مجال التوجه (مثل نماذج كובاي-وارن-كارتر [KWC] وهنري-ميلينثين-بلاب [HMP])، على نطاق واسع لمحاكاة نمو الحبيبات في المواد متعددة البلورات. تمثل هذه النماذج توجه البلورة كمتغير مجال مستمر ($\theta$) مقترن بمتغير ترتيب ($\phi$) يميز داخل الحبيبة عن حدود الحبيبة (GB).

القصور الجوهري:
تعاني نماذج مجال التوجه الحالية من خلل رياضي أساسي: فهي لا تستطيع إعادة إنتاج تناقص في الطاقة الحرة لحدود الحبيبات ($\gamma_{GB}$) مع زيادة زاوية عدم التوافق ($\Delta\theta$).

• في هذه النماذج، تكون $\gamma_{GB}$ دالة تزايدية رتيبة لـ $\Delta\theta$ بشكل صارم.
• هذا يمنع محاكاة المواد متعددة البلورات الواقعية حيث غالباً ما تُظهر $\gamma_{GB}$ نتوءات (قمماً محلية) عند زوايا "خاصة" معينة من عدم التوافق (مثل حدود التوأمة منخفضة الطاقة).
• وبناءً على ذلك، لا تستطيع هذه النماذج محاكاة الأنظمة المتناحية (مشكلة "رغوة الصابون") أو المواد ذات تضاريس طاقة حدود حبيبات غير رتيبة ومعقدة.

2. المنهجية

يقترح المؤلفون صياغة معدلة لنموذج KWC للتغلب على قيد الرتابة من خلال إدخال اللا-محلية (non-locality) في دالة الطاقة الحرة الخاصة بالنموذج.

أ. الاشتقاق النظري

• إثبات القصور: أثبت المؤلفون تحليلياً أن نماذج KWC وHMP القياسية تعاني من قصور، حيث أن مشتق طاقة حدود الحبيبة بالنسبة لعدم التوافق ($d\gamma_{GB}/d|\Delta\theta|$) تتحدد بواسطة معامل الاقتران $s$ ودالة الاقتران $g(\phi)$. وبما أن $g(\phi)$ يجب أن تزداد رتيباً لتوطين الحدود، فإن $d\gamma_{GB}/d|\Delta\theta|$ يجب أن تكون غير سالبة.
• الحل: للسماح لـ $\gamma_{GB}$ بالتناقص، يجب أن تعتمد معاملات الاقتران في دالة الطاقة الحرة على عدم التوافق اللا-محلي ($\Delta\theta$) بدلاً من مجرد التدرج المحلي ($\nabla\theta$).
• الطاقة الحرة المعدلة: تم تحديث دالة الطاقة الحرة لنموذج KWC لتصبح:
  $$F[\phi, \theta] = \int_{\Omega} \left[ V(\phi) + \frac{\alpha^2}{2}|\nabla\phi|^2 + s(\Delta\theta)g(\phi)|\nabla\theta| + \frac{\epsilon^2}{2}s(\Delta\theta)^2 h(\phi)|\nabla\theta|^2 \right] dV$$
  هنا، يصبح المعامل $s$ دالة في إجمالي عدم التوافق $\Delta\theta$ عبر الحدود.

ب. التنفيذ العددي

• حساب عدم التوافق اللا-محلي: بدلاً من تكامل $\nabla\theta$ على كامل النطاق، اقترح المؤلفون طريقة فعالة حسابياً:
  1. تحديد متجه العمودي المحلي لحدود الحبيبة ($\hat{n}$).
  2. استكمال (interpolation) مجال التوجه $\theta$ عند مسافة ثابتة $d$ في كلا الاتجاهين على طول $\hat{n}$ (أي عند $\vec{x} \pm d\hat{n}$).
  3. حساب $\Delta\theta = \theta(\vec{x} + d\hat{n}) - \theta(\vec{x} - d\hat{n})$.
• البارامترية (Parameterization): باستخدام العلاقة المستمدة بين $\gamma_{GB}(\Delta\theta)$ المدخلة والمعامل $s$، اشتق المؤلفون صيغة عكسية (المعادلة 39) لتحديد دالة $s(\Delta\theta)$ المطلوبة لأي ملف تعريف طاقة مرغوب.
• التحكم في الحركية (Mobility): لمنع "حركة المصعد" (elevator motion) غير الفيزيائية (دوران الحبيبة في الكتلة)، تم تعديل حركية التوجه $M_\theta$ بواسطة دالة $P(\epsilon|\nabla\theta|)$ تتباعد (diverge) عندما يقترب التدرج من الصفر، مما يؤدي فعلياً إلى تثبيت التوجه في الكتلة.

ج. الامتداد ثلاثي الأبعاد

تم توسيع النموذج ليشمل 3 أبعاد باستخدام الرباعيات الوحدة (unit quaternions) (متجهات رباعية الأبعاد) للتعامل مع درجات الحرية الثلاث للدوران. يتم حساب عدم التوافق باستخدام ضرب الرباعيات ($q q^{-1}$)، ويتم تكييف دالة الطاقة الحرة لتجميع تدرجات مكونات الرباعي.

3. المساهمات الرئيسية

1. الإثبات التحليلي: تقديم إثبات صارم أن نماذج مجال التوجه القياسية عاجزة رياضياً عن محاكاة تناقص طاقة حدود الحبيبات مع زيادة عدم التوافق.
2. الصياغة المبتكرة: تقديم اعتماد للمعاملات على زاوية عدم التوافق بشكل "لا-محلي"، مما يسمح بوظائف طاقة حدود حبيبات تعسفية.
3. خوارزمية اللا-محلية: تطوير خوارزمية استكمال لمسافة ثابتة لحساب عدم التوافق اللا-محلي بكفاءة، مما يتجنب التكلفة الحسابية العالية لطرق البحث التكرارية المستخدمة في المحاولات السابقة.
4. التعميم ثلاثي الأبعاد: اقتراح امتداد قوي للنموذج في ثلاثة أبعاد باستخدام الرباعيات والهارمونيات فوق الكروية (hyperspherical harmonics).

4. نتائج المحاكاة

قام المؤلفون بالتحقق من صحة النموذج باستخدام عمليات محاكاة أحادية وثنائية الأبعاد:

• ملفات طاقة تعسفية: نجح النموذج في تضمين دالة طاقة حدود حبيبات محددة (المعادلة 40) تتضمن نتوءاً حاداً عند $\Delta\theta = \pi/2$.
• الاستقرار: لم يتسبب وجود النتوء في عدم استقرار عددي. أعاد النموذج إنتاج ملف تعريف الطاقة المدخل بدقة (مع معايرة طفيفة لـ $\epsilon > 0$).
• العرض المحدود: حافظت حدود الحبيبات على عرض محدود حيث يتلاشى $\nabla\theta$ تماماً خارج الواجهة، وهو ما يتوافق مع التوقعات الفيزيائية.
• الحركية: لوحظ أن حركية حدود الحبيبات تتناقص مع زيادة عدم التوافق، وتتباعد عندما $\Delta\theta \to 0$. نجح النموذج في كبح دوران الحبيبات في الكتلة مع السماح بهجرة الحدود.
• الحد المتناحي (Isotropic Limit): أشار المؤلفون إلى أنه يمكن للنموذج نظرياً محاكاة أنظمة "رغوة الصابون" المتناحية عن طريق جعل $\gamma_{GB}$ ثابتاً، وهي مهمة كانت مستحيلة سابقاً مع نماذج مجال التوجه القياسية.

5. الأهمية

يمثل هذا العمل تقدماً كبيراً في نمذجة المجال الطوري للمواد متعددة البلورات:

• الواقعية: يسد الفجوة بين نماذج مجال التوجه الفعالة حسابياً وبين الطبيعة المعقدة وغير المتناحية لطاقات حدود الحبيبات الحقيقية.
• تعدد الاستخدامات: من خلال السماح بـ $\gamma_{GB}(\Delta\theta)$ تعسفي، يمكن للنموذج الآن محاكاة المواد ذات الحدود منخفضة الطاقة المحددة (مثل التوائم، أو حدود المواقع المشتركة الخاصة)، وهي أمور بالغة الأهمية لفهم خصائص المواد مثل القوة ومقاومة التآكل.
• المعايرة (Benchmarking): القدرة على محاكة الأنظمة المتناحية تسمح بالمعايرة المباشرة مقابل نهج أخرى للمجال الطوري (مثل نماذج المعلمات المتعددة) دون العبء الحسابي لتتبع كل حبيبة على حدة.
• التطبيقات المستقبلية: يفتح هذا الإطار الباب لمحاكة التطور المجهري المعقد في السبائك الهندسية حيث تحدد أنواع معينة من الحدود أداء المادة.

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