JOREK simulations of the X-point radiator formation and its movement in ASDEX Upgrade

ملخص تقني: محاكاة جورييك (JOREK) لتشكل وديناميكا مشع النقطة السينية في جهاز أزدكسس أبجريد (ASDEX Upgrade)

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

المنهجية
تستخدم هذه الدراسة كود الديناميكا المغناطيسية المائعة (MHD) الموسع وغير الخطي JOREK، المقترن بإطار عمل جسيمات حركي كامل النوع (full-f) للأنواع المتعادلة والشوائب. تم إجراء المحاكات في هندسة محور تناظر ثنائي الأبعاد (2D) تتوافق مع تفريغ رقم 38773 في جهاز ASDEX Upgrade.

• نموذج هجين: يتم نمذجة بلازما الديوتيريوم الخلفية باستخدام نموذج MHD مختزل (درجة حرارة واحدة، يتضمن انجرافات $E \times B$ وgrad-$B$ والديا مغناطيسية). بينما تُعامل جسيمات الديوتيريوم المتعادلة وشوائب النيتروجين كجسيمات حركية باستخدام تتبع مونتي كارلو، المرتبط بتبادل اتجاهي مع مجالات MHD.
• العمليات الفيزيائية: يتضمن النموذج التفاعلات الذرية مثل التأين، وتبادل الشحنة، وإعادة الاتحاد الحجمي، والإشعاع الخطّي والمستمر الفعال (المُحتسب عبر OpenADAS). ويتضمن قسم خاص من المحاكاة (القسم 7) تصادمات الشوائب النيوكلاسيكية باستخدام مؤثر هوما (Homma operator) لنمذجة الاحتكاك بين الشوائب والبلازما الخلفية.
• مراحل المحاكاة: تتبع الدراسة إعداداً ثلاثي المراحل:
  1. المرحلة الأولى: نموذج MHD الأساسي فقط لإرساء البروفيلات الأولية.
  2. المرحلة الثانية: إدخال الجسيمات المتعادلة والشوائب الحركية لدفع النظام من مرحلة الدفن المتصل للمشتتات (attached divertors) إلى الانفصال الكامل وتشكل مشع النقطة السينية (XPR).
  3. المرحلة الثالثة: ضبط معدلات الحقن والحقن لتحقيق وضع شبه مستقر لمشع النقطة السينية، متبوعاً بتغيرات ديناميكية (زيادة أو نقصان معدلات الحقن) لدراسة حركة وفقدان مشع النقطة السينية.

المساهمات الرئيسية والنتائج

1. تشكل مشع النقطة السينية والشروط المسبقة:
   نجحت عمليات المحاكاة في إعادة إنتاج الانتقال من المشتتات المتصلة إلى الحالة المنفصلة مع تشكل مشع النقطة السينية. قبل تشكل الـ XPR، يلتقط النموذج ظاهرة كثافة الجانب عالي المجال العالية (HFSHD) وانفصال أهداف المشتت الداخلي والخارجي. ويحدث تشكل الـ XPR نتيجة تكاثف الإشعاع، حيث يؤدي إشعاع الشوائب وتأين المواد المتعادلة إلى خلق حجم بارد ($T_e \approx 1.5$ إلكترون فولت) وعالي الكثافة ($n_e \approx 3 \times 10^{20}$ متر$^{-3}$) فوق النقطة السينية.

2. خصائص الاستقرار شبه المستقر لـ XPR:
   في الحالة الثابتة المرجعية، يتم الحفاظ على مشع النقطة السينية عند ارتفاع حوالي 6.8 سم فوق النقطة السينية ($\rho_{pol} = 0.995$).

• البنية: يتكون الـ XPR من قلب ($T_e \approx 1-2$ إلكترون فولت) محاط بغلاف شعاعي يمتد حتى $T_e \approx 20$ إلكترون فولت.
• تركيز الشوائب: تركيز الشوائب في القلب يتراوح بين 10-15%، وهو ما يتوافق مع الملاحظات التجريبية ومحاكاة SOLPS-ITER.
• آليات النقل: يتشكل منحدر جهد فوق الـ XPR، مما يدفع انزياح $E \times B$ قطبياً ينقل البلازما والشوائب نحو المستوى المتوسط الداخلي. وفي غياب الاحتكاك التصادمي، يؤدي ذلك إلى غلاف شعاعي مستطيل قطبياً.

3. الاستجابة الديناميكية لتغيرات معدل الحقن:

• تحول MARFE (عند زيادة الحقن): تؤدي زيادة معدل حقن النيتروجين بمقدار خمسة أضعاف إلى تحريك مشع النقطة السينية رأسياً للأعلى. ومع اشتداد ذروة الإشعاع وتداخلها مع جبهة التأين، تنخفض درجة حرارة القلب إلى أقل من 1 إلكترون فولت، وتزداد عملية إعادة الاتحاد الحجمي بشكل كبير. يؤدي هذا إلى فقدان الكثافة وتشكّل "مارفي" (MARFE) غير مستقر بالقرب من المستوى المتوسط الداخلي وفوق النقطة السينية.
• التراجع والفقد (عند خفض الحقن): يؤدي خفض معدل الحقن إلى الصفر إلى تراجع مشع النقطة السينية للأسفل. إذ يعمل انخفاض كثافة الشوائب على إضعاف تبديد القدرة الإشعاعية، مما يسخن القلب ويؤدي لحركة جبهة التأين نحو الداخل. يُفقد الـ XPR نهائياً عندما يستنزف المحتوى المتعادل في القلب، مما يتسبب في انهيار الغلاف الشعاعي داخل طبقة التخلص من النفايات (Scrape-Off Layer) عالية المجال في اتجاه المصب. وتظهر هذه العملية سلوك هستيرية (hysteresis)، حيث يتقلص الـ XPR عند كثافة شوائب أقل من تلك المطلوبة لتشكله.

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

الأهمية والنطاق
توضح الورقة قدرة كود JOREK، المعزز بإطار عمل جسيمي حركي، على محاكاة الديناميكا الزمنية لنظام مشع النقطة السينية (XPR)، بما في ذلك تشكله، وصيانته المستقرة، وتحوله إلى حالات MARFE. توفر هذه الدراسة أساساً لفهم التفاعل بين الفيزياء الذرية (التأين، إعادة الاتحاد، الإشعاع) ونقل البلازما (انجرافات $E \times B$) في سيناريوهات XPR.

يذكر المؤلفون صراحةً أن هذه المحاكاة ثنائية الأبعاد تعمل كأساس لأعمال مستقبلية، وبالتحديد:

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

العمل لا يدعي حل تحديات استنفاد الطاقة في مقاييس المفاعلات بشكل نهائي، بل يوفر أداة حسابية مُحققة لاستكشاف الآليات الكامنة لاستقرار وحركة مشع النقطة السينية، وهي أمور بالغة الأهمية لتصميم مفاعلات اندماج مستقبلية.