On the mechanism of Pedestal Relaxation Events -- Insights gained by turbulence simulations with GRILLIX

ملخص تقني: آلية أحداث استرخاء القاعدة في تفريغات نمط I (I-mode)

بيان المشكلة
تتطلب مفاعلات الاندماج المستقبلية أنظمة تشغيل توفر حجزاً عالياً للطاقة دون تجاوز حدود المواد عند جدار الوعاء. ويعد نظام الحجز المحسن (I-mode) مرشحاً لهذه الظروف، حيث يتميز بحجز مرتفع، ومحتوى منخفض من الشوائب، وغياب أنماط (ELMs) من النوع الأول الكبيرة. ومع ذلك، يمكن لتفريغات نمط I، وخاصة تلك القريبة من انتقال I-H، أن تظهر "أحداث استرخاء القاعدة" (PREs). وبينما تتضمن هذه الأحداث قذفات طاقة دورية تشبه الـ ELMs، إلا أنها تحدث عند مستويات طاقة أقل بكثير (~1% من الطاقة المخزنة) ويُعتقد أنها مدفوعة بآلية مختلفة. وبخلاف الـ ELMs، التي تُحفز بواسطة أنماط (Peeling-Ballooning)، لا تزال الفيزياء الكامنة وراء الـ PREs أقل فهماً. وقد أشارت محاولات سابقة لمحاكاة الـ PREs باستخدام نماذج الموائع الدورانية (gyrofluid) في هندسات بسيطة إلى أن أنماط التمزق الدقيقة (MTMs) هي المرشحة لذلك، لكنها لاحظت انتقالاً إلى اضطراب شبيه بالتبادل (interchange-like turbulence) دون وجود نقل كهرومغناطيسي ملموس. علاوة على ذلك، تعد محاكاة منطقة الحافة ذات التصادم المنخفض في أجهزة الاندماج (حيث $\rho_{pol} > 0.9$) باستخدام النماذج المائعة أمراً صعباً بسبب انهيار الافتراضات التصادمية القياسية.

المنهجية
استخدم المؤلفون كود الاضطراب المائعي العالمي العابر للتصادم GRILLIX لمحاكاة تفريغ نمط I من مفاعل (ASDEX Upgrade) (الدفعة رقم 37980). استخدمت المحاكاة نموذج "براغينسكي" (Braginskii) كهرومغناطيسي لتقليل الانجراف، تم توسيعه ليشمل الأنظمة العابرة للتصادم، مع تضمين ما يلي:

• إغلاق لاندو-مائعي (Landau-fluid closure) لتدفقات الحرارة الموصلة الموازية لمعالجة انخفاض التصادمية.
• توسيع نيوكلاسيكي للزوجة الأيونية.
• نموذج غاز متعادل ثلاثي العزوم.
• حدود مصادر تكيفية للحفاظ على قيم الحدود المركزية ($n, T_e, T_i$) عند $\rho_{pol} = 0.91$.

استمرت المحاكاة لأكثر من 5 مللي ثانية من الوقت الفيزيائي. وللتحقق من النتائج، قام المؤلفون بما يلي:

1. مقارنة مخرجات المحاكاة (البروفيلات، النشاط المغناطيسي، مدة الحدث) مع البيانات التجريبية.
2. إجراء تحليل مفصل لتكافؤ الأنماط (mode parity) وعلاقات التشتت لتحديد عدم الاستقرار الدقيق.
3. تتبع تطور النظام في فضاء طول التدرج ($a/L_{Te}$ مقابل $a/L_n$) مقابل تقدير معدل النمو الخطي المستمد من تحليل هندسي صفي (slab geometry) مبسط.
4. إجراء محاكاة مقارنة باستخدام إغلاق براغينسكي القياسي مع محددات تدفق الحرارة لعزل تأثير إغلاق لاندو-المائعي.
5. إجراء تحليل استقرار خطي لكلا الإغلاقين لتحديد شروط نمو أنماط MTM.

النتائج الرئيسية

