The L-H transition in tokamaks: power threshold, density minimum and toroidal-field asymmetry

إليك ملخص تقني مفصل للورقة البحثية "انتقال L–H في التوكماك: عتبة القدرة، الحد الأدنى للكثافة، وعدم تماثل المجال الطوري" لـ دي لوكا وآخرون.

1. بيان المشكلة

يعد انتقال L–H (من النمط المنخفض إلى النمط العالي) في أجهزة التوكماك ظاهرة حرجة لتحقيق طاقة الاندماج التجاري، وتتميز بتثبيط مفاجئ للاضطراب الحوافي ومضاعفة زمن احتجاز الطاقة. على الرغم من مرور أكثر من 40 عاماً من البحث، لا تزال الآلية الفيزيائية الكامنة وراء هذه الظاهرة محل نقاش. تشمل الملاحظات التجريبية الرئيسية التي تفشل النماذج النظرية الحالية في تفسيرها في آن واحد ما يلي:

• عتبة القدرة غير الرتيبة: تعتمد القدرة المطلوبة للانتقال ($P_{LH}$) بشكل غير رتيب على كثافة البلازما، حيث تظهر "حداً أدنى للكثافة".
• عدم تماثل المجال الطوري: يحدث الانتقال عند عتبة قدرة أقل عندما يشير انزياح $\nabla B$ الأيوني من المركز نحو نقطة X (التكوين "المواتي") مقارنة بالاتجاه المعاكس ("غير المواتي").
• الافتقار إلى القدرة التنبؤية: بينما لاحظت المحاكاة العالمية المدفوعة بالتدفق حدوث تشعبات، إلا أنها تفتقر إلى القدرات التنبؤية القائمة على المبادئ الأولية والتي تتوافق مع قوانين القياس التجريبية عبر مختلف الأجهزة والتكوينات.

2. المنهجية

استخدم المؤلفون محاكاة ثلاثية الأبعاد، مدفوعة بالتدفق، وثنائية السوائل باستخدام كود GBS، الذي يحل معادلات براغينسكي المخفضة للانزياح.

• الهندسة: أُجريت المحاكاة في هندسة "أحادية النقطة منخفضة" (lower-single-null) ذات تصريف، وهي نموذجية لأجهزة التوكماك الحديثة.
• المعلمات: حاكت الإعدادات معلمات نمط L النموذجية (مثل TCV، AUG، ITER) مع تصادمية حافة عالية ($\lambda_e/L_\parallel \lesssim 0.1$)، مما يبرر استخدام النموذج المائع. شملت المعلمات الرئيسية $R/\rho_s = 500$، و $\nu = 5 \times 10^{-2} c_s/R$، و $\beta = 10^{-4}$.
• بروتوكول المحاكاة: بدءاً من حالة شبه مستقرة لنمط L، تم رفع القدرة المدخلة فجأة بمعامل قدره $\sim 1.4$ (رفع تدريجي للقدرة). أُجريت المحاكاة لكل من اتجاهات المجال الطوري "المواتية" ($\sigma_T = -1$) و"غير المواتية" ($\sigma_T = +1$).
• التحليل النظري: أجرى المؤلفون حسابات شبه خطية لتدفق الزخم المضطرب لاستخلاص شروط توليد التدفق المتوسط. كما قاموا بتحليل آليات كسر التماثل واشتقاق قوانين قياس قائمة على المبادئ الأولى لعتبة القدرة، والحد الأدنى للكثافة، والقدرة الدنيا.

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

أ. ديناميكيات المحاكاة

• التكوين المواتي: عند حقن القدرة، حدث تشعب سريع في النقل ($\sim 5 R/c_s$)، مما أدى إلى زيادة بمقدار $\sim 2\times$ في زمن احتجاز الطاقة ($\tau_E$). نتج ذلك عن التعمق التلقائي لآبار المجال الكهربائي الشعاعي ($E_r$)، مما أدى إلى تثبيط النقل المضطرب وتدريج ملف ضغط الحافة.
• التكوين غير المواتي: عند نفس مستوى القدرة، لم يحدث أي انتقال. ولم يتحقق الانتقال إلى نمط H إلا عند مستويات قدرة أعلى بكثير، مما يؤكد عدم تماثل المجال الطوري.
• الآلية: يتم دفع الانتقال بواسطة التأثيرات الكهرومغناطيسية. وتحديداً، تولد اضطرابات موجة الانزياح (DW) ذات $\beta$ المحدودة، والناتجة عن تأثيرات $\beta$ المحدودة، تدفق $E \times B$ مقصراً عبر تدفق الزخم المضطرب.
  • أدى إزالة التأثيرات الكهرومغناطيسية ($\beta=0$) أو إزالة مكون $E \times B$ الزوني إلى منع الانتقال.
  • ينشأ تدفق الزخم من التفاعل بين موجات الانزياح المقاومة والقصورية، والتي ينظمها التصادم ($\nu$) والقص المغناطيسي.

