A simple model of current ramp down in the ITER tokamak

إليك ملخص تقني مفصل للورقة البحثية بعنوان "نموذج بسيط لخفض تيار الاندفاع في توكاماك ITER" للباحث ر. فيتزباتريك.

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

يتطلب الإغلاق الآمن لتفريغ التوكاماك خفضاً محكوماً لتيار البلازما حلقي الشكل ($I_p$) للتخلص من الطاقة المغناطيسية القطبية الهائلة المخزنة في النظام (حوالي 2 جيجا جول في ITER). وتتمثل العقبة الرئيسية أمام الإغلاق الآمن في استثارة عدم الاستقرار الماكروسكوبي لـ ديناميكا الموائع المغناطيسية (MHD) مع تطور ملفات التيار ومعامل الأمان ($q$).

• أنماط عدم الاستقرار الرئيسية: أنماط الكينك (kink) الخارجية المثالية وأنماط التمزق (tearing modes).
• المخاوف المحددة: تم تحديد نمط التمزق الكلاسيكي $m=2/n=1$ باعتباره الأكثر خطورة. فإذا نما هذا النمط بشكل كبير بما يكفي، يمكنه أن "ينغلق" (lock) على وعاء الفراغ (بسبب التفاعل مع مجالات الخطأ والتيارات الدوامية)، مما يؤدي إلى اضطراب كبير (disruption).
• الفجوة المعرفية: اعتمدت الدراسات السابقة غالباً على مخططات استقرار $q_a$-$l_i$ (معامل الأمان عند الحافة مقابل الحث الداخلي المعياري) التي تستخدم ملفات تيار مصطنعة لا تعكس فيزياء التوكاماك الفعلية. وهناك حاجة لتقييم ذاتي الاستقرار أثناء مرحلة خفض الاندفاع الديناميكية.

2. المنهجية

يستخدم المؤلف نموذج هندسة أسطوانية مبسط لمحاكاة خفض تيار الاندفاع في ITER. وبينما يعد النموذج الهندسي الأسطواني تقرييباً (يتجاهل تشكيل البلازما ونقاط X)، إلا أنه يسمح بحساب سريع وذاتي الاستقرار للتطور والاستقرار.

المكونات الأساسية للنموذج:

• التطور الكهرومغناطيسي: يحل معادلة الحث المقترنة بقانون أوم ($E = \eta j$) وقانون أمبير. يتم حساب المقاومة النوعية ($\eta$) باستخدام صيغة سبيتزر، والتي تعتمد على درجة حرارة الإلكترونات ($T_e$).
• النقل الحراري: يحل معادلة النقل الحراري للإلكترونات بما في ذلك التسخين الأومي، والتسخين المساعد، والانتشار الحراري ($\chi_\perp$). يفترض النموذج كثافة إلكترونات ثابتة ($n_e$) ويتجاهل تيار بوتستراب (وهو أمر مبرر في مرحلة خفض الاندفاع المتأخرة حيث تكون التدرجات منخفضة).
• معادلة الدائرة: تربط البلازما بالملف المركزي عبر الحث المتبادل، مما يحكم تطور جهد الحلقة والتيار.
• ملفات الاندفاع:
  • التيار: يتم خفضه بشكل خطي تقريباً عبر زمن محدد (مثلاً 60 ثانية).
  • نصف القطر الصغير: يتم تقليله في الوقت نفسه (سحق البلازما ضد نقطة X) للحفاظ على قيود معامل الأمان عند الحافة.
• تحليل استقرار MHD:
  • أنماط الكينك المثالية: يتم تقييمها باستخدام معيار نيوكومب للاستقرار ($\Delta_{ideal}$).
    An أنماط التمزق: يتم تقييمها باستخدام مؤشر استقرار التمزق ($\Delta_{tear}$). ومن المهم أن النموذج يتضمن تأثيرات التثبيت الناتجة عن متوسط الانحناء المغناطيسي (نظرية Glasser-Greene-Johnson) لحساب عتبة حرجة ($\Delta_{crit}$). يكون النمط غير مستقر فقط إذا كان المؤشر الفعال $\Delta_{eff} = \Delta_{tear} - \Delta_{crit} > 0$.
  • انغلاق النمط (Mode Locking): يتم تطبيق معيار لتحديد ما إذا كان عرض الجزيرة المغناطيسية المشبع ($W_{sat}$) يتجاوز عرض الانغلاق الحرج ($W_{crit}$). إذا كان $W_{sat} > W_{crit}$، ينغلق النمط على الجدار، مما يؤدي إلى اضطراب.

3. المساهمات الرئيسية

