Phase field as a front propagation method for modeling grain growth in additive manufacturing

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

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

ينتج التصنيع المضاف (AM)، وتحديداً صهر سرير المسحوق (PBF-AM)، مكونات معدنية معقدة ذات بنيات مجهرية حساسة للغاية لمعايير العملية. يؤدي التصلب السريع المدفوع بالتدرجات الحرارية العالية وسرعات الواجهة المرتفعة إلى ظهور سمات غير متوازنة، مثل الحبيبات العمودية فوق المحورية (epitaxial columnar grains) ذات النسيج الليفي القوي $\langle 001 \rangle$، مما يسبب تباينًا خواصيًا ميكانيكيًا.

تواجه نماذج المحاكاة الحالية مقايضة بين التكلفة الحسابية والدقة الفيزيائية:

• الخلايا الآلية (Cellular Automata - CA): فعالة ولكنها تعتمد على قواعد تجريبية، وتفشل في التقاط تأثيرات الانحناء، أو إعادة توزيع المذاب، أو حركة الواجهة المتسقة ديناميكيًا حراريًا، خاصة في الظروف العابرة (مثل انصهار الثقب المفتاحي/keyhole melting).
• نموذج المجال الطوري الكامل (Full Phase-Field - PF): دقيق ديناميكيًا حراريًا وقادر على تمثيل مورفولوجيا التشعبات (dendritic morphology) وانتشار المذاب، ولكنه مكلف حسابيًا للغاية بالنسبة لأحجام حوض الانصهار بمقياس المليمتر المعتادة في التصنيع المضاف.
• النهج الهرمي (Hierarchical Approaches): يستخدم نماذج حرارية عيانية (macroscopic) لتوجيه محاكاة المقياس المتوسط، ولكنه غالبًا ما يتطلب افتراضات حول حركية الواجهة ولا يمكنه حل التفاصيل على مقياس التشعبات.

هناك حاجة ماسة لـ نموذج ميزوسكوبي (متوسط المقياس) فعال حسابيًا وتنبؤي يمكنه محاكاة نمو الحبيبات وتطور النسيج عبر طبقات وممرات متعددة في التصنيع المضاف دون الحاجة إلى حل كل فرع تشعبي فردي.

2. المنهجية

قام المؤلفون بتطوير نموذج غلاف الحبيبة ميزوسكوبي مقترن بـ طريقة انتشار جبهة المجال الطوري (PF). يربط هذا النهج بين المقياس المجهري (أطراف التشعبات) والمقياس المتوسط (أغلفة الحبيبات).

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

• تمثيل غلاف الحبيبة: بدلاً من حل الأذرع التشعبية الفردية، يتم تمثيل الهيكل التشعبي متعدد الفروع للحبيبة بالكامل بواسطة سطح "غلاف" سلس يربط أطراف التشعب النشطة. يُفترض أن الجزء الداخلي صلب تمامًا أو في حالة عجين (mushy).
• انتشار جبهة المجال الطوري: يتم التعامل مع سطح الغلاف كواجهة منتشرة ($\phi$) وليس كحد حاد. يخضع تطور الغلاف لمعادلة معدل قياسية مشتقة من المتطابقة:
  $$\dot{\phi} \approx \frac{\pi}{\eta} \sqrt{\phi(1-\phi)} v_{env}$$
  حيث $\eta$ هو عرض الواجهة و $v_{env}$ هو سرعة الغلاف.
• القانون الحركي: لا يتم تحديد سرعة الغلاف عبر حدود الشعيرية القياسية (التي يتم تقليلها لتجنب التأثيرات الزائفة) بل من خلال قانون حركي قائم على القابلية المجهرية (microscopic-solvability). يتم حساب سرعة طرف التشعب ($v_{tip}$) بناءً على إجمالي التبريد الفائق ($\Delta T_{total}$):
  $$v_{tip} = a (\Delta T_{total})^n$$
  ثم تُشتق سرعة الغلاف من سرعة الطرف، والتوجه البلوري، والناظم المحلي للواجهة.
• النموذج الحراري: يتم حل معادلة التوصيل الحراري العابرة، مع دمج:
  • مصدر حراري غاوسي متحرك (ليزر/شعاع إلكتروني).
  • تحرر الحرارة الكامنة داخل منطقة العجين (يتم التعامل معه عبر سعة حرارية فعالة $C_p^*$).
  • شروط حدودية تحاكي تدرجات درجات الحرارة الثابتة.

