تخيل أنك تحاول التنقل في مدينة شوارعها عبارة عن أشرطة مغناطيسية غير مرئية وملتوية، والسيارات فيها هي موجات من الطاقة تنطلق عبر الهواء. هذا هو عالم فيزياء البلازما، دراسة الغاز فائق السخونة والمشحون كهربائيًا الذي يغذي النجوم، والذي إذا تمكنا من إتقانه، فقد يوفر طاقة نظيفة غير محدودة للأرض. وللحفاظ على هذا الغاز الناري من صهر وعائه، يستخدم العلماء مغناطيسات قوية لإنشاء "قفص" يحجز البلازما في شكل محدد. ولكن لتسخين هذا الغاز إلى درجات الحرارة اللازمة للاندماج، يقومون بإطلاق موجات راديو فيه، تمامًا كما يسخن الميكروويف الطعام. ومع ذلك، فإن هذه الموجات لا تسير في خطوط مستقيمة فحسب؛ بل تنحني وتلتوي وتمتصها البلازما بطرق معقدة اعتمادًا على شكل القفص المغناطيسي وسرعة الجسيمات بداخلها. وللتنبؤ بمكان ذهاب الطاقة بدقة، يستخدم العلماء "تتبع الأشعة" (ray tracing)، وهو طريقة حاسوبية تعمل مثل نظام تحديد المواقع (GPS) لهذه الموجات غير المرئية، حيث تحسب مسارها خطوة بخطوة لضمان إصابتها للمكان الصحيح لتسخين الوقود بكفاءة.
هنا يأتي دور BORAY-3D، وهو برنامج جديد طوره الباحثان يوشوان وانغ وهواشينغ شي، يعمل كجهاز تحديد مواقع (GPS) فائق القدرة لموجات البلازما هذه. قبل ظهور هذه الأداة، كان على العلماء استخدام خرائط مختلفة لأنواع مختلفة من الموجات أو لأشكال مختلفة من الأقفاص المغناطيسية. كانت بعض البرامج قادرة على التعامل مع الالتواءات ثلاثية الأبعاد المعقدة لجهاز "الستيلاراتور" (وهو آلة تشبه شكل الدونات بمركز ملتوي) ولكن فقط للموجات عالية التردد. بينما استطاعت برامج أخرى التعامل مع نطاق واسع من الترددات ولكن فقط للأشكال المستديرة المثالية والبسيطة التي تميز جهاز "التوكاماك". ويعد BORAY-3D الأول الذي يجمع كل شيء في حزمة واحدة؛ إذ يمكنه تتبع الموجات من الترددات المنخفضة جدًا (مثل همهمة الراديو) وصولًا إلى الترددات العالية للغاية (مثل ضوء الليزر)، ويمكنه القيام بذلك في كل من الأشكال المتماثلة البسيطة والمناظر المغناطيسية ثلاثية الأبعاد شديدة التعقيد.
اختبر الباحثون هذا الكود الجديد من خلال مقارنته ببرامج موثوقة أخرى وبيانات من العالم الحقيقي من تجارب اندماج قائمة. لقد قاموا بمحاكاة كيفية انتقال الموجات عبر "توكاماك" يحتوي على "تموج" في مجاله المغناطيسي (ناتج عن الفجوات بين المغناطيسات)، ووجدوا أن نتائجهم تطابقت تمامًا مع الأكواد المعتمدة. كما اختبروه على جهاز "وندلشتاين 7-X" (W7-X)، وهو "ستيلاراتور" ضخم في ألمانيا يتميز بشكل مغناطيسي ملتوي للغاية. وفي هذه المحاكاة، أظهر BORAY-3D أن تجاهل التأثيرات "النسبية" — حيث تتحرك الجسيمات بسرعة كبيرة تجعلها تخضع لنظرية أينشتاين للنسبية — يمكن أن يؤدي إلى أخطاء جسيمة في التنبؤ بكمية الطاقة الممتصة. ومن خلال تضمين هذه التأثيرات، استطاع الكود التنبؤ بدقة بكيفية تسخين الموجات للبلازما وكيف ستتوهج البلازما عائدةً نحو أجهزة الكشف. وسواء كان الأمر يتعلق بمحاكاة موجات بتردد 13.56 ميجاهرتز في جهاز تجريبي صغير أو موجات بتردد 220 جيجاهرتز في "ستيلاراتور" عملاق، فإن BORAY-3D يثبت أن أداة واحدة مرنة يمكنها التعامل مع كامل طيف تحديات تسخين الاندماج، مما يساعد العلماء على تصميم مفاعلات اندماج أفضل وأكثر كفاءة للمستقبل.
ملخص تقني لـ BORAY-3D: كود تتبع الأشعة لتكوينات البلازما المغناطيسية ثلاثية الأبعاد
بيان المشكلة تعد موجات التردد الراديوي (RF) ضرورية لتسخين البلازما ودفع التيار في أجهزة الحصر المغناطيسي، حيث تمتد تردداتها من الميغاهيرتز (موجات دورة أيون، هليكون، والتردد المنخفض - lower-hybrid) إلى الغيغاهيرتز (موجات دورة الإلكترون - EC). وبينما تعد أكواد تتبع الأشعة أدوات قياسية لنمذجة انتشار الموجات تحت تقريب البصريات الهندسية، فإن الأدوات الحالية تفتقر غالباً إلى القدرة الموحدة للتعامل مع كامل طيف موجات التردد الراديوي عبر مختلف هندسات البلازما. وتحديداً، فإن العديد من الأكواد إما توفر تغطية واسعة للترددات ولكنها تقتل بحدود على التوازنات متماثلة المحور (axisymmetric)، أو توفر قدرات ثلاثية الأبعاد متخصصة فقط في تطبيقات دورة الإلكترون (EC). علاوة على ذلك، فإن التنفيذات الأصلية غالباً ما تغفل تأثيرات الامتصاص النسبية الكاملة وتعتمد على أنظمة إحداثيات قائمة على دالة التدفق (flux-function)، مما يحد من قابليتها للتطبيق في مناطق خطوط المجال المفتوحة أو التكوينات ثلاثية الأبعاد المعقدة مثل الستيلاراتورات (stellarators) والتوكاماك ذات تموج المجال الحلقي.
المنهجية قام المؤلفون بتطوير BORAY-3D، وهو كود جديد لتتبع الأشعة يدمج نمذجة واسعة النطاق للترددات، ومعالجة الهندسة التعسفية ثنائية وثلاثية الأبعاد، والفيزياء النسبية في إطار عمل واحد. تتضمن المنهجية ما يلي:
نظام الإحداثيات ومعادلات الأشعة: يوسع الكود صياغة BORAY متماثلة المحور عبر حل معادلات الأشعة في الإحداثيات الأسطوانية (r,ϕ,z). وبخلاف الأكواد القائمة على إحداثيات التدفق، يقبل BORAY-3D بيانات المجال المغناطيسي والكثافة ودرجة الحرارة مباشرة في (r,ϕ,z) دون قيود دالة التدفق. وهذا يسمح بالمعالجة الموحدة لكل من المناطق ذات خطوط المجال المغلقة والمفتوحة. كما يُسمح بتغير العدد النمطي الحلقي nϕ، مما يتيح محاكاة المجالات غير المتماثلة محورياً.
نماذج التشتت والامتصاص:
الترددات المنخفضة والمتوسطة: بالنسبة لموجات دورة الأيون (ICW)، والهليكون، والتردد المنخفض (LHW)، يستخدم الكود نموذج البلازما الساخنة غير النسبي الموروث من BORAY، والمبني على موتر التشتت الماكسويلي (Maxwellian dispersion tensor).
الترددات العالية (ECW/ECE): بالنسبة لموجات دورة الإلكترون، يقوم الكود بتنفيذ نموذج امتصاص ماكسويلي نسبي كامل. فهو يحسب معامل الامتصاص باستخدام تكامل نسبي ويستخدم موتر Krivenski–Orefice لاستقطاب الموجة.
النقل الإشعاعي لـ ECE: يتم ربط وحدة نقل إشعاعي عكسي بوحدة الامتصاص النسبي لحساب أطياف انبعاث دورة الإلكترون (ECE)، باستخدام قانون كيرشوف في حد رايلي-جينز.
التنفيذ العددي: يتم دمج معادلات الأشعة باستخدام مخطط رانج-كوتا من الدرجة الرابعة. يُستخدم علاقة تشتت البلازما الباردة لتتبع الأشعة، بينما يتم حل موتر تشتت البلازما الساخنة بشكل تكراري لتحديد الامتصاص على طول مسار الشعاع.
المساهمات الرئيسية حدد البحث ثلاث مساهمات رئيسية تميز BORAY-3D عن الأدبيات الموجودة:
الصلاحية واسعة النطاق للترددات: يعمل الكود ضمن إطار نظري واحد يغطي نطاق تردد من 13.56 ميغاهيرتز (هليكون) إلى 220 غيغاهيرتز (انبعاث EC)، مع التعامل مع موجات ICW وLHW وECW مع معالجات امتصاص حركية مناسبة لكل منها.
المعالجة الهندسية الموحدة: من خلال استخدام الإحداثيات الأسطوانية وتجنب قيود دالة التدفق، يمكن للكود محاكاة تكوينات مغناطيسية تعسفية ثنائية الأبعاد (متماثلة المحور) وثلاثية الأبعاد بالكامل (بما في ذلك الستيلاراتورات والتوكاماك ذات التموج)، ودمج كل من مناطق خطوط المجال المغلقة والمفتوحة بسلاسة.
فيزياء EC النسبية: إن تنفيذ الامتصاص الماكسويلي النسبي الكامل ونقل ECE الإشعاعي العكسي يزيل القيود غير النسبية لكود BORAY الأصلي، وهو أمر بالغ الأهمية للنمذجة الدقيقة للهارمونيات العليا ضعيفة الامتصاص وأطياف ECE.
النتائج والاختبارات المعيارية تم اختبار الكود بشكل منهجي مقابل الأكواد الراسخة والنتائج التجريبية المنشورة عبر سيناريوهات مختلفة:
تموج المجال الحلقي في التوكاماك: تمت مقارنته بـ GENRAY لموجة تردد منخفض قدرها 3.7 جيجاهرتز في توكاماك ذي تموج قوي في المجال الحلقي. اتفقت مسارات الأشعة ضمن 7.1 ملم، وتطابقت اتجاهات امتصاص القدرة، مع عزو الفروقات إلى اختلاف نماذج البلازما الساخنة الأساسية.
موجات ECW في الستيلاراتور (W7-X): تمت مقارنته بـ Raytrax وTRAVIS لأشعة O-mode وX-mode بتردد 140 جيجاهرتز. اتفقت المسارات ضمن 5 ملم. وسلطت الدراسة الضوء على أن النماذج غير النسبية تقلل بشكل كبير من تقدير الامتصاص وتنحرف عن المنحنيات النسبية، مما يؤكد ضرورة المعالجة النسبية.
موجات ECW في الستيلاراتور (HSX): تم التحقق من صحته مقابل النتائج المنشورة لتسخين X2 عند 28 جيجاهرتز. استعاد الكود اعتماد الكثافة لامتصاص المرور الأول والترسيب المركزي، رغم اختلاف القيم المطلقة قليلاً بسبب عدم توفر بيانات المجال القابلة للقراءة آلياً في المرجع.
موجات ICW الشبيهة بـ LHD: تمت محاكاة موجة سريعة بقدرة 50 ميغاهيرتز في تكوين LHD مُقاس. استعاد الكود الاتجاهات الفيزيائية الواردة في الأدبيات، بما في ذلك الارتفاع في العدد الموجي الموازي المعتمد على الكثافة والتحول من هيمنة امتصاص التريتيوم إلى امتصاص الإلكترونات.
هليكون MUSE: أظهر تطبيقاً لموجة هليكون بقدرة 13.56 ميغاهرتز في ستيلاراتور ذي مغناطيس دائم، مبيناً تأثيرات مسار ثلاثية الأبعاد متميزة بناءً على مستويات تناظر الإطلاق.
ECE في W7-X: حسب أطياف انبعاث X-mode (115–220 جيجاهرتز) للجانب منخفض المجال والجانب عالي المجال. تنبأ النموذج النسبي بشكل صحيح بالعمق البصري لـ X3، في حين قللت النماذج غير النسبية من تقديرها بمقدار عشرة أضعاف.
الأهمية والادعاءات يزعم المؤلفون أن BORAY-3D يوفر أداة متعددة الاستخدامات لدراسة انتشار موجات التردد الراديوي في البلازما المغناطيسية من خلال توحيد القدرات التي كانت مشتتة سابقاً عبر أكواد مختلفة. تكمن أهميته في قدرته على التعامل مع كامل طيف الترددات الهندسية المعقدة ثلاثية الأبعاد (بما في ذلك خطوط المجال المفتوحة وتأثيرات التموج) دون الحاجة إلى تحويلات إحداثيات التدفق. إن تضمين امتصاص EC النسبي الكامل ونمذجة ECE يعالج قيوداً محددة في الإصدارات السابقة، مما يسمح بتنبؤات أكثر دقة للتسخين عالي التردد والتشخيصات في أجهزة الستيلاراتور والتوكاماك الحديثة. يشير البحث إلى أنه بينما يعد تقريب البصريات الهندسية قوياً، فإن انعكاس الموجة عند حواف البلازما وتأثيرات الحيود لا تزال تتطلب نمذجة الموجة الكاملة (full-wave modeling) للتحقق الكمي في تلك الأنظمة المحددة. ومن المخطط مستقبلاً معالجة التخميد التصادمي، وانعكاس الجدران، ودفع التيار، والتوزيعات غير الماكسويلية.