ب. كسر التماثل وعدم التماثل

تحدد الورقة أن عدم تماثل المجال الطوري ينشأ من كسر تماثل عكس الزمن بسبب أربعة شروط متزامنة:

1. تصادمية محدودة ($\nu \neq 0$).
2. قص $E \times B$ محدود ($\Omega_E \neq 0$).
3. قص مغناطيسي محدود ($\hat{s} \neq 0$).
4. عدم تماثل هندسي (أعلى-أسفل) (هندسة المصرف/Divertor).

توضح النظرية شبه الخطية أن المعلم الحرج للانتقال، $\beta^*_\nu$ (المتعلق بالتصادمية)، يكون أقل في التكوين المواتي ($\beta^*_\nu \approx 0.07$) منه في التكوين غير المواتي ($\beta^*_\nu \approx 0.11$). يؤدي هذا إلى توقع عامل عدم تماثل في القدرة قدره $\sim 2$، وهو ما يتوافق مع الملاحظات التجريبية.

ج. قوانين القياس القائمة على المبادئ الأولى

اشتق المؤلفون قوانين قياس تحليلية لانتقال L–H بناءً على النظرية شبه الخطية:

• فرع الكثافة العالية: $P_{LH} \propto a R^{1.1} B^{0.6} n^{1.05}$ (الأسات التقريبية). وهذا يتوافق مع قياس ITPA التجريبي ($R^2 = 0.87$) بشكل أفضل من الملاءمة التجريبية القياسية ($R^2 = 0.83$).
• فرع الكثافة المنخفضة: تم اشتقاقه لكل من الأنظمة التصادمية والأنظمة المحدودة بالغمد (sheath-limited)، مما أظهر توافقاً جيداً مع بيانات TCV و AUG.
• الحد الأدنى للكثافة ($n_{min}$) والقدرة الدنيا ($P_{min}$): يحدث التقاطع بين الفرعين عندما تكون التأثيرات القصورية والمقاومة متقاربة ($b_\mu \sim 1$). تتوافق القياسات المشتقة لـ $n_{min}$ و $P_{min}$ بشكل وثيق مع البيانات التجريبية للأجهزة ذات الجدران المعدنية (جدار شبيه بـ ITER).

4. التنبؤات للأجهزة المستقبلية

باستخدام القياسات المشتقة من المبادئ الأولى، قدم المؤلفون تنبؤات لـ ITER و SPARC:

• ITER (مجال كامل، $I_p=15$ MA، $B=5.3$ T):
  • $n_{min}$ المتوقع $\approx 5.8 \times 10^{19} \, \text{m}^{-3}$.
  • $P_{min}$ المتوقع $\approx 44 \pm 11$ MW.
  • الاستنتاج: من المتوقع أن يحقق ITER نمط H ضمن حد التسخين المساعد المخطط له وهو 73 MW.
• SPARC ($I_p=8.7$ MA، $B=12.2$ T):
  • $n_{min}$ المتوقع $\approx 2.6 \times 10^{20} \, \text{m}^{-3}$.
  • $P_{min}$ المتوقع $\approx 27 \pm 7$ MW.
  • الاستنتاج: الهامش ضيق بالنسبة لسعة التسخين المخطط لها وهي 25 MW، مما يؤكد المخاوف التي أثيرت في الأدبيات السابقة ولكن تم استخلاصها هنا من المبادئ الأولى.

5. الأهمية

• حل مشكلة دامت 40 عاماً: توفر الورقة تفسيراً فيزيائياً موحداً لآلية انتقال L–H، وتحديداً بربطها باضطرابات موجة الانزياح الكهرومغناطيسية والتوليد التلقائي لقص $E \times B$.
• تفسير عدم التماثل: نجحت الورقة في تفسير عدم تماثل اتجاه المجال الطوري من خلال كسر تماثل عكس الزمن، وهي ميزة غالباً ما تغيب في النماذج السابقة.
• القدرة التنبؤية: قوانين القياس المستمدة ليست مجرد ملاءمات تجريبية بل هي مشتقة من المبادئ الأولى. وقد نجحت في إعادة إنتاج البيانات التجريبية عبر أجهزة متعددة (TCV، AUG، JET، إلخ) وتتنبأ بأداء الأجهزة في الخطوة التالية (ITER، SPARC) بدقة عالية.
• التحقق من صحة النماذج المائعة: يثبت هذا العمل أن النماذج المائعة (معادلات براغينسكي) كافية لالتقاط الفيزياء الأساسية لانتقال L–H في الأنظمة ذات التصادمية العالية ذات الصلة بالمفاعلات الحالية والمستقبلية.