Pre L-H Transition Radial Electric Field and Transport Validations of Edge and Scrape-off Layer Gyrokinetic Simulations at ASDEX Upgrade

ملخص تقني: التحقق من المجال الكهربائي الشعاعي والنقل في مرحلة ما قبل انتقال L-H عبر محاكاة الجايروكينيتية للحافة وطبقة التخلص من الحطام (SOL) في جهاز ASDEX Upgrade

بيان المشكلة
يعد التنبؤ الدقيق بالاضطراب والنقل في حافة البلازما وطبقة التخلص من الحطام (SOL) تحديًا حاسمًا للمفاعلات الاندماجية المستقبلية، حيث تحدد هذه المناطق مستويات احتجاز الطاقة والقدرة الاندماجية. وبينما يتم التحكم في الانتقال من نمط الاحتجاز المنخفض (L-mode) إلى نمط الاحتجاز العالي (H-mode) من خلال عمليات اضطراب وحمل الحافة، فإن الاقتران القوي بين حافة البلازما المحتجزة وطبقة الـ SOL، بالإضافة إلى وجود نقاط X المغناطيسية، يتحدى قابلية تطبيق نماذج النقل المحلية. لقد نجحت دراسات التحقق السابقة باستخدام محاكاة الجايروكينيت (GK) المحلية من نوع $\delta f$ (مثل استخدام كود GENE) في إعادة إنتاج مستويات النقل الاضطرابي، لكنها اعتمدت على ملفات تعريف التوازن (equilibrium profiles) والمجالات الكهربائية الشعاعية ($E_r$) المفروضة تجريبيًا كمدخلات. وبناءً على ذلك، لا تستطيع هذه النهج المحلية تقديم وصف تنبئي لديناميكيات الحافة أو تطور بئر الـ $E_r$، وهو أمر مركزي لتنظيم اضطراب الحافة وعتبة قدرة انتقال L-H.