إعداد المحاكاة

• البرمجيات: تم التنفيذ باستخدام مكتبة Open-Phase مفتوحة المصدر.
• النطاق: ثنائي الأبعاد ($3 \times 8.5$ مم) وثلاثي الأبعاد ($1.5 \times 1.2 \times 5$ مم).
• المادة: سبيكة ثنائية من النيكل والألومنيوم (تمثل السبائك الفائقة).
• التهيئة: تم تهيئة الركائز بحبيبات فورونوي (Voronoi) عشوائية؛ ويحدث النوى (nucleation) على عمق ~30–40 $\mu$m تحت السطح.

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

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

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

• محاكاة المسار الواحد ثنائية وثلاثية الأبعاد:
  • أعاد النموذج بنجاح إنتاج النمو فوق المحوري (epitaxial growth) حيث تتقدم الحبيبات من الركيزة إلى السائل المبرد فائقًا.
  • لوحظ النمو التنافسي: الحبيبات ذات التوجه البلوري المفضل تجاه التدرج الحراري ($\langle 001 \rangle$ في المعادن المكعبة) تفوقت على جيرانها غير الموجهة، مما أدى إلى نسيج عمودي قوي.
  • كشفت المحاكاة ثلاثية الأبعاد عن تعقيد إضافي بسبب المنافسة خارج المستوى وتأثيرات الانحناء، والتي لا تستطيع النماذج ثنائية الأبعاد التقاطها.
• الدراسة البارامترية (درجة حرارة الركيزة والحركية):
  • درجة حرارة الركيزة المنخفضة: أدت إلى زيادة التدرج الحراري ($G$)، مما نتج عنه منطقة تبريد فائق أضيق ولكنها أكثر كثافة أمام جبهة التصلب.
  • الاستقرار المطلق: أدت التدرجات الحرارية العالية والسرعات العالية للواجهة إلى كبح رفض المذاب، مما دفع الواجهة نحو نظام "الاستقرار المطلق" حيث يتم تخميد الاضطرابات المورفولوجية.
  • تأثير الـ CET: يتحكم حجم وشدة منطقة التبريد الفائق مباشرة في احتمالية تكوين نوى جديدة (CET). يحدد النموذج كيفية ضبط معايير العملية للتحكم في هذا التحول.
• البناء متعدد الطبقات:
  • أظهرت محاكاة ثلاث طبقات متتالية الاختيار التدريجي للنسيج.
  • تم تجاوز الحبيبات متساوية المحاور العشوائية الأولية تدريجيًا بواسطة الحبيبات العمودية التي تمتد عبر الطبقات.
  • بحلول الطبقة الثالثة، نجت مجموعة محدودة فقط من التوجهات المهيمنة، مما أدى إلى ظهور حبيبات عمودية طويلة ومستمرة تمتد عبر طبقات متعددة، وهو ما يتوافق مع الملاحظات التجريبية في PBF-AM.

5. الأهمية

يمثل هذا العمل بديلًا فعالًا حسابيًا لنماذج المجال الطوري التشعبية الكاملة، مما يجعل محاكاة نمو الحبيبات عبر أطوال كبيرة (مليمترات) وجداول زمنية (بناء متعدد الطبقات) ذات صلة بعمليات التصنيع المضاف الصناعية أمرًا ممكنًا.

• القدرة التنبؤية: يمكن للنموذج التنبؤ بتطور البنية المجهرية والنسيج بناءً على معايير العملية (قدرة الشعاع، سرعة المسح، التسخين المسبق)، مما يساعد في تحسين عمليات التصنيع المضاف لتقليل التباين أو تعزيز الهياكل متساوية المحاور.
• تصميم العمليات: من خلال قياس العلاقة بين التدرجات الحرارية، والمعاملات الحركية، ومنطقة التبريد الفائق، يوفر النموذج أداة لتصميم استراتيجيات التحكم في التحول من العمودي إلى متساوي المحاور (CET)، وهو أمر بالغ الأهمية لتحسين الخواص الميكانيكية لمكونات التصنيع المضاف.
• القابلية للتوسع: هذا النهج مناسب للدمج في أطر أكبر لتصميم وتحسين العمليات، مما يسد الفجوة بين فيزياء التصلب الأساسية والتطبيقات الصناعية.