🔬 materials science

Phase-field simulations of nucleation, growth, and coarsening of $\beta_1$ precipitates in Mg-Nd alloys

تقدم هذه الدراسة سير عمل لمحاكاة حقل الطور المعلمي يدمج البيانات التجريبية المستمدة من مسبار الذرة الطوبوغرافي والمجهر الإلكتروني النافذ لنمذجة تنوي ونمو وتخشين ترسبات $\beta_1$ في سبائك المغنيسيوم والنيوديميوم كمياً، مما يتيح تحسين عمليات التعتيق للوصول إلى كثافات مستهدفة من الرواسب.

المؤلفون الأصليون: Lingxia Shi (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States), Stephen DeWitt (Department of Materials Science and Engineering, University of Mich

نُشر 2026-09-29

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Lingxia Shi (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States), Stephen DeWitt (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States), David Montiel (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States), Qianying Shi (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States), John Allison (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States), Katsuyo Thornton (Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, United States, Department of Nuclear Engineering and Radiological Sciences, University of Michigan, Ann Arbor, MI, United States)

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

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

بدأ الباحثون بإنشاء عينات من سبيكة المغنيسيوم وتركها لـ"تعتيق" (aging) في درجة حرارة محددة لفترات زمنية مختلفة، تتراوح من ثلاث ساعات إلى تسع ساعات. وبعد كل فترة تعتيق، استخدموا مجاهر قوية لالتقاط صور عالية الدقة لداخل المعدن. ولأن عينات المعدن قد قُطعت إلى شرائح رقيقة جداً للمعاينة، فإن الصور لم تظهر سوى مقطع عرضي ثنائي الأبعاد للبلورات. ولفهم الواقع ثلاثي الأبعاد الحقيقي، استخدم الفريق طريقة رياضية لتقدير عدد البلورات الموجودة في الحجم الكامل للمعدن بناءً على ما شاهدوه في الصور المسطحة. كما قاموا بقياس التركيب الكيميائي للمعدن المحيط بالبلورات لمعرفة كيف يتغير تركيز النيوديميوم مع تشكل البلورات.

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

أظهرت النتائج أن النموذج الذي تم ضبطه يمكنه بنجاح التقاط دورة حياة البلورات بأكملها. ففي المراحل الأولى من التعتيق، أظهرت المحاكاة اندفاعاً سريعاً لظهور بلورات جديدة، مما أدى إلى ارتفاع حاد في عدد البلورات. ومع مرور الوقت، تنبأ النموذج بشكل صحيح بأن هذا العدد سيصل إلى ذروته ثم يبدأ في الانخفاض. يحدث هذا الانخفاض لأن البلورات الأكبر تنمو على حساب البلورات الأصغر، وهي عملية تُسمى "التخشين" (coarsening)، حيث يحاول المعدن تقليل إجمالي طاقة سطحه. كما أعادت المحاكاة إنتاج أشكال واتجاهات البلورات المحددة، حيث أظهرتها كصفائح مسطحة تصطف في ثلاثة اتجاهات متميزة، تماماً كما هو ملاحظ في المعدن الحقيقي. علاوة على ذلك، تتبع النموذج بدقة انخفاض تركيز النيوديمي_ في المعدن المحيط مع تشكل البلورات، مطابقةً بذلك القياسات التي أُخذت من العينات الفعلية.

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

بيان المشكلة
تعتمد استجابة التصلد بالترسيب لسبائك المغنيسيوم-العناصر الأرضية النادرة (Mg–RE)، ولا سيما أنظمة Mg–Nd، على التوزيع المكاني ومورفولوجيا الرواسب المتكونة أثناء التعتيق. وبينما ساهمت محاكاة حقل الطور في تعزيز فهم تطور الرواسب، لا يزال النمذجة التنبؤية تشكل تحدياً بسبب عدم اليقين أو نقص المعاملات الديناميكية الحرارية والحركية الرئيسية. ويتمثل أحد الصعوبات المحددة في مقارنة نتائج المحاكاة مع البيانات التجريبية: حيث توفر التجارب عادةً بيانات مقطعية ثنائية الأبعاد (2D) (مثل المجهر الإلكتروني النافذ، TEM)، في حين أن التمثيل الفيزيائي الدقيق يتطلب غالباً نمذجة ثلاثية الأبعاد (3D). ومع ذلك، فإن محاكاة حقل الطور ثلاثية الأبعاد الكاملة لأنظمة تحتوي على رواسب عالية التباين مثل طور $\beta_1$ في سبائك Mg–Nd تعد مكلفة حاسوبياً بشكل باهظ. علاوة على ذلك، لم تنجح سوى دراسات قليلة في نمذجة تطور كثافة عدد الرواسب عبر جميع مراحل التعتيق (النوية، النمو، والخشونة) لتحسين ظروف التعتيق من أجل الحصول على الخصائص الميكانيكية المرغوبة.

المنهجية
طوّر المؤلفون سير عمل متكامل تجريبي وحسابي لتعريف معاملات نموذج حقل طور ثنائي الأبعاد (2D) للنوية العشوائية، والنمو، وخشونة رواسب $\beta_1$ في سبائك Mg–Nd المعتقة عند 523 كلفن.

التوصيف التجريبي:
- كثافة عدد الرواسب: تم الحصول على صور مجهر إلكتروني نافذ (TEM) عالية الدقة عند أوقات تعتيق مختلفة (من 3 إلى 9 ساعات). ونظراً لأن الـ TEM يوفر مقاطع عرضية ثنائية الأبعاد، فقد طُبقت طريقة قياس ستيريولوجية (تصحيح Sonderegger) لتقدير كثافات العدد ثلاثية الأبعاد ( $N_V$ ) من الأعداد ثنائية الأبعاد وسُمك الشريحة. ثم حُوّلت هذه الكثافات ثلاثية الأبعاد إلى كثافات عدد "فعالة" ثنائية الأبعاد ( $N_A$ ) عن طريق الحفاظ على المسافة المميزة بين الجسيمات ( $\lambda$ )، بافتراض أن $N_A = N_V^{2/3}$ .
- تركيب المصفوفة: استُخدم مجهر تحليل الطور الأيوني (APT) لقياس تركيز النيوديميوم (Nd) في المصفوفة عند أوقات تعتيق مختلفة، مما وفر بيانات حول استنفاد المذاب.
نمذجة حقل الطور:
- تم تنفيذ نموذج حقل طور ثنائي الأبعاد باستخدام إطار عمل PRISMS-PF. يستخدم النموذج صيغة كيم-كيم-سوزوكي (KKS) لربط التركيب ومعاملات الترتيب الهيكلي.
- تمثل ثلاثة معاملات ترتيب متغيرات التوافق الثلاثة لطور $\beta_1$ ، مما يلتقط تشوه الشبكة وطاقة الإجهاد المرن.
- تم دمج طريقة بذر النواة العشوائية، القائمة على نظرية النوية الكلاسيكية، مع النموذج لمحاكاة النوية والنمو المتزامنين. يعتمد معدل النوية على فوق الإشباع المحلي.
سير عمل تحديد المعاملات:
- نظراً لعدم اليقين في معاملات النوية (عامل زيلدوفيتش، تردد الارتباط، طاقة الواجهة)، تم استخدام تحسين متعدد الخطوات.
- الخطوة 1: تم اشتقاق كثافات العدد الفعالة ثنائية الأبعاد من البيانات التجريبية ثلاثية الأبعاد.
- الخطوة 2: تم تحسين معاملات معدل النوية ( $\rho_1, \rho_2, \tau_{in}$ ) مبدئياً باستخدام نموذج فوق إشباع موحد مبسط لمطابقة اتجاهات كثافة العدد التجريبية.
- الخطوة 3: تم تقدير الانتشار بناءً على قيم من الأدبيات وتحسينه لمطابقة معدل الخشونة.
- الخطوة 4: تم تطبيق عامل تصحيح المقياس الزمني ( $\xi$ ) على المعاملات الحركية ( $\rho_1, D, \tau_{in}$ ) لضبط الاتفاق بين محاكاة حقل الط� الكامل والبيانات التجريبية، وذلك لمعالجة القيود في فرضية فوق الإشباع الموحد.

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

الاتفاق الكمي: نجح النموذج المعرّف في إعادة إنتاج تطور كثافة عدد الرواسب، حيث التقط الزيادة السريعة خلال مرحلة النوية والانخفاض اللاحق أثناء الخشونة. وأظهرت كثافات العدد المحاكات توافقاً معقولاً مع البيانات التجريبية ثنائية الأبعاد الفعالة.
التطور المجهري: أعادت المحاكاة إنتاج المورفولوجيا الشبيهة بالصفائح لترسبات $\beta_1$ واتجاهات نموها الستة المتميزة (ثلاثة متغيرات توافق مع متغيري دوران لكل منها)، وهو ما يتوافق مع الملاحظات التجريبية. حقق النموذج ذلك دون الحاجة إلى تباين صريح في طاقة الواجهة، معتمداً بدلاً من ذلك على الخصائص المرنة متباينة الخواص.
تركيب المصفوفة: توافق متوسط تركيز النيوديميوم (Nd) في المصفوفة الذي تمت محاكاته جيداً مع قياسات APT، حيث وقع ضمن حدود عدم اليقين الإحصائي لمعظم نقاط البيانات.
المراحل الحركية: التقط النموذج الانتقال من النوية إلى النمو ثم إلى الخشونة. وأظهر كيف يؤثر استنفاد المذاب المحلي على معدلات النمو، حيث تنمو الرواسب في المناطق الغنية بالمذاب بشكل أكبر بينما تتوقف أو تذوب الرواسب في المناطق المستنفدة.

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

يذكر المؤلفون أن هذا النموذج الذي تم التحقق من صحته يمكّن من التنبؤ بالعلاقات بين العملية والبنية، مما يسمح بتحسين أوقات التعتيق للوصول إلى كثافات عدد الرواسب المستهدفة. كما يقرون بالقيود، مشيرين إلى أن عمليات المحاكاة ثنائية الأبعاد لا يمكنها التقاط التفاعلات الانتشارية ثلاثية الأبعاد بالكامل، وأن تحسين فوق الإشباع الموحد الأولي قد يؤدي إلى مجموعات معاملات غير مثالية. ومع ذلك، يخلصون إلى أن سير العمل يوفر أساساً لتحسين ظروف التعتيق في سبائك Mg–Nd ويمكن توسيعه ليشمل أنظمة Mg–RE أخرى وظروف تعتيق مختلفة. ويقترح المؤلفون أن التحسينات المستقبلية يمكن أن تتضمن بيانات تجريبية ذات دقة زمنية أعلى ودمج طرق بايز (Bayesian) لاستكشاف فضاء المعاملات.

أبحاث مشابهة