First-Principles Explanation of the Drift Configuration Dependence of the Radial Electric Field and High-Confinement Access in Tokamaks

ملخص تقني: تفسير من المبادئ الأولى لاعتماد تكوين الانجراف على المجال الكهربائي الشعاعي والوصول إلى نظام الحجز العالي في التوكاماك

بيان المشكلة
لا يزال أصل الاختلاف الكبير في عتبة قدرة نظام الحجز العالي (H-mode) بين تكوينات انجراف أيونات $\nabla B$ المواتية (fav) وغير المواتية (unfav) في أجهزة التوكاماك غير محلول. تجريبياً، تتطلب التفريغات غير المواتية أكثر من ضعف عتبة القدرة للتفريغات المواتية، وهو اختلاف مرتبط بوجود بئر أعمق للمجال الكهربائي الشعاعي ($E_r$) بالقرب من الفاصل (separatrix) في التكوين المواتي قبل انتقال L-H. وبينما تم اقتراح تأثيرات نيوكلاسيكية (NC)، تشير التجارب الأخيرة في جهاز "أزديكس أبراوب" (ASD-Upgrade/AUG) إلى أن هذه التأثيرات لا يمكنها تفسير فروق $E_r$ الملحوظة بالكامل، مما يشير إلى دور حاسم للفيزياء غير النيوكلاسيكية، وتحديداً التفاعلات غير الخطية بين الاضطراب والتدفق المتوسط. تفشل النماذج المختزلة الحالية والمنهجيات الخطية في التقاط الاقتران الذاتي غير الخطي بين توازنات المقاطع (profiles)، و$E_r$، والتدفقات، والاضطراب في الهندسات ذات النقاط المفرغة (diverted geometries) الواقعية.

المنهجية
يقدم المؤلفون أول دراسة من المبادئ الأولى باستخدام نظرية الجايروكينتيك (GK) لمقارنة التفريغات المواتية وغير المواتية في جهاز AUG باستخدام النسخة الطيفية من كود GENE-X. تستخدم المحاكاة معادلة "فالوس-جايروكينتيك" (GK Vlasov) كاملة الـ $f$، وكهرومغناطيسية، وطويلة الطول الموجي، وتصادمية مع نقاط X.