• تقييم استقرار ذاتي الاستقرار: على عكس النهج التجريبي السابق، يقوم هذا النموذج بتطوير ملفات التيار، ودرجة الحرارة، ومعامل الأمان بشكل ديناميكي ومتزامن، مما يوفر رؤية أكثر واقعية للاستقرار خلال مرحلة خفض الاندفاع العابرة.
• تحديد اعتماد درجة الحرارة: توضح الدراسة صراحة أن استقرار نمط التمزق $m=2/n=1$ حساس للغاية لدرجة حرارة البلازما الأولية.
• نقد مخططات $q_a$-$l_i$: تقدم الورقة دليلاً على أن مخططات $q_a$-$l_i$ قد تكون مضللة؛ حيث يمكن لسيناريوهين لهما نفس قيم $q_a$ و $l_i$ الأولية أن يعطيا نتائج استقرار مختلفة تماماً بناءً على ملف درجة الحرارة وتاريخ التسخين.
• التنبؤ بانغلاق النمط: يدمج النموذج نموذجاً مبسطاً لانغلاق النمط للتنبؤ بعتبات الاضطراب بناءً على عرض الجزيرة، بدم معدلات النمو الخطي فقط.

4. النتائج

أجرى المؤلف أربع عمليات محاكاة بمتغيرات مختلفة:

• المحاكاة 1 (بلازما باردة، بدون تسخين مساعد):

  • الإعداد: $q_a = 3.3$، بدون تسخين مساعد ($f_{aux}=0$)، درجة حرارة أولية $T_e \approx 1$ كيلو إلكترون فولت.
  • النتيجة: كان نمط التمزق $2/1$ غير مستقر منذ البداية ($\Delta_{eff} > 0$). تجاوز عرض الجزيرة المشبع عتبة الانغلاق فوراً.
  • الاستنتاج: من المرجح أن يؤدي خفض اندفاع البلازما الباردة إلى اضطراب فوري.

• المح المحاكاة 2 (بلازما ساخنة، مع تسخين مساعد):

  • الإعداد: نفس إعداد المحاكاة 1، ولكن مع تسخين مساعد كبير ($f_{aux}=4$)، مما أدى إلى $T_e \approx 4$ كيلو إلكترون فولت في البداية.
  • النتيجة: كان النمط مستقراً في البداية بسبب استقرار الانحناء الأعلى عند درجات الحرارة المرتفعة. أصبح غير مستقر في منتصف الاندفاع، لكن عرض الجزيرة المشبع ظل أقل بكثير من عتبة الانغلاق.
  • الاستنتاج: عملية خفض الاندفاع لمدة 60 ثانية المخطط لها هي ممكنة تقنياً بشرط بقاء البلازما ساخنة بما يكفي (على سبيل المثال، عبر تسخين جسيمات ألفا أو التسخين المساعد) عند بدء الاندفاع.

• المحاكاة 3 (اندفاع سريع):

  • الإعداد: تم مضاعفة معدل الاندفاع مقارنة بالمحاكاة 2.
  • النتيجة: أصبح النمط غير مستقر، ونما عرض الجزيرة بشكل كبير، مقترباً من عتبة الانغلاق.
  • الاستنتاج: من المرجح أن يكون الوضع حميداً ولكنه محفوف بالمخاطر؛ فالنظام قريب من حافة عدم الاستقرار.

• المحاكاة 4 (اندفاع سريع جداً):

  • الإعداد: تم مضاعفة معدل الاندفاع أربع مرات مقارنة بالمحاكاة 2 (حوالي 15 ثانية إجمالاً).
  • النتيجة: اندفع النمط ليصبح غير مستقر بسرعة. تجاوز عرض الجزيرة المشبع عتبة الانغلاق، مما أدى إلى توقع انغلاق النمط وحدوث اضطراب.
  • الاستنتاج: يُتوقع أن تؤدي عمليات خفض الاندفاع الأسرع بكثير إلى حدوث اضطرابات.

5. الأهمية

• الجدوى التشغيلية لـ ITER: تؤكد الدراسة أن خفض التيار لمدة 60 ثانية المخطط له في ITER هو أمر ممكن تقنياً، بشرط الحفاظ على درجات حرارة بلازما عالية عند بدء الاندفاع. وهذا يدعم استراتيجية بدء خفض الاندفاع بينما لا تزال البلازما في وضع H-mode أو خاضعة لتسخين جسيمات ألفا.
• هوامش الأمان: تضع الدراسة تحذيراً واضحاً بأن محاولة تقصير وقت خفض الاندفاع بشكل كبير (لزيادة توافر الماكينة) تحمل مخاطر عالية في إطلاق أنماط تمزق $2/1$ التي تنغلق على الوعاء، مما يسبب اضطراباً.
• رؤية منهجية: تسلط الورقة الضوء على حدود مخططات الاستقرار الثابتة ($q_a$-$l_i$) للسيناريوهات الديناميكية، وتدعو إلى استخدام نماذج التطور الزمني التي تأخذ في الاعتبار الملفات الحرارية واستقرار الانحناء.
• الاتجاهات المستقبلية: يقترح المؤلف أنه بينما يوفر النموذج الأسطواني رؤى قيمة، يجب أن تدمج الأعمال المستقبلية أكواد التوازن الحلقي (مثل RAPTOR) وأكواد الاستقرار (مثل RDCON) لأخذ تشكيل البلازما ونقاط X في الاعتبار للحصول على تنبؤات أكثر دقة.