• ملاحظة الـ PREs: أعادت المحاكاة إنتاج ثلاث دفعات متقطعة تم تحديدها كـ PREs. طابقت هذه الأحداث الخصائص التجريبية، بما في ذلك النشاط المغناطيسي القوي، والمدة التي تبلغ حوالي 0.3 مللي ثانية، وتسطيحاً أساسياً لبروفيل درجة حرارة الإلكترونات بدلاً من بروفيل الكثافة.
• تحديد أنماط MTMs: أكد التحليل التفصيلي أن الـ PREs مدفوعة بأنماط التمزق الدقيقة (Micro-Tearing Modes). وتشمل الأدلة:
  • تكافؤ التمزق (Tearing Parity): أظهر جهد المتجه الموازي ($\tilde{A}_{\parallel}$) تكافؤاً زوجياً، بينما أظهر الجهد الكهروستاتيكي ($\tilde{\phi}$) تكافؤاً فردياً.
  • خصائص النقل: هيمنت على الأحداث عمليات نقل حرارة إلكترونية شعاعية كهرومغناطيسية (flutter مغناطيسي).
  • علاقة التشتت: طابقت علاقة التشتت المحسوبة تقديرات الحركة الخطية لـ MTMs، حيث أظهرت انتشاراً في اتجاه دياماجنيتيز الإلكترون مع ذروة عند $k_y\rho_s \approx 0.2$.
• آلية الدورة العالمية: رسم المؤلفون دورة PRE نوعية:
  1. يزداد تدرج درجة حرارة الإلكترونات ($a/L_{Te}$) مع مرور الوقت.
  2. بمجرد أن يتجاوز التدرج عتبة حرجة، ينمو نمط MTM (معدل النمو الخطي $\hat{\gamma} > 0$).
  3. يؤدي نمو MTM إلى إحداث نقل حراري إلكتروني شعاعي قوي، مما يؤدي إلى تسطيح تدرج درجة الحرارة بسرعة.
  4. يدخل النظام في منطقة مستقرة خطياً ($\hat{\gamma} < 0$)، مما يؤدي إلى تلاشي النمط.
  5. يزداد التدرج ببطء مرة أخرى، مما يعيد تشغيل الدورة.
• دور الإغلاق المائعي: وُجد أن إغلاق لاندو-المائعي ضروري. أدت محاكاة مقارنة باستخدام إغلاق براغينسكي (مع محددات التدفق) إلى حالة هادئة دون وجود PREs. وكشف التحليل الخطي أن إغلاق لاندو-المائعي يُدخل تأثيرات غير موضعية (عبر حد $k_{\parallel}/|k_{\parallel}|$) تؤدي إلى معدل نمو موجب لـ MTMs، في حين يتنبأ إغلاق براغينسكي باستقرار خطي.
• التباين في الطاقة المخزنة: على عكس التجارب حيث تنخفض الطاقة المخزنة بعد الـ PRE، أظهرت المحاكاة زيادة في الطاقة. ويعزو المؤلفون ذلك إلى الشروط الحدودية التكيفية: فمع زيادة النقل الشعاعي، يقوم المصدر التكيفي بضخ المزيد من طاقة التسخين للحفاظ على القيم المستهدفة في المركز، مما يعزز الطاقة المخزنة بشكل اصطناعي.
• قيود المقاومة النوعية: أشار المؤلفون إلى أنه تمت ملاحظة الـ PREs في المحاكاة في تفريغ لم تحدث فيه تجريبياً. ويفترضون أن هذا يعود إلى نموذج مقاومة سبتزر (Spitzer) المستخدم، والذي يقلل من المقاومة النوعية عند درجات الحرارة العالية. هذا التقليل في المقاومة يقلل من حد التخميد في معادلة معدل النمو، مما يؤدي إلى زعزعة استقرار MTMs بشكل اصطناعي.

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

يصرح المؤلفون بتواضع أنه بينما تتنبأ المحاكاة بحدوث الـ PREs في وقت أبكر مما لوحظ في التفريغ التجريبي المحدد (على الأرجح بسبب قيود نمذجة المقاومة)، فإن الصورة النوعية لدورة الـ PRE — المدفوعة بالتفاعل بين تدرجات درجة حرارة الإلكترونات ونمو MTM — هي صورة قوية. ويخلصون إلى أن فهم هذه الآلية، وتحديداً اعتمادها على تدرج درجة حرارة الإلكترونات، يمكن أن يساعد في الاستراتيجيات التجريبية لتخفيف آثار الـ PREs. ويقترح البحث أن العمل المستقبلي يجب أن يحسن نماذج المقاومة النوعية داخل الأكواد المائعية أو يستخدم أكواد أكثر دقة مثل (GENE-X) لمعالجة أنظمة التصادم المنخفض بدقة أكبر.