المنهجية
يعتمد هذا العمل على استراتيجية تحقق تدريجية باستخدام محاكاة جايروكينيتية (GK) أولية، كهرومغناطيسية، كاملة الـ $f$، وتصادمية لحافة وSOL، بما في ذلك هندسة نقطة X. تستخدم المحاكاة كود GENE-X لنمذجة تفريغ هيدروجيني مخصص (#38176) في جهاز التوكاماك ASDEX Upgrade (AUG).

• استراتيجية المحاكاة: بدلاً من محاكاة التطور الزمني البطيء لملفات تعريف البلازما تحت قدرة مفروضة، يقوم المؤلفون بإجراء محاكاة مستقلة لأربع شرائح زمنية تجريبية متتالية لنمط L ($t = 2.7, 3.8, 4.5, 4.8$ ثانية) تقترب من انتقال L-H. يتم تثبيت حدود الحافة الداخلية على الملفات التجريبية، بينما تتطور ملفات تعريف الحافة والـ SOL ذاتيًا استجابةً للاضطراب، والنقل، ومصادر الكثافة.
• نمذجة مصدر الكثافة: لمعالجة غياب مصادر كثافة الحافة في النماذج القياسية، تم إدخال مصدر كثافة حركي (kinetic density source) مُعدل ($S_n$) بالقرب من الفاصل (separatrix) لمحاكاة تأين الغاز المتعادل. تم ضبط هذا المصدر ليتوافق مع ملفات تعريف الكثافة التجريبية بالقرب من الفاصل عند الشريحة الزمنية الأخيرة ($t=4.8$ ثانية).
• الشروط الحدودية: تُطبق شروط ديريكليه (Dirichlet) الحدودية عند الحافة الداخلية وحدود الـ SOL. تُجرى المحمات حتى الوصول إلى حالة شبه مستقرة (QSS)، حيث تتوازن تدفقات الجسيمات والطاقة مع المصادر والمصبات.
• التحليل: تتحقق الدراسة من ملفات تعريف المستوى المتوسط الخارجي (OMP) (الكثافة، درجة الحرارة، $E_r$) مقابل القياسات التجريبية. كما تقوم بتفكيك المجال الكهربائي الشعاعي عبر تحليل توازن القوى، وتوصيف اضطراب الحافة عبر علاقات التشتت، وتحليل تدفقات نقل الجسيمات والحرارة.

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

1. التحقق من ملفات تعريف $E_r$: تظهر المحاكاة توافقًا نوعيًا وكميًا ممتازًا مع ملفات تعريف $E_r$ التجريبية، لا سيما عمق وعرض بئر $E_r$ السالب بالقرب من الفاصل. وقد تبين أن إدراج مصدر الكثة ($S_n$) عند $t=4.8$ ثانية ضروري لإعادة إنتاج تعمق بئر $E_r$ الملاحظ تجريبيًا، مما جعل القيم الدنيا المحاكات (حوالي $-15$ كيلو فولت/متر) قريبة من قياسات نمط H.
2. توازن القوى والتدفقات الاستقطابية: يكشف تفكيك توازن القوة أنه بينما تهيمن تدرجات ضغط الأيونات والدوران الطردي على منطقة الحافة ($\rho_{pol} \lesssim 0.98$)، فإن هيكل بئر $E_r$ محكوم بـ التدفقات الاستقطابية المدفوعة بالاضطراب. هذه التدفقات تنحرف بشكل كبير عن تنبؤات النيوكلاسيكية (NC). ويرتبط تعمق بئر $E_r$ مباشرة بتقوية هذه التدفقات الاستقطابية غير الخطية والمدفوعة بالاضطراب.
3. دور مصدر الكثافة: يعد إدخال مصدر كثافة الحافة أمرًا حاسمًا لإعادة إنتاج ملفات تعريف الكثافة ذات الصلة تجريبيًا وتقليل التدفقات الحرارية للإلكترونات والأيونات المبالغ فيها. يعمل المصدر كمصب حراري موضعي، مما يؤدي إلى تسطيح الملفات واستقرار الاضطراب، وهو ما يؤدي بدوره إلى تقليل النقل الاضطرابي.
4. توصيف الاضطراب: تم تحديد اضطرب الحافة كصراع بين موجات درفت الإلكترونات (eDWs) وأنماط الإلكترونات المحاصرة (TEMs)، حيث تسود موجات eDWs. أدى إدراج مصدر الكثافة إلى تضييق طيف الاضطراب ونقل ذروته إلى أعداد موجية أقل، وهو ما يتوافق مع هيمنة الـ eDWs التصادمية.
5. تنبؤات النقل:
   • بدون مصدر الكثافة: تبالغ المحاكاة بشكل منهجي في تقدير التدفقات الحرارية، وخاصة تدفق الحرارة الإلكترونية ($P_e$)، وتفشل في إعادة إنتاج ملفات تعريف الكة الصحيحة بالقرب من الفاصل.
   • مع مصدر الكثافة: يؤدي إدراج $S_n$ إلى تقليل كل من تدفقات الحرارة الأيونية ($P_i$) والإلكترونية ($P_e$)، مما يحقق توافقًا كميًا مع تقديرات ASTRA التجريبية (ضمن 20% لـ $P_i$).
   • المساهمة الديامغناطية: بالقرب من انتقال L-H، لا ينشأ تدفق الحرارة الأيوني ($Q_i$) من النقل الاضطرابي فحسب، بل ومن مساهمة ديامغناطية كبيرة (حوالي 30% من إجمالي القدرة الأيونية). يسلط هذا الض p الضوء على ضرورة المعالجة المتسقة للتدفقات الاضطرابية والديامغناطية، وهو ما لا يمكن للنماذج المحلية لـ "أنبوب التدفق" التقاطه.
   • القياس (Scaling): في غياب مصدر الكثافة، يتبع تدفق الحرارة والجسيمات الأيونية مقياس "جيرو-بوهيم" (Gyro-Bohm) الذي يهيمن عليه الحمل الناقل غير الخطي لـ $E \times B$.

الأهمية والادعاءات
يزعم البحث أن هذا العمل يشكل علامة فارقة نحو محاكاة جايروكينيتية أولية تنبؤية لعتبة قدرة انتقال L-H. فمن خلال الانتقال من النماذج المحلية إلى المحاكاة العالمية كاملة الـ $f$، توضح الدراسة أن:

• محاكاة اضطراب الحافة وSOL العالمية يمكن أن تعيد إنتاج مكونات فيزياء نمط L الحرجة، وتحديدًا هيكل $E_r$ وتدفقات الحرارة الأيونية ($Q_i$)، والتي يُعتقد أنها أساسية لتحفيز انتقال L-H.
• يتم التحكم في تطور بئر $E_r$ من خلال تفاعلات التدفق المتوسط والاضطراب غير الخطي، وهي آلية غائبة في المحاكاة التي يتم فيها فرض $E_r$ خارجيًا.
• يتطلب التنبؤ الدقيق بنقل الحافة معالجة متسقة للتدفقات الاضطرابية والديامغناطية، بالإضافة إلى نموذج دقيق لمصادر جسيمات الحافة (تأين الغاز المتعادل).
• توفر النتائج مجموعة بيانات قيمة لتطوير نماذج النقل المختزلة (مثل TGLF) من خلال تقديم أطياف التدفق الاضطرابي من محاكاة GK عالمية كاملة الـ $f$.

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