Multiscale Modeling of Zr Corrosion: Predicting Microstructures of Oxide Films

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

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

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

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

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

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

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

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