Modeling of plasma transport during edge-localized mode in tokamak using a kinetic Vlasov-Poisson code

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

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

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

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

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