• الأساس التجريبي: تهدف الدراسة إلى نمذجة تفريغات الدوتيريوم في نظام L-mode المتطابقة (رقم 36983 للمواتي و#37375 لغير المواتي) بكثافة إلكترونات متوسطة الخط $n_e \simeq 2.7 \times 10^{19} \, \text{m}^{-3}$، وقدرة تسخين إلكتروني حلقي (ECRH) تبلغ 600 كيلوواط، ومجال مغناطيسي حلقي $|B_\phi| = 2.5$ تسلا، وتيار بلازما $|I_p| = 0.8$ ميجا أمبير.
• إعداد المحاكاة: تم توليد التوازن غير المواتي عن طريق عكس اتجاه المجال المغناطيسي الحلقي ($B_\phi$) بالنسبة للحالة المواتية، مما يؤدي إلى عكس اتجاه انجراف الأيونات $\nabla B$ مع الحفاظ على هندسة مغناطيسية ومعلمات بلازما متطابقة. يحاكي مصدر الكثافة تأين النيوترونات بالقرب من الفاصل (separatrix). تعمل عمليات المحاكاة حتى الوصول إلى حالة شبه مستقرة (حوالي 0.5 ميجاواط عبر الفاصل)، مع جمع الإحصائيات على نافذة زمنية قدرها 0.1 مللي ثانية.
• التحليل: تتحقق الدراسة من النتائج مقابل بيانات تشتت دوپلر المرتد (DBS) ومطيافية إعادة تركيب تبادل الشحنة (CXRS). كما تقوم بتحليل تطور متوسط الطاقة الحركية المستقطبة لسطح التدفق (FSA) وإجمالي إجهاد السرعة $\Pi$ (وهو تعميم لـ Reynolds stress في نموذج full-f) لتوضيح الآليات التي تقود $E_r$ والتدفقات المستقطبة.

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

1. التحقق: تعيد محاكاة GENE-X إنتاج مقاطع الكثافة ($n_e$) ودرجات الحرارة ($T_i, T_e$) والسرعة الحلقية ($U_{\phi i}$) التجريبية باتفاق جيد. ومن الأهمية بمكان أن تؤكد المحاكاة أن هذه المقاطع غير حساسة تقريباً لاتجاه الانجراف، مما يعزل السرعة المستقطبة ($U_{\theta i}$) و$E_r$ كعوامل تمييز رئيسية.
2. المجال الكهربائي الشعاعي ($E_r$) والتدفقات المستقطبة: تعيد المحاكاة إنتاج بئر $E_r$ الأكثر عمقاً الملاحظ تجريبياً في التكوين المواتي. يرتبط هذا العمق مباشرة بسعة أكبر بكثير للتدفق الأيوني المستقطب ($U_{\theta i}$) في الحالة المواتية بالقرب من الفاصل ($\rho_{pol} \sim 0.995$). في الحالة المواتية، يتجاوز $U_{\theta i}$ التوقعات النيوكلاسيكية، مما يشير إلى تسارع قوي مدفوع بالاضطراب، بينما في الحالة غير المواتية، يظل $U_{\theta i}$ أقرب إلى المستويات النيوكلاسيكية.
3. آلية تسريع التدفق: يتم دفع بئر $E_r$ الأعمق في الحالة المواتية بواسطة تعزيز انتقال الطاقة غير الخطي من الاضطراب إلى التدفقات المستقطبة المتوسطة. ويحكم هذا الانتقال التدرج الشعاعي لإجمالي إجهاد السرعة $\Pi$.
   • التكوين المواتي: يخلق توجيه تدفق $E \times B$ بالنسبة لنقطة X على جانب المجال المنخفض (LFS) قصاً تدفقياً موضعياً يتوافق مع إمالة الدوامات الناتجة عن القص المغناطيسي (MSI). يعمل هذا التوافق على تعديل ميل الدوامات المضطربة وتوزيع $\Pi$ المستقطب شعاعياً، مما يولد تدرجاً شعرياً قوياً سالباً لـ $\Pi$. هذا التدرج يدفع تسارعاً كبيراً للتدفقات المستقطبة.
   • التكوين غير المواتي: يعارض قص التدفق إمالة الدوامات (MSI)، مما يؤدي إلى تدرج أضعف بكثير في $\Pi$ وبالتالي تسارع أضعف للتدفق.
4. شدة الاضطراب: على عكس تسارع التدفق، وُجد أن شدة الاضطراب (سعة تقلبات الجهد والكثافة) أعلى في التكوين غير المواتي.
5. كفاءة انتقال الطاقة: يشير العمل المبذول بواسطة الإجهاد ($\Pi_W$) إلى انتقال طاقة أكبر وسالب بشكل ملحوظ من الاضطراب إلى الطاقة الحركية المستقطبة في الحالة المواتية مقارنة بالحالة غير المواتية. في الحالة المواتية، يهيمن النقل المضطرب للطاقة الحركية المتوسطة داخل بئر $E_r$، بينما في الحالة غير المواتية، يكون هذا الانتقال في حده الأدنى.

الأهمية والادعاءات
يدعي البحث أنه يقدم أول تفسير جايروكينتيكي (GK) مُتحقق منه ومتسق ذاتياً من المبادئ الأولى لكيفية تحكم تكوين الانجراف في الديناميكيات غير الخطية للمقاطع، و$E_r$، والتدفقات، والاضطراب، مما يحدد عتبة قدرة نظام H-mode.

يخلص المؤلفون إلى أن الاختلاف في الوصول إلى H-mode لا يعود فقط إلى التأثيرات النيوكلاسيكية أو شدة الاضطراب البسيطة، بل ينشأ من كفاءة انتقال الطاقة غير الخطي بين الاضطراب والتدفقات المتوسطة.

• في التكوين المواتي، يعزز التوافق المحدد بين قص التدفق والهندسة المغناطيسية من التدرج الشعاعي لإجهاد السرعة $\Pi$، مما يؤدي إلى توليد تدفق مستقطب أقوى، وبئر $E_r$ أعمق، وقص $E \times B$ أقوى، مما يسهل انتقال L-H.
• في التكوين غير المواتي، يتم كبح هذه الآلية، مما يؤدي إلى بئر $E_r$ أقل عمقاً ومستويات اضطراب أعلى، مما يتطلب قدرة تسخين أعلى لتحفيز الانتقال.

تؤسس هذه الدراسة أن اعتماد عتبة قدرة H-mode (وI-mode) على تكوين الانجراف متجذر أساساً في التفاعل غير الخطي بين الاضطراب والتدفقات المتوسطة بوساطة تدرج إجهاد السرعة $\Pi$ الكامل (full-f). ويشير المؤلفون إلى أن الاتجاهات المماثلة الملاحظة في محاكاة TCV تشير إلى عالمية هذه النتائج.