A model of grain growth in UN integrating molecular dynamics, phase-field modeling, and uncertainty quantification
تؤسس هذه الدراسة أول إطار كمي لنمو الحبيبات لنيتريد اليورانيوم (UN) من خلال دمج الديناميكا الجزيئية، ونمذجة حقل الطور، وتقدير عدم اليقين لتحديد أن حركية حدود الحبيبات الذاتية هي العامل المهيمن الذي يحكم حركية النمو، في حين أن سحب المسام مهمل وتأثير طاقة حدود الحبيبات ضئيل.
المؤلفون الأصليون:Mohamed AbdulHameed, Fadel M. Nasr, Wen Jiang, Mahmoud Yaseen, Benjamin Beeler
في قلب المفاعل النووي، لا يظل الوقود ساكناً. فحتى أثناء احتراقه، تظل البلورات الصغيرة التي تشكل حبيبات الوقود في حالة تحرك مستمر، وتندمج، وتنمو لتصبح أكبر حجماً. هذه العملية، المعروفة باسم "نمو الحبيبات"، هي رقصة أساسية للذرات تحدد مدى طول عمر الوقود ومدى أمان تشغيله. لعقود من الزمن، فهم العلماء هذا السلوك جيداً في أنواع الوقود النووي الشائعة مثل ثاني أكسيد اليورانيوم، لكن بديلاً واعداً يسمى "مونونيتريد اليورانيوم" ظل لغزاً. هذه المادة، وهي سيراميك صلب مكون من اليورانيوم والنيتروجين، توفر موصلية حرارية وكثافة فائقتين، مما يجعلها مرشحاً رئيسياً للمفاعلات من الجيل القادم. ومع ذلك، بدون معرفة كيفية نمو وتفاعل بلوراتها الداخلية بدقة، لا يستطيع المهندسون الوثوق بها تماماً لتشغيل مفاعل لسنوات. كانت القطعة المفقودة في أحجية اللغز هي رسم خريطة واضحة لكيفية تحرك الحدود بين هذه البلورات وكيفية دفع الطاقة الداخلية للمادة لهذا التحرك.
قام فريق من الباحثين الآن ببناء تلك الخريطة، حيث أنشأوا نموذجاً رقمياً شاملاً يربط بين حركة الذرات الفردية وسلوك حبة الوقود بأكملها. لم يعتمدوا على طريقة واحدة، بل نسجوا معاً ثلاث مقاربات متميزة لمحاكاة حياة المادة. أولاً، استخدموا محاكاة حاسوبية قوية لمراقبة كيفية ترتيب الذرات لنفسها عند الحدود بين البلورات في درجات حرارة تتراوح من الصفر المطلق حتى 2000 درجة مئوية. سمح لهم ذلك بحساب الطاقة المخزنة في هذه الحدود، وهي المحرك الرئيسي للنمو. بعد ذلك، لجأوا إلى تجربة واقعية أجريت منذ سنوات على وقود مشابه، باستخدام نهج رياضي ذكي لإزالة تأثير المسام الصغيرة واستخراج السرعة الحقيقية لتحرك حدود البلورات. وأخيراً، دمجوا هذه النتائج في محاكاة واسعة النطاق راقبت آلاف الحبيبات الافتراضية وهي تنمو بمرور الوقت، مع اختبار مدى حساسية النتائج للتغيرات الصغيرة في بيانات المدخلات.
كشفت النتائج عن قصة بسيطة بشكل مفاجئ حول ما يحرك هذا النمو. وجد الباحثون أن الطاقة المخزنة في الحدود بين البلورات تظل ثابتة تقريباً عند درجات الحرارة المنخفضة، ولكنها تبدأ في الارتفاع مع زيادة الحرارة. والأهم من ذلك، اكتشفوا أن جيوب الفراغ الصغيرة، أو المسام، التي غالباً ما تعيق حركة الحدود في مواد أخرى، هي في الواقع غير ضارة في هذا الوقود تحديداً تحت الظروف التي درسوها. ومن خلال تحليل البيانات التجريبية الوحيدة المتاحة لهذه المادة، أثبتوا أن هذه المسام لا تبطئ عملية النمو بشكل كبير. سمح لهم ذلك بحساب السرعة الجوهرية لحدود البلورات، وهو رقم لم يكن معروفاً بدقة من قبل. لقد حددوا أن الحدود تتحرك بمعامل سرعة معين وحاجز طاقة يجب التغلب عليه لحدوث الحركة، مما يوفر أول قواعد كمية لكيفية تطور مونونيتريد اليورانيوم بمرور الوقت.
ولضمان متانة هذه القواعد، أجرى الفريق آلاف عمليات المحاكاة، مع تغيير الأرقام التي حسبوها قليلاً لمعرفة مدى تغير الحجم النهائي للحبيبات. ووجدوا أن السرعة التي تتحرك بها الحدود هي العامل الأكثر أهمية؛ فإذا كانت السرعة غير دقيقة قليلاً، فإن الحجم المتوقع للبلورات يتغير بشكل كبير. أما الطاقة المطلوبة لتحريك الحدود فهي العامل الثاني في الأهمية، بينما كان للطاقة المخزنة في الحدود نفسها تأثير طفيف جداً. يخبر هذا التسلسل الهرمي العلماء المستقبليين بالضبط أين يجب أن يركزوا جهودهم: للتنبؤ بحياة هذا الوقود بثقة عالية، يجب عليهم قياس سرعة حدود البلورات بدقة أكبر. وتؤكد الدراسة أن المادة تتبع نمط نمو يمكن التنبؤ به، حيث يستقر توزيع أحجام الحبيبات في النهاية في شكل مستقر ومتشابه ذاتياً، بغض النظر عن كيفية بدء الوقود.
يوفر هذا العمل أول إطار كامل لفهم كيفية تقادم مونونيتريد اليورانيوم تحت تأثير الحرارة. ومن خلال الجمع بين الحسابات على المستوى الذري والمحاكاة واسعة النطاق والاختبار الصارم لعدم اليقين، حول الباحثون فهماً غامضاً للمادة إلى أداة تنبؤية دقيقة. لقد أظهروا أن المادة تسلك سلوكاً مباشراً ويمكن التنبؤ به، مدفوعاً بحركة حدودها الداخلية بدلاً من أن تعيقها العيوب الداخلية. وبينما أُجريت الدراسة على نوع محدد من الوقود وتحت ظروف محددة، فإن الأساليب التي طوروها يمكن تطبيقها على مواد نووية متقدمة أخرى. إن المسار نحو الأمام واضح الآن: مع كتابة قواعد اللعبة أخيراً، فإن الخطوة التالية هي تحسين قياسات الأرقام الأكثر حرجاً، لضمان أنه عندما يُستخدم هذا الوقود في المستقبل، يمكن الوثوق بأدائه بيقين مطلق.
ملخص تقني: إطار عمل متعدد المقاييس لنمو الحبيبات في نيتريد اليورانيوم
بيان المشكلة يُعد نيتريد اليورانيوم (UN) وقوداً نووياً متقدماً واعداً نظراً لكثافته العالية من المواد الانشطارية وموصليته الحرارية العالية. ومع ذلك، فإن اعتماده على نطاق واسع يعوقه نقص الفهم الأساسي لتطور بنيته المجهرية، وتحديداً حركية نمو الحبيبات. وبينما تمت دراسة نمو الحبيبات في أنواع أخرى من وقود الأكتينيدات مثل ثاني أكسيد اليورانيوم (UO₂) وسيليسايد اليورانيوم (U₃Si₂)، فإن الحركية والآليات في نيتريد اليورانيوم (UN) لا تزال غير مستكشفة. وتوجد فجوة حرجة في المعاملات الكمية المطلوبة للمحاكاة على المقياس المتوسط: وتحديداً طاقة حدود الحبيبات (GB energy) كدالة لدرجة الحرارة والتوجه الخاطئ (misorientation)، وحركية حدود الحبيبات الجوهرية الخالية من الآثار التجريبية مثل سحب المسام (pore drag). علاوة على ذلك، لم يتم قياس مدى حساسية تنبؤات نمو الحبيبات لعدم اليقين في هذه المعاملات بشكل منهجي.
المنهجية طوّر المؤلفون إطار عمل متكامل متعدد المقاييس يجمع بين الديناميكا الجزيئية (MD)، ونمذجة المجال الطوري (phase-field modeling)، والتقدير المساعد للنماذج البديلة لتقدير عدم اليقين (surrogate-assisted UQ).
الديناميكا الجزيئية (MD): لتحديد طاقات حدود الحبيبات، استخدمت الدراسة جهود "تسييبليايف وكوتشيفسكي" (Tseplyaev and Kocevski) في برنامج LAMMPS. قام الباحثون ببناء 27 نموذج بلورة ثنائية ذات ميل متماثل مع زوايا توجه خاطئ تتراوح بين 0° و90° حول المحور [001]. ولكل حد، تم إجراء بحث دقيق عن هيكل الطاقة الدنيا باستخدام عمليات الترجمة للأجسام الصلبة وتكرارات حذف الذرات (طريقة سطح γ). أُجريت محاكاة الديناميكا الجزيئية من 0 كلفن إلى 2000 كلفن لتقييم الاعتماد على درجة الحرارة.
استخراج الحركية: عالجت الدراسة ندرة البيانات الخاصة بنيتريد اليورانيوم (UN) من خلال تحليل مجموعة البيانات الوحيدة المتاحة لنيتريدات الأكتينيدات (Ronchi and Sari, 1983)، والتي تضمنت حبيبات (U,Pu)N. تم تطبيق نموذج ميكانيكي لسحب المسام، يتضمن اقتراناً جزئياً بين المسام وحدود الحبيبات، لاستخراج حركية حدود الحبيبات الجوهرية (الخالية من السحب) من حركية النمو المكعبة التجريبية (n=3). تضمن ذلك إجراء انحدار للحركية الفعالة مقابل نموذج يأخذ في الاعتبار قوى سحب المسام لتحديد ما إذا كانت آثار السحب جوهرية.
نمذجة المجال الطوري: أُجريت محاكاة نمو الحبيبات عند درجة حرارة ثابتة باستخدام إطار عمل MOOSE. استخدم النموذج معادلة "ألن-كان" (Allen-Cahn) مع معاملات ترتيب غير محفوظة لتمثيل الحبيبات. أُجريت المحاكاة في نطاقات ثنائية الأبعاد (2500 نانومتر) عبر درجات حرارة تتراوح بين 1500-2000 كلفن، باستخدام طاقات حدود الحبيبات المستمدة من الديناميكا الجزيئية ومعاملات الحركية المستخرجة.
تقدير عدم اليقين وتحليل الحساسية: استُخدم نهج مدعوم بالنماذج البديلة لتحليل تأثير عدم اليقين في المدخلات (المعامل المسبق للحركية M0، وطاقة التنشيط QM، وطاقة حدود الحبيبات γ) على المخرج (متوسط حجم الحبيبات). تضمن ذلك:
توليد 330 عينة للمعاملات عبر أخذ عينات "لاتين هايبركيوبيك" (Latin Hypercube Sampling).
تقليل أبعاد مخرجات حجم الحبيبات المعتمدة على الزمن باستخدام تحليل المكونات الرئيسية (PCA).
تدريب نماذج انحدار "العملية الغاوسية" (GP) على المكونات الرئيسية.
حساب مؤشرات "سوبول" (Sobol) للتحليل الحساس عبر الزمن لتحديد مساهمة كل معامل في تباين المخرجات.
النتائج الرئيسية
طاقة حدود الحبيبات: أظهر جهد "تسييبليايف" قيم طاقة حدود حبيبات متوافقة مع قيم نظرية الكثافة الوظيفية (DFT)، بينما قلل جهد "كوتشيفسكي" الطاقات بشكل منهجي بسبب التنافر غير الطبيعي بين ذرات اليورانيوم (U-U) الذي يمنع انغلاق حدود الحبيبات. طاقة حدود الحبيبات المتوسطة شبه مستقلة عن درجة الحرارة تحت 1000 كلفن ولكنها تزدادة عند درجات الحرارة الأعلى. تم وضع ارتباط متعدد الحدود: γ(T)=1.382+3.279×10−6T+5.748×10−18T5 (جول/م²)، مع انحراف معياري يبلغ حوالي 30% عبر الـ 27 حداً.
حركية حدود الحبيبات: كشف تحليل سحب المسام الميكانيكي أن سحب المسام ضئيل في ظل الظروف التجريبية لمجموعة البيانات المرجعية، حيث كانت عوامل تقليل الحركية s≈0.93–0.99. وبناءً على ذلك، اعتُمدت الحركية الفعالة المستمدة من البيانات التجريبية كحركية جوهرية. وكانت معاملات "أرينيوس" الناتجة هي M0=2.05×10−15 م⁴/(جول·ثانية) و QM=0.89 إلكترون فولت.
محاكاة المجال الطوري: أكدت المحاكاة نمو الحبيبات الطبيعي المدفوع بالانحناء، والذي يتميز بتطور خطي لـ (D2−D02) مع الزمن. تقاربت توزيعات حجم الحبيبات نحو شكل يشبه "هيليرت" (Hillert-like)، بشكل مستقل عن البنية المجهرية الأولية. تشبع معامل الاختلاف (CV) لأحجام الحبيبات عبر الظروف الأولية المختلفة بشكل أسي، حيث أدت درجات الحرارة الأعلى إلى تسريع الوصول إلى التشابه الذاتي ولكن مع إنتاج توزيعات حالة مستقرة أكثر اتساعاً.
تحليل الحساسية: أظهر تحليل الحساسية العالمي تسلسلاً هرمياً واضحاً للتأثير. يهيمن المعامل المسبق للحركية M0 على تباين المخرجات في جميع الأوقات (بمساهمة تبلغ 48% من تباين الدرجة الأولى عند t=5000 ثانية)، يليه طاقة التنشيط QM (28%). أما طاقة حدود الحبيبات γ فتساهم بشكل ضئيل (~7%). كما كشف التحليل أنه بينما تكون تفاعلات المعاملات كبيرة في الأوقات المبكرة، إلا أنها تتضاءل مع اقتراب النظام من نظام النمو المكافئ التقاربي.
الأهمية والادعاءات يدعي البحث أنه وضع أول إطار عمل كمي لنمو الحبيبات في نيتريد اليورانيوم (UN). وتتمثل مساهماته الرئيسية في:
توفير طاقات حدود الحبيبات المعتمدة على درجة الحرارة لـ 27 حداً من التماثل المتماثل، والتحقق منها مقابل DFT.
تطوير إجراء صارم لاستخراج حركية حدود الحبيبات الخالية من السحب من بيانات تجريبية محدودة من خلال نمذجة سحب المسام صراحةً واستبعاد تأثيراتها.
تقديم أول معاملات لحركية حدود الحبيبات الجوهرية (M0 و QM) لنيتريد اليورانيوم.
تحديد، من خلال تحليل الحساسية عبر الزمن، أن تقليل عدم اليقين في المعامل المسبق للحركية M0 وطاقة التنشيط QM هو الهدف الأعلى أولوية للجهود التجريبية المستقبلية، حيث تهيمن هذه المعاملات على عدم اليقين التنبئي.
يشير المؤلفون إلى أنه بينما تُعامل معاملات الحركية الحالية كتقديرات لترتيب المقدار نظراً لاستخدام بيانات (U,Pu)N والتقريبات في انتشار السطح، فإن الإطار يوفر منهجية قابلة للنقل، قائمة على الفيزياء، وواعية بعدم اليقين، يمكن تطبيقها في أكواد أداء الوقود وتوسيعها لتشمل أنظمة وقود نووي متقدمة أخرى. وتذكر الدراسة صراحة أن الطبيعة ثنائية الأبعاد للمحاكاة تؤدي إلى معدلات نمو تبلغ نصف المعدلات المتوقعة في كتل ثلاثية الأبعاد، ولكن السلوكيات النوعية وترتيبات الحساسية تظل صالحة.