Minute-Scale High-Fidelity Gyrokinetic Simulations with Portability from Laptop to Supercomputer

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

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

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

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

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

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

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

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